Abstract Zonal and meridional surface temperature gradients are considered to be determinants of large-scale atmospheric circulation patterns. However, there has been limited investigation of these gradients as diagnostic aids. Here, the twentieth-century variability in the Northern Hemisphere equator-to-pole temperature gradient (EPG) and the ocean–land temperature contrast (OLC) is explored. A secular trend in decreasing EPG and OLC is noted. Decadal and interannual (ENSO-related) variations in the joint distribution of EPG and OLC are identified, hinting at multistable climate states that may be indigenous to the climate or due to changing boundary forcings. The NH circulation patterns for cases in the tails of the joint distribution of EPG and OLC are also seen to be different. Given this context, this paper extends past efforts to develop insights into jet stream dynamics using the Lorenz-1984 model, which is forced directly and only by EPG and OLC. The joint probability distribution of jet stream and eddy energy, conditional on EPG and OLC scenarios, is investigated. The scenarios correspond to (i) warmer versus colder climate conditions and (ii) polarized ENSO phases. The latter scenario involves the use of a heuristic ENSO model to drive the Lorenz-1984 model via a modulation of the EPG or the OLC. As with GCMs, the low-order model reveals that the response to El Niño forcing is not similar to an anthropogenic warming signature. The potential uses of EPG and OLC as macro-level indicators of climate change and variability and for comparing results across GCMs and observations are indicated.