Jilin Wu, Ruike Wu, Zhijie Xiao, Mengxi Zhang
We propose an omnibus nonparametric test for structural changes in the multivariate volatility matrix. The test aggregates bounded generalized quantile scores over a range of quantile levels and has a weighted leave-$q$-out $U$-statistic representation. Deleting nearby index pairs renders the centering effect induced by serial dependence asymptotically negligible. All quantities required for implementation, including the variance estimator used for standardization, are constructed under the null, without specifying volatility dynamics under the alternative. The standardized statistic converges to a standard normal distribution. We establish consistency against fixed alternatives that generate a positive integrated quantile-score signal and derive nontrivial local power against smooth departures and increasingly sharp transitions approaching multiple structural breaks. The bounded-score construction avoids the finite fourth- or eighth-moment conditions commonly imposed by least-squares and quasi-likelihood procedures, while aggregation across quantiles uses more distributional information than single-quantile methods. Monte Carlo results show satisfactory size and favorable power under heavy-tailed innovations, with competitive performance under Gaussian innovations. An application to the Fama--French three-factor model provides evidence against stability of the factor covariance matrix over the full sample and several economically relevant subsamples.