Description: Ocean bottom pressure (OBP) variability serves as a proxy of ocean mass variability, the knowledge of which is needed in geophysical applications. The question of how well it can be modeled by the present general ocean circulation models on time scales in excess of 1 day is addressed here by comparing the simulated OBP variability with the observed one. To this end, a new multiyear data set is used, obtained with an array of bottom pressure gauges deployed deeply along a transect across the Southern Ocean. We present a brief description of OBP data and show large‐scale correlations over several thousand kilometers at all time scales using daily and monthly averaged data. Annual and semiannual cycles are weak. Close to the Agulhas Retroflection, signals of up to 30 cm equivalent water height are detected. Further south, signals are mostly intermittent and noisy. It is shown that the models simulate consistent patterns of bottom pressure variability on monthly and longer scales except for areas with high mesoscale eddy activity, where high resolution is needed to capture the variability due to eddies. Furthermore, despite good agreement in the amplitude of variability, the in situ and simulated OBP show only modest correlation.

Description: Extreme climatic and geophysical events pose a threat to societies and have the capacity to cause significant damage and losses whenever they occur, both in their immediate aftermath and in the medium- to long-term. Their consequences can be amplified even further when more than one event affects the same geographical areas within a short time. Be it cascading hazards, in which one event triggers the next, or simply hazards that happen to occur simultaneously (“compounding” hazards), estimation of their cumulative consequences is challenging because the action of one event affects the exposure and vulnerability to the next one. While the efforts from the research community to develop multi-hazard perspectives have increased considerably in recent years, multiple remaining challenges require strongly-coordinated efforts across different disciplines and areas of expertise to tackle them with the most appropriate tools. With a multidisciplinary team of scientists from four different Helmholtz research centres in Germany, we have started working on the CASCO project (2022-2024), in which we will develop an integrated risk workflow for CAScading and COmpounding hazards in COastal urban areas by focusing on a series of events occurring around Mount Etna (Italy). The case-scenario starts with a strong earthquake that triggers a submarine collapse at the eastern flank of Mount Etna, an area already known to be unstable, and both the earthquake and the landslide trigger a tsunami that hits the coasts of Sicily and Calabria. Almost concomitantly, a heatwave or heavy rainfall happens to affect the same regions, further stressing the population that had been affected by the combined effects of the earthquake and tsunami. The project will be directed towards the modelling of the cascading earthquake, landslide and tsunami events, the compounding heatwave and rainfall, as well as their immediate impacts in terms of cumulative damage and casualties. Moreover, the medium- to long-term response in urban dynamics and the effect of these extreme events on the economic development of the affected populations will be explored. By focusing on a tangible scenario, CASCO will not only tackle the challenges associated with bringing together the whole risk chain (which will be valid beyond our case-study) but also produce outcomes that help increase awareness of such extreme events and the need for societies to develop suitable strategies to strengthen their resilience and improve their disaster response.

Description: In this article, the tsunami model TsunAWI (Alfred Wegener Institute) and its application for hindcasts, inundation studies, and the operation of the tsunami scenario repository for the Indonesian tsunami early warning system are presented. TsunAWI was developed in the framework of the German-Indonesian Tsunami Early Warning System (GITEWS) and simulates all stages of a tsunami from the origin and the propagation in the ocean to the arrival at the coast and the inundation on land. It solves the non-linear shallow water equations on an unstructured finite element grid that allows to change the resolution seamlessly between a coarse grid in the deep ocean and a fine representation of coastal structures. During the GITEWS project and the following maintenance phase, TsunAWI and a framework of pre- and postprocessing routines was developed step by step to provide fast computation of enhanced model physics and to deliver high quality results.

Description: Earth Rotation Parameters (ERP) are used to improve estimates of the ocean circulation and mass budget. GRACE data can be used for verification or for further improvements. The Finite Element Sea-ice Ocean Model (FESOM) is used to simulate weekly ocean circulation and mass variations. The FESOM model is a hydrostatic ocean circulation model with a fully non-linear free surface. It solves the hydrostatic primitive equations with volume (Boussinesq approximation) and mass (Greatbatch correction) conservation. Fresh water exchange with the atmosphere and land is modelled as mass flux. This flux is the weakest part of the mass budget as it is the difference of large and uncertain quantities: evaporation, precipitation and river runoff. All uncertainties included in these parameters are directly reflected in the model results. ERP help in closing the budget in a realistic manner. Our strategy is designed for testing parametric estimation on a weekly basis. First, Oceanographic Earth rotation parameters (OERP) are calculated by subtracting atmospheric and hydrologic estimates from observed ERP. They are compared to OERP derived from a global ocean circulation model. The difference can be inverted to diagnose a correction of the oceanic mass budget. Additionally mass variations measured by GRACE are used for verification. In a second step, the global mass correction parameter, derived by the inversion, is used to improve the fresh water budget of FESOM.

Description: Representation of the shelf areas in the global ocean model: key study, questions and perspectives The East Siberian Arctic Shelf (ESAS), consisting of the Laptev, East Siberian and Chukchi Seas, represents the shallowest and broadest shelf region of the entire World Ocean. It occupies a little more than 20% of the total area covered by the Arctic Ocean (AO) and represents a critical physical and biochemical gateway for exchange between AO and terrestrial zone with complex oceanographic and biogeochemical regime influenced by both seawaterof Pacific and Atlantic origins. The is a growing need for better quality estimations of circulation and dynamics on the shelf to answer major present and future scientific, ecosystem and societal issues, because of changing climate. It is a complex task as soon as the ESAS represents wide area with variety of regimes and there is still substantial uncertainty in their role and feedbacks with the wider climate system. Making progress on this is largely dependent on the accurate reproducing of the physical environment in the coupled coastal-open ocean system. We would like to propose modeling system that will help to answer questions on the ESAS observed and future trends and dynamics features across time and space scales tracing the signal through the system Estuaries-ESAS-AO in both upscaling and downscaling directions. To reach mentioned goals, we built a coastal branch of the finite volume version of the global sea ice-ocean model FESOM (Danilov et al., 2004; Danilov, 2012; Wang et al., 2014). FESOM is the first model worldwide that provides multi-resolution functionality to large-scale ocean modeling, allowing to bridge the gap between scales. This unique feature is crucial for high efficient coupling, as soon as the exchange zone can be resolved similarly (with the same resolution) by the global and local solutions. Additional strong side of the elaboration of the coastal branch for the existing global model is a possibility to organize flux treatment in a same manner, increasing efficiency of coupling. Danilov, S., Kivman, G., Schröter, J. (2004). A finite-element ocean model: principles and evaluation, Ocean Model., 6, 125–150. Danilov, S. (2012). Two finite-volume unstructured mesh models for large-scale ocean modeling. Ocean Modell., 47, 14–25. Wang, Q., Danilov, S., Sidorenko, D., Timmermann, R., Wekerle, C., Wang, X., Jung, T., and Schröter, J. (2014). The Finite Element Sea Ice-Ocean Model (FESOM) v.1.4: formulation of an ocean general circulation model, Geosci. Model Dev., 7, 663–693.

Description: This contribution focuses on two applications of the FESOM model family. On the one hand, recent runs with the finite volume code FESOM2 on large global meshes with regional focus are presented. FESOM's shallow water branch TsunAWI is the subject of the second part. TsunAWI, still based on finite elements, is used as an operational model in the Indonesia Tsunami Early Warning System (InaTEWS). InaTEWS derives tsunami forecasts in two different ways: from scenarios in a pre-computed database or from an on-the-fly simulation. The pre-computed scenarios are based on TsunAWI simulations with inundation on a triangular mesh with a resolution ranging from 20km in the deep ocean to 300m - 50m in coastal areas. The on-the-fly propagation model EasyWave (Andrey Babeyko, GFZ) solves the linear shallow water equations on a regular finite-difference grid with a resolution of about 1 km and the coast line as a vertical wall. EasyWave is used after a tsunami has been generated in an area not covered by the database or after seismic measurements show an earthquake mechanism not present in the database. As the numerical settings of both models are quite different, variations in the outputs are to be expected; nevertheless, the differences in the warning levels should not be too large for identical sources. In the current study, we systematically compare the warning products like estimated wave height and estimated time of arrival by the two approaches.

Description: Abstract: The Indonesia Tsunami Early Warning System delivers simulated tsunami forecasts in two different ways: either matching scenario(s) from a pre-computed database or running on-the-fly tsunami simulation. Recently, the database has been extended considerably taking into account additional source regions not covered in earlier stages of the system. In this contribution, we present the current status of the data base coverage as well as a study investigating the warning products obtained by the two modeling approaches. The pre-computed tsunami scenarios are based on the finite element model TsunAWI that employs a triangular mesh with resolution ranging from 20km in deep ocean to 300m in coastal areas and to as much as 50m in some highly resolved areas. TsunAWI solves the nonlinear shallow water equations and contains a wetting-drying inundation scheme. The on-the-fly propagation model easyWave solves the linear shallow water equations on a regular finite-difference grid with a resolution of about 1 km and utilizes several simple options to estimate coastal impact. This model is used for forecasting after a tsunami has been generated in an area not covered by the database or after a moment tensor solution shows an earthquake focal mechanism not present in the database. Since warning products like estimated wave height (EWH) and estimated time of arrival (ETA) along the coast are based on modeling results, it is crucial to compare the resulting forecasted warning levels obtained by the two approaches. Resolutions and numerical settings of both models are quite different, therefore variations in the resulting outputs are to be expected; nevertheless, the extent of differences in warning levels should not be too large for identical sources. In the present study, we systematically investigate differences in resulting warning products along InaTEWS forecast points facing the Sunda arc.  Whereas the finite-element mesh of TsunAWI covers the coast up to a terrain height of 50m and warning products have been pre-calculated directly in the forecast points, easyWave offers several options for their approximation including projections from offshore grid points or vertical wall. Differences and potential reasons for variations of warning products like the role of bathymetry resolution as well as the general approach for the assessment of EWH and ETA for different modeling frameworks are discussed.