Recent magnetic satellite missions facilitate new birth of large-scale geomagnetic field models and their applications to tectonics. Here, we directly compare the global geomagnetic field models NGDC-720 with the tectonics and seismicity in Mainland China and surroundings. It is found that the tectonics and seismicity in this area show remarkable correlation with the vertical component of lithospheric magnetic field (B z) calculated at an altitude of 200 km. Previous thought was that earthquakes are more likely to occur in zero B z belts or in obvious anomaly gradient belts. On the contrary, we find that more than half (53.2%) of the earthquakes occurred in areas with B z of − 5 to − 3 nT or in areas with a relatively small horizontal gradient of B z in the same time interval with the satellite data. The percentage seismic energy in these areas (− 5 nT < B z < − 3 nT) is even as high as 94.6%. We explain this unexpected result with a two equivalent source dipole model, arguing that the viscosity difference caused by the temperature gradient within the lithosphere likely accounts for the correlation between magnetic anomalies and seismicity.