High pressure subsea natural gas dehydration (NGD) units using ethylene glycol (MEG) absorption have been proposed. To expand the experimental database and assist design qualification, new vapor-liquid equilibrium (VLE) experimental data have been measured for a 20-component glycol-water-natural gas mixture at T = (288-323) K, p = (6.0, 12.5) MPa, and w(MEG,feed) = (90, >99.8) %. MEG, H2O, CO2, N-2, and alkane (methane to n- and i-pentane) phase distributions have been quantified. Experimental uncertainty ranges from +/- 2-42%, with the greatest uncertainty for the quantification of trace components. Experimental results are modeled using the Cubic-Plus-Association (CPA) equation of state. Overpredictions (similar to 9%) are observed for the water content of the vapor phase. CO2 is shown to have a large effect on y(MEG), leading to modeling deviations in the order of 65%. A relatively accurate prediction of the natural gas partition coefficients was observed for major components C-1-C-3 and CO2, with modeling errors ranging from 5% for methane to 10% for CO2. More significant deviations were observed for trace components, with the largest deviation of 73% N-2. The CPA model provides both satisfactory and conservative results suitable for use in NGD process designs. On the basis of this work, operation at subsea conditions would significantly improve dehydration capability.