Capillary pressure vs. saturation (P-C(S-L)) Curves are fundamental to understanding liquid water transport and flooding in PEM gas diffusion layers (GDLs). P-C(S-L) curves convolute the influence of GDL pore geometry and internal contact angles at the three-phase liquid/solid/gas boundary. Even simple GDL materials are a spatially non-uniform mixture of carbon fiber and binder, making a Gaussian distribution of contact angles likely, based on the Cassie-Baxter equation. For a given Gaussian contact angle distribution with mean (theta(Mean)) and standard deviation (sigma), a realistic P-C(S-L) curve can be computed using a bundle of capillaries model and GDL pore size distribution data. As expected, computed P-C(S-L) curves show that theta(Mean) sets the overall hydrophilic (theta(Mean) < 90 degrees) or hydrophobic (theta(Mean) > 90 degrees) character of the GDL (i.e., liquid saturation level at a given capillary pressure), and sigma affects the slope of the P-C(S-L) curve. The capillary bundle model also can be used with (theta(Mean), sigma) as unknown parameters that are best-fit to experimentally acquired P-C(S-L) and pore size distribution data to find (theta(Mean), sigma) values for actual GDL materials. To test this, pore size distribution data was acquired for Toray TGP-H-090 along with hysteretic liquid and gas intrusion capillary pressure curve data. High quality best-fits were found between the model and combined datasets, with GDL liquid intrusion showing fairly neutral internal surface wetting properties (theta(Mean) = 92 degrees and (sigma = 10 degrees) whereas gas intrusion displayed a hydrophilic character (theta(Mean) = 52 degrees and sigma = 8 degrees). External liquid advancing and receding contact angles were also measured on this same material and they also showed major hysteresis. The new methods described here open the door for better understanding of the link between GDL material processing and the wetting properties that affect flooding. (C) 2008 Elsevier B.V. All rights reserved.