Abstract Increments in financial markets have anomalous statistical properties including fat-tailed distributions and volatility clustering (i.e., the autocorrelation functions of return increments decay quickly but those of the squared increments decay slowly). One of the central questions in financial market analysis is whether the nature of the underlying stochastic process can be deduced from these statistical properties. We have shown previously that a class of variable diffusion processes has fat-tailed distributions. Here we show analytically that such models also exhibit volatility clustering. To our knowledge, this is the first case where clustering of volatility is proven analytically in a model. Our results are compatible with the viewpoint that variable diffusion processes are possible models for financial markets.