The best answer to your question is probably given by the Nobel prize committee itself in "The Prize in Economic Sciences 2003 - Advanced Information" document. You should read it in full. Below is an excerpt.
According to the committee:
Financial economists have long since known that volatility in returns tends to
cluster and that the marginal distributions of many asset returns are leptokurtic, which means that they have thicker tails than the density of the normal distribution with the same mean and variance. Even though the time clustering of returns was known to many researchers, returns were still modeled as independent and identically distributed over time. Examples include Mandelbrot (1963) and Mandelbrot and Taylor (1967) who used so-called stable Paretian distributions to characterize the distribution of returns. Robert Engle’s modelling of time-varying volatility by way of autoregressive conditional heteroskedasticity (ARCH) thus signified a genuine breakthrough.
Suppose at time $t$ we observe the stochastic vector $(u_t, \mathbb{x}_t^{'})$ where $u_t$ is a scalar and $\mathbb{x}_t$ may contain some lags of $u_t$. The predictive model for the variable $u_t$ is $$ u_t = \mathbb{E}(u_t | \mathbb{x}_t) + \epsilon_t \tag{3.1} \label{one}$$
It is typically assumed that $\mathbb{E}(u_t | \mathbb{x}_t)$ has a parametric form and $\mathbb{E}\epsilon_t = 0$, $\mathbb{E}\epsilon_t \epsilon_s = 0, \forall t \not = s$
When estimating the parameters in $\mathbb{E}(u_t | \mathbb{x}_t)$, it was typically assumed that the unconditional variance of the error term $\epsilon_t$ is constant or time-varying in an unknown fashion. Engle considered the alternative assumption that, while the unconditional error variance − if it exists − is constant, the conditional error variance is time-varying.
This revolutionary notion made it possible to explain systematic features in the movements of variance over time and, a fortiori, to estimate the parameters of conditional variance jointly with the parameters of the conditional mean. The literature is wholly devoid of earlier work with a similar idea.
Engle parameterized the conditional variance of $\epsilon_t$ in the model $\eqref{one}$ such that large positive or negative errors $\epsilon_t$ were likely to be followed by another large error of either sign and small errors by a small error of either sign. Nowdays this is usually refered to as volatility clustering. He assumed that $\epsilon_t$ can be decomposed as $\epsilon_t = z_t \sigma_t^{\frac{1}{2}}$ where ${z_t}$ is a sequence of iid random variables with zero mean and unit variance and where $$ \sigma_t = \operatorname{Var}(\epsilon_t | \mathcal{F}_t) = \alpha_0 + \sum_{j=1}^m \alpha_j \epsilon_{t-j}^{2} \tag{3.2}\label{two}$$
In $\eqref{two}$ $\epsilon_t = u_t - \mathbb{E}(u_t | \mathbb{x}_t)$, $\alpha_0 > 0$ and $\alpha_j \geq 0, j = 1, \cdots m$, the information set $\mathcal{F}_t = \sigma(\{\epsilon_{t - j}, j \geq 1\})$
Equation $\eqref{two}$ defines the ARCH model introduced in "Autoregressive conditional heteroscedasticity with estimates of the variance of United Kingdom inflation", where the conditional variance is a function of past values of squared errors.
In this classic paper, Engle developed the estimation theory for the ARCH model, gave conditions for the maximum likelihood estimators to be consistent and asymptotically normal, and presented a Lagrange multiplier test for the hypothesis of no ARCH (equal conditional and unconditional variance) in the errors $\epsilon_t$.