Atmospheric water vapor plays a prominent role in climate change and atmospheric, meteorological, and hydrological processes. Because of its high spatiotemporal variability, precise quantification of water vapor is challenging. This study investigates Integrated Water Vapor (IWV) variability for the period 1995–2010 at 118 globally distributed Global Positioning System (GPS) sites, using additional UV/VIS satellite retrievals by GOME, SCIAMACHY, and GOME-2 (denoted as GOMESCIA below), plus ERA-Interim reanalysis output. Apart from spatial representativeness differences, particularly at coastal and island sites, all three IWV datasets correlate well with the lowest mean correlation coefficient of 0.878 (averaged over all the sites) between GPS and GOMESCIA. We confirm the dominance of standard lognormal distribution of the IWV time series, which can be explained by the combination of a lower mode (dry season characterized by a standard lognormal distribution with a low median value) and an upper mode (wet season characterized by a reverse lognormal distribution with high median value) in European, Western American, and subtropical sites. Despite the relatively short length of the time series, we found a good consistency in the sign of the continental IWV trends, not only between the different datasets, but also compared to temperature and precipitation trends.