The spatio-temporal resolution of atmospheric forcing plays a key role in the accuracy of simulated storm surges with hydrodynamic numerical models. Here, we generate five hydrodynamic hindcasts of coastal storm surges along the European Atlantic and the Mediterranean Sea coasts, forced with atmospheric fields of varying temporal (hourly and daily) and spatial (0.25 to 2) resolution since 1940. Our results, that are validated with insitu tide gauge observations, show that storm surges obtained with daily forcing underestimate the magnitude of coastal extreme sea level events by up to 50% compared to hourly simulations and observations. Nevertheless, low-resolution simulations capture the temporal variability of storm surges, including strong episodes. Furthermore, taking advantage of the consistent set of coastal storm surge hindcasts, we demonstrate that storm surges forced with daily mean atmospheric fields, when bias corrected via quantile mapping, provide accurate values of daily maxima as calculated by a high-resolution hindcast. This transformation paves the way to obtain daily maxima storm surge estimates from low-resolution atmospheric fields, as those typically provided by large-scale and global climate models, at a lower computational cost.
Submerged topography in shallow waters is fundamental in the propagation and dissipation of ocean waves in the surf and swash zones. However, obtaining accurate bathymetric data in this region is challenging due to the high temporal and spatial environmental variability. The bottom boundary condition can directly affect the accuracy of numerical models used for shallow water simulations. In this study, the performance of the SWASH numerical model in describing wave runup in the swash zone is assessed using different bathymetric boundary conditions. The first method involves using data measured in the surf zone obtained by a Unmanned Aerial Vehicle (UAV), and analyzing it using the cBathy algorithm. The second method utilizes a regular bathymetric mesh generated from Dean’s equilibrium profile combined with beach topography data. The third method relies exclusively on interpolation methods using data from deep waters and beach profiles. This interpolation approach is the most used among SWASH users when detailed or updated surf zone bathymetry is unavailable. Based on the numerical simulations performed, not incorporating data from the surf zone resulted in a 4% increase in the runup estimated and approximately a 2% difference in identifying the swash zone position. The method to obtain bathymetry through the cBathy algorithm, used in this article, is cost-effective and can be used to reduce uncertainties in surf zone numerical simulations, induced by the lack of knowledge about the bottom conditions.
In this short communication, we report initial success in utilizing existing Explainable Artificial Intelligence (XAI) methodology to investigate an emerging precursor of the El Niño-Southern Oscillation (ENSO), manifest as sea surface temperature anomalies (SSTA) in the Western North Pacific (WNP), and its impact on enhancing ENSO prediction accuracy. Our analysis reveals that integrating WNP SSTA with established XAI techniques significantly increases the predictability of ENSO states. We found marked improvement in prediction accuracy, from a 60 % baseline to over 85 % for forecasting moderate warm, cold, and neutral ENSO states one year ahead. For higher magnitude events, precision surpasses 90 %. This work, intended as a follow-up to recent studies, underscores the potential of augmenting emerging XAI with additional SST signals to advance long-term climate forecasting capabilities.
A realistic representation of the Southern Ocean (SO) in climate models is critical for reliable global climate projections. However, many models are still facing severe biases in this region. Using a fully coupled global climate model at non-eddying (1/2) and strongly eddying (1/10) grid resolution in the SO, we investigate the effect of a 0.5 °C, 1.0 °C and 1.6 °C warmer than observed SO on i) the spin-up behaviour of the high-resolution simulation, and ii) the representation of main dynamical features, i.e., the Antarctic circumpolar current (ACC), the subpolar gyres, the overturning circulation and the Agulhas regime in a quasi-equilibrium state. The adjustment of SO dynamics and hydrography critically depends on the initial state and grid resolution. When initialised with an observed ocean state, only the non-eddying configuration quickly builds up a strong warm bias in the SO. The high-resolution configuration initialised with the biased non-eddying model state results in immense spurious open ocean deep convection, as the biased ocean state is not stable at eddying resolution, and thus causes an undesirable imprint on global circulation. The SO heat content also affects the large-scale dynamics in both low- and high-resolution configurations. A warmer SO is associated with a stronger Agulhas current and a temperature-driven reduction of the meridional density gradient at 45S to 65S and thus a weaker ACC. The eddying simulations have stronger subpolar gyres under warmer conditions while the response in the non-eddying simulations is inconsistent. In general, SO dynamics are more realistically represented in a mesoscale-resolving model at the cost of requiring an own spin-up.