Abstract：This paper processes three issues of annual variation data of lithospheric magnetic field from 2016 to 2019 in central North China. Then we grid the annual change data of ΔY element in each period，and use the node value to calculate and assign a value to the ΔY gradient space. We analyze the gradient spatial distribution of ΔY. The results show that there is a certain correspondence between high value areas and active faults. High value areas are mostly extended along one or more active faults or distributed around active faults. By using the polynomial regression method to decompose the magnetic field of lithosphere in certain scales to obtain regional background anomalies and local anomalies, the results show that there is a good correspondence between the distribution of local anomalies and active faults. There are great differences in the morphology of gradient high value areas and local anomalies in each period, but there are also some high value areas or anomaly significant areas with continuity, especially near the faults, which may be more closely related to tectonic stress.