Abstract. The distribution of pressure on the vertical seaward front of an ice shelf has been shown to cause downward bending of the shelf if the ice is assumed to have vertically uniform viscosity. Satellite lidar observations show that many shelf edges bend upward and that the amplitude of upward deflections depends systematically on ice shelf thickness. A simple analysis is presented showing that upward bending of shelf edges can result from vertical variations in ice viscosity that are consistent with field observations and laboratory measurements. Resultant vertical variations in horizontal stress produce an internal bending moment that can counter the bending moment due to the shelf-front water pressure. Assuming a linear profile of ice temperature with depth and an Arrhenius relation between temperature and strain rate allows derivation of an analytic expression for internal bending moments as a function of shelf surface temperatures, shelf thickness and ice rheologic parameters. The effect of a power-law relation between stress difference and strain rate can also be included analytically. The key ice rheologic parameter affecting shelf edge bending is the ratio of the activation energy, Q, and the power-law exponent, n. For cold ice surface temperatures and large values of Q/n, upward bending is expected, while for warm surface temperatures and small values of Q/n downward bending is expected. The amplitude of bending should scale with the ice shelf thickness to the power 3/2, and this is approximately consistent with a recent analysis of shelf edge deflections for the Ross Ice Shelf. These scaling relations should help guide fully two-dimensional numerical simulations of shelf bending.