The catalytic hydrogenation of aqueous bicarbonate was modelled utilizing a central composite experimental design. The results were compared with a dataset derived from a measurement procedure where univariate evaluations of the possible reaction conditions were applied. The results show good correlation between the measured and the modelled data. The optimum conditions to produce the highest possible formate concentration were identified by response surface methodology.