We present a robust fault diagnosis method for uncertain multiple input multiple output (MIMO) linear parameter Varying (LPV) parity equations. The fault detection methodology is based on checking whether measurements are inside the prediction bounds provided by the uncertain MIMO LPV parity equations. The proposed approach takes into account existing couplings between the different measured outputs. Modelling and prediction uncertainty bounds are computed using zonotopes. Also proposed is an identification algorithm that estimates model parameters and their uncertainty such that all measured data free of faults will be inside the predicted bounds. The fault isolation and estimation algorithm is based on the use of residual fault sensitivity. Finally, two case studies (one based on a water distribution network and the other on a four-tank system) illustrate the effectiveness of the proposed approach. (C) 2012 Elsevier Ltd. All rights reserved.