Earth system models (ESMs) integrate interactions of atmosphere, ocean, land, ice, and biosphere to estimate the state of regional and global climate under a variety of conditions. Past climate field reconstructions with deterministic ESMs through the assimilation of climate proxies need to consider the required high computations and model non-linearity. Our tests indicate that iterative schemes based on the Kalman filter and careful sensitivity analysis are adequate for approaching the problem.