Detection and attribution of upper-tropospheric warming over the tropical western Pacific

Both model and satellite observations suggest that the warming in the tropical upper troposphere is more pronounced than in the lower troposphere. This study describes the tropical upper-tropospheric-warming (UTW) in the ECMWF Reanalysis Interim dataset, and reveals radiative and non-radiative processes contributing to the long-term change in UTW. The maximum intensity of the UTW is found between 250 and 400 hPa with a warming rate about 3K/decade over the tropical Atlantic and western Pacific oceans. The warming is more prominent from boreal winter to early spring than the rest of the year. In the tropical western Pacific, where the warmest ocean water resides, the long-term change in the annual mean UTW is mainly contributed by the variations of atmospheric deep convection and the oceanic dynamical process and heat storage (Ocn) with the latter playing a more important role. The UTW caused by the deep convection is mainly through the latent heat release associated with enhanced upward motion and condensation. The UTW related to Ocn is mainly a result of enhanced upward long-wave radiation from the ocean surface due to the large heat storage in the thermocline. The variations in Ocn mainly contribute to the UTW from winter to early spring, while the changes in deep convection play a more essential role from summer to fall. Also shown by the differences in deep convection, a strengthening of the Walker circulation since the early 2000s is also found in the reanalysis, which is most apparent in boreal summer and fall over the tropical Pacific, facilitating the deep convection to dominate.


Introduction
In the context of global warming, a robust feature in the tropics is the larger upper-tropospheric warming (UTW) than the lower-tropospheric and surface warming, which was predicted by Hansen et al. (2002) and the IPCC (Intergovernmental Panel on Climate Change) TAR (Third Assessment Report) general circulation models (GCMs). The IPCC AR4 (Fourth Assessment Report) more explicitly predicted that such UTW reached its maximum at about 200 hPa and decreased to zero near the tropical tropopause, although the feature might have been overestimated (Song et al. 2014). Moreover, the IPCC AR5 (Fifth Assessment Report) also reviewed this tropical UTW in Chap. 12 by Collins et al. (2013), showing a similar pattern to the former prediction with various magnitudes under different scenarios. Aside from the IPCC GCM simulation, previous studies also conducted observational verifications. For example, Fu and Johanson (2005) applied satellite-borne Microwave Sounding Unit observations of 1987-2003 to demonstrate that the 1 3 tropical troposphere warmed faster than the surface and that the trends of tropical tropospheric temperature increased with height, which confirms the GCM prediction from IPCC TAR. Fu et al. (2011) further proved the prediction from IPCC AR4 and by analyzing the same data set as in Fu and Johanson (2005) from 1987 to 2010.
The tropical UTW is associated with the climate sensitivity due to its impact on climate feedback processes as water vapor, lapse rate, cloud feedback, and atmospheric dynamics (Colman 2001;Hartmann and Larson 2002;Bony et al. 2006). A tropospheric warming that is larger at the upper levels than at the low levels has a directly negative effect on the lapse rate, which reduces the lapse rate leading to a negative radiative feedback and further affects the water vapor feedback (Bony et al. 2006;Soden and Held 2006). Spatial and temporal variations of the globally zonal mean UTW have been extensively studied using observational evidences and numerical experiment results, such as in Thorne et al. (2007); Santer et al. (2008); Christy et al. (2007Christy et al. ( , 2010, and others. While the UTW has been gradually proved in observations and predicted by models, large inter-model spreads of uncertainty still exist, due mainly to the currently poor understanding of climate feedback processes (Colman 2001;Fu and Johanson 2005). Before discussing the impact of such UTW on climate feedback processes, the reverse impact (i.e., the formation of the UTW via the impact of climate feedback on it) deserves detailed analyses. Song et al. (2014) applied a coupled atmosphere-surface "climate feedback-response analysis method" (CFRAM) to quantify the contributions of climate feedbacks to global tropospheric warming in the NCAR CCSM3.0 (Community Climate System Model Version 3.0 of National Center for Atmospheric Research) model due to a doubling of CO 2 concentration, and concluded that convective process and water vapor feedback are the two major contributors to the tropical uppertropospheric temperature change. In this study, we will apply this CFRAM developed by Lu and Cai (2009) and Cai and Lu (2009) to elucidate the formation of the tropical UTW by quantifying the relative contributions of each feedback process based on reanalysis data. The methodology will be detailed in Sect. 2.
As Ramage (1968) demonstrated, the tropical land generates a much greater amount of heat for export than the equatorial oceanic regions and the troposphere over the Maritime Continent in boreal winter is probably the single greatest energy source for the extratropical atmospheric circulation. Thus, as a secular tropical heating source, the latent heat released from the deep convection over the tropical western Pacific dominates that over other tropical regions and plays an important role in modulating global climate. Located amid the monsoon system, the deep convection over this region strongly interacts with the monsoon circulation. For example, the latent heat released by the enhanced convection in this region caused by strengthened northerly winds along the east coastline during strong East Asian winter monsoon events in turn strengthens the winter monsoon by enhancing the upper branch of local meridional circulation and accelerating the upper-tropospheric East Asian westerly jet stream (Yang and Webster 1990;Yang et al. 2002;Chang et al. 2004). The energy generated from the released latent heat can also exert an impact on the climate over not only East Asia, but also North America and Europe via teleconnection patterns such as the Pacific-Japan pattern and Pacific-North American pattern (He and Wu 2014;Ding et al. 2011). The anomalous circulation of the Pacific-Japan pattern, in turn, can intensify the regional convective activity via enhancing evaporation and moisture convergence and dynamically inducing anomalous ascent (Kosaka and Nakamura 2010;Wu et al. 2014). As Gill (1980) illustrated a considerable east-west asymmetric circulation by solutions for heating concentrated in a region of finite extent, the intensity and center of tropical heating also indicate the zonal tropical circulation such as that related to the El Niño-Southern Oscillation. Thus, the variation of the tropical heating over the tropical western Pacific (i.e., the UTW) can exert major impacts on global climate variations. Thus, this study will be mainly focused on the formation of the UTW over the tropical western Pacific.
The major objectives of the present work are: (1) to describe the observed long-term change in the UTW over the tropical western Pacific; (2) to address which radiative or non-radiative process is the main contributor to the formation of the long-term change; and (3) to understand how the variation of the main processes contributes to the long-term change in UTW.
Data and methods will be described in Sect. 2. Section 3 discusses the spatial and temporal variations of UTW. Attribution analyses of long-term change in UTW based on annual and monthly mean data will be presented in Sects. 4 and 5, respectively. Section 6 summarizes the major findings of this study.

Data and methods
The primary data set used in this study is the European Centre for Medium-range Weather Forecasts Reanalysis Interim (ERA-Interim, Dee et al. 2011). It covers the period from 1979 to the present and has a horizontal resolution of 1° longitude × 1° latitude with 37 pressure levels in the vertical ranging from 1000 to 1 hPa. Monthly outgoing long-wave radiation (OLR) from National Oceanic and Atmospheric Administration with a horizontal resolution of 2.5° by 2.5° is also used in this study.
The major analysis tool is the CFRAM, same as the method applied in Song et al. (2014), which is based on 1 3 the total energy balance within an atmosphere-surface column. We consider the total energy balance equation Yang and Ren 2017) separately for one decade as 1985-1996 and the other decade as 2004-2013, take the difference (∆) between the two decades (i.e., 2004-2013 minus 1985-1996), and then obtain where S (R) is the vertical profile of the net convergence (divergence) of short-wave (long-wave) radiation flux within each layer. Above the surface, ΔQ non−radiative is the vertical profile of the convergence of total energy due to atmospheric turbulent, convective, and advective motions. At the surface, ΔQ non−radiative corresponds to the loss of energy due to surface heat fluxes, as well as the net energy convergence in the ocean mixed layer if the surface is over ocean. The elements of E t represent change in the rate of energy storage. All terms in Eq. (1) are in units of W/m 2 .
As in Li et al. (2017), we express ∆S and ∆R as the sum of partial radiative energy flux convergence/divergence perturbations due to individual radiative feedback processes: In Eq. (2), superscripts "solar", "wv", and "α" represent solar insolation, water vapor, and surface albedo, respectively. R T is the Planck feedback matrix. Substituting Eq.
(2) into Eq. (1), rearranging the terms, we obtain The variables required as input to the radiative transfer model Liou 1992, 1993) including solar insolation at the top of the atmosphere (TOA), air/surface temperatures, specific humidity, cloud amount, cloud liquid/ ice water content, surface albedo, and ozone mixing ratio are all obtained from the ERA-Interim. ΔQ ATD in Eq. (3) is zero at the surface layer, and ΔQ ATD = −Δ(S − R) atmos in the atmospheric layers, representing the vertical profile of the energy perturbation in the atmosphere associated with atmospheric motions (including turbulent, convective (i.e., the latent heat release and adiabatic warming/cooling), and large-scale motions) and heat storage anomalies (expected to be very small), which is named as atmospheric dynamical process (ATD). ΔQ surf is zero in atmosphere layers and ΔQ surf = −Δ(S − R) surf at the surface layer, representing the energy perturbation at the surface due to surface turbulent sensible ( ΔQ sshf ) and latent ( ΔQ slhf ) heat flux anomalies, the changes in the rate of heat storage (of land or oceans) and the oceanic circulation if the surface is over oceans. The sum of the oceanic dynamics and land/ocean heat storage term (Ocn) is estimated as the residual of the surface energy balance equation, Equation (3) allows us to express the vertical profile of the temperature difference between two specific decades or months in the atmosphere-surface column at a given horizontal grid as the sum of the vertical profiles of the partial temperature changes (PTCs) due to the changes in solar insolation, water vapor, cloud, surface albedo, atmospheric dynamics, and oceanic dynamics and ocean/land heat storage.
To quantify the relative contributions of each individual process to 3D temperature anomalies between 250 and 400 hPa over a certain region based on its spatial pattern and mean amplitude, we compute the pattern-amplitude projection ( PAP i ) using the partial temperature anomalies ( ΔT iz ) in layer Z associated with each individual process i. Following Park et al. (2012), we have where and λ are latitude and longitude, respectively, a is the mean radius of the earth, ∆T is the observed temperature anomalies, Z is the vertical level, and A represents the region over which spatially averaged values are computed.
By definition, the sum of PAP coefficients at a certain layer is one (K). The difference between an area-averaged PTC and a PAP is that the latter emphasizes the "spatial similarity" between the PTC distribution and the total temperature change distribution and is thus a better metric for assessing the overall contribution of a specific process to the temperature change over the region considered.

Characteristics of UTW
The pattern of tropospheric warming (TW) trends at 250-1000 hPa (Fig. 1a) indicates that the warming over the middle-and high-latitudes is stronger than that over the lower latitudes, and the Northern Hemisphere warms faster than the Southern Hemisphere. In the tropics, significant warming covers the tropical Atlantic, the tropical Indian Ocean, and the tropical western Pacific. The IPCC AR4 has demonstrated that in the context of global warming, the upper troposphere especially above 400 hPa warms faster than the surface. The TW trends at 250-400 hPa are shown in Fig. 1b, which presents that the spatial distribution over the tropics resembles that in Fig. 1a, while the UTW over the high latitudes is not as strong as the warming at 250-1000 hPa, indicating that the lower-tropospheric and surface warming mainly contributes to the total TW over the high latitudes. Over the eastern Indian Ocean, the Maritime Continent, and the tropical western Pacific, the troposphere is the warmest area in the tropics (see contours), which also reveals a relatively high warming rate. Thus, we will focus on these regions from 100 to 160°E as shown by the purple box (Fig. 1a, b) in the following analysis. The vertical distribution of TW weighted averaged over 20°S-20°N depicts a pronounced UTW between 200 and 500 hPa with a higher magnitude (> 2 K/decade) in 250-400 hPa over the whole tropical region except the eastern tropical Pacific (Fig. 1c). Compared with the magnitude and amplitude of the UTW, the surface warming is relatively weaker globally and even presents cooling trends in the eastern tropical Pacific. The meridional-vertical distribution of TW averaged over 100°-160°E shows a significant warming core over 30°S-30°N between 200 and 500 hPa, with the strongest over 250-400 hPa (Fig. 1d). Another pronounced warming area is located in the lower troposphere north of 60°N, mainly concentrated below 700 hPa, also known as the Arctic warming (e.g., Cai 2005;Graversen et al. 2008).
UTW over 20°S-20°N/100°-160°E shows an interdecadally strengthening trend with a stronger interannual variation before the late-1990s, exactly in 1998, and a weaker interannual variation since then (Fig. 2). This feature is depicted more clearly by the 5-yr smoothing result as shown by the solid black curve in Fig. 2. The warming trend of UTW seems more prominent since the late-1990s, although the global surface warming experienced a hiatus of debate (Guan et al. 2015;Hu et al. 2017). We use two periods, 1986-1995 and 2004-2013, to represent the anomalously cooling and warming decades to demonstrate the main contributors to the interdecadal warming (IDW) by analyzing the changes in the radiative and non-radiative (dynamic and thermodynamic) processes between these two periods.

Attribution analysis
Based on the CFRAM analysis, the contributions of individual processes to the change in the mean temperature over 250-400 hPa between 1986250-400 hPa between and 1995250-400 hPa between and 2004250-400 hPa between -2013 are quantitatively revealed by their PTCs. As the spatial patterns shown in Fig. 3, the change in each variable of solar insolation, surface albedo, ozone, CO 2 , and water vapor contribute to the IDW at a magnitude about 0.1 K. Specifically, the changes in solar insolation contribute positively (positive PTCs) south of 10°S and negatively (negative PTCs) north of 10°N, due mainly to its interdecadal variation. The variation of surface albedo provides a positive contribution to the IDW over the Maritime Continent and the tropical western Pacific. PTCs caused by changes in ozone are uniformly negative while that caused by the change in CO 2 is uniformly positive with a large contribution band between 10°S and 10°N. Also, the more (less) water vapor over land and the Observations are shown in (l) with values significant at the 95% confidence level dotted 1 3 eastern Indian Ocean (the western Pacific) contributes to the UTW positively (negatively). Moreover, the change in cloud as a radiative process contributes negatively north of 5°S due to more cloud and slightly positively south of 10°N resulting from less cloud. Overall, among all the radiative processes, cloud exerts the largest impact on the IDW. The change in ATD warms most of the analysis region especially over the ocean domain north of 0°N, except the cooling effect over land. By the function of upward longwave radiation from the surface, the change in Ocn warms the region, with the largest contribution over the ocean east of 120°E, aside from the cooling effect over the southern South China Sea (SCS) and northern Australia. The spatial pattern of PTCs due to the change in surface latent heat flux (SLHF) resembles the opposite feature that is due to Ocn but with an opposite sign. The PTCs induced by the change in surface sensible heat flux (SSHF) exhibit a pattern similar to those caused by SLHF but with a smaller magnitude. Figure 1k shows the total temperature change calculated by adding all the PTCs of each individual process, which resembles the observed feature shown in Fig. 1l.
Based on the spatial distribution of the PTCs of each process, we explicitly illustrate the relatively essential contributions by the variations of cloud, ATD, Ocn, and SLHF to the formation of UTW by analyzing their PAP coefficients, which quantify the relative contributions of each process from both the spatial pattern and magnitude. An analysis of PAP indicates that the changes in ATD and Ocn are the two largest positive contributors to the UTW over 20°S-20°N/100°-160°E, with CO 2 as a secondary contributor, while variations of cloud and SLHF are the two largest negative contributors, with ozone as a third one (Fig. 4). The negative contribution of variations of cloud is mainly related to the increase in the middle and high clouds (figure not shown). And the offline error is negligible. Thus, the long-term change in the annual mean UTW over the tropical western Pacific is mainly modulated positively by ATD and Ocn and negatively by cloud and SLHF.

Long-term change in oceanic processes
As discussed before, the UTW caused by the change in Ocn occurs mainly via upward long-wave radiation from the surface, which is associated with the oceanic thermal status and dynamical process (here, land heat content is negligible). The forcing due to the variation of Ocn (Fig. 5a) shares the same spatial distributions as the PTCs (Fig. 3h). Here, a positive value means that energy is transferred from oceans to the atmosphere, and vice versa. Compared with the long-term change in the heat content of the 0-200-m mixed layer, the enhanced forcing over the ocean surface east of the Maritime Continent is contributed mainly by the more heat storage in the mixed layer, which is a result of the enhanced easterly trade wind. The long-term change in the heat content in the mixed layer shows a strengthened zonal gradient, with a thicker thermocline in the tropical western Pacific and a thinner thermocline in the tropical eastern Pacific (figure not shown), and the impact of such an enhanced gradient on the atmospheric circulation will be discussed in the following sub-section. The change in the oceanic dynamics plays a key role in the southern SCS, working against the warming effect of enhanced heat content in the mixed layer, leading to a net cooling effect on the surface temperature. As shown by the contours in Fig. 5a, the long-term change in surface heat fluxes including SLHF and SSHF are well collocated with the distribution of forcing, as the anomalous downward transport over the negative forcing and the enhanced upward transport over the positive forcing, which is also induced by the enhanced surface trade wind. Moreover, the long-term change in Ocn exerts an impact on the atmospheric circulation via anomalous heat fluxes. Thus, the overall positive

Long-term changes in atmospheric processes
Both the horizontal and vertical temperature transport strengthens over the research domain (figure not shown), and the latter is one order larger than the former. Following analyses are mainly focused on the long-term change in the vertical motion. Accompanied by the enhanced deep convection (indicated by OLR shown as the red curve in the lower panel of Fig. 6), significantly stronger ascending motion occurs mainly over 20°S-5°N/110°-160°E. An enhanced zonal circulation is observed in the tropics, with the moreeastward upward branch located over the ocean northeast of Australia. As shown in Fig. 6a, anomalous upward motion  is mainly located between 130 and 160°E, while the anomalously strengthened deep convection presents between 90 and 150°E. The long-term change in diabatic heating from the latent heat release of condensation exerts a warming effect between 110 and 160°E with a center around 150°E. The UTW in 250-400 hPa caused by ATD is mainly contributed by diabatic heating, which works against the adiabatic cooling of enhanced upward motion. The strengthened upward branch of the Hadley Cell is located over 20°S-5°N and centered between 10 and 0°N (Fig. 6b), while enhanced rainfall occurs between 20 and 20°N. Although diabatic heating is also enhanced, the total effect of ATD is to cool the troposphere over 10°S-0°N, implying that adiabatic cooling plays a key role. However, accompanied by such enhanced deep convection over the southern part, anomalous downward motions dominate over the northern portion of the domain (between 0° and 20°N) with less diabatic warming released. Thus, the UTW at 250-400 hPa over 0°-20°N is a result of enhanced adiabatic warming and weaker diabatic cooling. As the warming in the north is more dominant than the cooling in the south, overall warming exists at 250-400 hPa for the whole region. Therefore, the UTW induced by the ATD at 250-400 hPa in the region of consideration is mainly contributed by the interdecadally enhanced deep convection over 20°S-5°N/110°-160°E.

Attribution analysis
A further analysis of the long-term variations based on monthly mean data is expected to reveal more detailed characteristics of the UTW. The vertically-weighted mean temperature over 250-400 hPa/20°S-20°N/100°-160°E between the 2 decades is in the range of about 0.24-0.75 K. As illustrated in Fig. 7, the UTW peaks from December to March and is relatively weak in boreal summer and fall except in July and September. Figure 8 depicts the PAP results for long-term variation based on the CFRAM analysis using monthly data, similar to the results shown in Fig. 4 using annual mean data. As relatively weaker contributors, ozone and CO 2 act as consistently negative and positive contributors all year round. Relatively larger contributors in each month are the same as those for the annual mean, which include Ocn, ATD, SLHF, and cloud. The two main negative contributors for the annual mean, SLHF and cloud, show various seasonal variations: SLHF dominates in January-to-June, except February, and cloud plays an essential role in July-to-December.The two prominent positive contributors, Ocn and ATD, also exhibit different seasonal variations. Comparison of the PAP coefficients indicates that Ocn plays a key positive role from January to July, while ATD is the main contributor from August to December, which is similar to the seasonal evolution of SLHF and cloud. Thus, the UTW with various magnitudes in each month is modulated by different main radiative and non-radiative processes. Overall, the main positive contributor to the UTW in boreal winter and spring is Ocn, while that to the UTW in boreal summer and fall is ATD. Next we will briefly focus on the long-term changes in Ocn and ATD based on the seasonal mean instead of the monthly mean.

Long-term changes in oceanic processes
In boreal winter (Fig. 9a), the spatial pattern of Ocn forcing resembles that of the annual mean in Fig. 5a, with a strongly positive forcing center over the tropical western Pacific and anomalous negative forcing over the southern SCS that is not so pronounced as that shown by the annual mean. In boreal spring (Fig. 9b), almost all areas show positive forcing with a weaker strength than that in winter, while negative anomalies still exist in the southern SCS. In boreal summer (Fig. 9c), negative anomalies expand from the Bay of Bengal to the western Pacific east of the SCS, and also appear from the eastern Indian Ocean to northwestern Australia. Moreover, the positive forcing anomaly center moves eastward over 160°-180°E. Such differences in the seasonal variations of forcing between the southern and northern SCS also indicate differences in properties of the two parts . In boreal fall, negative forcing anomalies dominate and cover nearly the entire SCS and eastern Indian Ocean while positive anomalies further weaken and move eastward. Thus, as positive forcing gradually weakens and later moves out of the tropical western Pacific and negative forcing gradually strengthens and expands over the SCS, the positive contribution of Ocn fades from winter to fall or even turns into a negative contribution in fall. Moreover, as the anomalous surface heat fluxes indicate, the impact of Ocn on the atmospheric circulation also gradually weakens from winter to fall.

Long-term changes in atmospheric processes
The PTC contributed by ATD is modulated by the long-term variations of atmospheric non-radiative processes, which can be elucidated by the variations of temperature advection and adiabatic heating rate. As indicated by the variation of zonal circulation (Fig. 10a, c, e, g), enhanced convection mainly occurs over 130°-160°E. Specifically, the anomalous upward branch is broadened and tilted westward with height, and easterly anomalies are much enhanced in winter (Fig. 10a) and spring (Fig. 10c), resulting in a strengthened zonal circulation over the Indian Ocean. The upward branch is significantly strengthened over 140°-160°E without a clear tilt in summer (Fig. 10e) and fall (Fig. 10g), indicating that a pronounced zonal circulation is located over the tropical central and eastern Pacific. With such an interdecadally strengthened upward motion, the adiabatic heating released by condensation also shows an enhanced tendency with the strongest centers over 120°-150°E in winter and spring and over 130°-160°E in summer and fall. As shown by the PTC, the increased strength of adiabatic heating is more apparent than that of adiabatic cooling caused by enhanced upward motion, which is even prominent in summer and fall. Accompanied by the seasonally anomalous circulation, deep convection also shows obvious long-term variations (red curves in the upper panels of Fig. 10). In winter, enhanced precipitation is located mainly over 90°-150°E, and the deep convection moves slightly westward. From spring to fall, deep convection shows a more eastward location, especially in summer and fall, which should be involved with the aforementioned strengthened zonal circulation over the tropical Pacific.
As the long-term change in meridional circulation shows (Fig. 10b, d, f, and h), enhanced vertical circulation is located mainly over 10°S-10°N. In winter, upward anomalies are concentrated between 5°S and 5°N, while significantly enhanced precipitation is located around 10°S and 10°N. As the adiabatic cooling is compensated by the diabatic warming that dominates in the region that we focus on, the PTC is positive in the meridional direction. From spring to summer, anomalous upward motion is pronounced south of the equator, and anomalous downward motion is located north of the equator. The negative PTC in the upper troposphere south of the equator is a result of stronger adiabatic cooling (than diabatic warming), which is reversal over north of the equator where the stronger adiabatic warming (than diabatic cooling) leads to positive PTC in the upper troposphere. Such anomalous circulation features are much more prominent in summer and fall. The meridional distribution of deep convection shows a complicated seasonal variation. In boreal winter and spring, significantly enhanced precipitation occurs around 10°S and 10°N, heavier than that near the equator. Enhanced precipitation is located mainly north of the equator in summer and south of the equator in fall. The aforementioned anomalous downward motion north of the equator weakens the convection in situ, which somehow indicates that the upward branch of meridional circulation is more southward in the later decade. Thus, in all seasons the enhanced vertical circulations over 10°S-0°/140°-160°E exert positive contributions to the UTW. Because of the relatively stronger adiabatic cooling and the enhanced westward advection in the upper troposphere in winter and spring, the positive impact of ATD is weaker than that of Ocn in these seasons. From summer to fall, both diabatic heating and adiabatic warming provide positive contributions so that ATD is the dominant contributor to the UTW in August-to-December and February.

Conclusions
The fact that in the tropics the tropospheric temperature trends are larger than the surface warming trends has been studied previously using model outputs and satellite and radiosonde observations. In this study, we revisit this phenomenon using reanalysis data and further illustrate the radiative and non-radiative contributions to the formation of UTW using a climate feedback-response analysis method.
In the tropics, the UTW is pronounced over 250-400 hPa except over the eastern Pacific. Since climatologically the warmest 250-400-hPa UTW is located over the tropical western Pacific, the present study focuses on the UTW over 20°S-20°N/100°-160°E, where a significant interdecadal variation occurred in the air thickness at 250-400 hPa in the late 1990s. Thus, we choose two periods, 1986-1995 and 2004-2013, to represent the anomalously cooling and warming decades, respectively. To analyze the main contributors to this long-term change in the warming, CFRAM analysis is applied to quantify the contributions of radiative and non-radiative (dynamic and thermodynamic) processes.
As revealed by the annual mean data, the changes in Ocn and ATD are the two largest positive contributors while the variations in clouds and SLHF are the two main negative contributors. The cooling effect of ozone and the warming effect of CO 2 are two secondary factors that contribute to UTW. The positive contribution by Ocn is via enhanced upward long-wave radiation due to the large heat storage in the mixed layer. The positive contribution by ATD is modulated mainly by the enhanced diabatic heating induced by the strong convection over the tropical western Pacific, specifically over 10°S-0°/140°-160°E, while the adiabatic warming north of the equator also provides a contribution.
This study further discusses the long-term change for each month. The UTW is strongest in boreal winter from December to March and relatively weaker in summer and fall. The main positive and negative contributors vary with seasons. The prominent feature is that Ocn plays a key role in the formation of UTW in boreal winter and spring while ATD exerts a main impact in boreal summer and fall. Among the negative contributors, SLHF significantly affects the UTW in January-to-June and cloud plays an essential role in July-to-December. The strong and positive contribution of Ocn in winter and spring is modulated by the enhanced ocean heat content in the mixed layer over the tropical western Pacific. This contribution gradually weakens and "leaves" the tropical western Pacific while the negative forcing centered over the SCS gradually strengthens and expands, and as a result the positive contribution of Ocn diminishes from winter to fall and even turns into a negative contribution in fall. ATD exerts positive contributions to the UTW in all seasons. In boreal winter and spring, adiabatic cooling is relatively stronger, and with enhanced westward advection in the upper troposphere the positive impact of ATD is weaker than that of Ocn. In boreal summer and fall, both diabatic heating and adiabatic warming make positive contributions and ATD is the dominant contributor of the UTW. Moreover, a strengthened zonal circulation in the tropics is also observed over both the Indian and the Pacific Oceans, and that over the Pacific is pronounced in summer and fall, facilitating the ATD to dominate.