We revisit clear-sky one-dimensional radiative–convective equilibrium (1D-RCE) and determine its equilibrium climate sensitivity to a CO2 doubling (ECS) and associated uncertainty. Our 1D-RCE model, named konrad, uses the Rapid Radiative Transfer Model for GCMs (RRTMG) to calculate radiative fluxes in the same way as in comprehensive climate models. The simulated radiative feedbacks are verified by a line-by-line radiative transfer model, with which we also investigate their spectral distribution. Changing the model configuration of konrad enables a clear separation between the water vapor and the lapse rate feedbacks, as well as the interaction between the two. We find that the radiative feedback and ECS are sensitive to the chosen relative humidity profile, resulting in an ECS range of 2.09–2.40 K. Using larger CO2 forcings we find that the radiative feedback changes up to 10% for surface temperatures of 291–299 K. Although the ECS is similar to previous studies, it arises from the compensation of a larger clear-sky forcing (4.7 W m−2) and more strongly negative feedbacks (−2.3 W m−2 K−1). The lapse rate feedback and the feedback from the interaction of lapse rate and humidity compensate each other, but the degree of compensation depends on the relative humidity profile. Additionally, the temperature profile is investigated in a warming climate. The temperature change at the convective top is half as large as at the surface, consistent with the proportionally higher anvil temperature hypothesis, as long as the humidity is consistently coupled to the temperature profile.
Equilibrium climate sensitivity (ECS), the change in surface temperature in response to a doubling of atmospheric CO2, is arguably one of the most important quantities when discussing climate change. Since the pioneering work by Manabe and Wetherald (1967), a hierarchy of models has been developed to simulate Earth’s reaction to an external forcing. But even for the most simple models in this hierarchy, such as radiative–convective equilibrium (RCE) under fixed relative humidity, the value of the ECS is not known with precision (Schlesinger 1986), which makes constraining ECS for more realistic scenarios a challenge that continues to this day (Stevens et al. 2016).
The simple framework of one-dimensional RCE is fundamental to the understanding of climate change. Such a framework forms the theoretical backdrop against which more complex models have been developed and compared, and has greatly influenced the climate debate since its very beginning (National Research Council 1979).
In this paper we report on efforts to formulate a precise RCE problem, estimate its ECS, and explore its sensitivity to its numerical formulation and the choice of the assumed relative humidity profile. We assume clear skies, fixed relative humidity, fixed ozone, and a saturated isentropic lapse rate. Not taking clouds into account is a deliberate choice, as their inclusion introduces a considerable number of additional degrees of freedom, which are best tackled separately. The precise formulation of the problem allows other models to be compared to our calculations, thereby establishing a benchmark.
Because it is not immediately clear which relative humidity profile is most representative of the atmosphere’s response to warming, or even the extent to which RCE is sensitive to the homogeneity assumption implied by the selection of a single relative humidity profile, we explore the sensitivity of our calculations to different profiles of relative humidity. We do so by constructing an ensemble of different humidity distributions, which allows us to explore how the response to greenhouse forcing depends on the assumed profile.
The clear-sky formulation allows us to address the question as to whether there is a state dependence in the simplest possible representation of the climate system, as has been argued in a number of studies. Meraner et al. (2013), for instance, argue that an enhanced water vapor feedback leads to a state dependence of the ECS, whereas Colman and McAvaney (2009) find a balancing of water vapor and lapse rate feedback in a warmer climate in a global climate model (GCM).
For our RCE simulations we use a radiation scheme that trades accuracy for computational efficiency. By using line-by-line radiative transfer simulations we evaluate the fidelity of these base calculations and pinpoint how different radiative feedbacks are distributed spectrally.
Besides its impact on ECS, water vapor also impacts the temperature profile and its evolution in a warming climate. Hartmann and Larson (2002) described its role as a “thermostat” for the temperature of the tops of the deepest clouds in the fixed anvil temperature (FAT) hypothesis. A more general approach that takes changes in static stability into account results in a proportionally higher anvil temperature (PHAT; Zelinka and Hartmann 2010). Both FAT and PHAT are based on the interplay of clear-sky radiation and water vapor. In our study, the RCE model is run with different vertical relative humidity distributions and forced with a wide range of CO2 concentrations to quantify its agreement with the FAT, or PHAT, hypothesis.
2. Problem formulation
The model is configured in a way to simulate the tropical atmosphere, which is to first order in radiative–convective equilibrium (Popke et al. 2013). Like the original Manabe and Wetherald (1967) model, the new one-dimensional RCE model (1D-RCE) is based on a radiation scheme that calculates radiative heating for each atmospheric layer (see section 2a). We apply this heating to determine a provisional temperature profile including a coupled surface temperature. Afterward this predicted profile is modified through a convective adjustment (see section 2c). If no convective adjustment is made, the model calculates a pure radiative equilibrium, which results in a convectively unstable thermal structure at lower levels of the atmosphere. Finally, the absolute humidity is recalculated based on the new temperature profile and the assumed relative humidity (see section 2b). Alternatively, the model implements the assumption of fixed absolute humidity, if this step is omitted. The steps are repeated iteratively until an equilibrium state is reached. The full model flow is illustrated in Fig. 1. The model, konrad, is developed under the MIT License and is available on github.com/atmtools/konrad. The simulations in this study have been performed using version 0.6.6 (Kluft and Dacie 2019).
a. Radiation (RRTMG)
For computing radiation, we use the Rapid Radiative Transfer Model for GCMs (RRTMG; Mlawer et al. 1997). RRTMG is used in various GCMs and has been compared to reference radiative transfer models (Pincus et al. 2015). As part of this study, we also compare RRTMG to a benchmark line-by-line code to both test our implementation of RRTMG and assess the sensitivity to the representation of radiative transfer. We use RRTMG through the Climate Modelling and Diagnostics Toolkit (CliMT) Python interface (Monteiro et al. 2018).
Simulations are performed under clear-sky conditions without a diurnal cycle. Thus the choice of the solar geometry attempts to account for these simplifications so as to yield a control simulation that is not too different from conditions representative of the present-day Earth. Just decreasing the solar constant is not sufficient as the optical path through the atmosphere would then be underestimated. Following Cronin (2014) we choose a solar irradiance of 510 W m−2 and a zenith angle of 47.88°. This results in an incoming shortwave flux of 342 W m−2 at the top of the atmosphere while also keeping the amount of atmospheric absorption in the right proportion. The latter is important for computing radiative heating rates. The configured incoming radiation is lower than observed values in the tropics to compensate for energy transport to the extratropics, which is missing in konrad (Popke et al. 2013). A description of the resultant temperature profile for a present-day atmospheric composition is given in section 3.
RRTMG uses the correlated-k method to allow efficient calculation of fluxes integrated over the full electromagnetic spectrum (Mlawer et al. 1997). The scheme is tuned by using output data from line-by-line models to identify significant parts of the spectrum. The atmospheric states covered by this optimization are based on present-day climate as well as modified conditions like a doubling of CO2 concentration. Errors (other than those associated with an incorrect implementation) can occur when the climate differs significantly from the underlying reference (e.g., strongly increased CO2 or temperatures that exceed precalculated lookup tables). We compared the longwave heating rates for different equilibrium states to the line-by-line radiative transfer model ARTS (Atmospheric Radiative Transfer Simulator; Buehler et al. 2018; Eriksson et al. 2011). The differences found are smaller than 0.1 K day−1 in the troposphere and smaller than 0.6 K day−1 above the tropopause, in agreement with a comparison by Pincus et al. (2015). In addition, a comparison between the radiative feedbacks calculated using RRTMG and ARTS is given in section 4e.
b. Atmospheric state
The atmosphere is discretized in mass space, and thus we adopt fixed pressure p coordinates. The pressure grid is constructed using linearly increasing step sizes in ln(p) and is approximately given by
with surface pressure ps (1000 hPa), pressure at model top pt (1 Pa), and level index i in the range of [0, N]. This ensures adequate resolution of the upper stratosphere while keeping a sufficient number of grid points in the troposphere. The latter is important to resolve changes of the convective top and the cold-point tropopause height more accurately.
We performed simulations with the number of vertical levels ranging from 100 to 1000. The equilibrium surface temperatures for the lowest and highest vertical resolution differ by about 0.3 K. The difference vanishes with increasing number of levels and only affects absolute temperatures; the climate sensitivity is almost independent of the vertical resolution. In the following, simulations are performed with 500 vertical levels to better resolve the tropopause region. Our model atmosphere includes prescribed vertical profiles of relative humidity and ozone as a function of p. In addition, several gases with constant volume mixing ratios are included (see Table 1). The trace gas concentrations are chosen to be consistent with the ones used for the Radiative–Convective Equilibrium Model Intercomparison Project (RCEMIP; Wing et al. 2018).
The ozone volume mixing ratio is a function of pressure p and is also chosen to follow the RCEMIP prescription:
where g1 = 3.6478 ppmv, g2 = 0.832 09, and g3 = 11.3515 hPa. The ozone profile is set in the beginning of the simulation and kept constant throughout the whole simulation. In a companion study, we investigate how changes in this profile affect the equilibrium state (Dacie et al. 2019).
As baseline for the troposphere, we prescribe the vertical relative humidity distribution as a function of p in the same way as Manabe and Wetherald [1967, their Eq. (3)]:
with relative humidity at the surface RHs (77%) and surface pressure ps (1000 hPa).
This profile differs from the observed mean state of humidity mostly through the lack of a secondary peak in the upper troposphere. Our motivation for choosing this profile was the comparison to the historical study as well as the simplicity of its specification. In sections 4d and 5b, we explore the consequences of this assumption by investigating how our results change for a vertically uniform relative humidity as well as a profile with a second upper-tropospheric humidity peak.
Relative humidity is defined with respect to saturation over water for temperatures above 0°C and with respect to saturation over ice for temperatures below −23°C. For intermediate temperatures the equilibrium pressure is computed as a combination of the values over water and ice according to the Integrated Forecasting System (IFS) documentation [ECMWF 2018, their Eq. (12.13)]. The equilibrium pressures are calculated using empirical formulations by Murphy and Koop (2005).
The absolute humidity in the stratosphere is kept constant at the volume mixing ratio found at the cold-point tropopause. This leads to a moistening the stratosphere if the tropopause warms with surface warming.
c. Convective adjustment
Following Manabe and Wetherald (1967), we perform an energy-conserving convective adjustment in the troposphere, which acts to cool the surface and warm the atmosphere. Our implementation of the convective adjustment differs from that described in Manabe and Wetherald (1967) and Liou (2002), and is described in detail in the appendix. Briefly, we adjust our temperature profile to the saturated isentropic lapse rate, thus taking into account potential energy changes and contribution from the enthalpy of vaporization associated with convection in a saturated column. Although the relative humidity profile is not saturated, use of the saturated isentropic lapse rate is justified by the assumption that atmospheric convection only occurs in saturated regions, which cover a very small fraction of the tropics but set the temperature profile for the whole tropics in accordance with the weak temperature gradient balance (Sobel and Bretherton 2000; Charney 1963).
The way the convective adjustment redistributes energy results in a distinct convective top. Below this convective top the atmosphere is in a state of radiative–convective equilibrium, and above in a state of radiative equilibrium. In a companion study we explore the consequences of relaxing this hard adjustment (Dacie et al. 2019).
The saturated isentropic lapse rate depends on the atmospheric temperature and pressure and is calculated following Bohren and Albrecht [1998, their Eq. (6.111)]:
with gravitational acceleration g, isobaric specific heat capacity cp, enthalpy of vaporization lυ, saturation mixing ratio ws, gas constants for the surrounding (dry) air Rd and water vapor Rυ, and temperature T (see Table 2 for values). As a simplification lυ is assumed to be constant, thus not varying with temperature.
The assumed lapse rate does not account for fusion enthalpies and thus neglects contributions from the ice phase. The saturated isentropic lapse rate allows a temperature feedback to be taken into account: In a warming climate, the saturated isentropic lapse rate gets less steep, because ascending saturated air parcels are moister and therefore release more latent heat. Hence, the upper troposphere will warm more than the surface (Manabe and Stouffer 1980, their section 5). The strength of this lapse rate feedback is quantified alongside other decomposed feedbacks in section 4c.
We assume a slab surface with an albedo of 0.2 and a heat capacity of 215 MJ m−2 K−1. It can be interpreted as being a well-mixed ocean with a depth of 50 m. The heat capacity damps the surface warming rate and prevents strong vertical temperature changes, which could otherwise occur in single time steps. The chosen total heat capacity of the surface is a significant tuning parameter for the time scales of the model. For our studies, however, only the equilibrium states are relevant, which we have verified as being independent of the chosen heat capacity. Chosen surface constants are given in Table 2.
3. Control climate
Figure 2 shows the equilibrium temperature profiles for different CO2 concentrations as well as a temperature climatology for comparison. The climatology is based on tropical ocean profiles (30°S–30°N) from ERA5 from January 2008 to May 2018.
Konrad captures the temperature structure of the tropical atmosphere as represented by ERA5. The failure of konrad to form a sharper temperature inversion at the tropopause is indicative of the absence of several processes, like overshooting convection and the Brewer–Dobson circulation, both of which act to cool the atmosphere in those heights through adiabatic cooling (Dacie et al. 2019). In addition, radiative effects of clouds may also play a role.
Overall the RCE framework allows us to simulate a tropical temperature profile that is in qualitative agreement with observations. The resulting temperature profile and surface temperature are used as reference for our sensitivity studies.
4. Radiative feedback
a. Comparison to historical setup
In a first experiment, we reproduce the simulations by Manabe and Wetherald (1967) with our RCE model. Konrad is run under different assumptions about the humidity and lapse rate to quantify their effects, especially on feedbacks. The different model configurations are forced with a doubling of CO2. An overview of the corresponding ECS estimates for konrad as well as historical studies for comparison is given in Table 3.
The most basic configuration is an RCE with fixed absolute humidity distribution and a constant lapse rate of 6.5 K km−1. This configuration has an ECS of 1.34 K, which is in very close agreement with Manabe and Wetherald (1967).
Next, calculations are performed in which the relative humidity is kept constant following the profile specified in Eq. (3). This introduces a positive water vapor feedback that increases the ECS by 1.31 K (+98%) to a value of 2.65 K. This is slightly larger than Manabe and Wetherald (1967); however, Manabe and Wetherald (1975) mention a decrease in sensitivity after incorporating a superior radiation model, consistent with the interpretation that the larger value that they obtain arises from differences in the treatment of the radiative transfer.
Finally, the constant lapse rate is replaced by a saturated isentropic lapse rate. This introduces a negative lapse rate feedback that reduces the surface warming by 0.56 K (−21%). The estimate for an RCE with combined water vapor–lapse rate feedback is 2.09 K, which is in very good agreement with previous clear-sky estimates.
We conclude that the historical ECS estimates are in good agreement with, but slightly higher than, our modern implementation of a 1D-RCE model. The remaining differences are likely attributed to discrepancies in the exact composition of the atmosphere or details in the radiative transfer calculations.
b. Climate radiative feedback
We study the sensitivity of konrad to an abrupt change in the CO2 concentration by a factor of 0.5, 2, 4, or 8. The instantaneous radiative forcing for a doubling of the CO2 concentration is 2.92 W m−2 for the reference case, which incidentally is in perfect agreement with the mean of the CMIP5 models found by Collins et al. (2006). However, for different model configurations used throughout this study we find instantaneous forcings between 2.75–3.15 W m−2, which we mainly attribute to differences in the temperature profiles (Huang et al. 2016).
In contrast, the effective radiative forcing, which includes stratospheric adjustment, is 4.73 W m−2, which is larger than the number of 3.7 W m−2 often cited in literature. Possible reasons for this deviation from global mean estimates are the different temperature structure (Huang et al. 2016), an overestimation of the stratospheric cooling due to a missing ozone–temperature feedback (Shepherd and Jonsson 2008), and a missing cloud masking effect.
The surface temperature changes for all simulations are given in Table 4. In addition, the climate feedback parameter
is given. It quantifies how a surface temperature change ΔTs feeds back to the radiative imbalance ΔFTOA.
The climate feedback parameter λ is determined using a method introduced by Gregory et al. (2004): we regress the radiative imbalance at the top of the atmosphere against the surface temperature change for every time step. The climate feedback parameter is defined as the fitted slope. We exclude the phase of stratospheric adjustment by only using time steps after the net radiative flux reaches its peak. The Gregory method is justified by the almost perfect linear relationship between the radiative imbalance and surface temperature change (see Fig. 3). A comparison of the regression line with the actual data points shows differences that are smaller than 0.05 W m−2.
Figure 3 shows the effective radiative forcing Feff (y intercept of the linear regression), the surface temperature change ΔTs (x intercept), and climate feedback parameter (slope of the dashed lines) for each simulation. For every doubling in the CO2 concentration there is an almost constant increase in the simulated surface warming. As a consequence, the total climate feedback parameter is almost constant at 2.33 W m−2 K−1. Differences in the temperature change for different CO2 doublings are mostly attributable to changes in the forcing [a similar increase in forcing with increasing CO2 is also reported by Gregory et al. (2015)]. Although the radiative transfer becomes increasingly unreliable as one moves away from reference concentrations, konrad runs stably for CO2 concentrations ranging from 0.25 to 128 times the present-day values and the feedbacks do not diverge substantially from those inferred by extrapolating from the conditions shown here.
c. Decomposed feedbacks
Furthermore, we want to quantify the magnitude of the different radiative feedbacks in konrad. A common approach is the radiative kernel method described by Soden et al. (2008). The method allows response patterns of the Earth system to be transferred into radiative feedbacks using offline radiative transfer simulations. By contrast, our model formulation and the linearity of the climate feedback allow us to decompose the radiative feedback using different model configurations that leave out specific processes:
The Planck feedback λPL is defined by the regression of radiative imbalance against surface temperature for an experiment in which the tropospheric lapse rate and the absolute humidity are fixed constant.
In the water vapor configuration the relative humidity profile is fixed. The water vapor feedback λWV is defined as the increase of λ between the Planck and water vapor configuration. Note that, in contrast to existing literature (Soden et al. 2008), this definition of the water vapor feedback does not consider changes in the temperature lapse rate.
In the lapse rate configuration the temperature profile is convectively adjusted to a saturated isentropic lapse rate that is calculated from the atmospheric state [Eq. (4)] but with fixed absolute humidity. The lapse rate feedback λLR is defined as the increase of λ between the Planck and lapse rate configuration.
In the reference configuration the atmosphere is adjusted to a coupled saturated isentropic lapse rate and the relative humidity profile is fixed. Here, the water vapor and lapse rate feedbacks act combined, with introduces a nonlinear feedback caused by the lapse rate–driven increase of upper-tropospheric humidity. We define the magnitude of this additional water vapor–lapse rate feedback λWV∧LR such that the individual feedbacks add up to the total radiative feedback:
In the definition of Soden et al. (2008) the water vapor feedback includes both λWV and λWV∧LR. The decomposed magnitudes of the Planck, water vapor, lapse rate, and water vapor–lapse rate feedback are shown in Table 4.
The nonconstancy (increase) of the forcing with successive CO2 doubling is accompanied by a decrease in the strength of the Planck feedback. In addition, there is a slight strengthening of the water vapor feedback due to increased absorption in the atmospheric window (Koll and Cronin 2018). The lapse rate feedback on the order of 1.9 W m−2 K−1 also increases for stronger forcings. When evaluating changes in the feedback parameter one has to keep in mind that different model configurations, including or neglecting different adjustment processes, have different climate sensitivities. As a result, even an eightfold increase in CO2 does not lead to a large temperature change for the Planck and lapse rate configurations, which are both run with fixed absolute humidity. Considering all adjustment processes in our standard configuration, the changes of the water vapor and the lapse rate feedback have compensating signs, so that their sum is almost constant (Colman and McAvaney 2009; Cronin and Wing 2017; Soden et al. 2008; Vial et al. 2013).
The additional water vapor–lapse rate feedback is 1.47 W m−2 K−1, close to the value of the water vapor feedback. The combined λWV∧LR decreases more strongly with the radiative forcing than both λWV and λLR. The decrease is on the same order of magnitude as the decrease of the Planck feedback, resulting in an almost constant total feedback (compensating signs). Through our feedback decomposition, we find the upper-tropospheric component of the water vapor feedback (λWV∧LR) to be essential for the compensation. For our simulations, the climate feedback parameter changes by 2.1% across simulations spanning four doublings of the atmospheric CO2 concentration (Table 4). Over this range the state (temperature) changes by 6.5 K about a working temperature near 295 K. Over the same temperature range, Meraner et al. (2013) estimate a change in sensitivity of 30% (for their constant 80% relative humidity profile), which they attribute to changing feedbacks. Konrad’s climate sensitivity changes by 26% for these four doublings, similar to what is reported by Meraner et al. (2013), but almost all of konrad’s change can be attributed to the forcing increasing with progressive doublings rather than from changes in feedback. The near constancy of konrad’s net feedback arises from a balancing of individual feedbacks, which each change more than their net (see Table 4). Experiments using different solar forcing, much larger changes in CO2 (both not shown) or different humidity profiles (see section 4d)—which allow us to push the model to much larger temperatures—eventually experience larger changes in sensitivity. This appears to be indicative of larger changes in feedbacks (rather than just forcings) and is being investigated in a further study.
In general, both the radiative forcing and feedbacks are larger than GCM estimates from radiative kernels. The water vapor and lapse rate feedbacks are almost twice as large in magnitude compared to global mean estimates by Soden et al. (2008; for GCMs) and Cronin and Wing (2017; for 3D-RCE). The discrepancy may, to some extent, arise from differences in the spatial sampling. Cronin and Wing (2017) find much larger water vapor feedbacks, in closer agreement with our results, for high surface temperatures in their small-domain 3D-RCE configuration (T. Cronin 2019, personal communication). Another reason may be related to missing cloud masking in our clear-sky configuration. In comprehensive Earth system models the feedback in the middle and lower troposphere might be covered by cloud layers. This interpretation is consistent with the much larger water vapor kernels in the clear as compared to the all-sky calculations (Soden et al. 2008, their Fig. 2).
d. Sensitivity to assumed water vapor distribution
The assumption of a relative humidity profile that is exponentially decaying with height is not realistic for the tropical atmosphere. Observations have shown that especially the upper troposphere is much more humid (Zelinka and Hartmann 2011) and Minschwaner and Dessler (2004) have found the upper-tropospheric humidity to be important for the climate sensitivity. In the absence of overlying clouds, the midtropospheric humidity also plays an important role in governing the radiative feedback (Soden et al. 2008; Vial et al. 2013). Therefore, we test how our results are affected by the choice of the vertical relative humidity distribution. Although there are more realistic models of the tropical relative humidity distribution (Romps 2014), we deliberately decide to choose models that are simple in their specification.
In the tropical atmosphere, cloud detrainment causes an upper-tropospheric humidity (UTH) peak at the convective top. We modify Eq. (3) by including a UTH peak that is mathematically described using a normal distribution in ln(p) space:
with pressure level pUTH and magnitude rUTH (75%) of the peak. The actual relative humidity distribution is given by the maximum of the Manabe model [Eq. (3)] and the UTH peak [Eq. (7)] for each pressure level:
The UTH peak is idealized in order to fulfill two requirements. First, a shift of the UTH peak does not change the relative humidity in the lower troposphere, which also has a significant impact on the climate sensitivity and could obscure possible UTH effects. Second, a shift of the UTH peak does not change its magnitude. Both requirements are necessary to perform comprehensible sensitivity studies. The UTH peak is coupled to the temperature profile by setting its location pUTH equal to the pressure at the convective top as one might expect to happen if convective detrainment levels correlate with the height of the convective top (Zelinka and Hartmann 2010).
Another approach is to assume a vertically uniform relative humidity distribution, which is often used when constructing atmospheric profiles for radiative feedback simulations (Meraner et al. 2013; Koll and Cronin 2018; Thuburn and Craig 2002). We have chosen a vertically uniform relative humidity of 40% to account for subsaturated subsidence regions, which cover a vast part of the tropics. The chosen distribution also allows us to perform stable simulations over a wide range of radiative forcings. Higher humidities above roughly 60% produce a runaway feedback in konrad. This is in accordance with Pierrehumbert (1995), who finds the tropical atmosphere to be close to a runaway greenhouse if the fraction of subsidence regions, which act as “radiator fins,” decreases. In reality, if the tropics alone were close to a runaway greenhouse state and they did warm more than the rest of the planet, transport to middle and high latitudes would also increase.
Figure 4 shows the modified relative humidity profiles with UTH peaks at 125, 170, and 225 hPa (green) alongside the Manabe (blue) and the vertically uniform (yellow) distribution. The relative humidity distribution is used to calculate the water vapor amount in the troposphere. In the stratosphere the absolute humidity is fixed to the volume mixing ratio at the cold-point tropopause.
Experiments were first performed with different surface humidities for each profile (not shown) and indicate that for a given form of the humidity distribution the feedback varies roughly linearly with the integrated water vapor. This approximate linearity justifies the representation of the tropical atmosphere using a mean humidity profile in a single-column model.
Differences in ECS as a function of the given humidity profile are up to 20%. Some of these differences can be explained by the aforementioned, and near-linear, dependence on the integrated water vapor implied by each profile, but the shape of the profile also is important. To explore these effects we decompose the radiative feedbacks (Fig. 5). The Planck feedback is on the order of −3.5 W m−2 K−1 for all humidity distributions. Small differences arise from the fact that the initial states for each relative humidity distribution differ. The weakening of the Planck feedback for stronger forcings (darker colors) due to increased CO2 absorption is a robust feature for all configurations.
The water vapor feedback is on the order of 1.7 W m−2 K−1 for all configurations, with the coupled UTH peak giving the strongest feedback of 1.8 W m−2 K−1. The vertically uniform distribution has a slightly lower feedback compared to the other distributions. This is caused by decreased water vapor absorption in the atmospheric window due to the rather low relative humidity of 40%. All configurations show a slight increase of the water vapor feedback for stronger radiative forcings, which can be explained by a closing of the atmospheric window between 800 and 1200 cm−1 (Koll and Cronin 2018).
The additional water vapor–lapse rate feedback results from an increase in upper-tropospheric humidity associated with the change in temperature lapse rate. The increase in humidity leads to an upward shift of the emission level to a colder temperature, leading to a positive feedback, which we find to be in the range of 1.47–1.71 W m−2 K−1. For the Manabe and coupled UTH distribution, the relative humidity decreases with height, which limits the upward shift of the emission level. As a result, the water vapor–lapse rate feedback decreases under warming. This leads, in combination with an increased Planck feedback, to an almost constant total feedback. For the vertically uniform distribution, however, the water vapor–lapse rate feedback is almost constant irrespective of the applied warming. As the Planck feedback is no longer balanced, the total radiative feedback λ tends to show a state dependence. The increase of the climate feedback between a doubling and an octupling of the CO2 concentration is 7.4% for a surface temperature change of about 8 K.
In summary, the strength of the individual radiative feedbacks depends on the absolute humidity the atmosphere contains and how that humidity is distributed. A vertically uniform humidity distribution prevents a strengthening of the additional water vapor–lapse rate feedback under increased radiative forcing. This leads, in combination with a decreasing Planck feedback, to a slightly larger state dependence of the total climate feedback than for the reference configuration.
e. Comparison to a line-by-line radiative transfer model
We perform line-by-line radiative transfer simulations to verify the radiative feedbacks calculated with RRTMG and to gain insight in their spectral distribution. ARTS is used to calculate the outgoing longwave radiation (OLR) spectrum EOLR(ν) on a frequency grid that covers 30 000 equidistant grid points between 1 and 3000 cm−1. Gas absorption was taken into account by using the HITRAN (high-resolution transmission) database (Gordon et al. 2017) and the MT_CKD model (named for its authors: Mlawer, Tobin, Clough, Kneizys, and Davies) for the continuum absorption of water vapor, CO2, and molecular nitrogen (Mlawer et al. 2012, version 2.52).
We simulate the OLR after the stratospheric adjustment as well as in equilibrium to calculate the spectral radiative feedback λν, which can be integrated to determine the radiative feedback:
Here, the spectral radiative feedback is calculated by dividing the difference of the simulated OLR spectra by the surface temperature change. As ARTS is only capable of simulating the longwave radiative feedback the shortwave component is used from RRTMG. The results are shown as diamonds alongside the RRTMG feedbacks (circles) in Fig. 5. The Planck and lapse rate feedbacks are almost the same for ARTS and RRTMG (differences ≤1%). We find that RRTMG slightly underestimates the water vapor feedback by about 3% and overestimates the total radiative feedback by about 5%. The differences seem to be systematic and slightly increase from 4% to 6% for a doubling and an octupling of CO2 respectively. The comparison to the line-by-line model is robust for all vertical humidity distributions regarded in our study.
In addition, the line-by-line simulations can be used to interpret the spectral fingerprint of different radiative feedbacks. Figure 6 shows the OLR spectrum EOLR for konrad in reference configuration as a qualitative baseline. In addition, the 10 cm−1 running mean is shown to smooth the fractal character of the spectrum. The spectrum can be qualitatively divided into three spectral regions that are of interest for the radiative feedbacks: first, an atmospheric window region between 800 and 1200 cm−1, in which the overall absorption is weak and limited to water vapor absorption in the lower troposphere; second, the optically thick CO2 band between 500 and 800 cm−1; and third, a less absorbent region below 500 cm−1 dominated by upper-tropospheric water vapor absorption.
Figure 7 shows the spectral radiative feedbacks for different humidity distributions. The results are smoothed using a 10 cm−1 running mean. In the top-left panel (labeled PL, for Planck), results for konrad in the Planck configuration are shown. The spectral radiative feedback is almost the same for all humidity distributions, with small differences mainly caused by differences in the initial states. As reference, the radiative feedback of a blackbody at surface temperature is shown. The blackbody curve envelops the atmospheric Planck feedback, which is significantly decreased by CO2 and water vapor absorption.
In the top-right panel, the spectral water vapor (WV) feedback is shown, which is positive throughout the whole spectrum. In the CO2 band between 500 and 800 cm−1 the water vapor feedback is close to zero (Ingram 2010). Differences in the vertical humidity distribution can be directly observed in the spectral space: The vertically uniform distribution is drier in the lower troposphere than the other distributions, which leads to a smaller radiative feedback in the atmospheric window region between 800 and 1200 cm−1 (Koll and Cronin 2018). In contrast, the radiative feedback below 500 cm−1 is increased due to higher values of upper-tropospheric humidity.
The middle-left panel shows the spectral lapse rate (LR) feedback. The lapse rate feedback is governed by the temperature profile, so the curves are similar for all humidity distributions. The strongest lapse rate feedback is located around 500 cm−1, where the emission level of the OLR is in the upper troposphere, which warms more than the surface. There is another spectral region near 1300 cm−1 where the emission level is similar, but it constitutes much less radiative feedback due to a smaller temperature dependence of the Planck function there.
The additional combined water vapor–lapse rate feedback (WV∧LR) is shown in the middle-right panel. Like the water vapor feedback, it is positive throughout the whole spectrum, but it is mainly dominant for wavenumbers below 500 cm−1, which corresponds to the upper troposphere. The vertically uniform relative humidity results in the strongest water vapor–lapse rate feedback, as a deepening of the troposphere is accompanied by a persistent moistening of the upper troposphere. This leads to the most striking difference in the total spectral radiative feedback, which is the sum of all other feedbacks (bottom panel): For wavenumbers below 500 cm−1 the radiative feedback for the vertically uniform distribution is decreased by a factor of 2 compared to the other humidity distributions.
It is worth noting that the division of the traditional water vapor feedback into two terms, WV and WV∧LR, leads to a deeper understanding of the tendency for intermodel differences between water vapor and lapse rate feedbacks to cancel, as noted by many GCM studies: It is really the intermodel difference in LR and WV∧LR feedbacks that tends to cancel, since both are driven by changes in the tropospheric lapse rate.
Furthermore, it is worth noting that the feedback in the far-infrared spectral region below 600 cm−1 dominates the individual LR, WV, and WV∧LR feedbacks. Due to compensating effects the total feedback is up to an order of magnitude smaller than the individual feedbacks (Ingram 2010; Koll and Cronin 2018). Nonetheless, the total feedback in this region is particularly sensitive to differences in the vertical humidity structure, as it affects the relative strength of the individual feedbacks.
5. PHAT mechanism
a. Historical studies and modern theory
In a next step the vertical temperature structure in and around the tropopause and its evolution under increased CO2 concentrations are investigated. Early computations of ECS using RCE first started to include clouds on fixed pressure levels (FAP) as this was the easiest approach computationally. It later became apparent that it might be more plausible to consider fixing high clouds at temperature, rather than pressure levels (Augustsson and Ramanathan 1977), an idea whose theoretical justification was formulated in terms of the fixed anvil temperature hypothesis (FAT) by Hartmann and Larson (2002). The chosen approach has a significant impact on ECS, leading to estimates ranging from 1.21 to 2.31 K (Cess 1976; Ramanathan and Coakley 1978). Zelinka and Hartmann (2010) have shown that, if anvil clouds are associated with the maximum radiatively driven divergence, changes in static stability result in proportionally higher anvil temperatures with increasing surface temperature (PHAT). Latest observations and numerical simulations have shown PHAT to be a more realistic description of changes to cloud-top temperatures (Cronin and Wing 2017; Kuang and Hartmann 2007; Zelinka and Hartmann 2010, 2011). As both FAT and PHAT are based on clear-sky radiation, one of them should be verifiable in an RCE framework. We use the temperature at the convective top as proxy for the cloud-top temperature and how it might change with surface warming (Fig. 8).
b. Sensitivity to assumed water vapor distribution
The convective-top temperature in our reference configuration (blue circles) increases 1.17 K for a doubling of the CO2 concentration, which is roughly half the observed surface temperature change of 2.09 K and a fifth of what FAP predicts. The modified humidity distributions with a coupled UTH have a slightly weaker rate of increase of convective-top temperature with surface warming. Still they are well described by the PHAT hypothesis. Simulations with larger CO2 perturbations show that this ratio is robust for a wide range of surface temperature changes, consistent with GCM findings by Zelinka and Hartmann (2010).
In a follow-up experiment, the UTH peak is prescribed at different pressure levels between 125 and 225 hPa (every 25 hPa). For these configurations, the present-day equilibrium surface temperatures range from 294.5 to 296.6 K, which is about 1–3 K warmer than the reference. Every model configuration is forced with consecutive CO2 doublings between 0.5 and 8 times the present-day concentration. The presence of a fixed UTH peak has a significant impact on the evolution of convective-top temperatures (Fig. 8): Increasing the humidity below the level of the convective top (light green triangles) increases the rate of change of convective-top temperature with surface warming. A UTH peak that is fixed at the initial convective-top pressure is close to the PHAT line, but has a slightly weaker rate of increase of convective-top temperature with surface warming.
Most striking are the results for the two peaks located above the convective top (dark triangles). The additional humidity above the convective top increases the radiative cooling leading to a FAT-like behavior. The highest UTH peak at 125 hPa even leads to negative changes in convective-top temperatures, although this simulation is not realistic, as the UTH peak is rather high, which results in a too-moist lower stratosphere.
We conclude that the PHAT hypothesis describes the change in convective-top temperatures well for the vertical humidity distributions considered in our study. When using a more complex vertical humidity distribution its structure has to be coupled to the temperature profile in order to not artificially distort the PHAT mechanism.
c. Climate sensitivity
The effect of the UTH peaks on ECS can be inferred from the x distance between values in Fig. 8. Coupling the UTH peak to the convective top results in ECS estimates between the Manabe and a vertically uniform relative humidity distribution. A fixed UTH peak at different prescribed heights results in a spread of ECS estimates from 2.28 to 2.77 K, or feedbacks of −2.24 to −1.79 W m−2 K−1 respectively. This should be considered when trying to set up “realistic” 1D-RCE simulations using observed relative humidity distributions: imposing a UTH peak that is not properly coupled to the temperature structure will impact the ECS estimate.
6. Summary and conclusions
We defined and performed benchmark calculations for a simple clear-sky RCE problem. The ECS for konrad in reference configuration with fixed relative humidity, fixed ozone, and saturated isentropic lapse rate is 2.09 K, slightly higher than the 1.95 K found by Manabe and Wetherald (1975).
Furthermore, we quantified the decomposed climate radiative feedbacks. The water vapor and lapse rate feedback are almost twice as large in magnitude as global mean estimates, but compensate each other. Possible reasons for their increase are differences in the model region and a missing cloud masking effect in our clear-sky model. We find that a substantial component of the combined water vapor and lapse rate feedbacks comes from the interaction of the two, which explains their degree of compensation in more complex models.
A comparison of different vertical relative humidity distributions results in ECS ranges from 2.09 to 2.40 K with corresponding radiative feedbacks from −2.34 to −2.03 W m−2 K−1. Changes in the radiative feedbacks with successive doublings are less than 10% for a range of surface temperatures between 291 and 299 K. The strength of the state dependence itself depends on the shape of the assumed humidity profile: A vertically uniform relative humidity distribution results in an almost constant combined water vapor–lapse rate feedback for different forcings. As the Planck feedback decreases with stronger CO2 forcings, the total feedback shows tendencies of a state dependency.
The results are robust for a range of CO2 forcing from 0.5 to 8 times the present-day concentration (i.e., four doublings). A comparison of the radiation scheme RRTMG with the line-by-line model ARTS shows differences in the total climate feedback of about 5%. The errors are of systematic nature and do not substantially increase for stronger forcings.
The change in convective-top temperatures is half as much as the surface temperature change and is consistent with the PHAT hypothesis. The tendency of konrad to verify PHAT rests on the assumption that the relative humidity profile is coupled to the temperature profile, and departures from PHAT become apparent if they become decoupled.
All results are limited to clear-sky radiation. Future work will explore the implications of extra degrees of freedom as introduced by clouds.
We thank Timothy W. Cronin, David M. Romps, and two anonymous reviewers for a great and engaged discussion of the manuscript. This article contains modified Copernicus Climate Change Service Information (ERA5) . We thank Manfred Brath for many fruitful discussions and for providing the setup for the ARTS simulations. The contribution of Stefan A. Buehler was supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy–EXC 177 ‘Climate System Analysis and Prediction’, Contract O1LK1502B and EXC 2037 ‘Climate, Climatic Change, and Society’, Contract 390683824, contributing to the Center for Earth System Research and Sustainability (CEN) of Universität Hamburg. Kerry Emanuel developed an RCE model with the name CONRAD in the 1990s. We thank him for happily allowing us to use the name konrad for our new model. We thank Robert Pincus for fruitful discussion. The model source is developed under the MIT License and available on github.com/atmtools/konrad.
Instantaneous Adjustment of the Troposphere
A convective adjustment is performed if the atmosphere is unstable to convection in the troposphere. The adjustment fixes the temperatures according to a specified lapse rate, given in kelvins per kilometer, which is converted to kelvins per pascal under the assumption of hydrostatic equilibrium.
The change in enthalpy of an atmospheric layer due to a temperature change (ΔT)atm is proportional to ρcp(ΔT)atmdz, where ρ is the density, cp is the heat capacity at constant pressure (also assumed to be constant with height), and dz is the thickness of the layer. This can be rewritten (assuming hydrostatic equilibrium) in pressure coordinates as (cp/g)(ΔT)atmdp, where g is the gravitational acceleration, assumed to be constant with height. The corresponding energy change of the surface is Δzρscs(ΔT)s, where Δz is the thickness of the surface layer, ρs the density, and cs the heat capacity. The surface layer is assumed to be at one temperature (it does not vary with depth). Then, the energy conservation equation required for the convective adjustment is
The integral runs over the whole atmosphere, but there is no convective adjustment above a certain pressure level, pc, called the convective top. The convective temperature change in the atmosphere (ΔT)atm is a function of pressure and is given by
where Trad is the temperature profile after radiative heating rates have been applied and Tcon is the convectively adjusted temperature profile:
Here, γp is the change in temperature of the convective profile with pressure (γp > 0), ps is the surface pressure, and Tcon,s is the convectively adjusted surface temperature. We assume that convection acts to cool the surface and warm the troposphere and thus do not allow convection to cool the atmosphere. The pressure level pc is defined as the highest atmospheric level (lowest pressure level) that satisfies Tcon(pc) ≥ Trad(pc) and no convective adjustment is applied above this level.
The equations above, or similar equations with a different treatment of the surface, hold for the convective adjustment of many previous RCE models, including Manabe and Wetherald (1967) and MacKay and Khalil (1991). However, our computational implementation differs. We start with by guessing a surface temperature, Tcon,s and calculate the corresponding convectively adjusted temperature profile according to Eq. (A3). Then we test how close this profile is to satisfying energy conservation [Eq. (A1)]. We update our surface temperature and repeat the procedure iteratively, until we find a surface temperature and corresponding profile which satisfy Eq. (A1).