<p>The quantile-crossing spectrum is the spectrum of quantile-crossing processes created from a time series by the indicator function that shows whether or not the time series lies above or below a given quantile at a given time. This bivariate function of frequency and quantile level provides a richer view of serial dependence than that offered by the ordinary spectrum. A new estimator is proposed in this paper for the quantile-crossing spectrum as a bivariate function of frequency and quantile level. The proposed estimator is derived from a method called spline autoregression (SAR). It jointly fits an autoregressive (AR) model to the quantile-crossing series across multiple quantiles, where the functional AR coefficients are represented as spline functions of the quantile level and penalized for their roughness. Numerical experiments show that when the underlying spectrum is smooth in quantile level the proposed method is able to produce more accurate estimates in comparison with the alternative that ignores the smoothness.</p>