This paper presents a methodology that uses site-specific topographic and cosmogenic 10Be data to perform multi-objective model optimisation of a coupled coastal evolution and cosmogenic radionuclide production model. Optimal parameter estimation of the coupled model minimises discrepancies between model simulations and measured data to reveal the most likely history of rock coast development. This new capability allows a time series of cliff retreat rates to be quantified for rock coast sites over millennial timescales. Without such methods, long-term cliff retreat cannot be understood well, as historical records only cover the past ∼150 years. This is the first study that has (1) applied a process-based coastal evolution model to quantify...