Abstract: According to the typical formula of group optical atmospheric refractivity and the empirical formula of water vapor density in ITU-R recommendation, a new formula for calculating the optical atmospheric refractivity is given. Besides, by referring to the hydrostatic component of Hopfield refractivity model and the empirical model of water vapor density in ITU-R recommendation, a method of modeling optical atmospheric refractivity profile model is developed by using the historical radiosonde data. Take Qingdao area as an example, a local optical atmospheric refractivity profile model is established, and the key parameters in the model are derived from the local historical radiosonde data from 1986 to 1995 and the reference standard atmospheres. The comparison shows that the established model is more consistent with the radiosonde data than those calculated by the general Hopfield model. This method of modeling optical atmospheric refractivity profile will be beneficial to improve the accuracy of optical atmospheric refraction error correction in the data processing of trajectory measurement.