In this paper, we introduce the implied volatility from Black-Scholes model and suggest a model for constructing implied volatility surfaces by using the two-dimensional cubic (bi-cubic) spline. In order to utilize a spline method, we acquire grid (knot) points. To this end, we first extract implied volatility curves weighted by trading contracts from market option data and calculate grid points from the extracted curves. At this time, we consider several conditions to avoid arbitrage opportunity. Then, we establish an implied volatility surface, making use of the two-dimensional cubic spline method with previously estimated grid points. The method is shown to satisfy several properties of the implied volatility surface (smile, skew, and flattening) as well as avoid the arbitrage opportunity caused by simple match with market data. To show the merits of our proposed method, we conduct simulations on market data of S&P500 index European options with reasonable and acceptable results.