Blast furnace ironmaking process monitoring is an important and challenging task. Due to the influence of hot blast stove switching and large fluctuations in the quality of raw materials, the measurements of ironmaking processes show obvious non‐stationary characteristics, and in addition, the observed data are also characterized by time‐series dynamic and non‐Gaussian characteristics. In this paper, a dynamic stationary subspace analysis method based on the Gaussian mixture model (DSSA–GMM) is proposed to address the difficulties in blast furnace ironmaking process monitoring. The time‐series dynamic relationship of the data is conducted by introducing a sliding time window. The Gaussian mixture model (GMM) is used to deal with the non‐Gaussian characteristics of the data, and the parameters of the GMMs are estimated using the expectation–maximization algorithm. The stationary projection matrix is obtained by optimizing the Kullback–Leibler (K–L) divergence between GMMs of different periods to realize the stationary subspace separation. Finally, the convex hull of the stationary subspace is established for fault detection, thus realizing the monitoring for non‐stationary and non‐Gaussian dynamic processes. The effectiveness of the DSSA–GMM method is verified by a numerical simulation and a dataset collected from an actual blast furnace ironmaking process.