We calculated the ground deformation due to differently shaped rotational ellipsoidal pressure sources beneath the Aira caldera using a FEM model. Deeper horizontal pressure sources seem to be consistent with the observed data, as well as that expected from the Mogi pressure source. However, the variation of residual sum of square is small to decide a certain oblateness value of horizontal pressure sources.