Kinetic Monte Carlo (KMC) simulations are carried out to directly study diffusion of benzene through thin (37-100 nm) NaX zeolite membranes under a gradient in chemical potential. Nonlinearities in adsorbate loading near the membrane boundaries are shown to arise from the difference in adsorbate density between the zeolite and adjacent fluid phase. Direct extraction of the transport diffusivity from gradient KMC simulations enables testing of the Darken approximation. This rigorous approach reveals limitations of the Darken approximation and, for the first time, the potentially complex nonunique functionality and multiplicity of the transport diffusivity for strongly interacting adsorbates. In the companion paper we explore these nonlinear interfacial effects in the context of permeation through both single-crystal and polycrystalline membranes.