The article will focus on vector autoregressive‐moving average (VARMA) models with time‐dependent coefficients (td) to represent general nonstationary time series, not necessarily Gaussian. The coefficients depend on time, possibly on the length of the series , hence the name tdVARMA for the models, but not necessarily on the rescaled time . As a consequence of the dependency on of the model, we need to consider array processes instead of stochastic processes. Under appropriate assumptions, it is shown that a Gaussian quasi‐maximum likelihood estimator is consistent in probability and asymptotically normal. The theoretical results are illustrated using three examples of bivariate processes, the first two with marginal heteroscedasticity. The first example is a tdVAR(1) process while the second example is a tdVMA(1) process. In these two cases, the finite‐sample behavior is checked via a Monte Carlo simulation study. The results are compatible with the asymptotic properties even for small . A third example shows the application of the tdVARMA models for a real time series.