Past climate conditions are reconstructed using indirect and incomplete geological, biological, and geochemical proxy data. We propose that such reconstructions are best obtained by statistical inversion of hierarchical models that represent how multi–proxy observations and calibration data are produced by variation of environmental conditions in time and/or space. These methods extract new information from traditional proxies and provide robust, comprehensive estimates of uncertainty.