The ocean thermohaline circulation is important for transports of heat and the carbon cycle. We present results from PMIP2 coupled atmosphere‐ocean simulations with four climate models that are also being used for future assessments. These models give very different glacial thermohaline circulations even with comparable circulations for present. An integrated approach using results from these simulations for Last Glacial Maximum (LGM) with proxies of the state of the glacial surface and deep Atlantic supports the interpretation from nutrient tracers that the boundary between North Atlantic Deep Water and Antarctic Bottom Water was much shallower during this period. There is less constraint from this integrated reconstruction regarding the strength of the LGM North Atlantic overturning circulation, although together they suggest that it was neither appreciably stronger nor weaker than modern. Two model simulations identify a role for sea ice in both hemispheres in driving the ocean response to glacial forcing.