Abstract One of the precise widely used global Zenith Hydrostatic Delay (ZHD) model is based on the gravity value at the centroid of the atmospheric column at the station of observation and gravity value at the centroid is constant in this model for a specific location throughout the year. However, as the content and extent of atmosphere varies temporally, its centroid and consequently gravity value at the centroid also varies. Apart from this, the actual atmospheric condition of different region is not alike. Therefore, there is a need to develop a regional mean gravity model and development of such model has been discussed in this paper. To obtain the mean gravity model, first a regional model of centroid height of atmospheric column was developed as a function of the surface pressure and temperature. It was developed by multiple regressions between estimated centroid of the atmosphere and surface pressure, surface temperature using radiosonde data of five radiosonde stations spread over the Indian subcontinent. Three years radiosonde data from 2006 to 2008 was used for each station. The root mean square error in estimating centroid of the atmospheric column is about ±326 m, which is negligible considering the variability of the atmosphere and its content. The centroid height model has been used to formulate the mean gravity model, considering uniform lapse rate in gravity with height. It is found that proposed mean gravity model provides temporal variation of mean gravity values at the centroid and thus matches with the reality. The interesting advantage of the developed model is that the model shows diurnal variation of mean gravity. The accuracy of ZHD has shown of the order of about 0.3mm using the developed regional mean gravity model. However, already developed ZHD model has shown a slight inferior result compared to the developed model. These models have shown accuracy of about 0.8mm and 0.6mm.