Beaches around the world continuously adjust to daily and seasonal changes in wave and tide conditions, which are themselves changing over longer timescales. Different approaches to predict multi-year shoreline evolution have been implemented; however, robust and reliable predictions of shoreline evolution are still problematic even in short-term scenarios (shorter than decadal). Here we show results of a modelling competition, where 19 numerical models (a mix of established shoreline models and machine learning techniques) were tested using data collected for tairua beach, new Zealand with 18 years of daily averaged alongshore shoreline position and beach rotation (orientation) data obtained from a camera system. in general, traditional shoreline models and machine learning techniques were able to reproduce shoreline changes during the calibration period (1999-2014) for normal conditions but some of the model struggled to predict extreme and fast oscillations. During the forecast period (unseen data, 2014-2017), both approaches showed a decrease in models' capability to predict the shoreline position. this was more evident for some of the machine learning algorithms. A model ensemble performed better than individual models and enables assessment of uncertainties in model architecture. Research-coordinated approaches (e.g., modelling competitions) can fuel advances in predictive capabilities and provide a forum for the discussion about the advantages/disadvantages of available models. Quantitative prediction of beach erosion and recovery is essential to planning resilient coastal communities with robust strategies to adapt to erosion hazards. Over the last decades, research efforts to understand and predict shoreline evolution have intensified as coastal erosion is likely to be exacerbated by climatic changes 1-5. The social and economic burden of changes in shoreline position are vast, which has inspired development of a growing variety of models based on different approaches and techniques; yet current models can fail (e.g. predicting erosion in accreting conditions). The challenge for shoreline models is, therefore, to provide reliable, robust and realistic predictions of change, with a reasonable computational cost, applicability to a broad variety of systems, and some quantifiable assessment of the uncertainties.

Sandy shorelines respond to oceanographic and geologic processes at a variety of temporal (e.g., seconds to decades) and spatial (e.g., surf zone to ocean basin wide) scales. Thus, forecasting short-and long-term shoreline change remains difficult. Accurate predictions of shoreline movement in response to sea-level rise, changing wave climates, and reduced sediment supplies are increasingly sought to support coastal management, often out of fears that many beaches may be severely impacted by climate change (Le Cozan-

Rising seas coupled with ever increasing coastal populations present the potential for significant social and economic loss in the 21st century. Relatively short records of the full multidimensional space contributing to total water level coastal flooding events (astronomic tides, sea level anomalies, storm surges, wave run-up, etc.) result in historical observations of only a small fraction of the possible range of conditions that could produce severe flooding. The Time-varying Emulator for Short-and Long-Term analysis of coastal flood hazard potential is presented here as a methodology capable of producing new iterations of the sea-state parameters associated with the present-day Pacific Ocean climate to simulate many synthetic extreme compound events. The emulator utilizes weather typing of fundamental climate drivers (sea surface temperatures, sea level pressures, etc.) to reduce complexity and produces new daily synoptic weather chronologies with an auto-regressive logistic model accounting for conditional dependencies on the El Niño Southern Oscillation, the Madden-Julian Oscillation, seasonality, and the prior two days of weather progression. Joint probabilities of sea-state parameters unique to simulated weather patterns are used to create new time series of the hypothetical components contributing to synthetic total water levels (swells from multiple directions coupled with water levels due to wind setup, temperature anomalies, and tides). The Time-varying Emulator for Short-and Long-Term analysis of coastal flood hazard potential reveals the importance of considering the multivariate nature of extreme coastal flooding, while progressing the ability to incorporate large-scale climate variability into site specific studies assessing hazards within the context of predicted climate change in the 21st century.Plain Language Summary Predicting extreme coastal flooding is a present-day societal need and will only become more relevant as mean water levels increase due to sea level rise. However, the number of processes contributing to such events is too high for relatively short observational records to have measured all of the constructive combinations of waves, surge, wind, and sea level anomalies. We present a framework designed to create hypothetical combinations of relevant flood hazard potential processes by simulating the climate and weather patterns that drive coastal flooding. Including large-scale oceanic and atmospheric patterns as the drivers of coastal hazards reveals the climate a coastal community is most vulnerable to, which will be increasingly more important to understand as the climate changes during the 21st century.