An international joint research project, entitled Interhemispheric Coupling Study by Observations and Modelling (ICSOM), is ongoing. In the late 2000s, an interesting form of interhemispheric coupling (IHC) was discovered: when warming occurs in the winter polar stratosphere, the upper mesosphere in the summer hemisphere also becomes warmer with a time lag of days. This IHC phenomenon is considered to be a coupling through processes in the middle atmosphere (i.e., stratosphere, mesosphere, and lower thermosphere). Several plausible mechanisms have been proposed so far, but they are still controversial. This is mainly because of the difficulty in observing and simulating gravity waves (GWs) at small scales, despite the important role they are known to play in middle atmosphere dynamics. In this project, by networking sparsely but globally distributed radars, mesospheric GWs have been simultaneously observed in seven boreal winters since 2015/16. We have succeeded in capturing five stratospheric sudden warming events and two polar vortex intensification events. This project also includes the development of a new data assimilation system to generate long‐term reanalysis data for the whole middle atmosphere, and simulations by a state‐of‐the‐art GW‐permitting general circulation model using the reanalysis data as initial values. By analyzing data from these observations, data assimilation, and model simulation, comprehensive studies to investigate the mechanism of IHC are planned. This paper provides an overview of ICSOM, but even initial results suggest that not only GWs but also large‐scale waves are important for the mechanism of the IHC.

Based on three years of data observed by Beijing MST radar between 2012 and 2014, this study gives the distribution characteristics of the turbulence energy dissipation rate ε and the vertical turbulence diffusion coefficient ${\mathrm{K}}_{\mathrm{z}}$ in the troposphere–lower stratosphere over Beijing based on the spectral width method. Additionally, from the perspective of atmospheric stability expressed by the Brunt–Väisälä frequency N and the gradient Richardson number Ri, the seasonal variations of turbulence parameters at different heights are analyzed. The results show that ε (${\mathrm{K}}_{\mathrm{z}}$) is distributed in the range of 10⁻⁴ to 10^−1.5 m² s⁻³ (10⁰ to 10¹ m² s⁻¹), with median values of 10⁻³ to 10^−2.5 m² s⁻³ (10^0.4 to 10^0.75 m² s⁻¹). In the height range of −1 km to about +1 km from the tropopause, the profiles of ε have no noticeable seasonal difference from October to next May, and ε increases sharply with the slope of linear regression ∇log(ε)/∇h = 0.25 (m² s⁻³) km⁻¹. The seasonal variation of turbulence parameters has noticeable differences at different atmospheric layers. Furthermore, the atmospheric static/dynamic stability and turbulence intensity are the influencing factors of turbulence parameters. From the tropopause to 14 km, the values of ε and N both have a peak from October to next May, and this peak layer disappears from July to August. Within the range of 2–4 km from the ground, ${\mathrm{K}}_{\mathrm{z}}$ and the frequency of Ri < 0.25 have a significant positive correlation, with a Pearson correlation coefficient greater than 0.5, and both have large value layers in March to May and September to October. The turbulence intensity above 15 km is likely the main factor affecting the turbulence dissipation rate and vertical turbulence diffusion coefficient.

The tropopause region is crucial for the stratosphere‐troposphere exchange (STE) and acts as an indicator of climate change. Double tropopauses (DTs) act to increase the STE process but their driving mechanisms remain an open question. The present assessment offers for the first time the linkage between tropospheric gravity waves (GWs) and DT events by exploring a global data set of multi‐year radiosonde measurements. In the extratropics, the occurrence frequency of DT events keeps a remarkably consistent spatial‐temporal structure with GW total energy. Under the DT scenario, GW total energy has increased by 37.67% compared to single tropopause events. Based on a statistical assessment, the upward propagating GWs throughout the second tropopause region can probably raise Kelvin‐Helmholtz instability or turbulence, leading permanent irregularities in thermodynamic structure, and consequently, increasing the likelihood of DT.

Based on the quality-controlled observational spectral width data of the Beijing mesosphere–stratosphere–troposphere (MST) radar in the altitudinal range of 3–19.8 km from 2012 to 2014, this paper analyses the relationship between the proportion of negative turbulent kinetic energy (N-TKE) and the horizontal wind speed and the vertical shear of horizontal wind domain and gives the distributional characteristics of atmospheric turbulence parameters obtained by using different calculation models. Three calculation models of the spectral width method were used in this study – namely the H model (Hocking, 1985), N-2D model (Nastrom, 1997) and D–H model (Dehghan and Hocking, 2011). The results showed that the proportion of N-TKE in the H model, N-2D model and D–H model increases with the horizontal wind speed u and/or the vertical shear of horizontal wind speed $\frac{\partial u}{\partial z}$ , and the maximum values are 60 %, 45 % and 35 %, respectively. When the$\left|\frac{\partial u}{\partial z}\right|$ is greater than 0.006 s⁻¹, the N-TKE of the H model increases sharply with $\left|\frac{\partial u}{\partial z}\right|$ ; the increase rate is about $\frac{\mathrm{20}\mathit{\%}}{\mathrm{0.002}{\mathrm{s}}^{-\mathrm{1}}}$ . For these three models, the results are similar except that the vertical shear of the horizontal wind speed is greater than 0.006 s⁻¹. When $\frac{\partial u}{\partial z}>\mathrm{0.006}$ s⁻¹, the proportion of N-TKE in the N-2D and H models increases with $\frac{\partial u}{\partial z}$ , while the proportion in the D–H model is less than 10 % and has slight variation. However, it is still necessary to consider the applicability of the N-2D model and D–H model in some weather processes with strong winds. The distributional characteristics with height of the turbulent kinetic energy dissipation rate ε and the vertical eddy diffusion coefficient K_z derived by the three models are consistent with previous studies. Still, there are differences in the values of turbulence parameters. Also, the range resolution of the radar has little effect on the differences in the range of turbulence parameters' values. The median values of ε in the H model, N-2D model and D–H model is 10^−3.2–10^−2.7, 10^−3.0–10^−2.6 and 10^−3.3–10^−2.8 m² s⁻³, respectively. The median values of K_z in these three models are 10^0.3–10^0.7, 10^0.4–10^0.7 and 10^0.1–10^0.5 m² s⁻¹.

The performance of general circulation models (GCMs) in simulating horizontal winds is important because the distribution and variation in horizontal winds are central to investigating atmospheric dynamic characteristics and processes. Also, horizontal wind data can be used to extract some of the required information on gravity waves, tides, and planetary waves. In this context, the present paper evaluates the capability of the Institute of Atmospheric Physics atmospheric general circulation model high-top version (IAP-HAGCM) in simulating the horizontal winds and tides of the troposphere and lower stratosphere by presenting a climatological and statistical comparison against observations of the powerful Beijing mesosphere–stratosphere–troposphere (MST) radar (39.78°N, 116.95°E) during 2012–2014. The results illustrated that the IAP-HAGCM can successfully reproduce the time–altitude distribution of the monthly mean zonal wind and diurnal tide amplitude, albeit with some underestimation. The mean correlation coefficients and root-mean-square error for the zonal (meridional) winds were 0.94 (0.73) and 6.60 m s⁻¹ (2.90 m s^–1), respectively. Additionally, the IAP-HAGCM can capture the temporal variation in both the zonal and meridional winds. It is worth noting that, compared with the seven coupled model intercomparison project phase 6 (CMIP6) models, the IAP-HAGCM performs better in meridional wind simulations below 15 km. However, there are discrepancies in altitudinal ranges with large wind velocities, such as the westerly jet, in the transition region of the troposphere and stratosphere, and in February, April, July, and September. It is suggested that model users should take advantage of the model’s simulation ability by combining this information regarding when and where it is optimal with their own research purposes. Moreover, the evaluation results in this paper can also serve as a reference for guiding improvements of the IAP-HAGCM.

An international joint research project, entitled Interhemispheric Coupling Study by Observations and Modelling (ICSOM), is ongoing. In the late 2000s, an interesting form of interhemispheric coupling (IHC) was discovered: when warming occurs in the winter polar stratosphere, the upper mesosphere in the summer hemisphere also becomes warmer with a time lag of days. This IHC phenomenon is considered to be a coupling through processes in the middle atmosphere (i.e., stratosphere, mesosphere, and lower thermosphere). Several plausible mechanisms have been proposed so far, but they are still controversial. This is mainly because of the difficulty in observing and simulating gravity waves (GWs) at small scales, despite the important role they are known to play in middle atmosphere dynamics. In this project, by networking sparsely but globally distributed radars, mesospheric GWs have been simultaneously observed in seven boreal winters since 2015/16. We have succeeded in capturing five stratospheric sudden warming events and two polar vortex intensification events. This project also includes the development of a new data assimilation system to generate long-term reanalysis data for the whole middle atmosphere, and simulations by a state-of-art GW-permitting general circulation model using reanalysis data as initial values. By analyzing data from these observations, data assimilation, and model simulation, comprehensive studies to investigate the mechanism of IHC are planned. This paper provides an overview of ICSOM, but even initial results suggest that not only gravity waves but also large-scale waves are important for the mechanism of the IHC.