Wind turbine power curve (WTPC) serves as an important tool for wind turbine condition monitoring and wind power forecasting. Due to complex environmental factors and technical issues of the wind turbines, there are many outliers and inconsistencies present in the recorded data, which cannot be removed through any pre-processing technique. However, the current WTPC models have limited ability to understand such complex relation between wind speed and wind power and have limited non-linear fitting ability, which limit their modelling accuracy. In this paper, the accuracy of the WTPC models is improved in two ways: first is by developing multivariate models and second is by proposing MARS as WTPC modeling technique. MARS is a regression-based flexible modeling technique that automatically models complex the nonlinearities in the data using spline functions. Experimental results show that by incorporating additional inputs the accuracy of the power curve estimation is significantly improved. Also by studying the error distribution it is proved that multivariate models successfully mitigate the adverse effect of hidden outliers, as their distribution has higher peaks and lesser standard deviation, which proves that the errors, are more converged to zero compared to the univariate models. Additionally, MARS with its superior non-linear fitting ability outperforms the compared methods in terms of the error metrics and ranks higher than regression trees and several other popular parametric and non-parametric methods. Finally, an outlier detection method is developed to remove the hidden outliers from the data using the error distribution of the modeled power curves.