We investigated the incorporation of hydrogen into zircon at 1650 and 1550 °C, and pressures of 2.5 and 1.5 GPa under water-saturated conditions in a piston-cylinder apparatus. Concentrations were determined by polarized Fourier transform infrared spectroscopy using the zircon absorption coefficient εi = 36 241 cm−2 per mol H2O/L and range from ~90 to 200 ppm H2O by weight. Crystals grown in the presence of Ti4+ or Th4+ do not differ significantly in their H2O content. We also synthesized zircons with various concentrations of Lu2O3 and Al2O3 to characterize changes in band positions and hydrogen concentrations related to coupled substitutions in zircon. Trivalent cations correlate in a nearly 1:1 molar fashion with hydrogen highlighting a potentially important coupled substitution in high water activity environments. Bands from undoped and doped zircons in the OH stretching region of the infrared spectrum show broad agreement when compared to spectra from natural samples. Heating experiments at 1 atm and 1000 °C produce a decrease in the integrated area; while some bands disappeared entirely, others are particularly stable with little decrease in integrated area after 128 h at 1000 °C. Results presented here help eliminate uncertainties that arose from Fourier transform infrared studies of natural zircons and provide further clarification for the origin of band positions in natural samples. In addition to the water activity of the crystallizing medium, the H2O content of natural grains will likely be significantly influenced by trivalent cation concentrations. In crustal zircons especially, trivalent atomic contents generally exceed those of phosphorus, meaning that hydrogen may be particularly important for trivalent cation charge compensation. An unanticipated result of this study was the development of a reasonably effective technique that produces relatively homogenous zircons doped with minor impurities. This technique could potentially be utilized in studies aimed at developing zircon standards, because it yields crystals that appear to be more homogenous than those produced by the flux method, and are generally free of inclusions.