A multi-dimensional, non-isothermal, two-phase membrane electrode assembly (MEA) numerical model is developed where the micro-structure of the porous layers is characterized by a mixed wettability pore size distribution (PSD). The PSD model is used to predict local water saturation based on gas and liquid pressures, and can be used to study the effect of varying pore size and wettability. The MEA model accounts for gas transport via molecular and Knudsen diffusion, liquid water transport, sorbed water transport by back-diffusion, electro- and thermo-osmosis, and heat generation and transport. Multi-step kinetic models are used to predict anode and cathode electrochemical reactions. Local transport losses are accounted for using a local transport resistance. The PSD model is used to predict capillary pressure vs. saturation and saturation vs. relative liquid permeability curves based on PSDs from several GDLs obtained using mercury intrusion porosimetry. The PSD-based MEA model electrochemical performance predictions are also compared to experimental data from the literature. Results show that the numerical model is able to capture the performance changes associated with varying temperate and the introduction of a micro-porous layer. (C) The Author(s) 2017. Published by ECS. All rights reserved.