The 2010 Chile Maule earthquake ranks as the sixth largest earthquake ever recorded by a seismograph. In this paper, based on the parallel finite element method, the heterogeneous ellipsoidal earth model including topography, Moho discontinuity is setup to investigate the co seismic effect resulted from the 2010 Maule earthquake. Based on the Coulomb stress change, seismic activity of surrounding faults and triggering of aftershocks are estimated. The results show that: the effects of Earth's layered structure and curvature on co-seismic horizontal and vertical displacement are significant in the case of 2010 Maule earthquake. The plane model may underestimate the far-field co-seismic displacement and strain near the equator, overestimate and underestimate the far-field displacement and strain near Antarctica, respectively. If the effects of laterally heterogeneous medium and topography are neglected, co-seismic stress drop will be significantly underestimated, the error of the near-field east-westward and north-southward stress will be up to 95. 4% and 90. 8%, respectively. The Coulomb stress change in Illapel region is equivalent to 15 years' tectonic stress accumulation, which means the Maule earthquake hasten the failure of 2015 Illapel earthquake. Coulomb stress changes on optimally orientated faults are estimated, 70. 9% of aftershocks occur in regions where the failure stress exceeds 10 kPa. The MW 6. 9 and M-W 7. 0 aftershocks near Pichilemu were triggered by the great Maule earthquake.
