The first 2000 m of the global thermohaline structure of the ocean are statistically decomposed into vertical thermohaline modes, using a multivariate functional principal component analysis (FPCA). This method is applied on the Monthly Isopycnal and Mixed-Layer Ocean Climatology (MIMOC). The first three modes account for 92% of the joint temperature and salinity (T-S) variance, which yields a surprisingly good reduction of dimensionality. The first mode (69% of the variance) is related to the thermocline depth and delineates the subtropical gyres. The second mode (18%) is mostly driven by salinity and mainly displays the asymmetry between the North Pacific and Atlantic basins and the salty circumpolar deep waters in the Southern Ocean. The third mode (5%) identifies the low-and high-salinity intermediate waters, covarying with the freshwater inputs of the upper ocean. The representation of the ocean in the space defined by the first three modes offers a simple visualization of the global thermohaline structure that strikingly emphasizes the role of the Southern Ocean in linking and distributing water masses to the other basins. The vertical thermohaline modes offer a convenient framework for model and observation data comparison. This is illustrated by projecting the repeated Pacific section P16 together with profiles from the Array for Real-Time Geostrophic Oceanography (ARGO) global array of profiling floats on the modes defined with the climatology MIMOC. These thermo-haline modes have a potential for water mass identification and robust analysis of heat and salt content.