Connecting Soft and Hard: An Integrating Role of Systems Dynamics in Tsunami Modeling and Simulation

: Modeling and simulation have been used to study tsunamis for several decades. We created a review to identify the software and methods used in the last decade of tsunami research. The systematic review was based on the PRISMA methodology. We analyzed 105 articles and identified 27 unique software and 45 unique methods. The reviewed articles can be divided into the following basic categories: exploring historical tsunamis based on tsunami deposits, modeling tsunamis in 3D space, identifying tsunami impacts, exploring relevant variables for tsunamis, creating tsunami impact maps, and comparing simulation results with real data. Based on the outcomes of this review, this study suggests and exemplifies the possibilities of system dynamics as a unifying methodology that can integrate modeling and simulation of most identified phenomena. Hence, it contributes to the development of tsunami modeling as a scientific discipline that can offer new ideas and highlight limitations or a building block for further research in the field of natural disasters.


Introduction
Tsunamis are potentially extremely dangerous natural hazards that occur at low frequencies and are hard to predict.The "potentiality" depends on whether the tsunami hits an uninhabited area.In this case, nothing much happens regarding direct danger to people.However, it is an extreme risk if a populated area is affected.The main issue is that the events leading up to the tsunami do not clearly indicate its origin until it is generated [1].With frequently scarce [2], it is necessary to develop a coherent framework that incorporates existing assumptions (i.e., the system's general model) and various methods for hazard and risk analysis in order to assess the consequences of these events on various layers of society.The underlying objective of COST Action AGITHAR (Accelerating Global Science in Tsunami Hazard and Risk Analysis) is to develop further, standardize, and document such a framework, and this document is one of the Action's outputs.Numerous sources can generate tsunamis in the form of long propagating waves.Tsunamis are a phenomenon primarily generated along convergent margins and caused by undersea earthquakes [3].Tsunamis, once created, travel at high speeds and cover a large area of water.In case of a large tsunami, when waves reach coastal areas, they inundate land for several kilometers.
As a result, casualties occur, as does damage to or destruction of infrastructure and built-up areas.More than 700 million people live in areas of extreme sea-level events, including tsunamis [4].
Research on tsunamis is multidisciplinary.Tsunamis can be studied from many different perspectives, namely the historical perspective, where tsunami deposits are examined, and the period when the tsunami arrived at the selected location is estimated, for example [5].Furthermore, tsunamis are investigated from the perspective of their impact on the environment or on individual inhabited areas.Here, in particular, photo documentation is used before and after the wave impact, e.g., [6].Another way of exploring tsunamis is their modeling, where tsunamis are usually presented as a three-dimensional environment, e.g., [7].Models are usually developed to investigate impacts on a selected environment or on a populated area, e.g., [8], or on a hypothetical tsunami barrier, e.g., [9].Another way of examining tsunamis is from a statistical perspective, either by searching for essential variables for the part of the tsunami being examined, e.g., [10], or from the statistical comparison of the impacts of a real tsunami with a modeled one, e.g., [11].
To date, no comprehensive systematic review has been published focusing on methods and software packages applied to explore the multidisciplinary phenomenon of the tsunami.An example of a systematic review thematically focused on tsunamis and similar events is provided by Palupi [12], who focuses on the psychological preparedness of coastal communities for a natural disaster.Fernandez et al. [13] focused on the mapping of scientific evidence on the mental health impacts of floods caused by extended periods of heavy rain in river catchments.The aim of this study is to contribute to both aforementioned topics through a deep analysis of existing modeling and simulation techniques and methods in the field of tsunami research and an outline of the possibility of further extension.The achievement of this aim is based on highlighting the importance of finding a tool that is not over-specialized and narrowly focused on particular aspects of tsunamis.While existing tools used in particular fields of study need to be applied to find answers to particular domain-related questions. the multidisciplinary nature of tsunamis requires applying a systemic approach at a higher level of complexity.The holistic view and complex structures with mutually interconnected parts represent the primary attributes of the systems approach.Hence, data or information acquired by particular tools, techniques, or modeling software can be used as inputs to more complex models, improving understanding and insights into the tsunami phenomenon.System dynamics is a commonly used methodology whose application can successfully lead to implementing the system perspective.The goal is to prove its applicability and potential added value in tsunami research.
This manuscript is structured as follows.Section 2 describes the applied methodology, which is based on a systematic review methodology, Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA).Furthermore, technical details related to developing a system dynamics model are presented.The next two sections present achieved results.Section 3 provides a bibliographic analysis based on the specific keywords of identified articles.Furthermore, the results of the analysis of different types of software and methods for tsunami research are presented.Section 4 describes the developed system dynamics model, which demonstrates how various aspects can be integrated by one methodological approach.Particular modules are introduced.The final section concludes the paper with an outline of existing limitations and potential research pathways.Finally, the appendices capture an enumeration of the individual software applications, their methods, and reasons for their use.Moreover, an overview of the relationships between the software and the scientific disciplines and technical details associated with the developed model are provided.

Literature Analysis
Collecting a data set that might be used for further analysis was the initial step in the study.To conduct a review, one can choose from two alternatives.The first option involves looking for relevant information resources in databases maintained by individual publishing houses (Elsevier, Springer, Wiley, etc.).The second approach uses databases that have selected journals indexed (Web of Science, Scopus, EBSCO, etc.).Both strategies have advantages and disadvantages.The former, for example, provides a data collection from a broader set of resources, but the latter works with publications whose quality is recognized by authority and a community.We found the latter more suitable for this study due to the absence of redundant records.Thus, a search in the Web of Science database for published papers was conducted.
Articles focused on tsunamis that used software for this activity represented the target of the search.The main criteria for searching articles were the language used, namely English; the year of publication, between 2010 and 2020; full-text availability (in order to conduct the content analysis); and the use of a software package for tsunami exploration.Initially, we intended to use the following search command for the intended systematic search: However, articles containing the term "tsunami" broadly did not refer to tsunamis in the true sense of the word.An example of this phrase might be the tsunami of obesity.Therefore, we decided to narrow the search.The additional concept "water" filtered out articles that were not focused on the topic of this systematic search.Thus, we used the following command: The next step was to search for full-text versions of the articles.This step does not follow the exact PRISMA methodology [14]; however, due to the search for software used in the tsunami field and also the lack of a list of software used in the tsunami field, we decided to skip the abstract screening step.Another reason for skipping this stage is that applied software applications are not always mentioned in abstracts.
Consequently, we searched for the word "software" in the full-text versions of the articles.We used Adobe Acrobat Reader DC [15] software for the search.We recorded the name of each software found and entered it into the list of used software.After we had searched all available articles, we started a new search using the aforementioned software, which included the software names found in the previous full-text search.In this step, we sorted the articles based on those in which some of the found software was used and those in which it was not.We then searched the articles in which tsunami software was used (if the article did not use one of the found software, it was omitted with Reason 1).We then excluded articles from this set of articles that did not have a specific method (if the article did not use one of the found methods, it was omitted with Reason 2).The whole process is shown in Figure 1.Eventually, we found 2015 articles, from which we included 1573 in the full-text search, and 105 articles passed the full-text screening.
As a part of the analytical section of the systematic review, we performed a bibliographic synthesis of keywords using the VOSviewer [16] tool.Seventeen articles were not included in this stage due to failure to include keywords in the articles.Linguistic word preparation was conducted in which, for instance, keyword run-up and run-up were unified under the single term run-up, or the keyword tsunami was excluded from the synthesis because of the focus of all the articles on this topic.Moreover, the exclusion was necessary as we tried to answer the following research questions: (1) What software tools and methods are used in tsunami research?and (2) How are tsunami research tools interconnected?Figure 1.PRISMA stages (own work; template adopted from [14]).
As a part of the analytical section of the systematic review, we performed a bibliographic synthesis of keywords using the VOSviewer [16] tool.Seventeen articles were not included in this stage due to failure to include keywords in the articles.Linguistic word preparation was conducted in which, for instance, keyword run-up and run-up were unified under the single term run-up, or the keyword tsunami was excluded from the synthesis because of the focus of all the articles on this topic.Moreover, the exclusion was necessary as we tried to answer the following research questions: (1) What software tools and methods are used in tsunami research?and (2) How are tsunami research tools interconnected?

Model Development
System dynamics represents a specific and original methodological approach to modeling and simulating various types of systems.The core concepts of systems thinking, such as interconnectedness, feedback, adaptive capacity/resilience, self-organization, and emergence [17], are applied in system dynamics to help people make better decisions when confronted with complex, dynamic systems.The field provides a philosophy and tools to model and analyze dynamic systems.
Differential and difference equations are traditionally used to represent change in dynamic systems.However, system dynamics provides an intuitive modeling language, which is typical for all applications.This makes system dynamics an ideal tool for multidisciplinary work [18] as it enables the integration of subsystems that are distinct in their fundamental essence, i.e., soft disciplines such as economics or psychology can be connected with hard disciplines such as physics or geology.Exploration of complex phenomena such as tsunamis represents an appropriate example.Therefore, a complex model of various aspects of a tsunami has been developed to demonstrate the benefits and added value of incorporating this methodology into this field of study.

Model Development
System dynamics represents a specific and original methodological approach to modeling and simulating various types of systems.The core concepts of systems thinking, such as interconnectedness, feedback, adaptive capacity/resilience, self-organization, and emergence [17], are applied in system dynamics to help people make better decisions when confronted with complex, dynamic systems.The field provides a philosophy and tools to model and analyze dynamic systems.
Differential and difference equations are traditionally used to represent change in dynamic systems.However, system dynamics provides an intuitive modeling language, which is typical for all applications.This makes system dynamics an ideal tool for multidisciplinary work [18] as it enables the integration of subsystems that are distinct in their fundamental essence, i.e., soft disciplines such as economics or psychology can be connected with hard disciplines such as physics or geology.Exploration of complex phenomena such as tsunamis represents an appropriate example.Therefore, a complex model of various aspects of a tsunami has been developed to demonstrate the benefits and added value of incorporating this methodology into this field of study.
This study presents an original model presenting possibilities of tsunami analysis.The process of model development was based on the following procedure:

•
Identification of key topics in the scientific literature based on the literature review;

Application of Methods and SW Tools in the Current Research of Tsunamis
Figure 2 reveals that single methods or techniques are used mostly separately in research.One cluster is connected with the consequences of tsunamis in the form of numerical simulations that are primarily applied in the analysis of technical issues, such as construction vulnerability, and environmental issues, such as the influence of the tsunami wave on vegetation.The second cluster includes numerical modeling as well.This context deals with the origin of a tsunami and the wave itself.However, based on the discrepancy between the results from Table 1 presented below, where the interdependence of software tools usage reached a higher level, we decided to augment the input data with the name of the tool used.Thus, we added the names of the tools used to the keywords already provided by the authors of the articles.Then, we created a new outcome of the bibliographic analysis with the same settings as for the previous one.The results are presented in Figure 3.It is clear that most of the specific tools used in the tsunami field are interrelated.In fact, there is only one huge cluster and a few solitary outliers.It reveals that the applied methods are integrated mainly by the usage of OpenFoam, Matlab, and Calib.To understand the previous failure events and current hazard potential of the Mauritania Slide Complex along the NW-African continental margin.To conduct a PTHA for the Sydney metropolitan area, addressing the significant gap in published PTHAs, including inundation for Australia.
Xie and Chu, 2019 [88] OpenFOAM VOF To investigate the hydrodynamic forces on coastal structures impacted by tsunamis, focusing on the effects of the wave Froude number, the sizes and shapes of the structures, and the initial conditions.To assess tsunami evacuation planning for the Bodrum district along Turkey's western coast, which faces significant tsunami risks due to the active seismicity in the Eastern Mediterranean Sea.To investigate optimal dike shapes to mitigate tsunami overflow, focusing on numerical simulations based on waveforms from the Great East Japan Earthquake.
Based on Tables 1 and A1 presented in Appendix A, two types of the most numerous tsunami-related models can be identified, namely, models that directly model tsunamis or parts of tsunamis and models aimed at identifying events in history based on tsunami deposit examination.Table 2 shows the selected ratio of software use to the scientific disciplines (categories in which articles are classified in WoS) in which the articles were published (the whole table is in Appendix A, Table A3).As the Pareto rule of thumb suggests, 80.1% of software use belongs to 8 (30.8%) individual software.The situation is similar in the scientific fields, where 80.3% of the software used falls into 10 (34.5%) scientific fields related to tsunamis.
As mentioned earlier and shown in Figure 3, the different tsunami tools are interrelated.This interdependence is also evident in the scientific fields in which the papers were published.OpenFOAM (used in 32% of total cases) is the most multidisciplinary software application.For example, the MATLAB tool (used in 13.1% of total cases) was expected to be applied in various domains because of the multidisciplinary use of the tool in general.Surprisingly, the multidisciplinary application was not expected in the case of the CALIB tool (used in 8.6% of total cases).However, the results revealed the contrary.Note: These are selected fields and tools.For an overall overview, see Appendix A, Table A3.The fields of each article are retrieved from WoS.The number of fields varies for each article.

SD Model of Tsunami-Related Phenomena
As it is clear from the literature analysis, various methods and tools are used in tsunami research.However, system dynamics tools, particular software applications, or methods are not used.The reason is apparent.Applied techniques, tools, and software packages focus on a specific domain or set of issues.They enable highly specialized modeling and simulation.The aim of this section is to demonstrate the suitability of the system dynamics methodology for modeling and simulation of various aspects of tsunamis and the mutual interconnection of soft and hard disciplines.Hence, a model of tsunami-oriented aspects has been developed.The model contains seven main subsystems (modules), which represent topics under study, namely the module describing the tsunami as a natural phenomenon (The formation of the Tsunami) and six other modules showing the spheres affected by the tsunami wave, namely the modules People, Buildings, Infrastructure, Finance, Defensive elements, and Environment.Apparently, modules can be extended by any aspect we need to investigate.This configuration is created only for demonstrative purposes.

The Formation of the Tsunami
Apparently, it all starts with the formation of a tsunami.Developed links are arbitrary to enable meaningful and performable analysis.Together with the module containing defensive elements, the tsunami formation represents the "hard" part of the model in which exact relations among variables based on physical laws and principles must be modeled.Other modules contain more or less ambiguous relationships due to their softness.Their settings can be determined not only based on established relationships but also by model creators based on their expertise, best practices, or rules.The module contains tsunami formation and propagation and is presented in the following Figures 4 and 5.
This module is grounded in a set of converters, which mainly describe the behavior of the wave and their values.The module is based on equations presented in Appendix B. Water depth can be considered as the starting point of the module description.It is expressed using a graphical function.We have modeled a change in the depth of the water column, which is a function of time during the wave propagation and starts at the default value of 10,000 m and ends at 0 m.The change in the bottom terrain during the wave can be adjusted as needed based, for instance, on bathymetry data.Further variables associated with water depth are reckoned using the formulas from Appendix B.
The existence of shallow water waves must be verified, as tsunamis usually spread as this type of wave.This is verified by a simple condition that returns 0 if the shallow water wave does not exist and deactivates the module.
The wave height represents another important value associated with a wave at sea.Although unnecessary, the flow was chosen for this value instead of the converter, as this will allow a better connection with the wave amplitude value for which the conveyor must be used.The basic formula is used to determine the wave height at sea, which expresses that the more the sea depth decreases (the land approaches), the more the wave height increases.This value will be used to obtain information about the wave at the coast and the rise of the wave.
Another formula is based on maintaining wave energy flow and applies to so-called breaking and non-breaking waves.We also determine the value of the run-up factor and the maximum run-up height, which leads to the calculation of the Imamura-lida magnitude scale.This scale includes twelve different grades based on seismic intensity [123].Other values that the maximum run-up height needs as input are wave height before shoaling and inclination of the coast.Wave shoaling is an effect that causes the surface waves that enter the shore to change their height.The height of these waves is needed for the model.The value was set to 8 m for this particular wave.The inclination of the coast is also a given value of 1.30176852 radians.
Sci 2024, 6, x FOR PEER REVIEW which represent topics under study, namely the module describing the tsunami a ural phenomenon (The formation of the Tsunami) and six other modules show spheres affected by the tsunami wave, namely the modules People, Buildings, Inf ture, Finance, Defensive elements, and Environment.Apparently, modules can tended by any aspect we need to investigate.This configuration is created only for strative purposes.The model uses another slope, namely the Inclination of the water surface.This is also set as a fixed value for the given wave as 0.00021025 radians.This value is then used to calculate the flow velocity of the run-up.We also need to calculate the value of inundation depth and the so-called Manning's roughness coefficient.The Manning's roughness coefficient is used in the Manning formula to calculate flow in open channels.In this case, the individual materials are entered in the fields, and a corresponding coefficient is added to them for easy material change.Now, for example, the material Earth Channel Weedy, a grassy earth channel with a coefficient of 0.030, is used [124].This coefficient is also used for other calculations, such as inundation distance and flood distance.Its calculation uses only two values: coefficient and wave height at the shore, which is set at 1.928 m.Wave height at the shore is also used to calculate the loss in wave height per meter of inundation distance and the loss of wave height per meter of flooding distance.The last calculation is the wave run-up, which uses the above values and the so-called experimental constant, with a value of 0.5.
Five sectors further extend this basic module.These sectors deal with the emergence of the tsunami and can be linked to the primary model.Specific sectors are seismic activity, atmospheric disturbances, submarine landslides, meteorite impact, and volcanism.

People
This module illustrates the effect of a tsunami on the number of people in different categories, as presented in Figure 6.It is run entirely in units of People.Nine levels are included.In the first Population in the state reservoir, at time zero, the number of all people in the area affected by the tsunami is zero.In fact, except for the Population in the state variable, all other reservoirs are set to zero during the simulation initiation.This module is activated once the module The formation of the Tsunami is activated.The outflows separate the population into two other reservoirs, namely the affected population and the unaffected population, in both of these reservoirs at time zero because when the tsunami has not yet arrived, there is no need to distinguish between these two categories.People are further divided into those who are injured, dead, or uninjured into the Injured people, Dead people, and Unharmed people reservoirs.Rescuers is another reservoir where we have a number of rescuers in the area.Rescuers are then also divided into uninjured, wounded, and dead through model outflows.The last two reservoirs concern paramedics.The first is Paramedics in the state, where we have the number of all medics available, and the second is Paramedics in the area, where we have the number of all medics in a given location affected by the wave.

Buildings
This module, with the structure presented in Figure 7, demonstrates how the tsunami damaged or completely destroyed buildings in the area, as well as how reparation is funded from the Finance module later on.This module is influenced by The formation of the Tsunami module as well as the Finance module and is entirely in units of Buildings.The module is constructed based on the works of Leone at al. [125] and Rossetto at al. [126].
There are four reservoirs in this module, namely Buildings, Damaged buildings, Repaired buildings, and Destructed buildings.The Buildings repository contains a number of different types of buildings that are located in a given area.Buildings flow into the Damaged buildings storehouse, which are damaged by the tsunami and are classified according to the type of damage.Buildings that were completely destroyed by the tsunami wave flow into the Destructed buildings reservoir.After a while, the buildings will start to be repaired, and from the Damaged buildings reservoir, the buildings will start to flow into the Repaired buildings reservoir according to the volume of funds flowing from the last reservoir in this model, namely Finances divided into areas.This stock comes from the Finance module.

Buildings
This module, with the structure presented in Figure 7, demonstrates how the tsunami damaged or completely destroyed buildings in the area, as well as how reparation is funded from the Finance module later on.This module is influenced by The formation of the Tsunami module as well as the Finance module and is entirely in units of Buildings.The module is constructed based on the works of Leone at al. [125] and Rossetto at al. [126].
There are four reservoirs in this module, namely Buildings, Damaged buildings, Repaired buildings, and Destructed buildings.The Buildings repository contains a number of different types of buildings that are located in a given area.Buildings flow into the Damaged buildings storehouse, which are damaged by the tsunami and are classified according to the type of damage.Buildings that were completely destroyed by the tsunami wave flow into the Destructed buildings reservoir.After a while, the buildings will start to be repaired, and from the Damaged buildings reservoir, the buildings will start to flow into  In three reservoirs, Buildings, Damaged buildings, and Repaired buildings, we have used arrays, namely "Types of buildings", where the buildings are divided according to their durability.In the Damaged buildings and Repaired buildings reservoirs, another array is used, namely "Damage", where we have types of building damage (initial values are in oval brackets, array Damage has zero initial values):

•
Array Types of building • Wood and metal plates (180); • Brick building (550); In three reservoirs, Buildings, Damaged buildings, and Repaired buildings, we have used arrays, namely "Types of buildings", where the buildings are divided according to their durability.In the Damaged buildings and Repaired buildings reservoirs, another array is used, namely "Damage", where we have types of building damage (initial values are in oval brackets, array Damage has zero initial values): • Array Types of building

Infrastructure
Based on the work of Ghobarah and Saatcioglu [127] and Valencia et al. [128], this module represents the damage of tsunamis to traffic communications (see Figure 8).Then, the wave damage to water pipes and poles with power cables is modeled here.The module also shows how communications and cable poles are repaired after the wave leaves with units of kilometers or meters.This module is influenced by the module The formation of the Tsunami and the Finance module.module also shows how communications and cable poles are repaired after the wave leaves with units of kilometers or meters.This module is influenced by the module The formation of the Tsunami and the Finance module.
There are three reservoirs in these modules: communications (Roads, Railways), water pipes (plastic, copper, metal), and power cables.Regarding communications, there are three repositories, namely Communications, Destroyed communications, and Repaired communications with units of kilometers.While initial values of Roads and Railways are set to 500 and 350, respectively, Destroyed communications and Repaired communications are set to 0 as the initial condition before the tsunami's arrival.Similar to communications, water pipes consist of water pipes and destroyed water pipes with an array applied to different types of pipes, namely plastic, copper, and metal.Initial values are set to 10,000, 20,000, and 15,000, respectively.Destroyed water pipes have the initial value set to 0. Power cables contain three stocks: Power cables, Destroyed power cables, and New poles and cables.A two-dimensional array is used here: Pole material and height.Materials considered are wood, concrete, and metal.Height is divided into groups lower than 5 m, from 5 to 10 m, and higher than 10 m.With respect to arrays, initial values are [30,000, 25,000, 15,000; 200,000, 15,000, 3000; 120,000, 30,000, 20,000].Initial values of Destroyed power cables and New power cables are set to 0 again.There are three reservoirs in these modules: communications (Roads, Railways), water pipes (plastic, copper, metal), and power cables.Regarding communications, there are three repositories, namely Communications, Destroyed communications, and Repaired communications with units of kilometers.While initial values of Roads and Railways are set to 500 and 350, respectively, Destroyed communications and Repaired communications are set to 0 as the initial condition before the tsunami's arrival.Similar to communications, water pipes consist of water pipes and destroyed water pipes with an array applied to different types of pipes, namely plastic, copper, and metal.Initial values are set to 10,000, 20,000, and 15,000, respectively.Destroyed water pipes have the initial value set to 0. Power cables contain three stocks: Power cables, Destroyed power cables, and New poles and cables.A two-dimensional array is used here: Pole material and height.Materials considered are wood, concrete, and metal.Height is divided into groups lower than 5 m, from 5 to 10 m, and higher than 10 m.With respect to arrays, initial values are [30,000, 25,000, 15,000; 200,000, 15,000, 3000; 120,000, 30,000, 20,000].Initial values of Destroyed power cables and New power cables are set to 0 again.

Finance
In this module, selected financial issues are modeled (see Figure 9).It demonstrates the amount of pledged financial resources and investments into the impacted area.Moreover, distribution to various areas is also included.The work of Heger and Neumayer [129] and Kweifio-Okai [130] is used to build particular model elements and relationships.There are three main modules interacting with Finance, namely Building, Environment, and People.In this module, selected financial issues are modeled (see Figure 9).It demonstrates the amount of pledged financial resources and investments into the impacted area.Moreover, distribution to various areas is also included.The work of Heger and Neumayer [129] and Kweifio-Okai [130] is used to build particular model elements and relationships.There are three main modules interacting with Finance, namely Building, Environment, and People.This module has four reservoirs: Financial support, Finances of the affected state, Finance for damages, and Finance divided into areas.The first reservoir contains a value representing all the funds that have been raised to support the affected area from sources other than the state, such as contributions from other states, charities, or individuals.The second reservoir contains the finances of the affected state.Finances for damages add funds from the state and funds from other sources.The last reservoir contains an array of This module has four reservoirs: Financial support, Finances of the affected state, Finance for damages, and Finance divided into areas.The first reservoir contains a value representing all the funds that have been raised to support the affected area from sources other than the state, such as contributions from other states, charities, or individuals.The second reservoir contains the finances of the affected state.Finances for damages add funds from the state and funds from other sources.The last reservoir contains an array of particular areas, which are Food, Infrastructure renewal, Water, Healthcare, Shelters, and Defense elements.

Defense Elements
The Defense element module consists of five sectors focusing on the most essential tools used in practice.All of these sectors overtake initial data from the module The formation of the Tsunami.This model is not directly linked to other modules.However, connection to other modules can be considered.For instance, an apparent connection can be made with the Finance module.Due to the illustrative nature of this model, further development in this direction is not applied.Although the modules are not interconnected, the financial side is also, at least, tentatively mentioned in the individual sections.

•
Self-elevating Seawalls development in this direction is not applied.Although the modules are not interconnected, the financial side is also, at least, tentatively mentioned in the individual sections. •

Self-elevating Seawalls
Figure 10 presents the structure of the variables associated with self-elevating seawalls.Based on experiments and numerical analysis in OpenFOAM, the equations provided in Appendix C were used for calculations.The hydraulic head between the sea and the port side, where the gate is located, is represented in Equation 19in Appendix C by the letter h, and this value is set at 2 m.The experimental coefficient value is set at 0.0354 to best determine the velocity in the so-called seepage flow, which is the flow of fluid (water) in the permeable layers, in our case, a rubble wall composed of stones with a diameter of 50 cm.
Next, Surge force per unit width of the wall is reckoned, while the density of seawater is The hydraulic head between the sea and the port side, where the gate is located, is represented in Equation 19in Appendix C by the letter h, and this value is set at 2 m.The experimental coefficient value is set at 0.0354 to best determine the velocity in the so-called seepage flow, which is the flow of fluid (water) in the permeable layers, in our case, a rubble wall composed of stones with a diameter of 50 cm.
Next, Surge force per unit width of the wall is reckoned, while the density of seawater is 1.0273 kg m 3 and the surge height is 30 m.The time series of water levels could be represented by the Gaussian distribution, which can reproduce the rapid rise and subsequent fall of the water level η.Time is determined in seconds, and rise time indicates the time required to reach the top of the water level.
Flow discharge via the gap was calculated by taking advantage of Torricelli's theorem.It is set as a condition because four different situations can occur, depending on the gap width and discharge coefficient setting, which must be lower than 1.The hydraulic head between the seaside and the port side of the gate is again represented by the letter h.The friction loss within the gate gap is set at 2 m.The water depth in the gate gap is set at 4.5 m.The last important value is the correction factor.
Friction loss within the gate gap can be calculated using the Darcy-Welsbach formula.This value should be small when the gate width is narrow.The condition below is based on Torricelli's theorem.In any case, the discharge coefficient was derived purely from numerical analysis.Therefore, The correction factor aims to adjust the formula by incorporating other energy-dissipating mechanisms that cannot fully reproduce the hydrodynamic model Takagi et al. [131].
The discharge coefficient balances tsunami momentum and energy loss, leading to reduced flow, friction, shrinkage, expansion, and swirling.The larger the coefficient, the more momentum prevails, facilitating the flow of water through the gap.The regression between the discharge coefficient and the hydraulic head was derived for w = 10, 15, 20, 30 cm as follows.
The model sets the gap width and discharge coefficient as an array.We can also determine whether the wall is functional or not based on the idea that the width of the gap in the floating breakwater system should be less than 3% of the total port breakwater section Nakashima et al. [132].If applied condition 19 is met, the breakwater is functional.However, no technical procedure has been proposed to estimate the tsunami influx.The flow through such a constriction is usually complicated, so the resulting flow pattern is not easily subjected to any analytical solution.
The price is only estimated based on the Kamaishi Protection Breakwater in Japan, which is the deepest in the world.The price for 1 m 2 = USD 8.954 was determined and based on data from a hypothetical wall in a model for which we know the height, height, and length.So first, we calculate the size of the seawall size = height × lenght.We can use the formula cost o f the wall = size o f the wall × 8.954, for calculation of the cost of the wall, a hypothetical price in U.S. dollars.

• Breakwaters
The breakwaters sector, whose structure is presented in Figure 11, requires a wave height as an input value.Data for this particular model are taken from Ooya Harbor.Hudson's formula [133,134] was used as a starting point during past tsunami events to analyze the stability of the defensive breakwater.According to this formula, the weight of the required armor is proportional to the wave height of the proposed incident.The density of armor has a value of about 133 × 1000 kg/m 3 .In contrast, the relative underwater density of armor is only 15 × 1000 kg/m 3 .The slope of the structure was set to 30 rad for this situation.It should be noted that using the Empirically determined damage coefficient value is exposed to wind waves that do not exceed the peak for rubble structures.Therefore, the way the value is used here is not what it was intended for (i.e., for very long periods of waves exceeding rubble structures and breakwaters).Nevertheless, the armor units will benefit from the coupling effect when they resist the forces acting on them due to the tsunami currents.Without any better action, Esteban et al. [135] suggested that the Empirically determined damage coefficient values could be used, although it is clear that further research is needed on this issue.Unlike formulas such as Van der Meer [136], Hudson's formula does not provide an indication of the degree of damage that can be expected as a result of an event.However, in order to quantify the strength of the structure, the armor damage was measured by a factor S similar to that used by Van der Meer.The ratio was defined as the ratio between the weight of the required armor and the actual weight of the armor.

Vertical deep barrier
The vertical deep barrier sector with the structure presented in Figure 12 requires the input values water depth and gravitational acceleration.Pressure waves from earthquakes and landslides rebound from stable walls and then propagate back into the ocean [137].High tsunamis develop only at depths of less than about 500 m or even 200 m.All calculations are based on the work of Levin and Nosov [138].An array is inserted in the amplitude or wave height for easy switching between sample values.We can also determine if the barrier is effective and efficient based on a simple condition   =   ℎℎ Unlike formulas such as Van der Meer [136], Hudson's formula does not provide an indication of the degree of damage that can be expected as a result of an event.However, in order to quantify the strength of the structure, the armor damage was measured by a factor S similar to that used by Van der Meer.The ratio was defined as the ratio between the weight of the required armor and the actual weight of the armor.

• Vertical deep barrier
The vertical deep barrier sector with the structure presented in Figure 12 requires the input values water depth and gravitational acceleration.Pressure waves from earthquakes and landslides rebound from stable walls and then propagate back into the ocean [137].High tsunamis develop only at depths of less than about 500 m or even 200 m.All calculations are based on the work of Levin and Nosov [138].Unlike formulas such as Van der Meer [136], Hudson's formula does not provide an indication of the degree of damage that can be expected as a result of an event.However, in order to quantify the strength of the structure, the armor damage was measured by a factor S similar to that used by Van der Meer.The ratio was defined as the ratio between the weight of the required armor and the actual weight of the armor.

Vertical deep barrier
The vertical deep barrier sector with the structure presented in Figure 12 requires the input values water depth and gravitational acceleration.Pressure waves from earthquakes and landslides rebound from stable walls and then propagate back into the ocean [137].High tsunamis develop only at depths of less than about 500 m or even 200 m.All calculations are based on the work of Levin and Nosov [138].An array is inserted in the amplitude or wave height for easy switching between sample values.We can also determine if the barrier is effective and efficient based on a simple condition   =   ℎℎ a   ℎℎ  1  0 .We can also roughly determine the price based on the tsunami barrier height.An array is inserted in the amplitude or wave height for easy switching between sample values.We can also determine if the barrier is effective and efficient based on a simple condition IF e f f iciency = (tsunami barrier height > amplitude or wave height ) THEN 1 ELSE 0. We can also roughly determine the price based on the tsunami barrier height.

• Other defense alternatives
Two additional Defense elements were developed only for analytical purposes; their integration into the model is not necessary.These alternatives are the Caisson breakwater (see Figure 13) and the influence of coastal vegetation (see Figure 14).

Environment
This module shows how the tsunami will affect nature on the coast and at sea.There is also a model of drinking water contamination or sea pollution by various wastes.This module is affected by the module The formation of the Tsunami.This module (see Figure 15) is developed based on the work of Srinivas and Nakagawa [139].
There are four reservoirs in this module, namely Dead animals, Contaminated water, Destroyed coastal and sea life, and Waste in the sea.The values of all reservoirs are set to zero at the beginning of the simulation.The meanings of single reservoirs are self-explanatory.There are two arrays applied.In Destroyed coastal and sea life, types of vegetation, such as coral reefs, marine plants, mangrove forests, and coastal vegetation in cubic meters are used.In Waste in the sea, an array of specific types of waste is applied, namely (general) waste, soil, and debris.Waste is measured in kilograms and indicates the number of different types of waste that enter the sea when the wave recedes.This module shows how the tsunami will affect nature on the coast and at sea.There is also a model of drinking water contamination or sea pollution by various wastes.This module is affected by the module The formation of the Tsunami.This module (see Figure 15) is developed based on the work of Srinivas and Nakagawa [139].

Model Simulation
Due to the main purpose of the model and its illustrative nature, simulations were not used for the analysis or for finding a solution for a specific case or event.The aim was to show the possibilities of system dynamics in connecting various disciplines into one multidisciplinary descriptive mechanism.Therefore, all tests usually executed on the available model, such as what-if or sensitivity analysis, were not conducted.However, the model successfully passed the robustness test under extreme conditions and structural and behavioral tests.
The initial model parametrization is based on values indicated in previous sections presenting particular modules and sectors.Further parametrization is based on data associated with a specific tsunami wave from 2004 in the Banda Aceh area of Indonesia.Further settings can be found in the model, which is included in the Appendices A-C associated with this manuscript.There is one principal issue associated with the simulation.The model contains processes that take place in various time units.Simple unification is not possible in one model.We can either develop two separate models based on specific time units or use multiple simulations with simultaneous switching off of selected sectors.However, the solution is feasible.The one-hour step was selected as the time unit for the simulations used for all modules except the Formation of tsunami and Defense elements.These two modules use the one-second step.The following figures demonstrate the dynamics of the selected variable in each module.Figure 16 shows the situation of the population 100 h after the start of the simulation.The population was first divided into affected and unaffected populations.The affected population became either injured, uninjured, or dead.Figure 17 shows the size of the destruction in terms of communications and water pipe damage.As determined by the model's initial parametrization, roads are being repaired in terms tens of hours while the reparation of water pipes is postponed.Figures 18 and 19 present the outputs of the sensitivity analysis.Figure 18 captures the barrier utility at different levels of barrier effectiveness and different levels of potential wave height at impact.The right part of the graph (shown in green and separated by the black line) captures usefulness greater than 91%.The black line (on the right side) defines the region of usefulness between 89% and 91%.The part of the graph between the black lines describes the condition where the barrier application is useful.The left part of the graph (to the left black line) shows the situation where the application of the barrier is counterproductive.Figure 19 captures the change in the utility of barrier application relative to the effectiveness of the barrier and the potential wave height.Further settings can be found in the model, which is included in the Appendixes A-C associated with this manuscript.There is one principal issue associated with the simulation.The model contains processes that take place in various time units.Simple unification is not possible in one model.We can either develop two separate models based on specific time units or use multiple simulations with simultaneous switching off of selected sectors.However, the solution is feasible.The one-hour step was selected as the time unit for the simulations used for all modules except the Formation of tsunami and Defense elements.These two modules use the one-second step.The following figures demonstrate the dynamics of the selected variable in each module.Figure 16 shows the situation of the population 100 h after the start of the simulation.The population was first divided into affected and unaffected populations.The affected population became either injured, uninjured, or dead.

Study Limitations
This study has several principal limitations, which can be clustered into two segments.The first segment contains issues related to the literature review and analysis that was conducted.The second one is associated with the developed model.As for the former, the first limitation is the incomplete identification of all relevant articles.We were not able to circumvent this limitation because of the use of the term "tsunami" even in a field where it is used in a context other than its original meaning, an example being "obesity tsunami" or "addiction tsunami", etc.Based on the use of this word, we had to add an auxiliary relevant keyword to remove redundant articles.This step may have caused not all relevant articles to be found.Another limitation is in the search for software used, where not all studies mention the name of the software used and/or do not use the word "software" next to the software name.We were able to partially circumvent this limitation by using a double-machine full-text search, wherein the first step, the retrieved articles were searched for the term "software", followed by recording the names of the reported software and then performing a full-text search for the names of the retrieved software.Another limitation was that not all articles were found in the full-text version.Additional steps may be taken in future replications or extensions of this study to obtain a comprehensive pool of articles.Moreover, additional resources and paper repositories are available.Searching their content, on the other hand, would result in more redundancy in acquired papers rather than new discoveries.Regarding the model itself, we need to highlight that the purpose of the presented model is to demonstrate the possibility of integrating various disciplines in the investigation of complex multidisciplinary phenomena.Thus, particular modules and subsystems may seem to be incomplete or need to be elaborated.Indeed, they are, and they do need.The main intention is to outline how tsunami research can be at least partially consolidated.Elaboration would lead to a higher density of interconnections, which does not have to support comprehensibility and meaningfulness.Although the presented study is comprehensive, the list of tools, techniques,

Study Limitations
This study has several principal limitations, which can be clustered into two segments.The first segment contains issues related to the literature review and analysis that was conducted.The second one is associated with the developed model.As for the former, the first limitation is the incomplete identification of all relevant articles.We were not able to circumvent this limitation because of the use of the term "tsunami" even in a field where it is used in a context other than its original meaning, an example being "obesity tsunami" or "addiction tsunami", etc.Based on the use of this word, we had to add an auxiliary relevant keyword to remove redundant articles.This step may have caused not all relevant articles to be found.Another limitation is in the search for software used, where not all studies mention the name of the software used and/or do not use the word "software" next to the software name.We were able to partially circumvent this limitation by using a double-machine full-text search, wherein the first step, the retrieved articles were searched for the term "software", followed by recording the names of the reported software and then performing a full-text search for the names of the retrieved software.Another limitation was that not all articles were found in the full-text version.Additional steps may be taken in future replications or extensions of this study to obtain a comprehensive pool of articles.Moreover, additional resources and paper repositories are available.Searching their content, on the other hand, would result in more redundancy in acquired papers rather than new discoveries.Regarding the model itself, we need to highlight that the purpose of the presented model is to demonstrate the possibility of integrating various disciplines in the investigation of complex multidisciplinary phenomena.Thus, particular modules and subsystems may seem to be incomplete or need to be elaborated.Indeed, they are, and they do need.The main intention is to outline how tsunami research can be at least partially consolidated.Elaboration would lead to a higher density of interconnections, which does not have to support comprehensibility and meaningfulness.Although the presented study is comprehensive, the list of tools, techniques, methods, or software toolkits cannot be fully exhausted.Various software applications are used for particular and specific aspects of tsunami-related science [140].For instance, the topic of the fragility of constructions due to the tsunami impact (e.g., [141][142][143][144][145][146][147][148]) can serve as an example, which is not explored much in this study.However, the goal of this paper was not to create a comprehensive picture of tsunami research but only the part where tsunami-related software and modeling methods are used to investigate tsunamis.In future research, we could compare approaches to tsunami research using software and "software-free" tsunami research.The most significant challenge is related to the time units used in the model.There are processes that are more convenient to simulate in seconds or minutes, while other processes take place in the order of hours or days.Unification in the system dynamics model can be achieved, for instance, by developing the simulation of two separate models.This represents the main research challenge in this domain.

Conclusions
With advances in computing power and software for modeling and simulation, the capabilities in the vast majority of scientific fields have advanced [149], and this is also true in the field of tsunami research [150].In our study, we identified the most widely used software and methods in the tsunami field over the last decade.We also identified the sectors in the tsunami field where software and methods are used.Furthermore, we have assigned each software and method their function for the analyzed article.The most significant software packages in the tsunami research we identified include OpenFOAM, CALIB, MATLAB, ArcGIS, and COMCOT.The methods we have identified belong to various fields of study and are mostly focused on the tsunami origin or its propagation, exploring historical tsunamis based on tsunami deposits, modeling tsunamis in 3D space, identifying tsunami impacts, exploring relevant variables for tsunamis, creating tsunami impact maps, and comparing simulation results with real data.Various methods are applied, such as the Accelerator Mass Spectrometry, Computational Fluid Dynamics, Constant Rate of Supply Method, Digital Elevation Models, Finite Area Method, Fuzzy C-Means Clustering, and probabilistic tsunami hazard analysis.The importance of finding a tool that would not be over-specialized and narrowly focused on particular aspects of tsunamis needs to be emphasized as the main added value of this study.It is apparent that the existing specialized tools or methods used in particular fields of study need to be applied to find answers to particular domain-related research questions.However, due to the multidisciplinary nature of tsunamis, the more complex systems approach needs to be applied.The holistic view and feedback structures with mutually interconnected parts represent its primary attributes.Data or information acquired by tools, techniques, or modeling software can be used as inputs to more complex models, which would improve understanding and insights into tsunami phenomena.No study used or mentioned the possibility of applying system dynamics as a methodological tool, which focuses on capturing change and behavior over time.This approach is applied in various disciplines.This study demonstrates that the application of system dynamics as a commonly used modeling and simulation methodology can be successfully used for the implementation of the system perspective in tsunami research.To apply sensitivity analysis and propagation of uncertainties of cross-sections; to conduct sensitivity and uncertainty analysis 2 The TSUNAMI-1D, TSUNAMI-2D, and TSUNAMI-3D analysis sequences compute the sensitivity of keff and reaction rates to energy-dependent cross-section data for each reaction of each nuclide in a system model.The one-dimensional (1D) transport calculations are performed with XSDRNPM, the two-dimensional (2D) transport calculations are performed using NEWT, and the three-dimensional (3D) calculations are performed with KENO V.a or KENO-VI.To compare results with presented numerical solver; to simulate tsunami 2 A finite-volume nonlinear shallow water equation (NSWE) solver built on the OP2 domain-specific language (DSL) for unstructured mesh computations.VOLNA-OP2 is unique among tsunami solvers in its support for several high-performance computing platforms: central processing units (CPUs), the Intel Xeon Phi, and graphics processing units (GPUs).

Figure 2 .
Figure 2. The interconnectedness of methods and techniques.Note: This figure shows the interconnectedness of the keywords mentioned in the analyzed articles.In the figure, there is one larger cluster mainly focused on modeling and many smaller clusters that refer to the specifics of different models (own work; software: VOSviewer [16]).

Figure 2 .
Figure 2. The interconnectedness of methods and techniques.Note: This figure shows the interconnectedness of the keywords mentioned in the analyzed articles.In the figure, there is one larger cluster mainly focused on modeling and many smaller clusters that refer to the specifics of different models (own work; software: VOSviewer [16]).

Figure 3 .
Figure 3.The interconnectedness of methods, tools, and SW applications with their tools.Note: This figure shows the interconnectedness of using the different tsunami modeling tools.Unlike the previous figure focusing on keywords, the individual tools are quite significantly interconnected (own work; software: VOSviewer [16]).

Figure 3 .
Figure 3.The interconnectedness of methods, tools, and SW applications with their tools.Note: This figure shows the interconnectedness of using the different tsunami modeling tools.Unlike the previous figure focusing on keywords, the individual tools are quite significantly interconnected (own work; software: VOSviewer [16]).

Figure 5 .
Figure 5. Part of SD model: tsunami formation, variants of formation.Note: This image shows diagram sections comprising other possible tsunami sources.Compared to Figure 4a, these are rather minor in terms of importance and frequency of occurrence (own work; software: Stella Professional).

Figure 4 .Figure 4 .
Figure 4. (a) Part of SD model: tsunami formation, seismic activity.Note: This part of the diagram shows the formation of a tsunami from seismic activity (own work; software: Stella Professional).(b) Part of SD model: tsunami formation.Note: This segment of the diagram focuses on wave formation The area dealing with wave height shows that it is affected by water depth, wave amplitude, and wave propagation speed.The section dedicated to wave run-up illustrates how run-up is influenced by ground slope, surface roughness, and the loss in wave height per meter of inundation distance.The part of the diagram concerning waves behind a barrier includes factors like barrier level and barrier efficiency.This segment of the diagram focuses on tsunami simulation and includes parameters such as assumed water depth, assumed tsunami height, and wave distribution after propagating over various terrains (own work; software: Stella Professional).

Figure 5 .
Figure 5. Part of SD model: tsunami formation, variants of formation.Note: This image shows diagram sections comprising other possible tsunami sources.Compared to Figure 4a, these are rather minor in terms of importance and frequency of occurrence (own work; software: Stella Professional).

Figure 5 .
Figure 5. Part of SD model: tsunami formation, variants of formation.Note: This image shows diagram sections comprising other possible tsunami sources.Compared to Figure 4a, these are rather minor in terms of importance and frequency of occurrence (own work; software: Stella Professional).
further divided into those who are injured, dead, or uninjured into the Injured people, Dead people, and Unharmed people reservoirs.Rescuers is another reservoir where we have a number of rescuers in the area.Rescuers are then also divided into uninjured, wounded, and dead through model outflows.The last two reservoirs concern paramedics.The first is Paramedics in the state, where we have the number of all medics available, and the second is Paramedics in the area, where we have the number of all medics in a given location affected by the wave.

Figure 6 .
Figure 6.Part of SD model: impact on population.Note: This figure captures how the affected population is influenced by the impact rate and population density in the affected area.Mortality rate and injury rate determine the number of deaths and injured people, respectively.The diagram also includes the role of paramedics, showing how their assistance is divided between the unaffected and affected areas.Finally, the environment, specifically contaminated water, plays a role in the health outcomes of both injured people and rescuers, contributing to fatal injuries and deaths (own work; software: Stella Professional).

Figure 6 .
Figure 6.Part of SD model: impact on population.Note: This figure captures how the affected population is influenced by the impact rate and population density in the affected area.Mortality rate and injury rate determine the number of deaths and injured people, respectively.The diagram also includes the role of paramedics, showing how their assistance is divided between the unaffected and affected areas.Finally, the environment, specifically contaminated water, plays a role in the health outcomes of both injured people and rescuers, contributing to fatal injuries and deaths (own work; software: Stella Professional).Sci 2024, 6, x FOR PEER REVIEW 21 of 47

Figure 7 .
Figure 7. Part of SD model: impact on buildings.Note: This segment of the diagram focuses on the impact of a tsunami on buildings, specifically considering the assumed maximum tsunami height barrier.The flow depth and intensity of destruction affect the average level of damage to buildings.Damaged buildings are then categorized based on the rate of destruction (own work; software: Stella Professional).

Figure 7 .
Figure 7. Part of SD model: impact on buildings.Note: This segment of the diagram focuses on the impact of a tsunami on buildings, specifically considering the assumed maximum tsunami height barrier.The flow depth and intensity of destruction affect the average level of damage to buildings.Damaged buildings are then categorized based on the rate of destruction (own work; software: Stella Professional).

Figure 8 .
Figure 8. Part of SD model: impact on infrastructure.Note: This figure shows the infrastructure that a wave can damage or destroy.Several types of materials for each part can suffer from the tsunami impact.For example, water pipes can be plastic, copper, or steel (own work; software: Stella Professional).

Figure 8 .
Figure 8. Part of SD model: impact on infrastructure.Note: This figure shows the infrastructure that a wave can damage or destroy.Several types of materials for each part can suffer from the tsunami impact.For example, water pipes can be plastic, copper, or steel (own work; software: Stella Professional).

Figure 9 .
Figure 9. Part of SD model: financial impact.Note: This figure shows the financial aspect of the model, namely the income side.It considers state support and international financial aid, which usually comes with a delay (own work; software: Stella Professional).

Figure 9 .
Figure 9. Part of SD model: financial impact.Note: This figure shows the financial aspect of the model, namely the income side.It considers state support and international financial aid, which usually comes with a delay (own work; software: Stella Professional).

Figure 10
Figure 10 presents the structure of the variables associated with self-elevating seawalls.Based on experiments and numerical analysis in OpenFOAM, the equations provided in Appendix C were used for calculations.

Figure 10 .
Figure 10.Part of SD model: self-elevating seawalls.Note: This figure exhibits the mechanism of self-elevating seawalls.This part of the model can be switched off as required.This option reflects the requirements of reality, where not every area can have this type of barrier (own work; software: Stella Professional).

Figure 10 .
Figure 10.Part of SD model: self-elevating seawalls.Note: This figure exhibits the mechanism of self-elevating seawalls.This part of the model can be switched off as required.This option reflects the requirements of reality, where not every area can have this type of barrier (own work; software: Stella Professional).

Figure 12 .
Figure 12.Part of SD model: vertical deep barrier (own work; software: Stella Professional).

Figure 12 .
Figure 12.Part of SD model: vertical deep barrier (own work; software: Stella Professional).

Figure 12 .
Figure 12.Part of SD model: vertical deep barrier (own work; software: Stella Professional).

Figure 14 .
Figure 14.Part of SD model: influence of coastal vegetation.Note: This figure shows the impact of the wave on the vegetation and vice versa.There are several types of vegetation in this part of the model, and different vegetation parameters are assumed (own work; software: Stella Professional).

Figure 14 .
Figure 14.Part of SD model: influence of coastal vegetation.Note: This figure shows the impact of the wave on the vegetation and vice versa.There are several types of vegetation in this part of the model, and different vegetation parameters are assumed (own work; software: Stella Professional).

Figure 15 .
Figure 15.Part of SD model: impact on fauna and flora.Note: This part of the model shows the impact of the wave on fauna and flora.The amount of waste moved into the sea and the amount of marine life destroyed are examples of the main indicators monitored here (own work; software: Stella Professional).There are four reservoirs in this module, namely Dead animals, Contaminated water, Destroyed coastal and sea life, and Waste in the sea.The values of all reservoirs are set to zero at the beginning of the simulation.The meanings of single reservoirs are self-explanatory.There are two arrays applied.In Destroyed coastal and sea life, types of vegetation, such as coral reefs, marine plants, mangrove forests, and coastal vegetation in cubic meters are used.In Waste in the sea, an array of specific types of waste is applied, namely (general) waste, soil, and debris.Waste is measured in kilograms and indicates the number

Figure 15 .
Figure 15.Part of SD model: impact on fauna and flora.Note: This part of the model shows the impact of the wave on fauna and flora.The amount of waste moved into the sea and the amount of marine life destroyed are examples of the main indicators monitored here (own work; software: Stella Professional).

Figure 17
shows the size of the destruction in terms of communications and water pipe damage.As determined by the model's initial parametrization, roads are being repaired in terms tens of hours while the reparation of water pipes is postponed.Figures18 and 19present the outputs of the sensitivity analysis.Figure18captures the barrier utility at different levels of barrier effectiveness and different levels of potential wave height at impact.The right part of the graph (shown in green and separated by the black line) captures usefulness greater than 91%.The black line (on the right side) defines the region of usefulness between 89% and 91%.The part of the graph between the black lines describes the condition where the barrier application is useful.The left part of the graph (to the left black line) shows the situation where the application of the barrier is counterproductive.Figure19captures the change in the utility of barrier application relative to the effectiveness of the barrier and the potential wave height.

Figure 16 .
Figure 16.Temporal dynamics of casualties and emergency response following a tsunami.Note: The vertical axis is monitored in units of People; all three variables can reach values between predefined minima and maxima (Injured people min = 0, max = 80,000; Dead people min = 0, max = 500; Paramedics in the affected area min = 300, max = 1100) (own work; software: Stella Professional).

Figure 16 .
Figure 16.Temporal dynamics of casualties and emergency response following a tsunami.Note: The vertical axis is monitored in units of People; all three variables can reach values between predefined minima and maxima (Injured people min = 0, max = 80,000; Dead people min = 0, max = 500; Paramedics in the affected area min = 300, max = 1100) (own work; software: Stella Professional).

Figure 18 .
Figure 18.Simulation of barrier effectiveness in wave height reduction (own work; software: MatLab).

Figure 18 .
Figure 18.Simulation of barrier effectiveness in wave height reduction (own work; software: MatLab).Figure 18. Simulation of barrier effectiveness in wave height reduction (own work; software: MatLab).

Figure 18 .
Figure 18.Simulation of barrier effectiveness in wave height reduction (own work; software: MatLab).Figure 18. Simulation of barrier effectiveness in wave height reduction (own work; software: MatLab).

Figure 19 .
Figure 19.Impact of barrier efficiency on wave height attenuation (own work; software: MatLab).

Figure 19 .
Figure 19.Impact of barrier efficiency on wave height attenuation (own work; software: MatLab).

Table 1 .
SW tools and research methods were applied to analyzed studies.

Table 1 .
Cont.To address the lessons learned from Storm Xynthia in February 2010, which caused a significant sea surge along the French Atlantic coast, flooding low-lying coastal areas, particularly urbanized regions.To assess the impact of sea-level rise and wave-driven flooding on seabird colonies in the Pacific, mainly focusing on a globally important seabird rookery at Midway Atoll in the subtropical Pacific.
To reflect on the progress of tsunami preparedness in a coastal community in Aceh, Indonesia, nearly two decades after the catastrophic 2004 Indian Ocean tsunami.Rauter et al., 2021 [97] OpenFOAM CFD, PISO, MULES To evaluate the performance of the multiphase Navier-Stokes equations implemented in OpenFOAM for simulating impulse wave generation by landslides.Nemati et al., 2023 [98] ArcGIS, MATLAB DEM To report the results of numerical simulations for a potential subaerial landslide on the coast of Orcas Island and the resultant tsunami waves in the southern Strait of Georgia.Guo and Lo, 2022 [99] OpenFOAM CFD, PIMPLE, VOF To investigate the hydrodynamics of a solitary wave passing a vertical cylinder over a viscous mud bed for the first time.Attili et al., 2021 [100] OpenFOAM FVM, PIMPLE, VOF To focus on the numerical modeling of landslide tsunamis impacting dams.To perform a probabilistic hazard analysis of a tsunami generated by a subaqueous volcanic explosion at Taal Lake, located on Luzon Island in the Philippines.Song et al., 2023 [102] OpenFOAM VOF To investigate the impacts of tsunami-like waves on coastal bridge decks with superelevation, addressing a gap in existing research that typically focuses on wave impacts on flat bridge decks.Elsheikh et al., 2022 [103] OpenFOAM CFD, FVM To investigate the hydrodynamics of turbulent bores that propagate on a horizontal plane, closely resembling dam-break waves and tsunami-like hydraulic bores.Rahuman et al., 2022 [104] ANSYS CFD To visualize and compare the fluid flow patterns around the Rhizophora mangrove species' stilt roots and the Avicennia mangrove species' pneumatophore roots in the Pichavaram mangrove forest.To simulate the formation of a weak layer in the mountainous slope leading to the Taan Fiord landslide and to analyze the triggering factors from a geotechnical engineering perspective.Madden et al., 2023 [111] GeoClaw DEM To leverage Google's Tensor Processing Unit to rapidly evaluate different tsunami risk mitigation strategies, making high-performance computing accessible to communities.Yuan et al., 2021 [112] COMCOT PTHA, NSWE, LSWE To conduct a PTHA for mainland China and Taiwan Island.Celikbas et al., 2023 [113] ArcGIS DEM

Table 1 .
SW tools and research methods were applied to analyzed studies.

Table 2 .
Application of SW tools in particular research domains.

Table A1 .
Cont.GAMIT, GLOBK, and TRACK form a comprehensive suite of programs for analyzing GNSS measurements primarily to study crustal deformation.The software has been developed by MIT, Scripps Institution of Oceanography, and Harvard University with support from the National Science Foundation.The software may be obtained without written agreement or royalty fee by universities and government agencies for any non-commercial purposes.It is a variant of Clawpack for geophysical flows.Clawpack is a collection of finite volume methods for linear and nonlinear hyperbolic systems of conservation laws.Clawpack employs high-resolution Godunov-type methods with limiters in a general framework applicable to many kinds of waves.Clawpack is written in Fortran and Python.GeoWave is a software library that connects the scalability of distributed computing frameworks and key-value stores with modern geospatial software to store, retrieve, and analyze massive geospatial datasets.Hazus Program provides standardized tools and data for estimating risk from earthquakes, floods, tsunamis, and hurricanes.Hazus models combine expertise from many disciplines to create actionable risk information that increases community resilience.Hazus software is distributed as a GIS-based desktop application with a growing collection of simplified open-source tools.Risk assessment resources from the Hazus program are always freely available and transparently developed.ICEM Surf is the industry-leading Curve and Surface explicit geometry modeling tool for defining, analyzing, and performing high-end visualization of complex free-form shape CAD surface models to the highest quality.Used in product design processes throughout automotive, aerospace, consumer goods, and press-tool design industries, providing solutions for direct surface modeling, refinement, reconstruction, and scan modeling.OpenFOAM is the free, open-source CFD software developed primarily by OpenCFD Ltd. since 2004.It has a large user base across most areas of engineering and science, from both commercial and academic organizations.OpenFOAM has an extensive range of features to solve anything from complex fluid flows involving chemical reactions, turbulence, and heat transfer to acoustics, solid mechanics, and electromagnetics.It fully supports the DICOM standard for easy integration in your workflow environment and an open platform for the development of processing tools.It offers advanced post-processing techniques in 2D and 3D, exclusive innovative techniques for 3D and 4D navigation, and a complete integration with any PACS.RStudio's mission is to create free and open-source software for data science, scientific research, and technical communication.We do this to enhance the production and consumption of knowledge by everyone, regardless of economic means, and to facilitate collaboration and reproducible research, both of which are critical to the integrity and efficacy of work in science, education, government, and industry.the different geometry assemblings (surface reflection and refraction, borehole crosshole and tomography, and combination of borehole and surface measurements).You may also have a look at a one-sided brochure for GPR, reflections seismics, refraction seismics, and borehole application.SPAD https://ia-data-analytics.com/data-mining-software/, accessed on 15 March 2023 To apply statistical test on data 1 reflexW https://www.sandmeier-geo.de/reflexw.html,accessed on 15 March 2023 To process Ground-Penetrating Radar data 1 he software covers the complete range of wave data (seismic, GPR, ultrasound) and

Table A1 .
Cont.Coheris Analytics SPAD is the only software dedicated to Data Mining and Predictive analysis that provides a totally graphical and intuitive interface with powerful features.

Table A2 .
Usage of methods.As a part of OpenFOAM; model of FLOW3D; to solve a three-dimensional Reynolds Averaged Navier-Stokes equations; to investigate solitary wave-induced vertical and horizontal forces on coastal bridges; to predict flow character and dynamic loading profile from an idealized tsunami impact on a coastal community; to simulate impulse wave generation and propagation; to simulate To combine a probability of occurrence of earthquake; to estimate a relative error; to evaluate an overall uncertainty in tsunami hazard;to generate a synthetic earthquake catalog; to sample a resulting posterior; to confirm robustness of the created index; to vary all uncertain input parameters have been randomly within the specified distribution

Table A2 .
Cont.As a based method for enchanted method; for wave propagation from off-to onshore; to simulate waves; to study the nature of flows for an extreme wave above and in the interior of gravel bedformsSWE Shallow Water EquationsAs a base for a model; for modeling a tsunami with small numbers of observation points in more physically realistic settings; to calculate a modelVOF Volume of FluidTo capture a free surface; to obtain sea-level or mudslide interface location; to figure out the role of vegetation of finite width in energy reduction of flood flow; to track the free surface between two fluids

Table A3 .
Relationship among tools and research domains.