Abstract Proxyâbased reconstructions suggest that equilibrium changes in global mean sea surface temperature (ÎGMSST) are nearly equivalent to changes in mean ocean temperature (ÎMOT) on glacialâinterglacial timescales over the past 900,000 years. However, the underlying mechanisms responsible for this relationship remain poorly understood. Here we use simulations from Paleoclimate Modeling Intercomparison Project Phase 3 and 4 (PMIP3/4) to investigate equilibrium ÎMOT and its linkage to sea surface temperature changes between the Last Glacial Maximum (LGM, âź21,000 years ago) and preâIndustrial. Results show that PMIP3/4 simulations generally underestimate proxyâbased ÎMOT. Regression analysis reveals that LGM MOT is strongly modulated by midâlatitude SST cooling, with the Southern Ocean having a greater influence compared to other oceanic regions, thus helping explain why models with similar ÎGMSSTs exhibit significantly different ÎMOTs. Additionally, we find a strong relationship between simulated Antarctic seaâice coverage and Southern Ocean SST changes, with implications for constraining seaâice reconstructions.