Oceanography The Official Magazine of
The Oceanography Society
Volume 36 Supplement 1

View Issue TOC
Volume 36, No. 1
Pages 74 - 75

OpenAccess

Assessment of Tsunami Hazard Along British Columbia Coastlines from Coseismic Sources

By Soroush Kouhi , Reza Amouzgar, Mark Rankin, Cassandra Bosma, and Adrienne Shumlich 
Jump to
Full text Citation References Copyright & Usage
Full Text

Recent major tsunamis have demonstrated their potential for widespread destruction and associated danger to human populations. Tsunamis generated from local and distant sources in the Pacific Ocean have impacted the western coast of Canada over time, and destructive tsunamis are expected to affect west coast communities in the future. Ocean Networks Canada (ONC) has undertaken several collaborative tsunami hazard assessment projects along the west coast of Canada (Figure 1) with the goals of (1) identifying and assessing tsunami hazards and risks using tsunami modeling, mapping, community engagement, and Indigenous knowledge; and (2) raising public awareness by producing a documentary film, offering educational programs to K–12 students, and developing tsunami signage and online story maps. This article focuses on tsunami hazard identification through modeling and drawing upon Indigenous knowledge.

 

FIGURE 1. Several study areas along the coastlines of British Columbia, Canada, where tsunami risk assessment was performed by Ocean Networks Canada. Seismic source data are shown as seafloor vertical displacement (m) for the Alaska-Aleutian and Cascadia subduction zones. Blue indicates the seafloor/topography subsidence and red indicates seafloor/topography uplift. > High res figure

 

Modeling Methodology

Modeling the generation and propagation of tsunami waves requires three essential elements: (1) ocean floor (bathymetry) and land (topography) elevation data, (2) seismic source data from local and distant earthquakes, and (3) a suitable numerical model.

Multiple bathymetric and topographic data sets were integrated to develop the high-resolution digital elevation models (DEMs) required for tsunami modeling, including lidar and single beam and multibeam bathymetric data. For the local tsunami source, the Cascadia subduction zone, Natural Resources Canada (NRCan) provided a recent rupture model corresponding to a magnitude 9 earthquake (Gao et al., 2018). This source resembles great megathrust earthquakes that occur roughly once every 300–500 years in this region, with the latest in 1700. While there is no written record of the impact of this earthquake along the eastern Pacific, it has been recorded in oral history as well as in the coastal and offshore stratigraphy from northern California to Vancouver Island. The distant tsunami source, the Alaska-Aleutian subduction zone, was the origin site for numerous significant earthquakes of magnitude 8 and more in the twentieth century. The largest of recent earthquakes, the 1964 Alaska earthquake, produced the largest instrumentally recorded tsunami waves to date on the British Columbia coast; this event provides a realistic proxy for similar large events that may be generated by the Alaska-Aleutian subduction zone. The numerical simulation of the 1964 tsunami in Figure 1 was based on recent earthquake information acquired from the University of Alaska Fairbanks (Suleimani and Freymueller, 2020).

We modeled the propagation and inundation of tsunamis induced by seismic events using FUNWAVE-TVD (Shi et al., 2012). The first step was to set up a series of nested grids to enable the model to resolve the complex physics as a tsunami propagates from the deep ocean into shoaling coastal regions. The nested grids increase resolution from the outer grid (2 arc-minutes) to the inner grid (~10 m). The wave elevations and horizontal water velocities are transferred from the coarser resolution to finer resolution grids at the model boundaries. Several parameters were also adjusted in the model based on the geographic location of each study area.

The next step was to establish tsunami model scenarios based on current-day sea level and projected future sea level rise (SLR). Recently, James et al. (2021) studied several averaged and enhanced SLR scenarios for the water bodies surrounding Canada, based on different pathways of greenhouse gas emissions through the twenty-first century. The worst-case scenario for 2100 predicts a maximum of up to 1.5 m of SLR for the Enhance Scenario at 2100 along British Columbia coasts. For our tsunami simulations, we lowered the underlying DEM by the local projected SLR value as a way to account for SLR.

 

Community and Indigenous Engagement

Indigenous communities living where tsunami risk is high have fostered disaster resilience for millennia. Highlighting the knowledge and experience of the Indigenous community is critical for understanding current and future risks and implementing mitigation plans. Using a framework of “two-eyed seeing,” where different ways of knowing are used simultaneously for the benefit of all, helps to improve the science by learning from community stories. For the Northwest Vancouver Island tsunami risk assessment study, through both research and interviews conducted with Elders and community members, multiple stories, teachings, and oral histories of past tsunamis, including but not limited to the 1700 and 1964 tsunamis, were collected and used to authenticate modeling results (Figure 2).

 

FIGURE 2. This image of a silkscreen print created in 1977 by artist A-nii-sa-put (Tim Paul) from Mowachaht-​Muchalaht First Nation depicts the earthquake foot legend. Along with other stories, this one highlights the knowledge and teachings that exist regarding earthquakes and tsunamis. > High res figure

 

Results

The modeling results demonstrate that tsunami waves generated by both local and distant earthquakes would considerably impact British Columbia coastal communities. The impact of the locally generated tsunami is more significant in Boundary Bay, Port Alberni, and Northwest Vancouver Island, while a tsunami originating in Alaska would have greater impacts in Prince Rupert and Haida Gwaii. Tsunami arrival times to the outer coast of Northwest Vancouver Island and Boundary Bay were estimated at about 20 minutes and 3 hours, respectively, after an earthquake in the Cascadia subduction zone. The Alaska tsunami, however, reaches Haida Gwaii and Prince Rupert about 2.5 and 3.5 hours, respectively, after the earthquake occurs.

Results also suggest that SLR does not significantly influence the tsunami wave amplitude offshore. However, propagation of these waves over land will lead to a greater extent of flooding as a result of SLR. Several overwater hazard graphics (Figure 3) and inundation maps (Figure 4) were created to represent the tsunami hazard information in each study area, including the tsunami wave amplitudes, current velocities, tsunami arrival times, and extent of flooding.

 

FIGURE 3. Maximum wave amplitudes are shown for a Cascadia-generated tsunami in Esperanza Inlet and Nootka Sound, as determined by the Northwest Vancouver Island Tsunami Risk Assessment Study. > High res figure

 

FIGURE 4. Inundation extent map for a Cascadia Tsunami in Tahsis, BC, Canada, were produced by the Northwest Hydraulic Consultant Ltd. in collaboration with Ocean Networks Canada for the Northwest Vancouver Island Tsunami Risk Assessment Project. > High res figure

 

Acknowledgments

The authors would like to thank Northwest Hydraulic Consultants Ltd., Ocean Networks Canada partner at tsunami risk assessment projects for Prince Rupert, Northwest Vancouver Island, and Haida Gwaii. The authors would also like to thank the funders and involved communities, including the North Coast Regional District, the city of Prince Rupert, the village of Daajing Giids, the village of Masset, the village of Port Clements, Strathcona Regional District, Ka:’yu:’k’t’h’/Che:k:tles7et’h’ First Nations, Nuchatlaht First Nation, Ehattesaht/Chinehkint First Nations, Quatsino First Nations, Mowachaht/Muchalaht First Nations, and Semiahmoo First Nation.
Citation

Kouhi, S., R. Amouzgar, M. Rankin, C. Bosma, and A. Shumlich. 2023. Assessment of tsunami hazard along British Columbia coastlines from coseismic sources. In Frontiers in Ocean Observing: Emerging Technologies for Understanding and Managing a Changing Ocean. E.S. Kappel, V. Cullen, M.J. Costello, L. Galgani, C. Gordó-Vilaseca, A. Govindarajan, S. Kouhi, C. Lavin, L. McCartin, J.D. Müller, B. Pirenne, T. Tanhua, Q. Zhao, and S. Zhao, eds, Oceanography 36(Supplement 1):74–75, https://doi.org/10.5670/oceanog.2023.s1.24.

References

Gao, D., K. Wang, T.L. Insua, M. Sypus, M. Riedel, and T. Sun. 2018. Defining megathrust tsunami source scenarios for northernmost Cascadia. Natural Hazards 94:445–469, https://doi.org/10.1007/s11069-018-3397-6.

James, T.S., C. Robin, J.A. Henton, and M. Craymer. 2021. Relative sea-level projections for Canada based on the IPCC Fifth Assessment Report and the NAD83v70VG national crustal velocity model. Geological Survey of Canada, Open File 8764, 1 .zip file, https://doi.org/​10.4095/327878.

Shi, F., J.T. Kirby, J.C. Harris, J.D. Geiman, and S.T. Grilli. 2012. A high-​order adaptive time-stepping TVD solver for Boussinesq modelling of breaking waves and coastal inundation. Ocean Modelling 43–44:36–51, https://doi.org/10.1016/j.ocemod.2011.12.004.

Suleimani, E., and J.T. Freymueller. 2020. Near-field modeling of the 1964 Alaska tsunami: The role of splay faults and horizontal displacements. Journal of Geophysical Research: Solid Earth 125(7), e2020JB019620, https://doi.org/10.1029/2020JB019620.

Copyright & Usage

This is an open access article made available under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/), which permits use, sharing, adaptation, distribution, and reproduction in any medium or format as long as users cite the materials appropriately, provide a link to the Creative Commons license, and indicate the changes that were made to the original content. Images, animations, videos, or other third-party material used in articles are included in the Creative Commons license unless indicated otherwise in a credit line to the material. If the material is not included in the article’s Creative Commons license, users will need to obtain permission directly from the license holder to reproduce the material.