This paper presents an input–output reduced-order aerodynamic model employing the kriging interpolation for laminar flow. The kriging interpolation is a regression method, particularly useful for large motions. In the targeted flow regime, the response includes unsteady vortex shedding at high angles of attack. The reduced-order model (ROM) predicts the aerodynamic response of a rigid airfoil to a () gust under a wide range of angles of attack, gust lengths, and gust intensities. We study two algorithms for sampling training simulations, equally spaced and Latin hypercube methods. Both methods result in a ROM with an accurate prediction of aerodynamic forces and moments, given sufficient training points, particularly along the gust length axis. We use this ROM in coupling with the structural dynamics equations to study the aeroelastic response of a one-degree-of-freedom elastic airfoil. The aeroelastic gust response predictions are accurate for moderately large pitch angle deformation, using only the ROM created from computational fluid dynamics simulation of the rigid airfoil.