Formation mechanism for the amplitude of interannual climate variability in subtropical northern hemisphere: relative contributions from the zonal asymmetric mean state and the interannual variability of SST

The Amplitude Interannual climate Variability (AIV) differs among the subtropical northern hemisphere, and the Western North Pacific (WNP) was claimed to exhibit the largest AIV. The robustness of the AIV pattern is investigated in this study with different atmospheric variables from multiple datasets. As consistently shown by the interannual variance patterns of precipitation and circulation, the AIV over subtropical northern hemisphere closely follows the mean state of precipitation, where higher (lower) AIV is located at moister (drier) regions. The largest AIV is seen over the broad area from South Asia to WNP, followed by a secondary local maximum over the Gulf of Mexico. To further investigate the formation mechanism for the AIV pattern, numerical simulations are performed by Community Atmosphere Model version 4 (CAM4). The zonal asymmetry of AIV is reduced if the interannual SST variability is removed, and it almost disappears if the zonal asymmetry of SST mean state is removed. The results suggest that the zonal asymmetric AIV pattern primarily originates from the zonal asymmetric SST mean state, and it is amplified by the interannual SST variability. The atmospheric convection–circulation feedback plays a key role in connecting the AIV with the mean state precipitation. In both observation and CAM4 simulations, stronger (weaker) convection–circulation feedback is seen in moister (drier) regions. By modulating the mean state precipitation and the associated intensity of convection–circulation feedback, the zonal asymmetric SST mean state accounts for the zonal asymmetry of AIV in the subtropical northern hemisphere.


Introduction
The Amplitude of Interannual climate Variability (AIV) over subtropical western North Pacific (WNP) was claimed to be the largest over the subtropical northern hemisphere (Lu 2001;Sui et al. 2007;Wu and Zhou 2008;Chung et al. 2011). The strong interannual variability of atmospheric circulation over WNP has a profound influence on East Asia climate (Zhong et al. 2010), and it plays a crucial role in connecting the tropical Sea Surface Temperature (SST) anomalies with East Asian climate (Zhang et al. 1999;Chang et al. 2000;Chen et al. 2015a). As is widely noted by previous studies, an anomalous anticyclone over WNP favors anomalous water vapor transport into Yangtze River valley (Zhang et al. 1999;Zhou and Yu 2005), leading to disastrous flood such as in the summer of 1998 (Ma et al. 2015).
Although it was claimed that the amplitude of AIV over WNP is the largest among subtropical northern hemisphere, the only metric adopted by previous studies is geopotential height, either at 850 hPa (Lu 2001) or at 500 hPa (Sui et al. 2007;Wu and Zhou 2008;Chung et al. 2011). As shown in Fig. 1a, the AIV of 850 hPa geopotential height (H850) is much lower in the tropics than in the mid-latitudes, especially it approximates to zero near the equator. Similar results are obtained if 500 hPa geopotential height is used (figure not shown). This contradicts with the common Abstract The Amplitude Interannual climate Variability (AIV) differs among the subtropical northern hemisphere, and the Western North Pacific (WNP) was claimed to exhibit the largest AIV. The robustness of the AIV pattern is investigated in this study with different atmospheric variables from multiple datasets. As consistently shown by the interannual variance patterns of precipitation and circulation, the AIV over subtropical northern hemisphere closely follows the mean state of precipitation, where higher (lower) AIV is located at moister (drier) regions. The largest AIV is seen over the broad area from South Asia to WNP, followed by a secondary local maximum over the Gulf of Mexico. To further investigate the formation mechanism for the AIV pattern, numerical simulations are performed by Community Atmosphere Model version 4 (CAM4). The zonal asymmetry of AIV is reduced if the interannual SST variability is removed, and it almost disappears if the zonal asymmetry of SST mean state is removed. The results suggest that the zonal asymmetric AIV pattern primarily originates from the zonal asymmetric SST mean state, and it is amplified by the interannual SST variability. The atmospheric convection-circulation feedback plays a key role in connecting the AIV with the mean state precipitation. In both observation and CAM4 simulations, stronger (weaker) convection-circulation feedback is seen in moister (drier) regions. By modulating the mean state precipitation and the knowledge that the strongest source of interannual climate variability is El Nino-Southern Oscillation (ENSO) which originates in the equatorial Pacific. It may not be suitable to measure the AIV by using geopotential height, and an intercomparison of multiple atmospheric variables is needed to check the robustness of the subtropical AIV pattern.
Great attention has been paid on how WNP climate variability is forced by tropical SST anomalies. Five key oceanic regions are identified to be responsible for the interannual climate variability over WNP. These five regions include the equatorial Pacific Ocean (Chang et al. 2000;Wang et al. 2000Wang et al. , 2013Xiang et al. 2013), tropical Indian Ocean Li et al. 2008;Xie et al. 2009;Wu et al. 2010), tropical Atlantic Ocean (Lu and Dong 2005;Rong et al. 2010;Hong et al. 2014;Chen et al. 2015b), the vicinity of the maritime continent (Lu et al. 2006;Sui et al. 2007;Wu et al. 2009;Chung et al. 2011), and WNP (Wang et al. 2000;Wu et al. 2010;Wang et al. 2013;Xiang et al. 2013). The detailed mechanisms on how the WNP climate is modulated by these five regions are reviewed in He et al. (2015a).
Although a good understanding has been obtained on how WNP climate is forced by SST, it is still not clear why the AIV over WNP is larger than other subtropical regions, The formation mechanism for the spatial pattern of AIV has not been well explained. According to previous studies, two factors residing in SST may be responsible for the pattern of AIV.
First, the pattern of SST mean state, especially its zonal asymmetric pattern. Convection-circulation feedback is a key process which amplifies the interannual climate variability over WNP, i.e., an initial circulation anomaly results in a rainfall anomaly, and the latent heating anomaly associated with the rainfall anomaly in turn enhances the circulation anomaly (Xie et al. 2009;Jin et al. 2013;Xiang et al. 2013;He et al. 2015b). The intensity of the convection-circulation feedback depends on the abundance of mean state precipitation (Huang and Xie 2015), and the mean state precipitation is regulated by the SST (Seager et al. 2003).
Second, the interannual SST variability. The amplitude of SST variability is not uniformly distributed but characterized by a distinct spatial pattern (Fig. 1b). Larger SST variability is not necessarily associated with larger local atmospheric variability, because atmospheric variability is impacted by not only local SST anomaly but also non-local SST anomalies. The complicated relationship between the interannual SST variability and the AIV of the atmosphere calls for in-depth study.
More and more studies begin to pay attention to the response of interannual climate variability to global warming (Lu and Fu 2010;Cai et al. 2014;Lee et al. 2014;Zhou et al. 2014;Chen et al. 2016), but the fundamental formation mechanism of the current AIV pattern is yet not well understood. This study aims at answering the following two questions: (1) Is it a robust feature that WNP has a larger AIV than other subtropical regions? (2) How is the AIV pattern in the subtropical northern hemisphere regulated by SST, and what are the relative contributions from the mean state and the interannual variability of SST? The first question can be addressed by observational analysis, and numerical simulations are needed to address the second question.
The rest of this article is organized as follows. The descriptions on the data, model and methods are presented in Sect. 2. Analyses on the observed interannual climate variability is documented in Sect. 3. The relative contributions from the mean state and the interannual variability of SST are demonstrated by model experiments in Sect. 4. Finally, a summary is presented in Sect. 5. The AIV (shading) and the mean state (contours) for multiple variables. a 850 hPa geopotential height (unit m), b SST (unit °C), c precipitation (unit mm day −1 ), d 500 hPa vertical velocity (unit Pa s −1 ). The AIV shown here is scaled by its domain average and is dimensionless. The pattern correlation coefficient between the AIV and the mean state for each variable is marked on the upper-right corner of each panel 2 Data, model and methods We focus on the boreal summers from 1979 to 2014 in the observational analyses. The SST dataset adopted in this study is the Extended Reconstructed SST version 3 (Smith et al. 2008). The precipitation data adopted in this study include the GPCP (Adler et al. 2003) and CMAP (Xie and Arkin 1997) datasets. The atmospheric geopotential height and vertical velocity are adopted from the datasets of ERAIM (Dee et al. 2011), NCEP1 (Kalnay et al. 1996), and NCEP2 (Kanamitsu et al. 2002). All of the above datasets are referred to as "observation" for brevity. Multiple atmospheric datasets are adopted to address the observational uncertainty. The monthly values for June, July and August are averaged into the seasonal mean in summer, referred to as "JJA".
The atmospheric general circulation model employed in this study is the Community Atmosphere Model version 4 (CAM4) (Neale et al. 2013), under a finite volume dynamic core with a horizontal resolution of 1.9° in latitude and 2.5° in longitude. In the control (CTL) simulation, CAM4 is forced by observed interannual varying SST from 1900 to 2000, following the design of Atmospheric Model Intercomparison Project (AMIP) simulation (Gates et al. 1999). The observational SST data used as boundary forcing of CAM4 is created by Hurrell et al. (2008). It is a merged dataset between Hadley Center sea ice and SST dataset version 1 (HadISST1) and optimum interpolation (OI) SST version 2. The model outputs from 1901 to 2000 are used for analyses.
Two sensitivity simulations are also performed with CAM4, by modifying the SST without other changes ( Table 1). The interannual variability component of the observational SST is extracted by applying an 8-year highpass Fourier filter. In the "No asymmetry" (NoAsy) simulation, the model is forced by idealized SST which is created by adding the interannual variability component to the zonal mean climatology of SST, so that the interannual SST variability is the same as CTL simulation but the zonal asymmetry of SST mean state is removed. In the "No interannual variability" (NoVar) simulation, the interannual variability of SST is removed, and the model is forced by the climatological annual cycle of SST. These two sensitivity experiments are also integrated for 101 years and the last 100 years are adopted for analyses. The difference between CTL and NoAsy simulations demonstrates the role of zonal asymmetric SST mean state, while the difference between CTL and NoVar simulations demonstrates the effect of interannual SST variability.
The amplitude of interannual variability is measured by the interannual variance, i.e., the variance of the interannual variability component of the time series. All the observational and model time series at all grid points are filtered by an 8-year high-pass Fourier filter before calculating the interannual variance. The decadal oscillations with a period greater than 8 years are removed and the filtered time series contain only interannual variability. The interannual variance of precipitation is investigated, since precipitation has direct societal impacts and its usage is recommended by previous studies (Cai et al. 2014;Huang and Xie 2015). The AIV of circulation is measured in terms of the interannual variance of 500 hPa vertical velocity (W500), in comparison with the interannual variance of 850 hPa geopotential height (H850).
Pattern correlation coefficient is adopted to evaluate the similarity between two fields. The statistic significance is tested based on Student's t distribution, where the degree of freedom is estimated following the method of Leith (1973). The auto-correlation for the spatial patterns of either precipitation or W500 ranges from 0.65 to 0.97. In this study, the number of grid points for pattern correlation is 3024, and a rigorous estimate of the degree of freedom is −0.5 × ln(0.97) × 3024 = 46. Therefore, a pattern correlation coefficient exceeding ±0.28 (±0.37) is considered significant at the 95 % (99 %) confidence level.

Amplitude of interannual climate variability in observation
Is the pattern of AIV robust among multiple meteorological variables? Figure 1c-d shows the patterns of interannual variance for precipitation and W500, respectively. It is clear that the AIVs are larger over the moister regions, including the Inter-Tropical Convergence Zone (ITCZ), the broad area from South Asia to WNP, and a narrow area over the Gulf of Mexico. The patterns of interannual variance for W500 and precipitation are similar to each other, with a pattern correlation coefficient of 0.63. In contrast, the interannual variance pattern of H850 (Fig. 1a) is much different from that of precipitation, with a low pattern correlation coefficient of −0.07. It is more reasonable to measure the AIV of atmospheric circulation in terms of W500 instead of H850 for two reasons. First, the AIV of W500 is a good indicator for the AIV of precipitation. Second, the interannual climate variability over tropics and mid-latitudes are both well captured by W500, whereas H850 fails to capture the interannual climate variability over the tropics. Is the pattern of AIV robust among multiple datasets? Since the major focus of this study is the subtropics, we show the zonal profile of 20°N-30°N averaged interannual variance for multiple datasets in Fig. 2. The high interannual precipitation variance from South Asia to WNP is agreed by both GPCP and CMAP datasets (Fig. 2a). A remarkable difference is that CMAP data shows a peak at about 90˚E over South Asia (the southern slope of Himalayas) but GPCP data does not. The local maximum over the Gulf of Mexico from 100°W to 60°W is also agreed by both GPCP and CMAP datasets. As seen in the pattern of interannual W500 variance (Fig. 2b), the ERAIM and NCEP2 datasets agree well with each other that the highest interannual W500 variance appears over South Asia instead of WNP. But in NCEP1 dataset, the interannual W500 variance is higher over WNP than South Asia. A comparison between Fig. 2a, b shows that the zonal asymmetric patterns of AIV are similar between precipitation and W500.
How is the pattern of AIV related to the mean state? By comparing the shading with the contours in Fig. 1c, it is clear that larger interannual precipitation variability is located where there is more abundant rainfall. The pattern correlation coefficient between the AIV and the mean state for precipitation is 0.70, which is statistically significant at the 99 % confidence level according to Student's t test. Similarly, the spatial pattern of AIV is also closely related to the mean state for W500 (Fig. 1d), with a pattern correlation coefficient of −0.44, also exceeding the 99 % confidence level. These evidences indicate that larger AIV is associated with ascending motion and moister mean state, while smaller AIV is associated with descending motion and drier mean state.
Based on the above observational evidences, there are observational uncertainty on whether the AIV is greater over South Asia or WNP. However, it is robust that the AIV closely follows the mean state precipitation, where larger (smaller) AIV is seen over moister (drier) regions. It is also robust that the AIV over the broad area from South Asia to WNP is greater than other subtropical regions, followed by a secondary peak over the Gulf of Mexico. The AIV is the lowest over the arid regions, especially over eastern part of Pacific and Atlantic Oceans and North Africa.
Why is the AIV greater in moister regions? As implied by previous studies, the intensity of convection-circulation feedback is greater over moister regions (Xie et al. 2009;Jin et al. 2013;Xiang et al. 2013), and stronger local convection-circulation feedback may intensify local AIV. To measure the intensity of the convection-circulation feedback, Fig. 3 shows the regression coefficient of precipitation onto local W500 at interannual time scale. The regression coefficient is greater over moister regions (the ITCZ, South Asia, WNP and western Atlantic Ocean), but smaller over drier regions (eastern subtropical oceans and North Africa). As indicated by the spatial pattern of regression coefficients, greater convection and precipitation anomalies are generated by a unit vertical velocity anomaly at moister regions (Cai et al. 2014;Huang and Xie 2015), and therefore stronger latent heating anomaly is generated by precipitation anomaly. The stronger latent heating anomaly in turn stimulates stronger circulation anomalies. Therefore, (a) (b) Fig. 2 The zonal profiles of AIV averaged within 20°N-30°N. a The AIV of precipitation derived from GPCP and CMAP datasets (unit mm 2 day −2 ). b The AIV of 500 hPa vertical velocity derived from ERAIM, NCEP1 and NCEP2 datasets (unit 10 −4 Pa 2 s −2 ). The location of the zonal belt within 20°N-30°N is shown as the red box in Fig. 1b   Fig. 3 The observational intensity of convection-circulation feedback measured by the regression coefficient of precipitation onto W500 for the 1979-2014 period [unit (mm day −1 ) (Pa s −1 ) −1 ]. The regression coefficients are multiplied by −1, and the regression coefficients significant at the 95 % confidence level are stippled the intensity of convection-circulation feedback explains why the AIV is larger over moister regions.
The above observational evidences suggest that the spatially differentiated intensity of convection-circulation feedback plays a role in linking the mean state and AIV patterns. However, interannual SST variability is also crucial for the interannual atmospheric variability. The contribution of SST variability cannot be diagnosed by comparing the observational AIV patterns between SST and atmospheric variables, since precipitation and circulation are impacted by not only local but also remote SST anomalies. To further illustrate the mechanism for the AIV pattern, model simulations are needed.

Relative contributions from mean state and variability of SST
The observed mean state and AIV are captured by the CTL simulation of CAM4. Consistent with the observed mean state (Fig. 1c), the large AIVs of precipitation over the ITCZ, South Asia, WNP and the Gulf of Mexico are well captured by CAM4 (Fig. 4a). A major discrepancy from the observation is the overestimated AIV over South Asia and the underestimated AIV around the Philippines. Therefore, the zonal asymmetry in CAM4 is slightly weaker within the subtropical Pacific but stronger than observation in the entire subtropical northern hemisphere. The close relationship between the mean state and the AIV of rainfall is captured, with a pattern correlation coefficient of 0.81, exceeding the 99 % confidence level. The mean state and AIV of W500 are also well captured (Fig. 4b). The correlation coefficients between AIV and mean state for W500 is −0.59, consistent with the observational feature that higher (lower) AIV of circulation is located at the ascending (descending) region. In all, CAM4 well captures the observed AIV, mean state, and the relationship between AIV and mean state, adding to our confidence in further analyses on the sensitivity simulations. The role of the zonal asymmetric SST mean state is demonstrated as the difference in interannual variances between CTL and NoAsy simulations (Fig. 5a). Over the subtropical northern hemisphere, the zonal asymmetric SST mean state amplifies the interannual precipitation variance over South Asia, WNP, and the Gulf of Mexico, but damps the AIV over the eastern North Pacific and centraleastern Atlantic Ocean. The AIV of precipitation at the equatorial oceans is damped by the zonal asymmetric mean state, the cause of which is out of the scope of this study. The modeling evidence shows that the zonal asymmetric SST mean state does make a substantial contribution to the zonal asymmetric pattern of AIV over subtropical northern hemisphere.
The role of interannual SST variability is demonstrated as the difference between CTL and NoVar simulations (Fig. 5b). It is clear that the SST variability amplifies the interannual precipitation variance over the deep tropics, especially along the ITCZ, the Arab Sea and the Bay of Bengal. At the subtropics north of 20°N, the SST variability still amplifies the AIV of precipitation but its effect is  5 The differences of the AIV of precipitation between the CTL simulation and the sensitivity simulations (unit mm 2 day −2 ). a CTL simulation minus NoAsy simulation. b CTL simulation minus NoVar simulation. The difference in interannual variance statistically significant at the 95 % confidence level is stippled weaker. To focus on the subtropics, the zonal profiles for the 20°N-30°N averaged AIV and mean state are shown in Fig. 6.
The CTL simulation is characterized by large AIV from South Asia to WNP, with a secondary peak over tropical Atlantic (black line in Fig. 6a). The CTL simulation more resembles CMAP data that the highest peak is seen over South Asia rather than WNP, which is different from GPCP data. Compared with CTL simulation, the highest AIV of precipitation in NoAsy simulation slightly shifts eastward from 72.5°E to 85°E (red line in Fig. 6a). Over the broad subtropical regions other than the South Asia, the zonal difference of the AIV almost disappears. The AIV decreases substantially over WNP and the Gulf of Mexico, but increases over eastern Pacific and eastern Atlantic Oceans. In the NoVar simulation (blue line in Fig. 6a), the AIV in precipitation is also characterized by zonal asymmetric pattern but the zonal asymmetry is weaker than the CTL simulation. By regressing the profile of CTL simulation onto NoVar simulation, it is found that if a factor of 1.58 is multiplied on the zonal profile of NoVar simulation (orange line in Fig. 6a), it most closely fits to the CTL simulation. This indicates the role of the SST variability is to amplify the AIV of precipitation by 58 % across the subtropics, which magnifies the zonal asymmetry.
How is the AIV connected to the mean state precipitation in the simulations? As shown in Fig. 6b, the mean state of precipitation in NoVar simulation fits well with CTL simulation, suggesting the interannual variability of SST has no impact on the mean state of rainfall over subtropical northern hemisphere. In the NoAsy simulation, the zonal asymmetry of the mean state precipitation substantially weakens. This is consistent with Seager et al. (2003) that the zonal asymmetric SST amplifies the subtropical anticyclones which favors moist (dry) condition over western (eastern) subtropical oceans. In a moister atmosphere with more mean state precipitation, greater precipitation anomaly can be stimulated by a unit circulation anomaly (Cai et al. 2014;Huang and Xie 2015). Therefore, the weakened asymmetry in the mean state rainfall explains the weakened asymmetry of the AIV across the subtropical basins in NoAsy simulation. Compared with WNP, no substantial decrease of mean state precipitation is seen over South Asia Fig. 6b, since the South Asian area shown in Fig. 6b is dominated by land (see Fig. 1b), and land surface temperature is not artificially modified in NoAsy simulation.
How is the AIV of W500 impacted by the mean state and variability of SST? Does it share similar pattern with precipitation? As shown in the difference between CTL and NoAsy simulations (Fig. 7a), the AIV of W500 increases significantly at the 95 % confidence level over the Bay of Bengal, WNP and the Gulf of Mexico, but decrease significantly over eastern North Pacific and central-eastern North Atlantic. The zonal asymmetry of the SST mean state regulates the zonal asymmetric AIV pattern of not only precipitation but also circulation. As shown in the difference between CTL and NoVar simulations (Fig. 7b), the interannual SST variability significantly amplifies the AIV of W500 over the deep tropics, especially equatorial central Pacific. Compared with the deep tropics, the effect of SST variability is weaker at north of 20°N.
The 20°N-30°N averaged zonal profiles for the interannual W500 variance are displayed in Fig. 8a. Compared with CTL simulation, the zonal asymmetry of AIV of W500 weakens in the oceanic sectors in NoAsy simulation. However, in the NoVar simulation, the zonal asymmetric pattern as seen in CTL simulation is still clear despite of smaller amplitude. Regression of the zonal profiles between CTL and NoVar simulations suggest that, the zonal profile in NoVar simulation well fits the CTL simulation if the former is multiplied by a factor of 1.39 (orange line in Fig. 8a). This suggests the interannual SST variability amplifies the AIV of W500 for 39 % over the subtropical northern hemisphere, and therefore amplifies the existing zonal asymmetry. This value of amplification is slightly smaller than that for the precipitation (Fig. 6a), possibly suggesting the direct role of SST variability on precipitation and its indirect role on W500 through changing the precipitation. In all, the zonal asymmetric AIV in the subtropics is determined by the zonal asymmetry of the SST mean state, whereas the variability of SST plays a secondary role in amplifying the zonal asymmetry.
To explain how the AIV of W500 is regulated by SST, the 20°N-30°N averaged regression coefficients of precipitation onto W500 at interannual timescale are shown in Fig. 8b. In the CTL simulation, the convection-circulation feedback is strong over South Asia, WNP and the Gulf of Mexico, but weak in the eastern Pacific and eastern Atlantic. In the NoAsy simulation, the zonal asymmetry in the intensity of convection-circulation feedback weakens substantially, and little difference is seen between the western part and eastern part of the oceans. This is caused by the weakened zonal asymmetry of the mean state precipitation (see Fig. 6b). In the NoVar simulation, the intensity of convection-circulation feedback closely resembles that of the CTL simulation, consistent with the high consistency of precipitation mean state between these two simulations (see Fig. 6b). The above evidences suggest that intensity of convection-circulation feedback highly depends on the precipitation mean state, and the spatially differentiated intensities of convection-circulation feedback are responsible for the AIV pattern of atmospheric circulation.
Based on the above numerical simulations, it can be concluded that the zonal asymmetry in the AIV over the subtropical northern hemisphere is determined by the zonal asymmetry of SST mean state, while the SST variability plays a secondary role in amplifying the zonal asymmetry. Although the SST variability has substantial effect on AIV in the equatorial regions such as ITCZ, its effect on the subtropical northern hemisphere is limited. The relative importance of the mean state and the variability of SST differs between tropics and subtropics.

Conclusion and discussion
The amplitude of interannual climate variability is characterized by a distinct spatial pattern over the subtropical northern hemisphere. The robustness of its pattern is evaluated in terms of different meteorological variables in multiple datasets. The formation mechanism for the AIV pattern is investigated through sensitivity simulations performed by CAM4. By removing the zonal asymmetric mean state and the interannual variability of SST, the relative roles of zonal asymmetric mean state and the interannual variability of SST are demonstrated. The major conclusions are summarized as follows.
1. It is a robust feature that the interannual climate variability is the strongest over the broad area from South Asia to WNP over the subtropical northern hemisphere, followed by a secondary local peak over the Gulf of Mexico. The AIV pattern of precipitation is consistent with that of circulation in terms of W500. The pattern of AIV depends strongly on the mean state pre- The zonal profile of the 20°N-30°N averaged regression coefficient of precipitation onto W500 multiplied by −1 [unit 10 2 (mm day −1 )·(Pa s −1 ) −1 ]. In both (a) and (b), the black, red and blue lines are for the CTL, NoAsy, and NoVar simulations, respectively. In (a), the orange line stands for the AIV in the NoVar simulation multiplied by a factor of 1.39 cipitation, and higher (lower) AIV is seen over moister (drier) regions. Our results show that the W500 is more suitable to measure the AIV of circulation compared with geopotential height, since the former well captures the tropical variability and the connection of circulation to precipitation. 2. The zonal asymmetric pattern of AIV over the subtropical northern hemisphere is dominated by the zonal asymmetry of SST mean state, and it is amplified by the interannual variability of SST. The zonal asymmetry of the AIV disappears if the zonal asymmetric SST mean state is removed, and it is weakened if the SST variability is removed. The interannual variability of SST plays a role in amplifying the AIV of the precipitation (by 58 %) and circulation (by 39 %), and it therefore amplifies the zonal asymmetry of the AIV. 3. Convection-circulation feedback dependent on the mean state is the key physical process connecting the mean state with AIV. Over the subtropical northern hemisphere, the zonal asymmetric pattern of AIV is well explained by the zonal asymmetric intensity of convection-circulation feedback, where stronger AIV is associated with stronger local convection-circulation feedback. In both observation and model simulations, the intensity of convection-circulation feedback is stronger (weaker) over moister (drier) regions. By modulating the zonal distribution of the mean state precipitation, the zonal asymmetric SST mean state results in the zonal asymmetric intensity of convection-circulation feedback, and accounts for the AIV pattern over the subtropical northern hemisphere.