Water vapor vertical profiles are important in numerical weather prediction, moisture transport, and vertical flux calculation. This study presents the Multi-Axis Differential Optical Absorption Spectroscopy (MAX-DOAS) retrieval algorithm for water vapor vertical profiles and the retrieved results are validated with corresponding independent datasets under clear sky. The retrieved Vertical Column Densities (VCDs) and surface concentrations are validated with the Aerosol Robotic Network (AERONET) and National Climatic Data Centre (NCDC) datasets, achieving good correlation coefficients (R) of 0.922 and 0.876, respectively. The retrieved vertical profiles agree well with weekly balloon-borne radiosonde measurements. Furthermore, the retrieved water vapor concentrations at different altitudes (100–2000 m) are validated with the corresponding European Centre for Medium-range Weather Forecasts (ECMWF) ERA-interim datasets, achieving a correlation coefficient (R) varying from 0.695 to 0.857. The total error budgets for the surface concentrations and VCDs are 31% and 38%, respectively. Finally, the retrieval performance of the MAX-DOAS algorithm under different aerosol loads is evaluated. High aerosol loads obstruct the retrieval of surface concentrations and VCDs, with surface concentrations more liable to severe interference from such aerosol loads. To summarize, the feasibility of detecting water vapor profiles using MAX-DOAS under clear sky is confirmed in this work.