A two-dimensional mathematical model integrating fluid flow and solutes mass transport with a particle-based approach for crystal nucleation and growth was developed to study gypsum (CaSO4 center dot 2H(2)O) scaling in membrane feed channels. The model is able to describe the generally observed trends for performance decline due to scaling: permeate flux deterioration when constant pressure is applied and trans-membrane pressure increase when a constant permeate flux is set. Model results for precipitation in the feed channel show a gradual increase of precipitated gypsum in the axial direction. The zones around spacer-membrane contacts are most prone to scaling even when the feed water is under-saturated. Simulations highlight the importance of both thermodynamic and kinetic effects on gypsum precipitation for various salinities. Despite having the same degree of saturation of the feed water, variable feed stream compositions result in different amounts of precipitate formed. The numerical results suggest that the degree of saturation of feed water (DSin) may not be a reliable indicator for precipitate formation in membrane feed channels. Results support the need for considering the dynamics of local flow and solute patterns, combined with ionic activity calculations for understanding mineral precipitation in spiral-wound membrane modules used in desalination. (C) 2013 Elsevier B.V. All rights reserved.