Multistage membrane distillation (MD) systems can have significantly higher efficiencies than their single stage counterparts. However, multistage MD system design has received limited attention. In this paper, the performance of a multistage vacuum membrane distillation (MSVMD) which is thermodynamically similar to a multi-stage flash distillation (MSF) is evaluated for desalination, brine concentration, and produced water reclamation applications. A wide range of solution concentrations were accurately modeled by implementing Pitzer's equations for NaCl-solution properties. The viability of MSVMD use for zero liquid discharge (ZLD) applications is investigated, by considering discharge salinities close to NaCl saturation conditions. Energy efficiency (gained output ratio or GOR), second law efficiency, and the specific membrane area were used to quantify the performance of the system. At high salinities, the increased boiling point elevation of the feed stream resulted in lower fluxes, larger heating requirements and lower GOR values. The second law efficiency, however, is higher under these conditions since the least heat for separation increases faster than the system's specific energy consumption with increase in salinity. Under high salinity conditions, the relative significance of irreversible losses is lower. Results indicate that MSVMD systems can be as efficient as a conventional MSF system, while using reasonable membrane areas and for a wide range of feed salinities. Given MD's advantages over MSF such as lower capital requirement and scalability, MSVMD can be an attractive alternative to conventional thermal desalination systems. (C) 2015 Elsevier B.V. All rights reserved.