  Statistical fluctuations in the transverse distribution of sources in
relativistic heavy ion collisions and an asymmetric emission profile associated
with the wounded nucleons lead to rapidity dependence of the reaction plane.
The size of this effect is estimated for the gold-gold collisions at the
highest RHIC energy ($\sqrt{s_{NN}}=200 \mathrm{GeV}$) in the Glauber model
using GLISSANDO. The hydrodynamical evolution of the resulting torqued fireball
is carried out in 3+1 dimensions considering a perfect fluid. Hadronization is
simulated using THERMINATOR including non-flow contribution coming from
resonance decays. Some experimental measures that can be used to detect the
torque effect are proposed.
