Roles of vertical distributions of atmospheric transient eddy dynamical forcing and diabatic heating in midlatitude unstable air–sea interaction

Atmospheric transient eddy dynamical forcing (TEDF)-driven midlatitude unstable air–sea interaction has recently been recognized as a crucial positive feedback for the maintenance of the extratropical decadal variabilities. Our recent theoretical work (Chen et al., Clim Dyn https://doi.org/10.1007/s00382-020-05405-0, 2020) has characterized such an interaction through building an analytical midlatitude barotropic atmospheric model coupled to a simplified upper oceanic model. This study extends the analytical model to including a two-layer quasi-geostrophic baroclinic atmospheric model and then identifies the roles of vertical distributions of atmospheric TEDF and diabatic heating in midlatitude unstable air–sea interaction. It is found that midlatitude air–sea coupling with more realistic vertical profiles of atmospheric TEDF and diabatic heating destabilizes oceanic Rossby wave modes over the entire range of zonal wavelengths, in which the most unstable coupled mode features an equivalent barotropic atmospheric low (high) pressure over a cold (warm) oceanic surface. Spatial structure and period of the most unstable mode are more consistent with the observation than those from in previous model. Although either TEDF or diabatic heating alone can lead to a destabilized coupled mode, the former makes a dominant contribution to the instability. The increase of low-layer TEDF stimulates the instability more effectively if the TEDF in upper layer is larger than in lower layer, while the TEDF in either high or low layers can individually cause the instability. The surface heating always destabilizes the air–sea interaction, while the mid-level heating always decays the coupled mode. The results of this study further confirm the TEDF-driven positive feedback mechanism in midlatitude air–sea interaction proposed by recent observational and numerical experiment studies.


Introduction
Observational studies have revealed that there is a significant decadal variability in the midlatitude North Pacific ocean-atmosphere system (Trenberth 1990;Graham et al. 1994;Minobe 1997;Mantua et al. 1997;Enfield and Mestas-Nunez 1999;Zhu and Yang 2003;Fang et al. 2006). The decadal mode of North Pacific SST anomalies, also known as the Pacific Decadal Oscillation (PDO), is well correlated with the atmospheric circulation anomalies that exhibit an equivalent barotropic structure in the vertical direction, with geopotential lows (highs) above cold (warm) water (Kushnir et al. 2002;Namias and Cayan 1981;Fang and Yang 2016). Such an equivalent barotropic cold/trough (warm/ ridge) structure is the typical feature of midlatitude climate variabilities on decadal time scale (Cayan 1992;Deser and Blackman 1993).
Usually, to generate an oscillation, a positive feedback mechanism is needed to make the initial perturbation to develop, and a delayed negative feedback mechanism is also needed for the phase transition. For PDO, the negative feedback mechanisms are attributed to be those slow upperocean adjustment processes involving the oceanic gyre (Latif and Barnett 1994;Fang et al. 2006), the subduction (Gu and Philander 1997), and the oceanic Rossby wave propagation (Qiu et al. 2007), in which the decadal time scale of PDO is determined by the adjustment period, while the positive feedback mechanism is considered to be the local large-scale ocean-atmosphere interaction in the middle latitudes (Latif andBarnett 1994, 1996;Latif 1999;Miller and Schneider 2000;Robertson 1996;Zorita and Frankignoul 1997). Though a number of observations, theoretical analyses and GCM simulations have provided evidence for the active feedback of extratropical SST on large-scale atmospheric circulation McWilliams 1997, 1998;Neelin and Weng 1999;Latif and Barnett 1996;Liu and Wu 2004;Zhong et al. 2008;Zhu et al. 2008;Fang and Yang 2011), the mechanism responsible for the midlatitude unstable air-sea interaction remains unclear, which is primarily because the processes by which the midlatitude sea surface temperature (SST) anomaly affects the seasonal-mean atmospheric anomaly have not been fully understood.
Different from the tropical atmospheric circulation that is mainly driven by the SST-related convective diabatic heating, the midlatitude atmospheric circulation is both thermally-and eddy-driven. Since the atmospheric stratification in the middle latitudes is stable, the SST-induced diabatic heating is relatively weak and mainly confined to the lower troposphere. On the other hand, the midlatitude atmosphere is strongly baroclinic, and the synoptic transient eddies are active, particularly over the midlatitude oceanic frontal zones, featuring storm tracks (Ren et al. 2010;Chu et al. 2013;Liu et al. 2014;Nakamura et al. 2004;Small et al. 2008;Wang et al. 2017). The transient eddies can systematically transport heat and momentum, driving and maintaining mean atmospheric circulations in the middle latitudes via the nonlinear eddy-to-mean flow feedback (Ren et al. 2010;Xiang and Yang 2012;Zhang et al. 2012;Nie et al. 2013Nie et al. , 2014. From the perspective of quasi-geostrophic potential vorticity (QGPV) dynamics, the seasonal-mean potential vorticity (PV) sources include the diabatic heating, the transient eddy heating, and the transient eddy vorticity forcing (Fang and Yang 2016). The latter two items caused by the convergence of the heat and vorticity transport by transient eddies, respectively, can be considered as indirect forcing sources of atmospheric circulation. Therefore, the midlatitude SST anomalies can affect the atmosphere by changing both the diabatic heating and the transient eddy forcing. Fang and Yang (2016) identified the features of those PV sources associated with PDO and quantitatively diagnosed their effects on the atmospheric anomalies. Based on their analyses, a positive feedback mechanism for midlatitude unstable air-sea interaction in the North Pacific was hypothesized as follows. An initial midlatitude surface westerly anomaly accompanied by intensified Aleutian low tends to force a negative SST anomaly by increasing upward surface heat fluxes and driving southward Ekman current anomaly. Then, the basin-scale SST cooling tends to increase the meridional SST gradient in its southern flank, thus enhancing the subtropical oceanic front. In adjustment to the enhanced oceanic front, the low-level atmospheric baroclinicity tends to be strengthened, and more transient eddies are generated, resulting in enhanced meridional thermal and momentum transportations. Although all the diabatic heating, transient eddy heating and transient eddy vorticity forcing are increased during these processes, only the atmospheric response to the transient eddy vorticity forcing tends to be barotropic. The vorticity forcing that dominates the total atmospheric forcing tends to produce an equivalent barotropic atmospheric low pressure that intensifies the initial anomalies of the midlatitude surface westerly and Aleutian low. Therefore, the midlatitude air-sea interaction, in which the oceanic front and the atmospheric transient eddy are the indispensable ingredients, can provide a positive feedback mechanism for the development and maintenance of the observed decadal anomalies in the midlatitude North Pacific ocean-atmosphere system. This hypothesis has been confirmed by the later observational and modeling studies (Wang et al. 2017(Wang et al. , 2019Tao et al. 2020;Zhang et al. 2020).
Since atmospheric transient eddy dynamical forcingdriven midlatitude unstable air-sea interaction is recognized as a crucial positive feedback for the maintenance of the extratropical decadal variabilities, our recent theoretical work (Chen et al. 2020, hereinafter CFY2020) has characterized such an interaction through building an analytical midlatitude barotropic atmospheric model coupled to a simplified upper oceanic model. In the analytical coupled model, the atmospheric component is governed by the barotropic QGPV equation with transient eddy vorticity forcing which is parameterized to be linearly proportional to the second-order meridional derivative of SST anomaly in terms of observational analyses. It is found that the oceanic Rossby wave stimulated by air-sea coupling can be unstable within a wide range of wavelengths. The configuration structure of the most unstable mode is similar to that of the observed decadal variability. Although both surface diabatic heating and transient eddy vorticity forcing can cause unstable mode, the latter plays a dominant role. The results provide a theoretical support for the eddy-driven midlatitude unstable air-sea interaction hypothesis.
Note that a simplified assumption of the barotropic atmosphere is used in CFY2020. In reality, both diabatic heating and transient eddy vorticity forcing for the atmosphere have vertical variations. Studies using linear baroclinic model show that the atmospheric response is sensitive to the vertical structure of the idealized heat source. In response to the low-level heating, the atmosphere tends to produce a barotropic low-pressure response, while to the mid-level heat source, the atmospheric response tends to be baroclinic (Qiu et al. 2014;Hoskins and Karoly 1981;Fang and Yang 2011). The influence of vertical distribution of the transient eddy vorticity forcing on the atmospheric response remains to be clarified.
In the present study, the atmospheric component of the analytical midlatitude coupled air-sea model in CFY2020 is extended to a two-layer baroclinic quasi-geostrophic model, and the oceanic component remains to be a 1.5-layer upper ocean model. The vertical structures of diabatic heating and transient eddy vorticity forcing are both taken into account to further explore their relative contributions to midlatitude unstable air-sea interaction. The rest of the paper is organized as follows. Description of the analytical coupled model is presented in Sect. 2. The dynamical features of midlatitude air-sea coupled mode and the relative contributions of diabatic heating and transient eddy vorticity forcing with different vertical profiles are examined in Sects. 3 and 4, respectively. A comparison of the results from current model with those from previous model is presented in Sect. 5. Influences of other processes and factors such as the oceanic adjustment processes, the air-sea coupling intensity and the background atmospheric flow on the unstable coupled mode are discussed in Sect. 6. The final section is devoted to conclusions and discussion.

An extended analytical midlatitude
coupled air-sea model Figure 1 shows a schematic diagram of the extended analytical coupled air-sea model used in this study. The atmospheric component is a two-layer baroclinic quasi-geostrophic model with a shallow boundary layer. The oceanic component remains to be a quasi-geostrophic upper-ocean model, as in CFY2020, which consists of a constant-depth mixed layer, a thin entrainment layer, and a thermocline layer. Both the atmosphere and the ocean are governed by their respective quasi-geostrophic potential vorticity (QGPV) equations on a midlatitude beta plane. For the atmosphere, the diabatic heating can occur both at the top of boundary layer (surface heating, Q s ) and between the two layers (mid-level heating, Q ). The transient eddy vorticity forcing ( F eddy ) exists in both layers. Since the adjustment time scale of the atmosphere is much shorter than that of the ocean, the atmospheric fluctuations are assumed to be a steady-state response to the forcing sources. Under these considerations, seasonal mean QGPV equations for the two atmospheric layers can be written as, where the overbar denotes seasonal mean, q is the potential vorticity defined by x the Jacobian operator, L a = (N a H a )∕(2f ) the atmospheric baroclinic Rossby radius of deformation with N a the atmospheric Brunt-Väisälä buoyancy frequency and H a the depth of the free troposphere. In this model, seasonal mean atmospheric anomalies are driven by PV sources due to diabatic heating [the second terms on the right hand of Eqs. (1) and (2)] and due to transient eddy vorticity forcing ( F eddy ). The PV sources due to diabatic heating is proportional to the vertical gradient of diabatic heating, while F eddy = −∇ ⋅ ⃗ V � h � , representing the seasonal mean convergence of vorticity transport by synoptic transient eddies (denoted by primes). r∇ 2 represents the damping effect that is proportional to the relative vorticity, and r is the diffusive coefficient with a characteristic time scale taken as r −1 = 5 days (Pedlosky 1970).
The atmospheric basic state is assumed to be a horizontally uniform zonal flow with vertical shear. Then, seasonal mean atmospheric streamfunction can be expressed as a sum of a part determined by the basic flow (climatological mean) and a perturbation (seasonal anomaly) like 1 = 1 − U 1 y , 2 = 2 − U 2 y . Further, atmospheric perturbations can be decomposed into barotropic and baroclinic components defined as Then the barotropic and baroclinic PV perturbation equations can be linearized about the basic state, respectively, as.
where ̂ = 2 , and ∼ = 2Ũ∕L 2 a . Accordingly, the right-hand terms in Eqs. (7) and (8) represent their seasonal anomalies. It should be noted that the surface diabatic heating and the atmospheric transient eddy vorticity forcing in each layer drive the barotropic and baroclinic components together, whereas the mid-level heating only acts on the baroclinic component.
As assumed in CFY2020, the anomalous transient eddy vorticity forcing is parameterized to be proportional to the second-order meridional derivative of SST anomaly. Thus here the forcing in the upper layer is specified to be equal to that in the barotropic model (CFY2020), i.e., F eddy1 = 2 T 1 ∕ y 2 , with T 1 representing the sea surface temperature anomaly. Meanwhile, following Shutts (1987) and Fang and Yang (2011), the mid-level heating is prescribed to have the form of Newtonian relaxation of the atmospheric potential temperature perturbation at middle level a to a radiative-convective equilibrium temperature anomaly * a , given by where a can be converted into ∼ using the thermal wind relation, ao is a typical value of a , a is the atmospheric air-sea heat exchange parameter, and * a is set to be the function of SST anomaly. For simplicity, it is assumed that * a = T 1 .
Based on the observational analysis on vertical distributions of the anomalous diabatic heating and transient eddy vorticity forcing over the North Pacific region ( Fig. 2b and c, quoted from CFY2020), the surface heating is estimated to be twice as large as that of the mid-level heating ( Q s = 1 Q , 1 = 2 ), and the transient eddy vorticity forcing in the lower atmosphere is estimated to be 0.2 times of that in the upper atmosphere ( F eddy2 = 2 F eddy1 , 2 = 0.2 ). Then, Eqs. (7) and (8) can be rewritten as The oceanic component of the coupled model keeps the same as in CFY2020 and Fang and Yang (2011). The upperlayer motion represented by the first baroclinic Rossby mode is controlled by a linear QGPV equation with a rest basic state driven by the curl of the surface wind stress. The evolution of SST is determined by physical processes including horizontal advection, vertical entrainment and heat exchange at air-sea interface. The governing equations for the ocean are given by where o is the upper-layer oceanic geostrophic streamfunction, L o the oceanic baroclinic Rossby radius of deformation. c 1 , c 2 , c 3 and c 4 are all positive coefficients relevant to advection and entrainment processes, and o the oceanic air-sea heat exchange parameter. Equations (10)-(13) constitute the final closed equations for the coupled air-sea model. Compared with the coupled model used in Fang and Yang (2011), the role of the atmospheric transient eddy dynamical feedback is introduced in the present model. Also, the coupled model with barotropic atmosphere used in CFY2020 can be regarded as a special case of the present model.

Dynamical features of the midlatitude unstable coupled modes
Atmospheric and oceanic perturbations are assumed to be limited within a channel of [0, π ] in y direction and satisfy the standard plane wave solution in x direction with the form where is the complex frequency, k and l are the horizontal wavenumbers in x and y directions, respectively. Substituting (14) into Eqs. (10)-(13), we obtain the following dispersion relations of the coupled modes  Oceanic specific heat at constant pressure Mean zonal wind for atmospheric layer 1 17m s −1 U 2 Mean zonal wind for atmospheric layer 2 8m s −1 c 4 Mean oceanic meridional temperature gradient  The two solutions of (15) correspond to two coupled modes caused by midaltitude air-sea interaction. The real part of represents the frequency of the coupled mode and the imaginary part the growth rate of the amplitude of coupled mode. The coupled mode is unstable if the imaginary part of is positive. For calculating , values of all parameters are taken to be the same as in CFY2020 (see Table 1), except for the basic zonal wind speeds which are set to be U 1 = 17m s −1 and U 2 = 8m s −1 for upper and lower atmosphere, respectively. The frequency and growth rate of the two coupled modes can be calculated as a function of zonal wavelength. The frequency curve of the first mode (the blue line in Fig. 3a) is very close to that of the free oceanic Rossby wave (the red line), indicating that the first coupled mode ( + ) is essentially a coupled oceanic Rossby wave mode modified by air-sea interaction. The coupled oceanic Rossby wave mode propagates westward and is unstable over the whole range of zonal wavelengths (Fig. 3b). The other coupled mode ( − ) corresponds to an eastward-propagating decaying mode since its growth rate is always negative within the entire range of zonal wavelengths (Figure not shown). Because the focus in this study is mainly on the unstable air-sea interaction, the damping mode will no longer be discussed below.
The most unstable coupled Rossby wave mode occurs at the wavelength around 9900 km (Fig. 3b), with a frequency of 0.78 per year (Fig. 3a) corresponding to a period at about 8 years, which is close to that of the observed decadal variability (~ 10 years, Qiu et al. 2007). Figure 4 shows the spatial phase relationships between the SST anomaly (shading) and the upper-and lower-layer atmospheric streamfunction anomalies as well as the upper-oceanic streamfunction anomalies (contours) for the most unstable coupled mode. The upper-and lower-layer atmospheric streamfunction anomalies are both in phase with the SST anomaly, indicating that the atmospheric response to the SST anomaly is equivalent barotropic with highs (lows) over warm (cold) water ( Fig. 4a and b). In addition, the wind-driven upper oceanic streamfunction anomaly has a slight eastward phase shift relative to SST and atmospheric streamfunction anomalies, indicating the importance of the meridional advection process in determining the SST variation. Such a phase relationship between atmospheric and oceanic anomalies is quite consistent with that of the observed decadal variability ( Fig. 5c and d). These results suggest that midlatitude air-sea interaction through atmospheric thermal and transient eddy dynamical feedback can produce an unstable coupled oceanic Rossby wave mode on decadal time scale with equivalent barotropic cold trough/warm ridge structure.
Note that, according to the form of plane wave solutions used above in this study, all the variables are required to be in phase in y direction. Thus the meridional phase shift as shown in Fig. 5c and d cannot be obtained. If the form of plane wave solutions is assumed as ∝ e i(kx+ly− t) , the

Role of diabatic heating with different vertical profiles
In this coupled model, diabatic heating and atmospheric transient eddy vorticity forcing are the two PV sources that jointly affect the seasonal-mean atmospheric anomaly. To examine the relative contributions of different PV sources in destabilizing the coupled oceanic Rossby wave mode, the transient eddy vorticity forcing term is firstly removed by setting = 0 in Eqs. (10) and (11), and only diabatic heating is kept. Three heating profiles are investigated, respectively: surface heating only ( Q s ≠ 0, Q = 0 ), mid-level heating only ( Q s = 0, Q ≠ 0 ), and both surface and mid-level heating ( Q s , Q ≠ 0). The dispersion relation curves (frequencies) of all the three cases are almost the same as the free oceanic Rossby wave, shown in Fig. 6a, c and e. This indicates that the frequency of the oceanic Rossby wave has little been affected by the air-sea coupling. However, the instability appears to be sensitive to the heating profiles. When the surface heating is only considered, the coupled oceanic Rossby mode can destabilize over the entire range of wavelengths ( Fig. 6b) with growth rate that has a significant reduction compared with the standard full PV sources case (Fig. 3b), and the most unstable mode shifts to a shorter wavelength at around 8500 km. When the mid-level heating is only considered, the coupled oceanic Rossby mode becomes a weak damping  Fig. 4, but for the surface diabatic heating only (left panels) and the mid-level diabatic heating only (right panels) mode. By comparing the phase relationships of atmospheric and oceanic variables for the two cases (Fig. 7), we can find that the atmospheric anomalies for the two cases are nearly out of phase, relative to the same negative SST anomalies. The atmospheric response to the surface cooling induced by cold SST anomaly shows an equivalent barotropic low above cold water anomaly (Fig. 7a and c), similar to the situation in the standard case but with an eastward phase shift relative to SST anomalies. The surface low pressure further sustains the cyclonic oceanic flow and enhances the cold SST anomaly by transporting cold water (Fig. 7e), thus forming a positive air-sea feedback, indicating that the coupled Rossby wave mode is a growing mode. However, in the case of the negative SST anomaly-induced mid-level cooling, the upper-layer atmospheric streamfunction exhibits a high anomaly response just over the negative SST anomaly (Fig. 7b), while the high anomaly of the lower-layer atmospheric streamfunction leads that of the upper layer by nearly ∕2 phase (Fig. 7d). The surface high pressure tends to weaken the cyclonic oceanic flow and thus reduce the cold SST anomaly. In this case, the coupled Rossby wave mode decays gradually.
When the surface and mid-level heating are both included, the results are very similar to that of the surface heating only case (Fig. 6f), indicating that the contribution of diabatic heating to the unstable air-sea coupled mode is dominated by the surface heating. Mid-level heating that tends to cause baroclinic atmospheric response does not cause any unstable coupled mode while surface heating always generates the instability. Figure 10a shows the variation of unstable growth rate with zonal wavelength and surface heating intensity. Here the mid-level heating is fixed and the surface heating is modified continuously by changing 1 . It can be seen that, corresponding to the fastest-growing mode, the stronger the surface heating is (the larger 1 is), the more unstable the coupled mode emerges.

Role of transient eddy vorticity forcing with different vertical profiles
By setting a = 0 in Eqs. (10) and (11), the diabatic heating is neglected. Atmospheric transient eddy vorticity forcing becomes the only PV source of the atmosphere. Also, three vertical profiles of the transient eddy vorticity forcing are discussed: the forcing occurring only in the lower troposphere ( F eddy2 ≠ 0, F eddy1 = 0 ), only in the upper troposphere ( F eddy2 = 0, F eddy1 ≠ 0 ), and in both layers ( F eddy2 , F eddy1 ≠ 0 ). In all the three cases, the oceanic Rossby wave mode is unstable over the entire range of wavelengths (Fig. 8). However, in the lower-layer forcing-only case, the instability of the coupled mode is very weak. Growth rate of the most unstable mode is only 0.02 per year, and corresponding wavelength and period are 8800 km and 7.67 years, respectively ( Fig. 8b and e). In the upper-layer forcing-only case, growth rate of the most unstable mode is greatly increased, corresponding to a longer wavelength of 9900 km and a longer period of 8.17 years ( Fig. 8a and d). The spatial phase relationship between the atmospheric streamfunction and SST anomalies ( Fig. 9a and c) is much closer to that in the standard case ( Fig. 3a and b). When the transient eddy vorticity forcings  Fig. 6, but for the transient eddy vorticity forcing in the upper layer only (left panels), lower layer only (mid panels), and both layers (right panels), respectively. Note that the green lines in the lower panels indicate the associated growth rate when the transient eddy vorticity forcing in the lower layer is increased to be equal to that in the upper layer at both layers are considered, growth rate is similar to that in the upper-layer forcing-only case, but the most unstable growth rate is increased slightly with a zonal wavelength of around 10,200 km and a period of 8.34 years ( Fig. 8c and f). From the above results, it can be concluded that the upperlayer transient eddy vorticity forcing plays a dominant role in producing unstable air-sea interaction mode.
It should be mentioned that, in the present study, the lower-layer transient eddy vorticity forcing is set to be 0.2 times of that of upper layer ( 2 = 0.2 ) according to the observational analysis. If we assume the magnitudes of the forcing at both layers are equal ( 2 = 1 ), which is equivalent to the setting in the barotropic atmospheric model (CFY2020), then the unstable growth rate is enhanced greatly no matter whether the upper-layer forcing is considered or not (as shown by the green lines in Fig. 8e and f). Under this setting, the contribution of individual lower-layer vorticity forcing on the unstable mode has overtaken that of individual upper-layer forcing ( Fig. 8d and e), suggesting that the lowerlayer vorticity forcing causes instability more effectively. This is also illustrated by Fig. 10b, in which the upper-layer forcing is fixed and the lower-layer forcing is changed continuously. With the increase of low-layer forcing, the unstable growth rate increases, and the fastest-growing mode moves to longer wavelengths slightly. This result can be explained by Eqs. (7) and (8). According to the relationship between transient eddy vorticity forcing and SST anomalies, a cold SST anomaly causes positive vorticity forcing in both layers. The lower-layer forcing F eddy2 makes the same contribution to the barotropic component of the atmosphere as the upperlayer forcing F eddy1 does. But for the baroclinic component, F eddy2 tends to decrease the baroclinic component by offsetting the contribution of F eddy1 , indicating that the lower-layer atmospheric response of geopotential low anomaly becomes stronger which strengthens the cold SST anomaly more effectively through driving cyclonic oceanic flow. Although either diabatic heating or transient eddy vorticity forcing can lead to a destabilized oceanic Rossby wave mode with similar equivalent barotropic structure of atmospheric anomaly, the latter has a larger contribution to the unstable growth rate (Figs. 3b, 6f, and 8f). This demonstrates again that atmospheric transient eddy dynamical feedback plays a crucial role in unstable midlatitude air-sea interaction, as pointed out by previous studies (Fang and Yang 2016;CFY2020).

Comparison with the results of coupled model with barotropic atmosphere
In CFY2020, since the atmosphere is assumed to be barotropic, only the surface diabatic heating that occurs at the top of boundary layer and the transient eddy vorticity forcing with no vertical variation can have effect. Nevertheless, as expressed by Eqs. (1) and (2)   Compared with the results of CFY2020 in which the atmosphere is assumed to be barotropic, the period of the most unstable mode in this study becomes longer, and the associated SST and atmospheric streamfunction anomalies are exactly in phase (Fig. 4), which is more consistent with the observation (Fig. 5c and d), instead of a slight phase shift in x direction in the barotropic atmospheric model ( Fig. 5a  and b). Another striking difference is on the growth rate for the most unstable mode which is largely reduced (Fig. 11a) in the present model when the vertical profiles of diabatic heating and transient eddy vorticity forcing are taken to be close to the observed. One reason is due to the offset effect by the mid-level diabatic heating (Fig. 11b), and other reason is due to the smaller lower-layer transient eddy vorticity forcing (Fig. 11c). As mentioned in Sect. 3, because the vertical structures of atmospheric PV sources and basic flow are more realistic in the present baroclinic model, the associated period and spatial configuration patterns of the atmospheric and oceanic anomalies for the unstable coupled mode coincide with the observation much better.

Relative roles of oceanic entrainment and advection processes
As expressed in Eq. (12), the oceanic processes determining SST evolution include horizontal advection, vertical entrainment and air-sea flux exchange, which represent the oceanic adjustment to the atmospheric forcing. Relative roles of different oceanic processes in destabilizing coupled oceanic Rossby wave mode are examined below. Since air-sea flux exchange just stimulates very weak damping mode (figure not shown), main processes emphasized here are entrainment and advection. Figure 12 shows the growth rate as a function of wavelength for the coupled oceanic Rossby wave mode for three cases: full processes, entrainment-only, and advection-only, respectively. It can be found that both entrainment and advection have contributions to the unstable coupled mode. When the wavelength is shorter than 10,000 km, total unstable growth rate is contributed primarily from advection, while within the wavelength range beyond 10,000 km, it is contributed jointly from advection and entrainment.

Influence of coupling strength and basic zonal flows
The unstable coupled mode may be influenced by two kinds of critical parameters, coupling coefficients ( a and ) determining the coupling strength between ocean and atmosphere through thermal and dynamical processes, and background zonal wind speeds in the two atmospheric layers ( U 1 and U 2 ). In the standard case, values of these parameters are set as: = 11.2s −2 m 2 K −1 , a = 8 × 10 −7 s −1 , U 1 = 17m s −1 , and U 2 = 8m s −1 . The influence of air-sea coupling strength and basic zonal flows on unstable air-sea interaction is estimated by changing those parameters. Figure 13 shows the variation of growth rate for the most unstable Rossby wave mode with zonal wavelength versus dynamical coupling coefficient ( ), thermal coupling coefficient ( a ), upper-layer basic zonal wind speed ( U 1 ), lower-layer basic zonal wind speed ( U 2 ), and vertical wind shear ( Ũ = U 1 − U 2 ), respectively. The growth rate for the standard case is marked by " × " in the figure. It can be seen from Fig. 13a and b that the coupled oceanic Rossby wave is always unstable within the entire range of air-sea coupling parameters. With an increase of coupling strength, no matter dynamical coupling or thermal coupling, the instability of the coupled mode is enhanced. However, compared with the thermal coupling (Fig. 13b), the influence of transient eddy dynamical feedback intensity is much more prominent, because the amplitude of growth rate is much larger and the corresponding unstable wavelength range is much Fig. 12 Growth rate as a function of wavelength for the coupled oceanic Rossby wave mode with SST evolution determined by a full oceanic processes, b entrainment process-only, and c advection process-only. Note that a is the same as Fig. 3b wider (Fig. 13a). When the air-sea coupling is strengthened, the SST anomaly can force the atmosphere more strongly, especially through transient eddy vorticity forcing, thus the coupled mode becomes more unstable.
The intensity and the vertical shear of basic zonal wind are also crucial factors determining the unstable growth. When basic zonal wind speed in the upper layer is larger than that in the lower layer ( U 1 > 8m s −1 ), the unstable growth rate decreases with increase of U 1 (Fig. 13c). Conversely, when basic zonal wind speed in the lower layer is less than that in the upper layer ( U 2 < 17m s −1 ), the unstable growth rate increases with increase of U 2 (Fig. 13d). In other words, the larger the vertical wind shear, the weaker the instability (Fig. 13e). When basic zonal wind speeds in both layers are equal, as in the barotropic model, the unstable growth rate reaches its maximum (Fig. 13e). That is because the vertical wind shear tends to weaken the barotropic atmospheric response by vorticity advection and thus is not conducive to the unstable coupled mode.

Conclusions and discussion
Midlatitude air-sea interaction is one of the crucial sources for decadal climate variabilities. However, mechanism responsible for the midlatitude air-sea interaction remains unclear primarily due to a lack of deep understanding on how the midlatitude SST can affect the atmosphere. Fig. 13 As in Fig. 10, but for a the dynamical coupling coefficient ( ), b thermal coupling coefficient ( a ), c the upperlayer background zonal wind speed ( U 1 ), d the lower-layer background zonal wind speed ( U 2 ), and e the vertical wind shear ( U 1 − U 2 ), respectively Generally, it is believed that the midlatitude SST may affect the atmosphere through diabatic heating, as in the tropics. But such a thermally-driving mechanism fails to explain observed equivalent barotropic structure of decadal atmospheric anomalies.
The atmosphere features strong baroclinicity in the midlatitudes where synoptic eddies accompanied with jet stream develop vigorously, and heat and momentum transports by those transient eddies can in turn feedback onto the midlatitude mean flow. The eddy-to-mean flow feedback has been found to be a key approach for midlatitude SST to affect the atmosphere by a number of recent studies (Wang et al. 2017(Wang et al. , 2019Tao et al. 2020). Atmospheric transient eddy dynamical forcing (TEDF)-driven midlatitude unstable air-sea interaction has recently been recognized as a crucial positive feedback for the maintenance of the extratropical decadal variabilities (Fang and Yang 2016;Chen et al. 2020). Our recent theoretical work by Chen et al. (2020) has characterized such an interaction through building an analytical midlatitude barotropic atmospheric model coupled to a simplified upper oceanic model.
As a theoretical development, this study extends the analytical coupled air-sea model of Chen et al. (2020) to including a two-layer quasi-geostrophic baroclinic atmospheric model and then identifies the roles of vertical distributions of atmospheric TEDF and diabatic heating in midlatitude unstable air-sea interaction. In the extended coupled model, both the atmospheric and oceanic perturbations satisfy quasi-geostrophic dynamics and are governed by their respective linearized QGPV equations on a midlatitude beta plane. A linearized mixed layer temperature equation is considered to describe the SST evolution. The ocean is driven by the anomalous wind stress of the atmosphere and in turn forces the atmosphere through diabatic heating and atmospheric transient eddy vorticity forcing. The diabatic heating is specified at the top of atmospheric boundary layer and at the mid-level of two-layer atmosphere, which is linearly dependent of SST anomaly. The transient eddy vorticity forcing in each layer is parameterized to be proportional to the second-order meridional derivative of SST anomaly. Based on observation, the surface heating is set to be twice of the mid-level heating and the transient eddy vorticity forcing in the lower layer is 0.2 times of that in the upper layer for the standard case.
Analytical solutions of the model show that the midlatitude air-sea coupling with more realistic vertical profiles of atmospheric TEDF and diabatic heating destabilizes the oceanic Rossby wave mode that propagates westward over the entire range of zonal wavelengths, in which the most unstable or the fastest-growing coupled mode occurs at the wavelength of around 9900 km, with a period of nearly 8 years. For this mode, atmospheric streamfunction anomalies in the two layers are both in phase with SST anomaly, with atmospheric lows (highs) over cold (warm) water. Correspondingly, the oceanic streamfuction anomalies also exhibit a cyclonic (anticyclonic) anomaly which is slightly shifted to east of the cold (warm) SST anomaly. Thus the unstable coupled mode exhibits an equivalent barotropic cold trough (warm ridge) structure, quite similar to the observed feature of decadal climate variability. Figure 14 illustrates a schematic diagram of the phase relationships among atmospheric and oceanic variables for the most unstable coupled oceanic Rossby wave mode in the midlatitudes. If there is an initial cold (warm) SST anomaly in the north (south), it induces a strong surface cooling (warming) and a weak mid-level cooling (warming), leading to a positive (negative) PV source in both atmospheric layers. Meanwhile, the above SST anomaly pattern tends to give rise to a larger meridional SST gradient and thus a stronger low-level atmospheric baroclinicity. Atmospheric synoptic eddies become more active, causing positive (negative) transient eddy vorticity forcing (also, PV source) in the north (south) for both layers. Therefore, through both diabatic heating and transient eddy dynamical forcing, the initial cold (warm) SST anomaly generates an in-phase equivalent barotropic atmospheric low (high) anomaly over the SST anomaly. The surface low (high) pressure anomaly tends to drive the upper oceanic Rossby wave and generate a cyclonic (anticyclonic) upper oceanic flow anomaly with a slight eastward phase shift relative to the SST anomaly, which thus further amplify the initial cold (warm) SST anomaly mainly through oceanic advection process. Consequently, the amplitude of the oceanic Rossby wave grows due to the air-sea coupling.
Relative roles of diabatic heating and transient eddy vorticity forcing with different vertical profiles in destabilizing the coupled oceanic Rossby wave are analyzed. A SST anomaly can induce in-phase surface and mid-level diabatic heating as well as an opposite-phase transient eddy vorticity forcing in both layers of the atmosphere. The surface heating favors the barotropic atmospheric response, while the mid-level heating always excites the baroclinic atmospheric response. In the observation, the surface heating is larger than the mid-level heating, the contribution of diabatic heating on the unstable coupled mode is dominated by the surface heating. On the other hand, the transient eddy vorticity forcing at each layer can force barotropic atmospheric response and induce unstable coupled mode, while the lower-layer forcing can affect the unstable mode more effectively by weakening the baroclinic response. However, in the standard case, since the transient eddy vorticity forcing at the lower layer is much smaller than that in the upper layer, the upper-layer transient eddy vorticity forcing makes more prominent contribution to the unstable growth rate. Although either diabatic heat forcing or transient eddy vorticity forcing can generate unstable coupled oceanic Rossby wave alone, the transient eddy dynamical feedback plays a much more important role.
Besides the SST-induced forcing on the atmosphere, the oceanic adjustment processes driven by the atmospheric windstress anomaly, which mainly include advection and entrainment, can also affect the unstable coupled mode. When SST is adjusted by the advection process alone, the coupled oceanic Rossby wave with zonal wavelength shorter than 10,000 km tends to be unstable. When the entrainment process is considered alone, the unstable coupled mode occurs within a wider range of wavelengths (> 10,000 km). Moreover, the influences of air-sea coupling strength and background zonal wind speed on the unstable coupled mode are also examined. When air-sea coupling (either dynamical coupling or thermal coupling) is strengthened, the instability of the coupled mode is enhanced. But the effect of dynamical coupling through transient eddy vorticity forcing is more significant. For the basic zonal wind speed, its barotropic part favors the growth of coupled mode, whereas its baroclinic part determined by the vertical wind shear tends to decay the coupled mode.
Compared with the results derived from CFY2020 in which the atmospheric component is barotropic, the spatial configuration structure and period of the unstable coupled mode in this study are more consistent with the observation, although the growth rate of the unstable coupled mode becomes smaller due to the mid-level atmospheric heating and the TEDF with vertical varation which both tend to contribute to baroclinic atmospheric response and check the instability.
It should also be mentioned that there are three PV sources in the atmospheric QGPV equation due to diabatic heating, transient eddy heating and transient eddy vorticity forcing, respectively (Fang and Yang 2016). The PV source due to transient eddy heating is proportional to the vertical gradient of transient eddy heating ( Q eddy ) that represents the convergence/divergence of heat flux transported by transient eddies. Transient eddy heating that has the same dimension as diabatic heating can be seen as an indirect heat source of the atmosphere, reflecting the transient eddy thermal forcing on the atmosphere. The magnitude of transient eddy heating is comparable with that of diabatic heating, and it is mainly confined to the middle of the troposphere (Fig. 2c in CFY2020). Hence in the barotropic model, it vanishes after vertical integration. In the present study, it is not explicitly expressed in the coupled model. However, since it is equivalent to a mid-level heating, the increased mid-level heating or decreased surface heating (i.e., decreased 1 in Fig. 10a) in this study can be equivalent to the case that the transient eddy heating is considered in the coupled model. At this time, the growth rate of the most unstable coupled mode is decreased a little due to the enhancement of baroclinic atmospheric response induced by transient eddy heating.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.