Atheoretical ionospheric model is developed in order to study the variation of ionospheric density profile during geomagnetic storms, considering the disturbances of thermospheric temperature, composition, wind and electric field. The results show that the increase of N
2/O is the most efficient way to decrease the electron density as well as the peak density N
mF
2, but the up and down movement of peak height (h
mF
2) is mainly caused by neutral wind or electric field disturbance. The temperature change has little effect on the profile variation. It is possible to distinguish these different physical processes using f
0F
2 and h
mmF
2 parameters in practical measurement.