The Impact of Cloud Radiative Effects on the Tropical Tropopause Layer Temperatures

A single-column radiative-convective model (RCM) is a useful tool to investigate the physical processes that determine the tropical tropopause layer (TTL) temperature structures. Previous studies on the TTL using the RCMs, however, omitted the cloud radiative effects. In this study, we examine the impact of cloud radiative effects on the simulated TTL temperatures using an RCM. We derive the cloud radiative effects based on satellite observations, which show heating rates in the troposphere but cooling rates in the stratosphere. We find that the cloud radiative effect warms the TTL by as much as 2 K but cools the lower stratosphere by as much as −1.5 K, resulting in a thicker TTL. With (without) considering cloud radiative effects, we obtain a convection top of ≈167 hPa (≈150 hPa) with a temperature of ≈213 K (≈209 K), and a cold point at ≈87 hPa (≈94 hPa) with a temperature of ≈204 K (≈204 K). Therefore, the cloud radiative effects widen the TTL by both lowering the convection-top height and enhancing the cold-point height. We also examine the impact of TTL cirrus radiative effects on the RCM-simulated temperatures. We find that the TTL cirrus warms the TTL with a maximum temperature increase of ≈1.3 K near 110 hPa.


Introduction
The tropical tropopause layer (TTL), ranging from ≈14.5 km to 18.5 km, is a transition layer between the troposphere and stratosphere and has both tropospheric and stratospheric characteristics [1][2][3][4][5][6].Since air enters the stratosphere preferentially through upwelling in the TTL, the TTL exerts important controls over the stratospheric composition and global climate [6][7][8][9].For example, the water vapor amount in the tropical lower stratosphere is largely determined by the cold TTL temperatures through the dehydration that dries the upwelling air [6,[10][11][12][13].In view of the importance of the TTL for the climate system, it is important to understand the physical processes that govern the TTL temperature structures, which include dynamics on a vast range of scales, radiation, and clouds.
A single-column radiative-convective model (RCM) is a useful tool as an essentially correct first approximation to investigate the physical processes that determine the TTL structure and its sensitivities to constituents and dynamical cooling associated with large-scale ascent [14][15][16].In the RCM framework, the TTL is defined as the layer between the cold-point tropopause (i.e., the coldest point in the temperature profile) and the top of the convection.The RCMs are indeed able to capture the fact that the cold-point tropopause is several kilometers above the convectively adjusted region [14,17,18], indicating a cold point decoupled from convection due to radiative processes.
Thuburn and Craig [14] suggested that a localized heating in the 15 µm CO 2 band, owing to the sharp curvature in the temperature profile near the convection top, forces apart the cold point and convection top.Birner [15] examined a stratospheric radiative equilibrium solution by prescribing a tropospheric temperature profile.He found that the separation of the top of the convection and the tropopause in the tropics (i.e., the existence of the TTL) is almost entirely due to stratospheric dynamics.Recently, Birner and Charlesworth [16] investigated the role of water vapor and ozone radiative heating on the detailed temperature structure of the TTL based on the RCM simulations.It was shown that adiabatic cooling due to upwelling reduces the radiative control by water vapor, resulting in a negative transport-radiation feedback while the radiative control by ozone is enhanced due to upwelling, i.e., a positive transport-radiation feedback.All previous studies on the TTL using the RCMs, however, omitted the cloud radiative effects [14][15][16].
Cloud radiative effects play an important role in the TTL radiative energy budget, largely through the anvil clouds associated with the convective outflow between about 10 and 15 km, and thin cirrus clouds driven by the large-scale vertical motion between ≈14.5 to 17 km [19][20][21][22].Here we examine the impact of cloud radiative effects on the simulated TTL temperatures using an RCM.Consistent with the previous studies [14,[16][17][18], a cold-point tropopause is well separated from the convectively adjusted regions in our RCM simulations even without considering the dynamic cooling associated with the large-scale ascent.In contrast to Thuburn and Craig [14], however, we find that the TTL does not owe its existence to the CO 2 radiative effects.We further show that the cloud radiative effect leads to a lower top of the convectively adjusted region but a higher cold-point tropopause, resulting in an increase in the TTL thickness.Section 2 presents the cloud radiative effects derived from observations and Section 3 describes the RCM used in this study and compares the simulations with Thuburn and Craig [14].The impact of cloud radiative effects on simulated TTL temperatures is shown in Section 4, and a summary is given in Section 5.

Cloud Radiative Effects in Tropics
We derived the cloud radiative effect, which is defined as the difference between radiative heating rate profiles in whole sky and those in clear sky, for each 2.5 • × 2.5 • grid box over the tropics and for each month from June 2006 to May 2007 based on observations.Collocated cloud fields from satellite observations including Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO), CloudSat, and a Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the Aqua satellite were used.CALIPSO lidar can detect optically-thin cirrus clouds that prevail in the TTL but the lidar signal becomes attenuated when cloud optical depth is larger than ≈3.0-4.0 [23,24].We used the CALIPSO level 2 version 3 5 km cloud and aerosol layer products.The clouds included in the aerosol layer product as so-called "stratospheric features" were considered here by separating clouds from aerosols following Tseng and Fu [25].For cloud layers that the CALIPSO lidar could penetrate, we used cloud top and cloud base heights from CALIPSO.Cloud optical depth was retrieved following Yang et al. [21] but using an updated effective lidar ratio of 18 sr for ice clouds.Within a given cloud layer, a homogenous vertical distribution of the cloud extinction coefficient and a generalized particle size of 30 µm were assumed [21].Given a cloud optical depth, the derived cloud radiative effect was not sensitive to the assumed ice particle size.The collocated MODIS (version 6) and CALIPSO cloud optical depths agreed well by examining cloud columns that CALIPSO could penetrate with an average difference of less than 20%.For optically thick clouds that CALIPSO could not penetrate, we used the collocated CloudSat cloud geometric heights and vertical profiles of ice and liquid water content scaled by the MODIS cloud optical depth.We employed the CloudSat 2B-CWC-RO dataset [26] and MODIS level 2 MYD 06 cloud product version 6 [27].Noting that the MODIS gave the total cloud optical depth for the whole column, the optical depth for the opaque clouds was obtained by subtracting the total non-opaque cloud optical depth detected using CALIPSO from the MODIS observation as , where τ MODIS is the MODIS observed total column cloud optical depth and ∑ non−opaque τ CALIPSO is the summation of the cloud optical depths retrieved from CALIPSO for the non-opaque cloud layers above the opaque cloud layers.We then used τ opaque to scale the profiles of cloud water content from CloudSat.Here the vertical distributions of cloud particle sizes were derived based on cloud water content according to Heymsfield et al. [28] for ice clouds and Fueglistaler and Fu [29] for water clouds.
To derive the cloud radiative heating based on observations, monthly mean vertical profiles of temperature, water vapor mixing ratio, and ozone mixing ratio were obtained from the ERA-Interim gridded atmospheric reanalysis [30].However, the coarse vertical resolution around the tropical tropopause had limited the ability of the temperature profile from the ERA-Interim to resolve a realistic cold-point tropopause temperature, which is essential to obtain accurate radiative heating rates in the TTL where the radiative relaxation time is long [31].The COSMIC uses the GPS Radio Occultation (RO) limb-sounding technique for sounding the Earth's atmosphere from satellites [32,33].It has the advantages of high accuracy in the upper troposphere and lower stratosphere [34][35][36] with a high vertical resolution [37].COSMIC temperature profiles from June 2006 to May 2007 were obtained from the COSMIC Data Analysis and Archive Center (CDAAC).For each grid box and month, the ERA-Interim and COSMIC temperature profiles were blended as: where weighting function W linearly increases from 0 to 1 within each layer, T ERA and T COSMIC are the ERA-Interim and COSMIC temperature profile, respectively, and T blend-in is the blended temperature (see Figure 1 as an example).For water vapor and ozone profiles from ERA-Interim, we applied a cubic spline interpolation to avoid abrupt jumps near the cold-point tropopause.All atmospheric profiles were interpolated into a finer vertical resolution of ≈100 meter in the TTL following Yang et al. [21].
Atmosphere 2018, 9, x FOR PEER REVIEW 3 of 13 profiles of cloud water content from CloudSat.Here the vertical distributions of cloud particle sizes were derived based on cloud water content according to Heymsfield et al. [28] for ice clouds and Fueglistaler and Fu [29] for water clouds.
To derive the cloud radiative heating based on observations, monthly mean vertical profiles of temperature, water vapor mixing ratio, and ozone mixing ratio were obtained from the ERA-Interim gridded atmospheric reanalysis [30].However, the coarse vertical resolution around the tropical tropopause had limited the ability of the temperature profile from the ERA-Interim to resolve a realistic cold-point tropopause temperature, which is essential to obtain accurate radiative heating rates in the TTL where the radiative relaxation time is long [31].The COSMIC uses the GPS Radio Occultation (RO) limb-sounding technique for sounding the Earth's atmosphere from satellites [32,33].It has the advantages of high accuracy in the upper troposphere and lower stratosphere [34][35][36] with a high vertical resolution [37].COSMIC temperature profiles from June 2006 to May 2007 were obtained from the COSMIC Data Analysis and Archive Center (CDAAC).For each grid box and month, the ERA-Interim and COSMIC temperature profiles were blended as: where weighting function W linearly increases from 0 to 1 within each layer, TERA and TCOSMIC are the ERA-Interim and COSMIC temperature profile, respectively, and Tblend-in is the blended temperature (see Figure 1 as an example).For water vapor and ozone profiles from ERA-Interim, we applied a cubic spline interpolation to avoid abrupt jumps near the cold-point tropopause.All atmospheric profiles were interpolated into a finer vertical resolution of ≈100 meter in the TTL following Yang et al. [21].The NASA Langley Fu-Liou radiation model [38][39][40][41][42][43][44][45] was used to calculate radiative heating rates.A delta-four stream approach for the radiative transfer was used for the shortwave while a delta-two and four stream combination was used for the longwave [39,41,46].The correlated kdistribution method was used to parameterize the non-gray gaseous absorption by H2O, CO2, O3, N2O, and CH4 [38].A treatment of CFCs was included in the 8-12 μm window region [43].The H2O The NASA Langley Fu-Liou radiation model [38][39][40][41][42][43][44][45] was used to calculate radiative heating rates.A delta-four stream approach for the radiative transfer was used for the shortwave while a delta-two and four stream combination was used for the longwave [39,41,46].The correlated k-distribution method was used to parameterize the non-gray gaseous absorption by H 2 O, CO 2 , O 3 , N 2 O, and CH 4 [38].A treatment of CFCs was included in the 8-12 µm window region [43].The H 2 O continuum absorption was included in the whole thermal infrared spectra following the CKD 2.4 [47].For ice clouds, the single scattering properties, including the extinction coefficient, single-scattering albedo, and asymmetry factor, were parameterized with ice water content and generalized effective particle size [40,42].For water clouds, the single scattering properties were parameterized based on Mie calculations using liquid water content and mean effective radius.Monthly and latitudinal dependence of solar insolation and the diurnal variation of solar radiation were taken into account explicitly.
We calculated the cloud radiative effects for all cloud profiles detected using the CALIPSO combined with CloudSat and MODIS for a given month and over a given 2.5 • × 2.5 • grid box.We then averaged them to obtain the monthly mean cloud radiative effect over each grid box, from which we derived the annual mean cloud radiative effects averaged over the tropics.
Figure 2 (left) shows the cloud radiative effects averaged over 20 • N-20 • S (solid line) and western Pacific (15 • N-15 • S, 120 • E-180 • E) (dashed line).The cloud radiative effects were positive throughout the troposphere except near the surface.The three maxima in the cloud radiative effects near 800-900 hPa, 550 hPa, and 200-300 hPa were in line with the trimodal structure of tropical clouds [48][49][50].The maximum heating rates near 200-300 hPa were 0.5 K/day over the tropics and 0.95 K/day over the western Pacific.The cooling effects of clouds in the lower stratosphere above ≈90 hPa level were also seen in Fueglistaler and Fu [29].The cloud radiative effects over the western Pacific were about twice of those averaged over 20 • N-20 • S. Figure 2 (right) shows the cloud radiative effects of TTL cirrus, which was derived as the difference of the radiative heating rates for the whole skies with and without the TTL cirrus.The TTL cirrus is clouds with cloud base height higher than 14.5 km [51,52].The TTL cirrus introduces heating in the TTL, with a maximum heating of 0.09 K/day (0.17 K/day) near 110 hPa over the tropics (western Pacific).The impact of the cloud radiative effects shown in Figure 2 on the TTL temperatures are presented in Section 4.
continuum absorption was included in the whole thermal infrared spectra following the CKD 2.4 [47].For ice clouds, the single scattering properties, including the extinction coefficient, singlescattering albedo, and asymmetry factor, were parameterized with ice water content and generalized effective particle size [40,42].For water clouds, the single scattering properties were parameterized based on Mie calculations using liquid water content and mean effective radius.Monthly and latitudinal dependence of solar insolation and the diurnal variation of solar radiation were taken into account explicitly.
We calculated the cloud radiative effects for all cloud profiles detected using the CALIPSO combined with CloudSat and MODIS for a given month and over a given 2.5° × 2.5° grid box.We then averaged them to obtain the monthly mean cloud radiative effect over each grid box, from which we derived the annual mean cloud radiative effects averaged over the tropics.
Figure 2 (left) shows the cloud radiative effects averaged over 20° N-20° S (solid line) and western Pacific (15° N-15° S, 120° E-180° E) (dashed line).The cloud radiative effects were positive throughout the troposphere except near the surface.The three maxima in the cloud radiative effects near 800-900 hPa, 550 hPa, and 200-300 hPa were in line with the trimodal structure of tropical clouds [48][49][50].The maximum heating rates near 200-300 hPa were 0.5 K/day over the tropics and 0.95 K/day over the western Pacific.The cooling effects of clouds in the lower stratosphere above ≈90 hPa level were also seen in Fueglistaler and Fu [29].The cloud radiative effects over the western Pacific were about twice of those averaged over 20° N-20° S. Figure 2 (right) shows the cloud radiative effects of TTL cirrus, which was derived as the difference of the radiative heating rates for the whole skies with and without the TTL cirrus.The TTL cirrus is clouds with cloud base height higher than 14.5 km [51,52].The TTL cirrus introduces heating in the TTL, with a maximum heating of 0.09 K/day (0.17 K/day) near 110 hPa over the tropics (western Pacific).The impact of the cloud radiative effects shown in Figure 2 on the TTL temperatures are presented in Section 4.

Radiative-Convective Model and Simulations
The radiative-convective model used here closely follows that in Thuburn and Craig [14] for a comparison purpose but using different radiation schemes.Thuburn and Craig [14] employed the broadband radiation scheme developed by Morcrette [53] and Zhong and Haigh [54] and Zhong et al. [55].Here we use the NASA Langley Fu-Liou radiation model to calculate the radiative heating rates.We also used the rapid radiative transfer model (RRTM) [56,57].Since they produce very similar results, we only show results from the RCM with the NASA Langley-Fu-Liou radiation scheme.The vertical domain of the RCM was from 1013 hPa to 0.01 hPa.An equal logarithmic spacing was used for 601 levels.The radiation heating rates were calculated for the entire domain to adjust the temperature.We only analyzed the results below 10 hPa when the simulated temperatures reached equilibrium for all altitudes below a 10 hPa level.We use a solar constant of 1361 Wm −2 , a mean solar zenith angle of 50.5 • , and a daytime fraction of 0.5.
The ozone profile used was that found in Yang et al. [58], Yang et al. [21], and Lin et al. [22], taken from the SHADOZ observations [59], with data above 0.2 hPa obtained from the HALOE [4,60].Similar results were obtained by using a standard tropical ozone profile [61].The CO 2 , CH 4 , and N 2 O were assumed to be well-mixed with a concentration of 356 ppmv, 1.7 ppmv, and 0.28 ppmv, respectively.Following Thuburn and Craig [14], the water vapor profile was determined at each time step with a fixed relative humidity (RH) of 50%, subject to the constraint that the water vapor mixing ratio may not increase with altitude.Thus, there was a constant RH in the troposphere but a constant mixing ratio that begins from the temperature minimum level and extends through the stratosphere.
Following Thuburn and Craig [14], the surface temperature Ts and the tropospheric lapse rate γ are specified as 300 K and 6.5 K/km, respectively, in the RCM.At each time step, if the temperature at altitude z fell below Ts − γz, then it was reset to equal Ts − γz.The model was integrated in time from an arbitrary initial condition until equilibrium was reached.At each time step, we calculated the water vapor profile, updated temperatures with the calculated radiative heating rates, and then applied the convective adjustment.The temperature profile in the convectively adjusted region was thus fixed by the imposed parameters while the top of the convectively adjusted region and the radiative equilibrium profile above were determined by the RCM [14].
In our control simulations using an RH of 50% and CO 2 concentration of 356 ppmv, we obtained a convection top of ≈150 hPa with a temperature of ≈209 K and a cold point at ≈95 hPa with a temperature of ≈204 K.For comparison, Thuburn and Craig [14] obtained the convection top at 145 hPa and at a temperature of 211.4 K, and the cold point at 105 hPa and a temperature of 206.9 K.The small differences may be caused by the different ozone profiles used and the different radiation models used, as well as the differences in some details in the RCMs.As in Thuburn and Craig [14], our RCM in radiative equilibrium produced a TTL layer that separated the cold point from the top of the convectively-adjusted region.
For a further comparison with results from Thuburn and Craig [14], we repeated their sensitivity study to RH by changing the RH from 1% to 90% (keeping all other parameters the same as in the control experiment) (Figure 3).Our results were very similar to those from Thuburn and Craig [14] (see their Figure 7), if not identical.Figure 3 shows that the cold point height was insensitive to the RH but the convection-top height increased with the RH.Both temperatures at the cold point and top of the convective region decreased with the RH.We also carried out the sensitivity runs to carbon dioxide concentration following Thuburn and Craig [14].Here we changed the CO2 mixing ratio from the 356 ppmv value by multiplying a factor ranging from 0 to 10. Figure 4 shows the height (Figure 4a) and temperature (Figure 4b) of the cold point and the convection top versus the carbon dioxide factor (f).The symbol "X" represents the equilibrium states for f = 0. Different from Thuburn and Craig [14], the separation of the cold point and convection top (i.e., the existence of the TTL) was always seen regardless of the CO2 concentrations.For f = 0 (i.e., a zero CO2 concentration), the separation was still about 0.7 km and 1.1 K (Figure 4).It is worth noting that we also obtained a separation of 0.7 km and 1.0 K using the RRTM in the RCM.Therefore, our results indicated that in contrast to Thuburn and Craig [14], the separation of the cold point and the convective top was not due to the CO2 radiative effects.The dependences of the height and temperature of the cold point and the convection top on the carbon dioxide concentration (Figure 4) were also somewhat different from those shown in Thuburn and Craig [14] (see their Figure 5).For example, the convection-top height was insensitive to the CO2 factor in our simulations, while it decreased with the CO2 factor in Thuburn and Craig [14].When f < 0.2, the coldpoint height increased with f in our simulations, while the dependence from Thuburn and Craig [14] was not monotonic.In addition, the temperature of the cold point and the convection top was less sensitive to the CO2 factor in our simulations as compared with Thuburn and Craig [14], who showed a steady increase with the CO2 mixing ratio for f < 0.2.We believe that the differences between the present results and those from Thuburn and Craig [14] are largely associated with the different radiation models used.We have more confidence in our results since we obtained similar results We also carried out the sensitivity runs to carbon dioxide concentration following Thuburn and Craig [14].Here we changed the CO 2 mixing ratio from the 356 ppmv value by multiplying a factor ranging from 0 to 10. Figure 4 shows the height (Figure 4a) and temperature (Figure 4b) of the cold point and the convection top versus the carbon dioxide factor (f ).The symbol "X" represents the equilibrium states for f = 0. Different from Thuburn and Craig [14], the separation of the cold point and convection top (i.e., the existence of the TTL) was always seen regardless of the CO 2 concentrations.For f = 0 (i.e., a zero CO 2 concentration), the separation was still about 0.7 km and 1.1 K (Figure 4).It is worth noting that we also obtained a separation of 0.7 km and 1.0 K using the RRTM in the RCM.Therefore, our results indicated that in contrast to Thuburn and Craig [14], the separation of the cold point and the convective top was not due to the CO 2 radiative effects.The dependences of the height and temperature of the cold point and the convection top on the carbon dioxide concentration (Figure 4) were also somewhat different from those shown in Thuburn and Craig [14] (see their Figure 5).For example, the convection-top height was insensitive to the CO 2 factor in our simulations, while it decreased with the CO 2 factor in Thuburn and Craig [14].When f < 0.2, the cold-point height increased with f in our simulations, while the dependence from Thuburn and Craig [14] was not monotonic.In addition, the temperature of the cold point and the convection top was less sensitive to the CO 2 factor in our simulations as compared with Thuburn and Craig [14], who showed a steady increase with the CO 2 mixing ratio for f < 0.2.We believe that the differences between the present results and those from Thuburn and Craig [14] are largely associated with the different radiation models used.We have more confidence in our results since we obtained similar results using two different modern radiation models.It is interesting to notice that both our simulations and those from Thuburn and Craig [14] show that the height and temperature of the cold point and the convection top had little dependence on the carbon dioxide concentration when f > 0.2.Furthermore, these two studies show similar separations in both height and temperatures, which are about 2.5 km and 4.5 K, respectively, when f > 0.2.
using two different modern radiation models.It is interesting to notice that both our simulations and those from Thuburn and Craig [14] show that the height and temperature of the cold point and the convection top had little dependence on the carbon dioxide concentration when f > 0.2.Furthermore, these two studies show similar separations in both height and temperatures, which are about 2.5 km and 4.5 K, respectively, when f > 0.2.

Impact of Cloud Radiative Effects
The cloud radiative effects derived from observations (Figure 2) were used to examine their impact on the TTL temperature structures in the RCM simulations.The cloud radiative effects were thus considered in our radiative-convective model as a fixed heating rate profile.That is, at each time step, the total radiative heating rates used to update the temperatures were that from the radiation model for the clear-sky condition plus the fixed cloud radiative effects.
Figure 5 (left) shows the temperature profiles simulated from the RCM with (red) and without (blue) considering the cloud radiative effects and Figure 5 (right) shows their differences.The cloud radiative effect warmed the TTL (by as much as 2 K near 130 hPa) but cooled the lower stratosphere (by as much as −1.5 K), which led to a thicker TTL.With (without) considering cloud radiative effects, we obtained a convection top of ≈167 hPa (≈150 hPa) with a temperature of ≈213 K (≈209 K), and a cold point at ≈87 hPa (≈94 hPa) with a temperature of ≈204 K (≈204 K).Therefore, the cloud radiative

Impact of Cloud Radiative Effects
The cloud radiative effects derived from observations (Figure 2) were used to examine their impact on the TTL temperature structures in the RCM simulations.The cloud radiative effects were thus considered in our radiative-convective model as a fixed heating rate profile.That is, at each time step, the total radiative heating rates used to update the temperatures were that from the radiation model for the clear-sky condition plus the fixed cloud radiative effects.
Figure 5 (left) shows the temperature profiles simulated from the RCM with (red) and without (blue) considering the cloud radiative effects and Figure 5 (right) shows their differences.The cloud radiative effect warmed the TTL (by as much as 2 K near 130 hPa) but cooled the lower stratosphere (by as much as −1.5 K), which led to a thicker TTL.With (without) considering cloud radiative effects, we obtained a convection top of ≈167 hPa (≈150 hPa) with a temperature of ≈213 K (≈209 K), and a cold point at ≈87 hPa (≈94 hPa) with a temperature of ≈204 K (≈204 K).Therefore, the cloud radiative effects deepened the TTL by both lowering the convection-top height and enhancing the cold point height.
effects deepened the TTL by both lowering the convection-top height and enhancing the cold point height.Figure 6 is the same as Figure 5 but using the cloud radiative effect over the western Pacific.Because of a significantly larger cloud radiative effect over the western Pacific than the tropical mean (Figure 2), the TTL was warmed by more than 4 K near 140 hPa but cooled by more than 2 K in the lower stratosphere (Figure 6 right).Using the cloud radiative effect over the western Pacific, we obtained a convection top of ≈250 hPa with a temperature of ≈230 K, and a cold point at ≈82 hPa with a temperature of ≈204 K. Figure 6 is the same as Figure 5 but using the cloud radiative effect over the western Pacific.Because of a significantly larger cloud radiative effect over the western Pacific than the tropical mean (Figure 2), the TTL was warmed by more than 4 K near 140 hPa but cooled by more than 2 K in the lower stratosphere (Figure 6 right).Using the cloud radiative effect over the western Pacific, we obtained a convection top of ≈250 hPa with a temperature of ≈230 K, and a cold point at ≈82 hPa with a temperature of ≈204 K.It should be noted that the cloud radiative effects were closely coupled with the convection, which was not considered in our RCM by using a fixed cloud radiative effect profile.For a tropical region with strong convection (e.g., over the western Pacific), we expected a higher convective top.However, this effect was partly compensated by the larger cloud radiative effect that led to a lower convective top (Figure 6).Thus, the cloud radiative effect had a negative feedback to the convection top.The relative role of these two effects could not be addressed using the RCM but can be investigated with a cloud-resolving model that explicitly considers the coupling of convection, clouds, and radiation by prescribing Ts but using various large-scale forcings to control the strength of the convection.
With a weak Coriolis force in the tropics, the direct response to a diabatic heating anomaly was mainly the vertical advection instead of a local temperature change.The temperature response necessary to balance the cloud radiative effects, which was through a modification of longwave radiative heating, would be observed over the entire domain [62].Thus, it was reasonable to use a column model to determine the tropical average temperature response to the tropical average radiative heating perturbation.However, the interpretation of the response to the cloud radiative effects only over the western Pacific (Figure 6) should be done with caution in addition to a sensitivity of the response to different cloud radiative effects.
We also examined the impact of TTL cirrus radiative effects on the RCM-simulated temperatures.The TTL cirrus warmed the TTL with a maximum temperature increase of ≈1.3 K (≈2.5 K) near 110 hPa over the tropics (the western Pacific).Most of the impact of cloud radiative effects on temperatures near 110 hPa was caused by the TTL cirrus clouds.It should be noted that the cloud radiative effects were closely coupled with the convection, which was not considered in our RCM by using a fixed cloud radiative effect profile.For a tropical region with strong convection (e.g., over the western Pacific), we expected a higher convective top.However, this effect was partly compensated by the larger cloud radiative effect that led to a lower convective top (Figure 6).Thus, the cloud radiative effect had a negative feedback to the convection top.The relative role of these two effects could not be addressed using the RCM but can be investigated with a cloud-resolving model that explicitly considers the coupling of convection, clouds, and radiation by prescribing Ts but using various large-scale forcings to control the strength of the convection.
With a weak Coriolis force in the tropics, the direct response to a diabatic heating anomaly was mainly the vertical advection instead of a local temperature change.The temperature response necessary to balance the cloud radiative effects, which was through a modification of longwave radiative heating, would be observed over the entire domain [62].Thus, it was reasonable to use a column model to determine the tropical average temperature response to the tropical average radiative heating perturbation.However, the interpretation of the response to the cloud radiative effects only over the western Pacific (Figure 6) should be done with caution in addition to a sensitivity of the response to different cloud radiative effects.
We also examined the impact of TTL cirrus radiative effects on the RCM-simulated temperatures.The TTL cirrus warmed the TTL with a maximum temperature increase of ≈1.3 K (≈2.5 K) near 110 hPa over the tropics (the western Pacific).Most of the impact of cloud radiative effects on temperatures near 110 hPa was caused by the TTL cirrus clouds.

Summary and Conclusions
A single-column radiative-convective model (RCM) was used to investigate the physical processes that determine the tropical tropopause layer (TTL) temperature structures.Consistent with previous studies, our RCM simulated the cold-point tropopause that is several kilometers above the convectively adjusted region, indicating a cold point decoupled from convection due to radiative processes.In contrast to a previous study, however, we found that the TTL did not owe its existence to the CO 2 radiative effects.
Cloud radiative effects played an important role in the TTL radiative energy budget.However, previous studies on the TTL using the RCMs all omitted the cloud radiative effects.In this study, we examined the impact of cloud radiative effects on the simulated TTL temperatures using an RCM.We derived the cloud radiative effects based on satellite observations, which showed heating rates in the troposphere but cooling rates in the stratosphere.
We found that the cloud radiative effect warmed the TTL by as much as 2 K but cooled the lower stratosphere by as much as −1.5 K, resulting in a thicker TTL.With (without) considering cloud radiative effects, we obtained a convection top of ≈167 hPa (≈150 hPa) with a temperature of ≈213 K (≈209 K), and a cold point at ≈87 hPa (≈94 hPa) with a temperature of ≈204 K (≈204 K).Therefore, the cloud radiative effects widened the TTL by both lowering the convection-top height and enhancing the cold-point height.We also examined the impact of TTL cirrus radiative effects on the RCM-simulated temperatures.We found that the TTL cirrus warmed the TTL with a maximum temperature increase of ≈1.3 K near 110 hPa.

Figure 1 .
Figure 1.An example of the merged temperature profile (black) from the ERA-Interim (red) and COSMIC (blue) temperature profiles.

Figure 1 .
Figure 1.An example of the merged temperature profile (black) from the ERA-Interim (red) and COSMIC (blue) temperature profiles.

Figure 2 .
Figure 2. (Left) Cloud radiative effect (CRE) over the tropics (solid) and western Pacific (dashed).(Right) The same as left but for the TTL cirrus cloud radiative effect.

Figure 2 .
Figure 2. (Left) Cloud radiative effect (CRE) over the tropics (solid) and western Pacific (dashed).(Right) The same as left but for the TTL cirrus cloud radiative effect.

Figure 3 .
Figure 3. (a) Dependence of the convection-top height (solid) and the cold-point height (dashed) on the relative humidity.(b) Dependence of the convection-top temperature (solid) and the cold-point temperature (dashed) on the relative humidity.

Figure 3 .
Figure 3. (a) Dependence of the convection-top height (solid) and the cold-point height (dashed) on the relative humidity.(b) Dependence of the convection-top temperature (solid) and the cold-point temperature (dashed) on the relative humidity.

Figure 4 .
Figure 4. (a) Dependence of the convection-top height (solid) and the cold-point height (dashed) on the carbon dioxide scaling factor.(b) Dependence of the convection-top temperature (solid) and the cold-point temperature (dashed) on the carbon dioxide scaling factor.The "X" symbols represent the states for f = 0 (i.e., zero CO2 mixing ratio).

Figure 4 .
Figure 4. (a) Dependence of the convection-top height (solid) and the cold-point height (dashed) on the carbon dioxide scaling factor.(b) Dependence of the convection-top temperature (solid) and the cold-point temperature (dashed) on the carbon dioxide scaling factor.The "X" symbols represent the states for f = 0 (i.e., zero CO 2 mixing ratio).

Figure 5 .
Figure 5. (Left) Equilibrium temperature profiles from the radiative-convective model with (red) and without (blue) cloud radiative effects over the tropics (20° N-20° S).The square symbols indicate the top of convection.Circles indicate the cold point.(Right) Impact of CRE at the equilibrium temperature.

Figure 5 .
Figure 5. (Left) Equilibrium temperature profiles from the radiative-convective model with (red) and without (blue) cloud radiative effects over the tropics (20 • N-20 • S).The square symbols indicate the top of convection.Circles indicate the cold point.(Right) Impact of CRE at the equilibrium temperature.

Figure 6 .
Figure 6.As in Figure 5 but over the western Pacific.

Figure 6 .
Figure 6.As in Figure 5 but over the western Pacific.