<p>In simulations of the climate during the Last Glacial Maximum (LGM), we employ two different isotope-enabled atmospheric general circulation models (NCAR iCAM3 and MPI ECHAM6-wiso) and use simulated (by coupled climate models) as well as reconstructed (from a new global climatology of the ocean surface duing the LGM, GLOMAP) surface conditions.</p><p>The resulting atmospheric fields reflect the more pronounced structure and gradients in the reconstructions, for example, the precipitation is more depleted in oxygen-18 in the high latitudes and more enriched in low latitudes, especially in the tropical convective regions over the maritime continent in the equatorial Pacific and Indian Oceans and over the equatorial Atlantic Ocean. Furthermore, at the sites of ice cores and speleothems, the model-data fit improves in terms of the coefficients of determination and root-mean square errors.</p><p>In additional sensitivity experiments, we also use the climatologies by Annan and Hargreaves (2013) and Tierney et al. (2020) and consider the impact of changes in reconstructed sea-ice extent and the global-mean sea-surface temperature.</p><p>Our findings imply that the correct simulation or reconstruction of patterns and gradients in sea-surface conditions are crucial for a successful comparison to oxygen-isotope data from ice cores and speleothems.</p>