This paper proposed a hybrid intelligent process model, based on finite-element method (FEM) and Gaussian process regression (GPR), for electrical discharge machining (EDM) process. A model of single-spark EDM process has been constructed based on FEM method, considering the latent heat, variable heat distribution coefficient of cathode (f c ), and plasma flushing efficiency (PFE), to predict material removal rate (MRR) and surface roughness (Ra). This model was validated using reported analytical and experimental results. Then, a GPR model was proposed to establish relationship between input process parameters (pulse current, pulse duration, and discharge voltage) and the process responses (MRR and Ra) for EDM process. The GPR model was trained, tested, and tuned using the