The position and strength of the Hadley cell circulation determine the habitable zones in the tropics, yet our understanding of and ability to predict changes in the circulation is limited. One potentially important source of uncertainty is the dependence of the Hadley cell on turbulent drag. Here, the sensitivity of the Hadley cell and associated features such as the intertropical convergence zone to variations in the magnitude of the turbulent drag is explored with an atmospheric general circulation model in aquaplanet configuration. The tropical circulation and precipitation, and extratropical features such as the polar jet stream, displayed a strong sensitivity to the strength of the parameterized turbulent drag, with distinct low- or high-drag regimes. However, the response of the meridional heat transport produced a surprising departure from previous expectations: with greater drag, simulations exhibited less heat transport than low-drag simulations, which is in the opposite sense to that from Held and Hou. This may be due to the energetic constraints in the present model framework. When exposed to a uniform global warming, the response of the ITCZ precipitation depends strongly on the choice of drag, whereas most simulations exhibit a poleward expansion of the subtropics.  

The NASA Atmospheric Tomography (ATom) mission built a photochemical climatology of air parcels based on in situ measurements with the NASA DC-8 aircraft along objectively planned profiling transects through the middle of the Pacific and Atlantic oceans. In this paper we present and analyze a data set of 10 s (2 km) merged and gap-filled observations of the key reactive species driving the chemical budgets of O₃ and CH₄ (O₃, CH₄, CO, H₂O, HCHO, H₂O₂, CH₃OOH, C₂H₆, higher alkanes, alkenes, aromatics, NO_x, HNO₃, HNO₄, peroxyacetyl nitrate, and other organic nitrates), consisting of 146 494 distinct air parcels from ATom deployments 1 through 4. Six models calculated the O₃ and CH₄ photochemical tendencies from this modeling data stream for ATom 1. We find that 80 %–90 % of the total reactivity lies in the top 50 % of the parcels and 25 %–35 % in the top 10 %, supporting previous model-only studies that tropospheric chemistry is driven by a fraction of all the air. Surprisingly, the probability densities of species and reactivities averaged on a model scale (100 km) differ only slightly from the 2 km ATom 10 s data, indicating that much of the heterogeneity in tropospheric chemistry can be captured with current global chemistry models. Comparing the ATom reactivities over the tropical oceans with climatological statistics from six global chemistry models, we find generally good agreement with the reactivity rates for O₃ and CH₄. Models distinctly underestimate O₃ production below 2 km relative to the mid-troposphere, and this can be traced to lower NO_x levels than observed. Attaching photochemical reactivities to measurements of chemical species allows for a richer, yet more constrained-to-what-matters, set of metrics for model evaluation.

The NASA Atmospheric Tomography (ATom) mission built a photochemical climatology of air parcels based on in situ measurements with the NASA DC-8 aircraft along objectively planned profiling transects through the middle of the Pacific and Atlantic oceans. In this paper we present and analyze a data set of 10 s (21(m) merged and gap-filled observations of the key reactive species driving the chemical budgets of O-3 and CH4 (O-3, CH4, CO, H2O, HCHO, H2O2, CH3OOH, C2H6, higher alkanes, alkenes, aromatics, NOx, HNO3, HNO4, peroxyacetyl nitrate, other organic nitrates), consisting of 146 494 distinct air parcels from ATom deployments 1 through 4. Six models calculated the O-3 and CH4 photochemical tendencies from this modeling data stream for ATom 1. We find that 80 %-90 % of the total reactivity lies in the top 50 % of the parcels and 25 %-35 % in the top 10 %, supporting previous model-only studies that tropospheric chemistry is driven by a fraction of all the air. In other words, accurate simulation of the least reactive 50 % of the troposphere is unimportant for global budgets. Surprisingly, the probability densities of species and reactivities averaged on a model scale (100 km) differ only slightly from the 2 km ATom data, indicating that much of the heterogeneity in tropospheric chemistry can be captured with current global chemistry models. Comparing the ATom reactivities over the tropical oceans with climatological statistics from six global chemistry models, we find excellent agreement with the loss of O-3 and CH4 but sharp disagreement with production of O-3. The models sharply underestimate O-3 production below 4 km in both Pacific and Atlantic basins, and this can be traced to lower NOx levels than observed. Attaching photochemical reactivities to measurements of chemical species allows for a richer, yet more constrained-towhat-matters, set of metrics for model evaluation.

In this paper we introduce a Bayesian framework, which is explicit about prior assumptions, for using model ensembles and observations together to constrain future climate change. The emergent constraint approach has seen broad application in recent years, including studies constraining the equilibrium climate sensitivity (ECS) using the Last Glacial Maximum (LGM) and the mid-Pliocene Warm Period (mPWP). Most of these studies were based on ordinary least squares (OLS) fits between a variable of the climate state, such as tropical temperature, and climate sensitivity. Using our Bayesian method, and considering the LGM and mPWP separately, we obtain values of ECS of 2.7K (0.6-5.2, 5th-95th percentiles) using the PMIP2, PMIP3, and PMIP4 datasets for the LGM and 2.3K (0.5-4.4) with the PlioMIP1 and PlioMIP2 datasets for the mPWP. Restricting the ensembles to include only the most recent version of each model, we obtain 2.7K (0.7-5.2) using the LGM and 2.3K (0.4-4.5) using the mPWP. An advantage of the Bayesian framework is that it is possible to combine the two periods assuming they are independent, whereby we obtain a tighter constraint of 2.5K (0.8-4.0) using the restricted ensemble. We have explored the sensitivity to our assumptions in the method, including considering structural uncertainty, and in the choice of models, and this leads to 95% probability of climate sensitivity mostly below 5K and only exceeding 6K in a single and most uncertain case assuming a large structural uncertainty. The approach is compared with other approaches based on OLS, a Kalman filter method, and an alternative Bayesian method. An interesting implication of this work is that OLS-based emergent constraints on ECS generate tighter uncertainty estimates, in particular at the lower end, an artefact due to a flatter regression line in the case of lack of correlation. Although some fundamental challenges related to the use of emergent constraints remain, this paper provides a step towards a better foundation for their potential use in future probabilistic estimations of climate sensitivity.

The Earth's equilibrium climate sensitivity (ECS) to a doubling of atmospheric CO2, along with the transient climate response (TCR) and greenhouse gas emissions pathways, determines the amount of future warming. Coupled climate models have in the past been important tools to estimate and understand ECS. ECS estimated from Coupled Model Intercomparison Project Phase 5 (CMIP5) models lies between 2.0 and 4.7 K (mean of 3.2 K), whereas in the latest CMIP6 the spread has increased to 1.8-5.5 K (mean of 3.7 K), with 5 out of 25 models exceeding 5 K. It is thus pertinent to understand the causes underlying this shift. Here we compare the CMIP5 and CMIP6 model ensembles and find a systematic shift between CMIP eras to be unexplained as a process of random sampling from modeled forcing and feedback distributions. Instead, shortwave feedbacks shift towards more positive values, in particular over the Southern Ocean, driving the shift towards larger ECS values in many of the models. These results suggest that changes in model treatment of mixed-phase cloud processes and changes to Antarctic sea ice representation are likely causes of the shift towards larger ECS. Somewhat surprisingly, CMIP6 models exhibit less historical warming than CMIP5 models, despite an increase in TCR between CMIP eras (mean TCR increased from 1.7 to 1.9 K). The evolution of the warming suggests, however, that several of the CMIP6 models apply too strong aerosol cooling, resulting in too weak mid-20th century warming compared to the instrumental record.