A sparse grid method is developed to build a surrogate model of the RZWQM2. The global optimization algorithm QPSO is used for parameter estimation. Calibrating the surrogate gives comparable results as the calibration of the RZWQM2. The surrogate global optimization greatly improves computational efficiency.