The simulation of dynamically coupled ice sheet, ice stream, and ice shelf-systems poses a challenge to most numerical ice sheet models. Here we review present ice sheet model limitations targeting a broader audience within Earth Sciences, also those with no specific background in numerical modeling, in order to facilitate cross-disciplinary communication between especially paleoglaciologists, marine and terrestrial geologists, and numerical modelers. The 'zero order' (Shallow Ice Approximation, SIA)-, 'higher order'-, and 'full Stokes' ice sheet models are described conceptually and complemented by an outline of their derivations. We demonstrate that higher order models are required to simulate coupled ice sheet-ice shelf and ice sheet-ice stream systems, in particular if the results are aimed to complement spatial ice flow reconstructions based on higher resolution geological and geophysical data. The zero order SIA model limitations in capturing ice stream behavior are here illustrated by conceptual simulations of a glaciation on Svalbard. The limitations are obvious from the equations comprising a zero order model. However, under certain circumstances, simulation results may falsely give the impression that ice streams indeed are simulated with a zero order SIA model. (C) 2011 Elsevier Ltd. All rights reserved.