## 5 Conclusions  
Zonal revised of spherical cap harmonic analysis(R-SCHA) was used to model the lithospheric magnetic field of the Tibetan Plateau at satellite altitude. The accuracy and stability of this model were tested and evaluated in the horizontal and radial directions. Furthermore, the relationship between magnetic anomalies and geological structures was analyzed. One of the key factors for modeling is refining satellite magnetic data. A complete data preprocessing procedure was designed, which involved data screening, model correction, track filtering, and hierarchical gridding. Reliability is another important factor of a model. Thus, the model must not generate significant shock effects in the horizontal and vertical directions but still achieve a certain spatial resolution. Considering that the difference between the latitude and longitude over the Tibetan Plateau is large, a small cap partition modeling can be used in a particular truncation level to provide good spatial resolution while ensuring accuracy. Therefore, the Tibetan Plateau was divided into two independent parts and the half-angle of the cap (approximately 1\({}^{\circ}\)) was moderately amplified to eliminate the influence of the boundary effect. Meanwhile, hierarchical gridding was conducted for satellite data distributed at a height of 250-510 km to guarantee the balance of each layer of data. The small semi angle of cap improved the fitting accuracy of the model throughout the region, but the difficulty of this method lies at the joint of each subregion model. The present paper used iterative inversion to create a method for seamless splicing between adjacent spherical caps. Moreover, model error and reliability analyses played important roles in the modeling process. The difference between the model value and the satellite observation value were initially calculated before the RMS of the model was quantitatively analyzed at the heights of the CHAMP and Swarm satellites. The fitting error of the model is small, and the RMS of each geomagnetic component is better than 0.3 nT. Moreover, the stability of the model in horizontal and vertical directions was examined by introducing the EMM2015 (16-720 order) high-accuracy lithospheric magnetic field model. The model does not have any significant shock effects in the horizontal direction. The model can provide a stable space continuation in the radial direction. The analysis indicates that the R-SCHA algorithm can accurately describe the middle- to long-scale lithospheric magnetic field in the Tibetan Plateau at the satellite altitude layer. The separation of internal and external field sources remains a problem in the modeling process. The R-SCHA algorithm does not consider the joint estimation of the external source field. Instead, it facilitates the separation of an external field source through the model and filtering method in the data preprocessing steps. This method may introduce new errors in the process of model fitting. Therefore, Combining Swarm constellation differential gradient measurement is necessary to effectively suppress the interference of the long wavelength signal that remains in the external field and the main magnetic field. However, a map of a high-resolution lithospheric magnetic field over the Tibetan Plateau cannot be developed because of the lack of support for ground magnetic survey and aeromagnetic data. Although a good relationship between the distribution of magnetic anomalies and deep crustal tectonic structures at 300 km satellite altitude was identified, the model fails to extract abundant shallow geological structures information. Future research will be conducted on a downward continuation algorithm for geomagnetic anomaly data and multi-scale modeling will be integrated to obtain a high-resolution and high-accuracy lithospheric magnetic field. This endeavor will provide technical support for the study on the crustal structures of the Tibetan Plateau and related dynamic mechanisms.