A local Gaussian process regression approach is presented, which allows to model nonlinearities of internal combustion engines more accurate than global Gaussian process regression. By building smaller models, the prediction of local system behavior improves significantly. In order to predict a value, the algorithm chooses the nearest training points. The number of chosen training points depends on the intensity of estimated nonlinearity. After determining the training points, a model is built, the prediction performed and the model discarded. The approach is demonstrated with a benchmark system and air charge test bed measurements. The measurements are taken from a turbocharged SI gasoline engine with both variable inlet valve lift and variable inlet and exhaust valve opening angle. The results show how local Gaussian process regression outmatches global Gaussian process regression concerning model quality and nonlinearities in particular.