The stoichiometric ratio and flow channel geometry play a vital role in the performance of high temperature proton exchange membrane (HT-PEM) fuel cells. Because of the high cost of experiments or simulations, most analyses and optimization of the stoichiometric ratio and flow channel geometry are limited to several points in the entire design domain. In this study, an analysis and optimization method for HT-PEM fuel cells based on the surrogate model was proposed. Surrogate models were constructed using some of the available budgets of samples to analyze and optimize the entire design domain. With this method, it was indicated that the effect of the cathode stoichiometric ratio is more significant to the cell performance than the anode stoichiometric ratio and there are significant nonlinear interactions among the flow channel geometry parameters. At the fixed operating voltage, the flow channel geometry with the maximum current density and that with the maximum real power were obtained. Compared with the base design, the designs obtained by the surrogate model improve the current density and real power by 10.54% and 3.93%, respectively. Thus, this analysis and optimization method is demonstrated to be helpful and deserves attention in future research. (C) 2020 Hydrogen Energy Publications LLC. Published by Elsevier Ltd. All rights reserved.