Improving the accuracy of numerical weather predictions remains a challenging task. The absence of sufficiently detailed temporal and spatial real-time in-situ measurements poses a critical gap regarding the proper representation of atmospheric moisture fields, such as water vapor distribution, which are highly imperative for improving weather predictions accuracy. The estimated amount of the total vertically integrated water vapor (IWV), which can be derived from the attenuation of global positioning systems (GPS) signals, can support various atmospheric models at global, regional, and local scales. Currently, several existing atmospheric numerical models can estimate the IWV amount. However, they do not provide accurate results compared with in-situ measurements such as radiosondes. Here, we present a new strategy for assimilating 2D IWV regional maps estimations, derived from combined GPS and METEOSAT satellite imagery data, to improve Weather Research and Forecast (WRF) model predictions accuracy in Israel and surrounding areas. As opposed to previous studies, which used point measurements of IWV in the assimilation procedure, in the current study, we assimilate quasi-continuous 2D GPS IWV maps, combined with METEOSAT-11 data. Using the suggested methodology, our results indicate an improvement of more than 30% in the root mean square error (RMSE) of WRF forecasts after assimilation relative standalone WRF, when both are compared to the radiosonde measured data near the Mediterranean coast. Moreover, significant improvements along the Jordan Rift Valley and Dead Sea Valley areas are obtained when compared to 2D IWV regional maps estimations. Improvements in these areas suggest the impact of the assimilated high resolution IWV maps, with initialization times which coincide with the Mediterranean Sea Breeze propagation from the coastline to highland stations, as the distance to the Mediterranean Sea shore, along with other features, dictates its arrival times.