Fast Krylov subspace methods for geostatistical inverse problems