## 1 Introduction  
The lithospheric magnetic field accounts for approximately 4% of the total energy of the earth's magnetic field [1]. It originates from the crust and upper mantle of magnetic rocks and reflects a 3D spatial distribution of rock magnetism. Given its different rock magnetic environments, temperature pressure conditions, magnetic carriers, and tectonic evolutions, the lithospheric magnetic field carries abundant and complex geological structures and tectonic movement information [2, 3]. Therefore, analyzing the magnetic distribution frame in the depth of the crust and the top of the upper mantle using a high-resolution and high-accuracy lithospheric magnetic field model can reveal important information regarding regional tectonics and their evolution, the division of rock masses, and the heat flow of the lithosphere [4, 5, 6, 7, 8]. The accuracy of lithospheric magnetic field models and spatial resolutions improved significantly with the recent accumulation of Oersted, CHAMP, Swarm, and other satellite magnetic data, as well as the continuous improvement of inversion algorithms and separation technology of the various source fields [9]. A series of representative lithospheric magnetic field models was successively released internationally. The minimum spatial wavelength reaches 130 km when satellite observations are only used for the recently released HDGM-2016 model ([[http://www.ngdc.noaa.gov/geomag/HDGM/index.html](http://www.ngdc.noaa.gov/geomag/HDGM/index.html)]([http://www.ngdc.noaa.gov/geomag/HDGM/index.html](http://www.ngdc.noaa.gov/geomag/HDGM/index.html))). The minimum spatial wavelength is approximately 300 km when the spherical harmonic function expansion is 133 degrees for the POMME-10 model ([[http://geomag.org/models/pomme10.html](http://geomag.org/models/pomme10.html)]([http://geomag.org/models/pomme10.html](http://geomag.org/models/pomme10.html))). Unlike ground and aero magnetic surveys, satellite magnetic survey is characterized by high coverage and repetition rate; this approach can be used to observe the change of the earth's magnetic field in real time [10; 11]. The altitude of satellite observation data is relatively high. Thus, observation data can effectively suppress the interference of surface environments and shallow geological bodies. Furthermore, the medium/long-wavelength signal of the lithospheric magnetic field was extracted and used in studies on tectonic structures in the deep crust and the top of the upper mantle, as well as the depth of the Curie surface and heat flow anomaly [12; 13; 14; 15]. [PERSON] et al. calculated the magnetization intensity and the thickness of a global large-scale lithospheric magnetic field using a revised spherical cap harmonic model to analyze the spatial power spectrum of the lithosphere magnetic field [16]. [PERSON] et al. estimated the depth of magnetic sources by analyzing the radial power spectrum of aeromagnetic anomaly by assuming that magnetization has a fractal distribution; a fractal crustal magnetization model was used to map the depth of the Curie. An obvious corresponding relationship between the characteristics of the bottom of magnetic sources and the prominent heat flow anomalies was identified. The study provided new information regarding geological structures and heat flow anomalies [17]. [PERSON] et al. (2013) found that the Z component of the satellite magnetic anomalies was consistent with the geological features of the old craton in China. Moreover, the relationship between the geological detail features and the radial component of magnetic anomalies was thoroughly explained by integrating ground magnetic survey data [18]. [PERSON] used satellite magnetic data to estimate the heat flow underneath the Antarctic ice sheet. With the enhancement of the resolution of the satellite magnetic field models, the method is believed to be capable of identifying the detailed features of heat flow and may be utilized to eventually discover volcanoes that have not yet been discovered beneath the glacier [12].