Past and projected changes in global hydroclimate in Earth system models have been examined. The Budyko framework that relates the partitioning of precipitation into evaporation to a location’s aridity has been modified to account for the effect of interannual terrestrial water storage and compared to traditional methods. The new formulation better fits climate model data over most of the globe. Old and new formulations are used to quantify changes in the spatial patterns of hydroclimate based locally on year-to-year variations water and energy cycle variables. Focus is on multi-model median responses to changing climate. The changes in hydroclimate from preindustrial to recent historical (1965-2014) conditions often have different patterns and characteristics than changes due only to increasing CO2. For simulations with gradually increasing CO2, differing model treatments of vegetation are found specifically to have categorically different impacts on hydroclimate, particularly altering the relationship between aridity and the fraction of precipitation contributing to evaporation in models that predict vegetation changes. Models that predict vegetation phenology have consistently different responses to increasing CO2 than models that do not. Dynamic vegetation models show more widespread but less consistent differences than other models, perhaps reflecting their less mature state. Nevertheless, there is clearly sensitivity to vegetation that illustrates the importance of including the representation of biospheric shifts in Earth system models.

Earth’s hydrological cycle is expected to intensify in response to global warming, with a ‘wet-gets-wetter, dry-gets-drier’ response anticipated. The subtropics (~15-30°N/S) are predicted to become drier, yet proxy evidence from past warm climates suggests these regions may be characterised by wetter conditions. Here we use an integrated data-modelling approach to reconstruct global- and regional-scale rainfall patterns during the early Eocene (~48-56 million years ago), with an emphasis on the subtropics. Model-derived precipitation–evaporation (P–E) estimates in the tropics (0-15° N/S) and high latitudes (>60° N/S) are positive and increase in response to higher temperatures, whereas model-derived P–E estimates in the subtropics (15-30° N/S) are negative and decrease in response to higher temperatures. This is consistent with a ‘wet-gets-wetter, dry-gets-drier’ response. However, some DeepMIP model simulations predict increasing – rather than decreasing – subtropical precipitation at higher temperatures (e.g., CESM, GFDL). Using moisture budget diagnostics we find that the models with higher subtropical precipitation are characterised by a reduction in the strength of subtropical moisture circulation due to weaker meridional temperature gradients. These model simulations (e.g., CESM, GFDL) agree more closely with various proxy-derived climate metrics and imply a reduction in the strength of subtropical moisture circulation during the early Eocene. Although this was insufficient to induce subtropical wetting, if the meridional temperature was weaker than suggested by the DeepMIP models, this may have led to wetter subtropics. This highlights the important role of the meridional temperature gradient when predicting past (and future) rainfall patterns.