A series of hydrothermal diamond anvil cell experiments was conducted to constrain the equilibrium distribution of molecular H2 between H2O-saturated sodium aluminosilicate melts and H2O at elevated temperatures (600–800 °C) and pressures (317–1265 MPa). The distribution of H2 between the silicate liquid and the aqueous fluid was achieved through real-time monitoring of the H-H stretching vibration under in situ conditions using Raman vibrational spectroscopy. Results show that the solubility of H2 in silicate melts saturated with H2O decreases as the temperature increases, with control exerted by the mole fraction of H2O in the melt. The dissolution of H2 in the hydrous silicate melts appears to follow Henrian behavior, resembling that of an inert, neutral non-polar species. To express species solubility as a function of temperature (T in K) an empirical equation was developed:ln Km/f = 11.4 (±1.3) *1000 / T (K) − 18.7 (±1.1)where Km/f is the equilibrium constant for the reaction H2(g) = H2(melt). This equation was derived by integrating data from the current and prior experimental studies that include silicate melts with varying H2O saturation levels. It should be deemed applicable within the temperature range of 600–1450 °C and pressures ranging from 0.3 to 3 GPa. The implications are extended into developing an understanding of the H partitioning between H2-rich atmospheres blanketing magma oceans in the early history of planetary bodies. For example, transferring H from primordial atmospheric envelopes to the interior of rocky exoplanets may be less efficient than previously believed, which should be considered in models of volatile retention. Experimental data also suggest that minimal amounts of solar nebula H2 are likely to dissolve in the molten surface of primitive objects in the protoplanetary disk (∼10−5 to 10−9 mol faction H2 in melt), contradicting the highly reducing conditions observed in chondrule mineral compositions.