The Mars Atmosphere-Ice Coupler MAIC-2 is a simple, latitudinal model, which consists of a set of parameterisations for the surface temperature, the atmospheric water transport and the surface mass balance (condensation minus evaporation) of water ice. It is driven directly by the orbital parameters obliquity, eccentricity and solar longitude (Ls) of perihelion. Surface temperature is described by the Local Insolation Temperature (LIT) scheme, which uses a daily and latitude-dependent radiation balance. The evaporation rate of water is calculated by an expression for free convection, driven by density differences between water vapor and ambient air, the condensation rate follows from the assumption that any water vapour which exceeds the local saturation pressure condenses instantly, and atmospheric transport of water vapour is approximated by instantaneous mixing. Glacial flow of ice deposits is neglected. Simulations with constant orbital parameters show that low obliquities favour deposition of ice in high latitudes and vice versa. A transient scenario driven by a computed history of orbital parameters over the last 10 million years produces essentially monotonically growing polar ice deposits during the most recent 4 million years, and a very good agreement with the observed present-day polar layered deposits. The thick polar deposits sometimes continue in thin ice deposits which extend far into the mid latitudes, which confirms the idea of "ice ages" at high obliquity.

The role of water ice in the solar system is reviewed from a fluid-dynamical point of view. On Earth and Mars, water ice forms ice sheets, ice caps and glaciers at the surface, which show glacial flow under their own weight. By contrast, water ice is a major constituent of the bulk volume of the icy satellites in the outer solar system, and ice flow can occur as thermal convection. The rheology of polycrystalline aggregates of ordinary, hexagonal ice Ih is described by a power law, different forms of which are discussed. The temperature dependence of the ice viscosity follows an Arrhenius law. Therefore, the flow of ice in a planetary environment constitutes a thermo-mechanically coupled problem; its model equations are obtained by inserting the flow law and the thermodynamic material equations in the balance laws of mass, momentum and energy. As an example of gravity-driven flow, the polar caps of Mars are discussed. For the north-polar cap, large-scale flow velocities of the order of 0.1...1 mm/a are likely, locally enhanced by a factor ten or more in the vicinity of surface scarps/troughs. By contrast, the colder south-polar cap is expected to be almost stagnant. Tidally heated convection is discussed for the example of the icy crust of Europa, where a two-dimensional model predicts the formation of an upper, conductive lid and a lower, convective layer with flow velocities of the order of 100 mm/a. Very little is known about the fluid-dynamical relevance of high-pressure phases of water ice as well as ices made up of other materials.