Geothermal heat flux (GHF) is one of the key thermal boundary conditions for ice-sheet models. We assess the sensitivity of the Lambert-Amery glacial system in East Antarctica to four different GHF datasets using a regional ice-sheet model. A control solution of the regional model is initialised by minimising the misfit to observations through an optimisation process. The Lambert-Amery glacial system simulation contains temperate ice up to 150 m thick and has an average basal melt of 1.3 mm a−1, with maximum basal melting of 504 mm a−1. The simulations which use a relatively high GHF compared to the control solution increase the volume and area of temperate ice, which causes higher surface velocities at higher elevations, which leads to the advance of the grounding line. The grounding line advance leads to changes in the local flow configuration, which dominates the changes within the glacial system. To investigate the difference in spatial patterns within the geothermal datasets, they were scaled to have the same median value. These scaled GHF simulations showed that the ice flow was most sensitive to the spatial variation in the underlying GHF near the ice divides and on the edges of the ice streams.