Flooding of the membrane electrode assembly (MEA) and dehydrating of the polymer electrolyte membrane have been the key problems to be solved for polymer electrolyte membrane fuel cells (PEMFCs). So far, almost no papers published have focused on studies of the liquid water flux through differently structured gas diffusion layers (GDLs). For gas diffusion layers including structures of uniform porosity, changes in porosity (GDL with microporous layer (MPL)) and gradient change porosity, using a one-dimensional model, the liquid saturation distribution is analyzed based on the assumption of a fixed liquid water flux through the GDL. And then the liquid water flux through the GDL is calculated based on the assumption of a fixed liquid saturation difference between the interfaces of the catalyst layer/GDL and the GDL/gas channel. Our results show that under steady-state conditions, the liquid water flux through the GDL increases as contact angle and porosity increase and as the GDL thickness decreases. When a MPL is placed between the catalyst layer and the GDL, the liquid saturation is redistributed across the MPL and GDL. This improves the liquid water draining performance. The liquid water flux through the GDL increases as the MPL porosity increases and the MPL thickness decreases. When the total thickness of the GDL and MPL is kept constant and when the MPL is thinned to 3 μm, the liquid water flux increases considerably, i.e. flooding of MEA is difficult. A GDL with a gradient of porosity is more favorable for liquid water discharge from catalyst layer into the gas channel; for the GDLs with the same equivalent porosity, the larger the gradient is, the more easily the liquid water is discharged. Of the computed cases, a GDL with a linear porosity 0.4x + 0.4 is the best.