Surface observations and CHAMP measurement data are employed to develop a three-dimensional surface spline (3DSS) model of mainland China. The magnetic field distribution at the satellite level is then demonstrated using the model obtained. The results of this model are compared and verified by deriving the corresponding two (2DTY) and three-dimensional (3DTY) Taylor polynomial models. Issues such as the removal of disruptive geomagnetic fields, the data gap between the surface and satellite levels, and boundary effects are carefully considered during modeling. We then focus on evaluating the modeling effect of the satellite data. Ten satellite points not involved in the modeling procedure are selected, and the residuals, absolute change rates, and RMSEs of these points are calculated. Results show that the distribution of the magnetic field determined by the 3DSS model is highly consistent with that obtained from the IGRF12 model. Expect for component Y, the absolute change rates of other components are less than 0.5%. Specifically, the RMSE of Y of 3DSS is nearly 60% lower than those of 3DTY and 2DTY; the RMSE of other components of the former are also over 90% lower than those of the latter. This finding implies that the 3DSS model has good performance for modeling satellite data and its results are reliable. Moreover, the modeling effect of 3DTY is better than that of 2DTY.