We show that the geometric lifting of the RSK correspondence introduced by A.N. Kirillov (2001) is volume preserving with respect to a natural product measure on its domain, and that the integrand in Givental's integral formula for GL(n,R)-Whittaker functions arises naturally in this context. Apart from providing further evidence that Whittaker functions are the natural analogue of Schur polynomials in this setting, our results also provide a new `combinatorial' framework for the study of random polymers. When the input matrix consists of random inverse gamma distributed weights, the probability distribution of a polymer partition function constructed from these weights can be written down explicitly in terms of Whittaker functions. Next we restrict the geometric RSK mapping to symmetric matrices and show that the volume preserving property continues to hold. We determine the probability law of the polymer partition function with inverse gamma weights that are constrained to be symmetric about the main diagonal, with an additional factor on the main diagonal. The third combinatorial mapping studied is a variant of the geometric RSK mapping for triangular arrays, which is again showed to be volume preserving. This leads to a formula for the probability distribution of a polymer model whose paths are constrained to stay below the diagonal. We also show that the analogues of the Cauchy-Littlewood identity in the setting of this paper are equivalent to a collection of Whittaker integral identities conjectured by Bump (1989) and Bump and Friedberg (1990) and proved by Stade (2001, 2002). Our approach leads to new `combinatorial' proofs and generalizations of these identities, with some restrictions on the parameters.