Multivariate functional data are increasingly common in various applications. This paper introduces a covariate-regulated sparse subspace learning (CSSL) model to address the challenges of complex, time-varying cross-correlation in such data, with applications in process monitoring and fault isolation. The efficacy of the CSSL model is demonstrated through simulations and a case study using SCADA data from wind turbines.