Abstract. Observational evidence, including offshore moraines and sediment cores confirm that at the Last Glacial maximum (LGM) the Greenland ice sheet (GrIS) grew to a significantly larger spatial extent than seen at present, grounding into Baffin Bay and to the continental shelf break. Given this larger spatial extent and it is close proximity to the neighboring Laurentide (LIS) and Innuitian Ice sheet (IIS), it is likely these ice sheets will have had a strong non-local influence on the spatial and 15 temporal behaviour of the GrIS. Most previous paleo ice sheet modelling simulations recreated an ice sheet that either did not extend out onto the continental shelf; or utilized a simplified marine ice parametersiation and therefore did not fully include ice shelf dynamics, and or the sensitivity of the GrIS to this non-local signal from the surrounding ice sheets.In this paper, we investigated the evolution of the GrIS over the two most recent glacial-interglacial 20 cycles (240 kyr BP to present day), using the ice sheet-ice shelf model, IMAU-ICE and investigated the influence of the LIS and IIS via an offline relative sea level (RSL) forcing generated by a GIA model. This RSL forcing controlled via changes in the water depth below the developing ice shelves, the spatial and temporal pattern of sub-ice shelf melting, which was parametrised in relation to changes in water depth. 25In the suite of simulations, the GrIS at the glacial maximums coalesced with the IIS to the north, expanded to the continental shelf break to the south west but remained too restricted to the north east.In terms of an ice-volume equivalent sea level contribution, at the Last Interglacial (LIG) and LGM the ice sheet added 1.46m and -2.59m to the budget respectively. The estimated lowering of the sea level by the Greenland contribution is considerably more (1.26 m) than most previous studies indicated 30 whereas the contribution to the LIG high stand is lower (0.7 m). The spatial and temporal behaviour of the northern margin was highly variable in all simulations, controlled by the sub surface melt (SSM), which was dictated by the RSL forcing and the glacial history of the IIS and LIS. In contrast, the southwestern part of the ice sheet was insensitive to these forcing's, with a uniform response in an all simulations controlled by the surface air temperature (SAT) forcing, derived from ice cores. 35Clim. Past Discuss., https://doi