We improved the surface spline method to make the centers of the base functions fixed. This improvement makes it meaningful to analyze the time series of coefficients. Based on the revised surface spline (RSS) method and 2382 geomagnetic data measured in mainland China, we built Chinese Geomagnetic Reference Field 2020 (CGRF2020) to describe the spatial distribution of the geomagnetic field in most of China and surrounding areas. The root mean square residual is 0.1318 degrees for D, 0.0857 degrees for I, 93.70 nT for F, 73.94 nT for X, 65.02 nT for Y, and 97.82 nT for Z. The values of F, X, Y, and Z are about 20 nT less than those by the SCH method and IGRF13. We also preliminarily analyzed the correlation between geomagnetic anomalies and geological structures.