babak tehranirad - Academia.edu (original) (raw)
Papers by babak tehranirad
Coastal Structures and Solutions to Coastal Disasters 2015, 2017
Tsunamis are able to move large amounts of sediment during their inundation process. However, pre... more Tsunamis are able to move large amounts of sediment during their inundation process. However, previous studies conducted to assess tsunami hazard have not considered tsunamiinduced sediment transport in their simulations, although it is possible that morphological adjustments during tsunami inundation increase the levels of hazard. In this paper, we present a model which is able to simulate sediment processes caused by tsunami action. The model is validated by comparing results to an experiment of sediment transport under breaking solitary waves. Finally, the model is applied to the simulation of sediment transport over a typical barrier island geometry, and we study the resulting changes in tsunami runup on the mainland behind the barrier, as compared to fixed bathymetry/topography conditions. We show that considering morphological changes during tsunami inundation modeling can increase the runup on the mainland behind the barrier. We conclude that assuming a fixed bathymetry when modeling tsunami coastal hazard may lead to underestimating runup and inundation values and cause errors in estimating the tsunami inundation line.
Coastal Structures and Solutions to Coastal Disasters 2015, 2017
ABSTRACT Landslide generated tsunamis pose a significant hazard over a wide range of spatial scal... more ABSTRACT Landslide generated tsunamis pose a significant hazard over a wide range of spatial scales. At large scale, tsunamis resulting from failures along the continental shelf margins represent a significant component of coastal tsunami hazard. At smaller scales, tsunamis resulting from subaerial or submarine slope failures represent a significant hazard in many steep-walled fjords. A relatively well known tsunami-genic landslide event occured in Kitimat, British Columbia in 1975 (Prior et al., 1983). Multiple large waves associated with the tsunami caused severe damage to infrastructure. Earlier attempts at slide reconstruction and tsunami modeling have only yielded partial agreement with anecdotal observations of wave heights, draw down and run-up (Murty, 1979; Skvortsov and Bornhold, 2007). In this paper, we consider several models for landslide generation of dispersive tsunami waves. In each case, we model the main water column using the 3-D, nonhydrostatic model NHWAVE, described originally by Ma et al (2012). NHWAVE solves the RANS equations in a surface and terrain following σ coordinate system. In perfect fluid applications, NHWAVE has been shown to provide accurate modeling of dispersive surface wave propagation using as few as three vertical σ levels, making it comparable in cost to high-order Boussinesq models. A recent reformulation of terms representing diffusion mechanisms in the model have led to improved descriptions of wave breaking, turbulence and residual flow generation in periodic and transient wave cases (Derakhti et al, 2015a,b). A variety of tsunami generation approaches are described here. Ma et al (2012) described the generation of a tsunami by a non-deforming slide by specifying a given ground motion in the usual kinematic bottom boundary condition for NHWAVE. Resulting tsunami wave response was validated against laboratory measurements by Enet and Grilli (2007), who studied the motion of a freely sliding mass on a plane slope. Experimental and numerical results clearly demonstrate the importance of frequency dispersion in characterizing the studied event. The resulting model has been used to study landslide tsunami generation in a number of settings, including slope failures along the Great Bahama Bank carbonate platform, flanking the Straits of Florida, and dependencies on terminal outrun velocity, slope inclination and failure volume (Schnyder et al., 2013). At an opposite extreme of model complexity, Ma et al (2013) have described a landslide mass as a suspended sediment load accelerating under its own weight. The resulting model requires an extensive increase in vertical resolution in order to resolve the mechanics of the evolving sediment mass, rendering it computationally expensive in comparison to the 2012 model. The model was extensively validated against detailed laboratory measurements of velocity profiles and sediment concentration in gravity currents (Garcia, 1993), and was further used to examine the dependence of generated tsunami waves on a range of parameters.
a b s t r a c t We derive fully nonlinear, weakly dispersive model equations for propagation of s... more a b s t r a c t We derive fully nonlinear, weakly dispersive model equations for propagation of surface gravity waves in a shallow, homogeneous ocean of variable depth on the surface of a rotating sphere. A numerical model is developed for the weakly nonlinear version of the model based on a combined finite-volume and finite-difference method with a fourth-order MUSCL-TVD scheme in space and a third-order SSP Runge–Kutta scheme in time. In the context of tsunami generation and propagation over trans-oceanic distances, a scaling analysis reveals that the importance of frequency dispersion increases with a decrease of the source width, while the effect of the Coriolis force increases with an increase of the source width. A sen-sitivity analysis to dispersive and Coriolis effects is carried out using the numerical model in a series of numerical experiments in an idealized ocean using Gaussian and di-polar sources with different source sizes. A simulation of the Tohoku 2011 tsunami is u...
The Cartesian version of the fully non-linear Boussinesq model FUNWAVE-TVD has been benchmarked i... more The Cartesian version of the fully non-linear Boussinesq model FUNWAVE-TVD has been benchmarked in Tehranirad et al. (2011) for tsunami application using PMEL-135 benchmarks provided by Synolakis et al. (2007). This report presents results from the spherical version of FUNWAVE-TVD (Kirby et al., 2012) using the same set of benchmark tests. The results presented here represent testing of Version 2.0 of the code, and will be updated online at
Recent observations of the coastal impact of large tsunamis (e.g., Indian Ocean 2004; Tohoku 2011... more Recent observations of the coastal impact of large tsunamis (e.g., Indian Ocean 2004; Tohoku 2011) and related numerical and theoretical works have made it increasingly clear that tsunami waves arrive nearshore as a series of long waves (so-called N-waves) with, often, the superposition of undular bores around each crest. Such wave trains are much more complex and very much in contrast with the solitary wave paradigm which for a long time was the accepted idealization of tsunami waves in both experimental and numerical work. The dissipation associated with these breaking bores can be very large, particularly over a wide and shallow continental shelf such as along the east coast of North America, particularly for the shorter waves associated with tsunamis generated by Submarine Mass Failures (SMFs). In this paper, we perform numerical simulations of tsunami coastal impact in the context of both idealized laboratory experiments and several tsunami case studies. We attempt to clarify t...
Recent observations of the coastal impact of large tsunamis (e.g., Indian Ocean 2004; Tohoku 2011... more Recent observations of the coastal impact of large tsunamis (e.g., Indian Ocean 2004; Tohoku 2011) and related numerical and theoretical works have made it increasingly clear that tsunami waves arrive nearshore as a series of long waves (so-called N-waves) with, often, the superposition of undular bores around each crest. Such wave trains are much more complex and very much in contrast with the solitary wave paradigm which for a long time was the accepted idealization of tsunami waves in both experimental and numerical work. The dissipation associated with these breaking bores can be very large, particularly over a wide and shallow continental shelf such as along the east coast of North America, particularly for the shorter waves associated with tsunamis generated by Submarine Mass Failures (SMFs). In this paper, we perform numerical simulations of tsunami coastal impact in the context of both idealized laboratory experiments and several tsunami case studies. We attempt to clarify t...
Pure and Applied Geophysics, 2020
We describe the development of a time-resolved model for tsunami-induced morphology change, based... more We describe the development of a time-resolved model for tsunami-induced morphology change, based on an existing Boussinesq model for weakly-dispersive free-surface waves coupled to a depth-integrated model for sediment concentration and an equation for bed level change. The model allows for spatial variability in bottom friction coefficients and accounts for the presence of non-erodible beds or finite depths of available sediment. The model is verified using one laboratory data set and against field observations of morphology change in the Crescent City, CA harbor during the 2011 Tohoku tsunami event. The model is then applied to a hypothetical example illustrating the impact of tsunami inundation on an undeveloped barrier island.
Scientific reports, Nov 4, 2016
Submarine slope failures are a likely cause for tsunami generation along the East Coast of the Un... more Submarine slope failures are a likely cause for tsunami generation along the East Coast of the United States. Among potential source areas for such tsunamis are submarine landslides and margin collapses of Bahamian platforms. Numerical models of past events, which have been identified using high-resolution multibeam bathymetric data, reveal possible tsunami impact on Bimini, the Florida Keys, and northern Cuba. Tsunamis caused by slope failures with terminal landslide velocity of 20 ms(-1) will either dissipate while traveling through the Straits of Florida, or generate a maximum wave of 1.5 m at the Florida coast. Modeling a worst-case scenario with a calculated terminal landslide velocity generates a wave of 4.5 m height. The modeled margin collapse in southwestern Great Bahama Bank potentially has a high impact on northern Cuba, with wave heights between 3.3 to 9.5 m depending on the collapse velocity. The short distance and travel time from the source areas to densely populated ...
Pure and Applied Geophysics, 2015
Your article is protected by copyright and all rights are held exclusively by Springer Basel. Thi... more Your article is protected by copyright and all rights are held exclusively by Springer Basel. This e-offprint is for personal use only and shall not be self-archived in electronic repositories. If you wish to self-archive your article, please use the accepted manuscript version for posting on your own website. You may further deposit the accepted manuscript version in any repository, provided it is only made publicly available 12 months after official publication or later and provided acknowledgement is given to the original source of publication and a link is inserted to the published article on Springer's website. The link must be accompanied by the following text: "The final publication is available at link.springer.com".
Natural Hazards, 2014
With support from the US National Tsunami Hazard Mitigation Program (NTHMP), the authors have bee... more With support from the US National Tsunami Hazard Mitigation Program (NTHMP), the authors have been developing tsunami inundation maps for the upper US East Coast (USEC), using high-resolution numerical modeling. These maps are envelopes of maximum elevations, velocity, or momentum flux, caused by the probable maximum tsunamis identified in the Atlantic oceanic basin, including from far-field coseismic or volcanic sources, and near-field Submarine mass failures (SMFs); the latter are the object of this work. Despite clear field evidence of past large-scale SMFs within our area of interest, such as the Currituck slide complex, their magnitude, pre-failed geometry, volume, and mode of rupture are poorly known. A screening analysis based on the Monte Carlo simulations (MCS) identified areas for possible tsunamigenic SMF sources along the USEC, indicating an increased level of tsunami hazard north of Virginia, potentially
Coastal Engineering Proceedings, 2012
Recent observations of the coastal impact of large tsunamis (e.g., Indian Ocean 2004; Tohoku 2011... more Recent observations of the coastal impact of large tsunamis (e.g., Indian Ocean 2004; Tohoku 2011) and related numerical and theoretical works have made it increasingly clear that tsunami waves arrive nearshore as a series of long waves (so-called N-waves) with, often, the superposition of undular bores around each crest. Such wave trains are much more complex and very much in contrast with the solitary wave paradigm which for a long time was the accepted idealization of tsunami waves in both experimental and numerical work. The dissipation associated with these breaking bores can be very large, particularly over a wide and shallow continental shelf such as along the east coast of North America, particularly for the shorter waves associated with tsunamis generated by Submarine Mass Failures (SMFs). In this paper, we perform numerical simulations of tsunami coastal impact in the context of both idealized laboratory experiments and several tsunami case studies. We attempt to clarify t...
Ocean Modelling, 2013
We derive fully nonlinear, weakly dispersive model equations for propagation of surface gravity w... more We derive fully nonlinear, weakly dispersive model equations for propagation of surface gravity waves in a shallow, homogeneous ocean of variable depth on the surface of a rotating sphere. A numerical model is developed for the weakly nonlinear version of the model based on a combined finite-volume and finitedifference method with a fourth-order MUSCL-TVD scheme in space and a third-order SSP Runge-Kutta scheme in time. In the context of tsunami generation and propagation over trans-oceanic distances, a scaling analysis reveals that the importance of frequency dispersion increases with a decrease of the source width, while the effect of the Coriolis force increases with an increase of the source width. A sensitivity analysis to dispersive and Coriolis effects is carried out using the numerical model in a series of numerical experiments in an idealized ocean using Gaussian and di-polar sources with different source sizes. A simulation of the Tohoku 2011 tsunami is used to illustrate the effects of dispersive and Coriolis effects at large distances from the source region.
Ocean Modelling, 2012
We present a high-order adaptive time-stepping TVD solver for the fully nonlinear Boussinesq mode... more We present a high-order adaptive time-stepping TVD solver for the fully nonlinear Boussinesq model of Chen (2006), extended to include moving reference level as in Kennedy et al. (2001). The equations are reorganized in order to facilitate high-order Runge-Kutta time-stepping and a TVD type scheme with a Riemann solver. Wave breaking is modeled by locally switching to the nonlinear shallow water equations when the Froude number exceeds a certain threshold. The moving shoreline boundary condition is implemented using the wetting-drying algorithm with the adjusted wave speed of the Riemann solver. The code is parallelized using the Message Passing Interface (MPI) with non-blocking communication. Model validations show good performance in modeling wave shoaling, breaking, wave runup and waveaveraged nearshore circulation.
Coastal Structures and Solutions to Coastal Disasters 2015, 2017
Tsunamis are able to move large amounts of sediment during their inundation process. However, pre... more Tsunamis are able to move large amounts of sediment during their inundation process. However, previous studies conducted to assess tsunami hazard have not considered tsunamiinduced sediment transport in their simulations, although it is possible that morphological adjustments during tsunami inundation increase the levels of hazard. In this paper, we present a model which is able to simulate sediment processes caused by tsunami action. The model is validated by comparing results to an experiment of sediment transport under breaking solitary waves. Finally, the model is applied to the simulation of sediment transport over a typical barrier island geometry, and we study the resulting changes in tsunami runup on the mainland behind the barrier, as compared to fixed bathymetry/topography conditions. We show that considering morphological changes during tsunami inundation modeling can increase the runup on the mainland behind the barrier. We conclude that assuming a fixed bathymetry when modeling tsunami coastal hazard may lead to underestimating runup and inundation values and cause errors in estimating the tsunami inundation line.
Coastal Structures and Solutions to Coastal Disasters 2015, 2017
ABSTRACT Landslide generated tsunamis pose a significant hazard over a wide range of spatial scal... more ABSTRACT Landslide generated tsunamis pose a significant hazard over a wide range of spatial scales. At large scale, tsunamis resulting from failures along the continental shelf margins represent a significant component of coastal tsunami hazard. At smaller scales, tsunamis resulting from subaerial or submarine slope failures represent a significant hazard in many steep-walled fjords. A relatively well known tsunami-genic landslide event occured in Kitimat, British Columbia in 1975 (Prior et al., 1983). Multiple large waves associated with the tsunami caused severe damage to infrastructure. Earlier attempts at slide reconstruction and tsunami modeling have only yielded partial agreement with anecdotal observations of wave heights, draw down and run-up (Murty, 1979; Skvortsov and Bornhold, 2007). In this paper, we consider several models for landslide generation of dispersive tsunami waves. In each case, we model the main water column using the 3-D, nonhydrostatic model NHWAVE, described originally by Ma et al (2012). NHWAVE solves the RANS equations in a surface and terrain following σ coordinate system. In perfect fluid applications, NHWAVE has been shown to provide accurate modeling of dispersive surface wave propagation using as few as three vertical σ levels, making it comparable in cost to high-order Boussinesq models. A recent reformulation of terms representing diffusion mechanisms in the model have led to improved descriptions of wave breaking, turbulence and residual flow generation in periodic and transient wave cases (Derakhti et al, 2015a,b). A variety of tsunami generation approaches are described here. Ma et al (2012) described the generation of a tsunami by a non-deforming slide by specifying a given ground motion in the usual kinematic bottom boundary condition for NHWAVE. Resulting tsunami wave response was validated against laboratory measurements by Enet and Grilli (2007), who studied the motion of a freely sliding mass on a plane slope. Experimental and numerical results clearly demonstrate the importance of frequency dispersion in characterizing the studied event. The resulting model has been used to study landslide tsunami generation in a number of settings, including slope failures along the Great Bahama Bank carbonate platform, flanking the Straits of Florida, and dependencies on terminal outrun velocity, slope inclination and failure volume (Schnyder et al., 2013). At an opposite extreme of model complexity, Ma et al (2013) have described a landslide mass as a suspended sediment load accelerating under its own weight. The resulting model requires an extensive increase in vertical resolution in order to resolve the mechanics of the evolving sediment mass, rendering it computationally expensive in comparison to the 2012 model. The model was extensively validated against detailed laboratory measurements of velocity profiles and sediment concentration in gravity currents (Garcia, 1993), and was further used to examine the dependence of generated tsunami waves on a range of parameters.
a b s t r a c t We derive fully nonlinear, weakly dispersive model equations for propagation of s... more a b s t r a c t We derive fully nonlinear, weakly dispersive model equations for propagation of surface gravity waves in a shallow, homogeneous ocean of variable depth on the surface of a rotating sphere. A numerical model is developed for the weakly nonlinear version of the model based on a combined finite-volume and finite-difference method with a fourth-order MUSCL-TVD scheme in space and a third-order SSP Runge–Kutta scheme in time. In the context of tsunami generation and propagation over trans-oceanic distances, a scaling analysis reveals that the importance of frequency dispersion increases with a decrease of the source width, while the effect of the Coriolis force increases with an increase of the source width. A sen-sitivity analysis to dispersive and Coriolis effects is carried out using the numerical model in a series of numerical experiments in an idealized ocean using Gaussian and di-polar sources with different source sizes. A simulation of the Tohoku 2011 tsunami is u...
The Cartesian version of the fully non-linear Boussinesq model FUNWAVE-TVD has been benchmarked i... more The Cartesian version of the fully non-linear Boussinesq model FUNWAVE-TVD has been benchmarked in Tehranirad et al. (2011) for tsunami application using PMEL-135 benchmarks provided by Synolakis et al. (2007). This report presents results from the spherical version of FUNWAVE-TVD (Kirby et al., 2012) using the same set of benchmark tests. The results presented here represent testing of Version 2.0 of the code, and will be updated online at
Recent observations of the coastal impact of large tsunamis (e.g., Indian Ocean 2004; Tohoku 2011... more Recent observations of the coastal impact of large tsunamis (e.g., Indian Ocean 2004; Tohoku 2011) and related numerical and theoretical works have made it increasingly clear that tsunami waves arrive nearshore as a series of long waves (so-called N-waves) with, often, the superposition of undular bores around each crest. Such wave trains are much more complex and very much in contrast with the solitary wave paradigm which for a long time was the accepted idealization of tsunami waves in both experimental and numerical work. The dissipation associated with these breaking bores can be very large, particularly over a wide and shallow continental shelf such as along the east coast of North America, particularly for the shorter waves associated with tsunamis generated by Submarine Mass Failures (SMFs). In this paper, we perform numerical simulations of tsunami coastal impact in the context of both idealized laboratory experiments and several tsunami case studies. We attempt to clarify t...
Recent observations of the coastal impact of large tsunamis (e.g., Indian Ocean 2004; Tohoku 2011... more Recent observations of the coastal impact of large tsunamis (e.g., Indian Ocean 2004; Tohoku 2011) and related numerical and theoretical works have made it increasingly clear that tsunami waves arrive nearshore as a series of long waves (so-called N-waves) with, often, the superposition of undular bores around each crest. Such wave trains are much more complex and very much in contrast with the solitary wave paradigm which for a long time was the accepted idealization of tsunami waves in both experimental and numerical work. The dissipation associated with these breaking bores can be very large, particularly over a wide and shallow continental shelf such as along the east coast of North America, particularly for the shorter waves associated with tsunamis generated by Submarine Mass Failures (SMFs). In this paper, we perform numerical simulations of tsunami coastal impact in the context of both idealized laboratory experiments and several tsunami case studies. We attempt to clarify t...
Pure and Applied Geophysics, 2020
We describe the development of a time-resolved model for tsunami-induced morphology change, based... more We describe the development of a time-resolved model for tsunami-induced morphology change, based on an existing Boussinesq model for weakly-dispersive free-surface waves coupled to a depth-integrated model for sediment concentration and an equation for bed level change. The model allows for spatial variability in bottom friction coefficients and accounts for the presence of non-erodible beds or finite depths of available sediment. The model is verified using one laboratory data set and against field observations of morphology change in the Crescent City, CA harbor during the 2011 Tohoku tsunami event. The model is then applied to a hypothetical example illustrating the impact of tsunami inundation on an undeveloped barrier island.
Scientific reports, Nov 4, 2016
Submarine slope failures are a likely cause for tsunami generation along the East Coast of the Un... more Submarine slope failures are a likely cause for tsunami generation along the East Coast of the United States. Among potential source areas for such tsunamis are submarine landslides and margin collapses of Bahamian platforms. Numerical models of past events, which have been identified using high-resolution multibeam bathymetric data, reveal possible tsunami impact on Bimini, the Florida Keys, and northern Cuba. Tsunamis caused by slope failures with terminal landslide velocity of 20 ms(-1) will either dissipate while traveling through the Straits of Florida, or generate a maximum wave of 1.5 m at the Florida coast. Modeling a worst-case scenario with a calculated terminal landslide velocity generates a wave of 4.5 m height. The modeled margin collapse in southwestern Great Bahama Bank potentially has a high impact on northern Cuba, with wave heights between 3.3 to 9.5 m depending on the collapse velocity. The short distance and travel time from the source areas to densely populated ...
Pure and Applied Geophysics, 2015
Your article is protected by copyright and all rights are held exclusively by Springer Basel. Thi... more Your article is protected by copyright and all rights are held exclusively by Springer Basel. This e-offprint is for personal use only and shall not be self-archived in electronic repositories. If you wish to self-archive your article, please use the accepted manuscript version for posting on your own website. You may further deposit the accepted manuscript version in any repository, provided it is only made publicly available 12 months after official publication or later and provided acknowledgement is given to the original source of publication and a link is inserted to the published article on Springer's website. The link must be accompanied by the following text: "The final publication is available at link.springer.com".
Natural Hazards, 2014
With support from the US National Tsunami Hazard Mitigation Program (NTHMP), the authors have bee... more With support from the US National Tsunami Hazard Mitigation Program (NTHMP), the authors have been developing tsunami inundation maps for the upper US East Coast (USEC), using high-resolution numerical modeling. These maps are envelopes of maximum elevations, velocity, or momentum flux, caused by the probable maximum tsunamis identified in the Atlantic oceanic basin, including from far-field coseismic or volcanic sources, and near-field Submarine mass failures (SMFs); the latter are the object of this work. Despite clear field evidence of past large-scale SMFs within our area of interest, such as the Currituck slide complex, their magnitude, pre-failed geometry, volume, and mode of rupture are poorly known. A screening analysis based on the Monte Carlo simulations (MCS) identified areas for possible tsunamigenic SMF sources along the USEC, indicating an increased level of tsunami hazard north of Virginia, potentially
Coastal Engineering Proceedings, 2012
Recent observations of the coastal impact of large tsunamis (e.g., Indian Ocean 2004; Tohoku 2011... more Recent observations of the coastal impact of large tsunamis (e.g., Indian Ocean 2004; Tohoku 2011) and related numerical and theoretical works have made it increasingly clear that tsunami waves arrive nearshore as a series of long waves (so-called N-waves) with, often, the superposition of undular bores around each crest. Such wave trains are much more complex and very much in contrast with the solitary wave paradigm which for a long time was the accepted idealization of tsunami waves in both experimental and numerical work. The dissipation associated with these breaking bores can be very large, particularly over a wide and shallow continental shelf such as along the east coast of North America, particularly for the shorter waves associated with tsunamis generated by Submarine Mass Failures (SMFs). In this paper, we perform numerical simulations of tsunami coastal impact in the context of both idealized laboratory experiments and several tsunami case studies. We attempt to clarify t...
Ocean Modelling, 2013
We derive fully nonlinear, weakly dispersive model equations for propagation of surface gravity w... more We derive fully nonlinear, weakly dispersive model equations for propagation of surface gravity waves in a shallow, homogeneous ocean of variable depth on the surface of a rotating sphere. A numerical model is developed for the weakly nonlinear version of the model based on a combined finite-volume and finitedifference method with a fourth-order MUSCL-TVD scheme in space and a third-order SSP Runge-Kutta scheme in time. In the context of tsunami generation and propagation over trans-oceanic distances, a scaling analysis reveals that the importance of frequency dispersion increases with a decrease of the source width, while the effect of the Coriolis force increases with an increase of the source width. A sensitivity analysis to dispersive and Coriolis effects is carried out using the numerical model in a series of numerical experiments in an idealized ocean using Gaussian and di-polar sources with different source sizes. A simulation of the Tohoku 2011 tsunami is used to illustrate the effects of dispersive and Coriolis effects at large distances from the source region.
Ocean Modelling, 2012
We present a high-order adaptive time-stepping TVD solver for the fully nonlinear Boussinesq mode... more We present a high-order adaptive time-stepping TVD solver for the fully nonlinear Boussinesq model of Chen (2006), extended to include moving reference level as in Kennedy et al. (2001). The equations are reorganized in order to facilitate high-order Runge-Kutta time-stepping and a TVD type scheme with a Riemann solver. Wave breaking is modeled by locally switching to the nonlinear shallow water equations when the Froude number exceeds a certain threshold. The moving shoreline boundary condition is implemented using the wetting-drying algorithm with the adjusted wave speed of the Riemann solver. The code is parallelized using the Message Passing Interface (MPI) with non-blocking communication. Model validations show good performance in modeling wave shoaling, breaking, wave runup and waveaveraged nearshore circulation.