Certainly not going to argue with the theory, I'm too early in a physics program to dispute any of the ideas.
However, from my Thermal Physics class, a hard drive appears to be a large system of paramagnets. The individual bits can either be 0, 1, or in terms of the little magnets, up or down. This is a 2 state system, known as a paramagnet.
Now, from the book we use, the magnetization of the system is given by (THG sucks for equations):
M = N (mu) * tanh((mu)B/(kT))
where N is the number of paramagnets, mu is the magnet dipole moment (of a single paragamagnet), B is the external magnetic field applied (writing arm, for example), k is Boltzmann's constant, and T is the temperature.
Now, if you're not familiar with the hyperbolic tangent function, it has a range of -1 to 1. As the argument goes to infinity, the tanh function reaches a value of 1. As the arguemet goes to 0, though, the tanh function goes to 0.
Now, as T increases, the argument will be lower, so M is lower. Then when the data is written and it cools back down, T is lower so the argument is higher and M is higher.
Perhaps someone else can answer that.