The long-term climate variations of the Quaternary were primarily influenced by concurrent changes in Earth's orbit, greenhouse gases, and ice sheets. However, because climate changes over the coming century will largely be driven by changes in greenhouse gases alone, it is important to better understand the separate contributions of each of these forcings in the past. To investigate this, idealized equilibrium simulations are conducted in which the climate is driven by separate changes in obliquity, precession, CO 2 , and ice sheets. To test the linearity of past climate change, anomalies from these single-forcing experiments are scaled and summed to compute linear reconstructions of past climate, which are then compared to mid-Holocene and last glacial maximum (LGM) snapshot simulations, where all forcings are applied together, as well as proxy climate records. This comparison shows that much of the climate response may be approximated as a linear response to forcings, while some features, such as modeled changes in sea ice and Atlantic meridional overturning circulation (AMOC), appear to be heavily influenced by nonlinearities. In regions where the linear reconstructions replicate the full-forcing experiments well, this analysis can help identify how each forcing contributes to the climate response. Monsoons at the mid-Holocene respond strongly to precession, while LGM monsoons are heavily influenced by the altered greenhouse gases and ice sheets. Contrary to previous studies, ice sheets produce pronounced tropical cooling at the LGM. Compared to proxy temperature records, the linear reconstructions replicate long-term changes well and also show which climate variations are not easily explained as direct responses to long-term forcings.