Abstract. Until now, the lack of time-continuous, terrestrial paleoenvironmental data from the Pleistocene Arctic has made model simulations of past interglacials difficult to assess. Here, we compare climate simulations of four warm interglacials at Marine Isotope Stages (MISs) 1 (9 ka), 5e (127 ka), 11c (409 ka) and 31 (1072 ka) with new proxy climate data recovered from Lake El'gygytgyn, NE Russia. Climate reconstructions of the mean temperature of the warmest month (MTWM) indicate conditions up to 0.4, 2.1, 0.5 and 3.1 • C warmer than today during MIS 1, 5e, 11c and 31, respectively. While the climate model captures much of the observed warming during each interglacial, largely in response to boreal summer (JJA) orbital forcing, the extraordinary warmth of MIS 11c compared to the other interglacials in the Lake El'gygytgyn temperature proxy reconstructions remains difficult to explain. To deconvolve the contribution of multiple influences on interglacial warming at Lake El'gygytgyn, we isolated the influence of vegetation, sea ice and circum-Arctic land ice feedbacks on the modeled climate of the Beringian interior. Simulations accounting for climate-vegetation-land-surface feedbacks during all four interglacials show expanding boreal forest cover with increasing summer insolation intensity. A deglaciated Greenland is shown to have a minimal effect on northeast Asian temperature during the warmth of stages 11c and 31 (Melles et al., 2012). A prescribed enhancement of oceanic heat transport into the Arctic Ocean does have some effect on Lake El'gygytgyn's regional climate, but the exceptional warmth of MIS l1c remains enigmatic compared to the modest orbital and greenhouse gas forcing during that interglacial.