Stress changes due to earthquake rupture can disturb geomagnetic fields significantly. In order to investigate the impact of the 2021 MW 7.3 Maduo earthquake on geomagnetic fields, a piezomagnetic model is constructed based on the coseismic slip to calculate the static coseismic piezomagnetic anomalies (PMs). The PMs are considerable in near-field. However, the PMs are negligible in regions tens of kilometers from the fault rupture. The PMs of our model are consistent with those of other strike-slip earthquakes, indicating that our piezomagnetic model is reasonable. The east component of observed coseismic geomagnetic changes and calculated PMs on a geomagnetic repeat station located about 6 km from fault trace are +4.8 ± 2.2 nanotesla and +4.3 nanotesla, respectively. It seems that the piezomagnetic model can explain the observed data. The PMs are up to 10 nanotesla in the near-field with the initial magnetization of 3 A/m and stress sensitivity of 2 × 10−3 MPa−1. Consequently, considerable coseismic geomagnetic changes that are above error could be observed along the fault, especially at locations with geometrical complexities.