Multiscale Quantiﬁcation of Tsunami Hazard Exposure in a Paciﬁc Small Island Developing State: The Case of Samoa

: This study presents a scenario-based approach for identifying and comparing tsunami exposure across different sociopolitical scales. In Samoa, a country with a high threat to local tsunamis, we apply scenarios for the 2009 South Paciﬁc tsunami inundation at different grid resolutions (50 and 10 m) to quantify building and road exposure at the national, district and village levels. We show that while the coarser 50 m model is adequate for use in the rapid identiﬁcation of exposure at the national and district levels, it can overestimate exposure by up to three times more at the village level. Overestimation typically occurs in areas characterized by ﬂat, low-lying, gentle-rising terrain. Overall, a 35% increase in buildings exposed to the 50 m model is observed compared with the 10 m scenario on southeast Upolu island. Similarly, a 31% increase in road exposure is observed for the 50 m scenario. These observations are discussed within the context of tsunami evacuation planning and logistics. Notwithstanding the variability in exposure, a precautionary approach leads us to conclude that while higher-resolution models are recommended where available data and/or ﬁnancial resources permit, the absence of such datasets should not preclude the use of coarser hazard datasets in risk assessments. Finer-resolution models provide more credence in detailed local-level exposure evaluation. While the results of this study are speciﬁc to the Samoan context, the results can be applied to the multiscale assessment of tsunami risk exposure in similar hazard contexts.


Introduction
Tsunami hazard and risk assessments provide crucial decision-support tools to inform disaster resilience in coastal locations, particularly in the context of preparedness, evacuation and land-use planning [1][2][3][4][5][6]. Indeed, with global tsunami-related casualties and damages estimated to be one hundred times greater between 1998 and 2017 compared to the preceding 20 year period, due largely to the 2004 Indian Ocean and 2011 Tohoku-oki tsunamis [7], the need for hazard risk information in tsunami-prone regions is essential (e.g., [8,9]).
In the Small Island Developing States (SIDS) of the Pacific, tsunami resilience planning is guided by the Framework for Resilient Development in the Pacific (FRDP) [10]. Aligned with the Sendai Framework for Disaster Risk Reduction 2015-2030 (SFDRR) [3], the FRDP provides recommendations on integrative approaches to address climate change and disaster risk management within Pacific Island country contexts. This includes enhancing adaptative capacity and resilience to tsunami hazards.
A significant proportion of SIDS populations and business hubs in the central Pacific are in coastal low-lying areas which are vulnerable to the impacts of tsunamis [10]. This was clearly demonstrated by the 2009 South Pacific tsunami (SPT) which impacted the Independent State of Samoa (herein Samoa) [11][12][13][14], causing widespread loss of life and infrastructure damage along the southern shores of Upolu and Savaii islands [15][16][17][18][19] ( Figure 1). While planning tools to support tsunami resilience in Samoa have been developed, such as identifying evacuation and hazard mitigation zones (e.g., [20,21]), these were purposefully designed to be conservative and derived using relatively coarse datasets. Consequently, they may not necessarily be representative of onshore hazard risk characteristics at local village levels.
GeoHazards 2021, 2, x FOR PEER REVIEW 2 of 18 Aligned with the Sendai Framework for Disaster Risk Reduction 2015-2030 (SFDRR) [3], the FRDP provides recommendations on integrative approaches to address climate change and disaster risk management within Pacific Island country contexts. This includes enhancing adaptative capacity and resilience to tsunami hazards. A significant proportion of SIDS populations and business hubs in the central Pacific are in coastal low-lying areas which are vulnerable to the impacts of tsunamis [10]. This was clearly demonstrated by the 2009 South Pacific tsunami (SPT) which impacted the Independent State of Samoa (herein Samoa) [11][12][13][14], causing widespread loss of life and infrastructure damage along the southern shores of Upolu and Savaii islands [15][16][17][18][19] (Figure 1). While planning tools to support tsunami resilience in Samoa have been developed, such as identifying evacuation and hazard mitigation zones (e.g., [20,21]), these were purposefully designed to be conservative and derived using relatively coarse datasets. Consequently, they may not necessarily be representative of onshore hazard risk characteristics at local village levels.  [20]. Red zone denotes areas at potential threat to small tsunami which pose only a marine threat. Orange zone encompasses all areas at potential threat to tsunami inundation from regional/distant earthquakes which are not locally felt. Yellow zone denotes all areas at potential threat to tsunami inundation following a long and strongly felt earthquake. Districts located along southeast Upolu experienced the greatest losses in terms of casualties, damage to property and the built environment from the 2009 tsunami. Districts cited in text: 1 = Palauli West; 2 = Palauli I le Falefa; 3 = Satupaitea; 4 = Palauli East; 5 = Faasaleleaga-I; 6 = Aiga I le Tai; 7 = Aana Alofi-III; 8 = Falelatai/Samatau; 9 = Lefaga/Faleseela; 10 = Siumu; 11 = Falealili; 12 = Lotofaga; 13 = Lepa; 14 = Aleipata Itupa I Luga; 15 = Aleipata Itupa I Lalo. Villages cited in text: a = Saleapaga; b = Lalomanu; c = Satitoa; d = Malaela. Star = Apia city.  [20]. Red zone denotes areas at potential threat to small tsunami which pose only a marine threat. Orange zone encompasses all areas at potential threat to tsunami inundation from regional/distant earthquakes which are not locally felt. Yellow zone denotes all areas at potential threat to tsunami inundation following a long and strongly felt earthquake. Districts located along southeast Upolu experienced the greatest losses in terms of casualties, damage to property and the built environment from the 2009 tsunami. Districts cited in text: 1 = Palauli West; 2 = Palauli I le Falefa; 3 = Satupaitea; 4 = Palauli East; 5 = Faasaleleaga-I; 6 = Aiga I le Tai; 7 = Aana Alofi-III; 8 = Falelatai/Samatau; 9 = Lefaga/Faleseela; 10 = Siumu; 11 = Falealili; 12 = Lotofaga; 13 = Lepa; 14 = Aleipata Itupa I Luga; 15 = Aleipata Itupa I Lalo. Villages cited in text: a = Saleapaga; b = Lalomanu; c = Satitoa; d = Malaela. Star = Apia city.
Similarly, detailed exposure data representing built environments is limited in Pacific SIDS. The distribution of built elements, particularly buildings and roads for example, provide useful indicators for population clusters (e.g., [19]). Building and road exposure to the physical characteristics of tsunamis (e.g., flow depth, extent, distance from coast, shielding) can support the identification of vulnerable clusters or 'risk hotspots' (e.g., [8]), whereby mitigation and risk reduction strategies could be targeted (e.g., identifying compromised structures and routes where evacuation might be impeded).
In general, there is a dearth in country specific tsunami hazard and exposure information to inform coastal risk management and preparedness planning in Pacific SIDS. However, tsunami inundation and runup characteristics along with exposure and vulnerability data are key requirements in the risk quantification process (e.g., [22]). This would enable the comparison of risk characteristics between communities and the identification of areas where prevention and/or adaptation interventions should be prioritized.
The aim of this study is to quantify and assess the influence which tsunami inundation model resolution has on multiscale exposure distributions of built elements at risk. We apply representations of the 2009 SPT hydrodynamic model developed in [14], as well as digital representations of buildings and roads to quantify risk exposure distributions across multiple sociopolitical levels. The key objectives are to: (1) apply a scenario-based approach to assess building and road exposure to model 2009 SPT inundation at 50 and 10 m grid resolution as developed by [14]; and (2) to demonstrate the influence of inundation model resolution on exposure distributions at national, district, and village levels, and how these can support preparedness planning. We discuss the suitability of model resolution in rapid tsunami risk exposure analysis as well as the implications of our findings within the context of tsunami preparedness and evacuation planning.

Study Location and Tsunami Hazard Context
Samoa is located between 13 • -15 • S and 171 • -173 • W and has a population of approximately 200,000 spread across 41 districts comprising 310 villages [23,24]. The capital, Apia, comprises approximately 18% of the population and is located in north Upolu island, itself containing close to 75% of the total population. Salelologa in the Faasaleleaga-I district is the central business area on Savaii island to the west.
Given that approximately 70% of Samoa's total population and human environments are within coastal zones, the threats to inundation-related hazards, which includes tsunami, are extreme [25,26]. This was emphasized by the 2009 SPT which originated from a Mw 8.1 earthquake sequence involving complex near-simultaneous normal and reverse faulting at the Northern Tonga Trench (NTT) [11,12,27]. The 2009 SPT represents a maximum likely local event and was the second most devastating disaster in this region over the preceding century in terms of loss of life and trauma, i.e., second only to the 1918 influenza pandemic disaster [28]. While Samoa is known to have experienced far-field tsunamis such as that from the 1960 Valdivia event [29], the overall threat to safety from far-field tsunamis is much lower compared with local events which can impact the shores of Samoa in less than 20 min (e.g., [14,30]).
Recent studies have revealed an apparent long-term local tsunami hazard in Samoa associated with the NTT source region [31]. Indeed, the most recent predecessor to the 2009 SPT was the 1917 Mw 8.3 Samoa-Tonga earthquake and tsunami, which is known to have inundated the southern coasts of Savaii and Upolu [28,32] (Figure 1). Detailed historical records on the built environment's damage and loss of life from the 1917 event are scant (e.g., [28,31]). Available sedimentary and anecdotal evidence suggests that inundation was more geographically widespread than the 2009 SPT, extending west from Falealupo on Savaii eastward to Lalomanu on Upolu (e.g., [31]). However, there does not appear to be evidence of impacts comparable to the magnitude of the devastation experienced due to the 2009 event, which was focused on southeast Upolu. Nevertheless, both the 2009 and 1917 events reinforce the high hazard potential posed by tsunamis sourced at the NTT, providing analogies to inform and assist with planning.
In 2011, tsunami evacuation zones were developed for Samoa [20], providing baselines to identify built features that are located in areas which could be inundated in a tsunami event and which might require a rapid evacuation response. These zones are characterized as red, orange and yellow and denote areas at potential threat to varying levels of hazard risk. For example, the yellow zone denotes the maximum area at potential threat to tsunami inundation after a strongly-felt (or 'local') earthquake, while the orange zone encompasses areas believed to be at risk following regional/distant earthquakes which are not locally felt. On the other hand, the red zone denotes areas at potential threat to small tsunami which might only pose a marine threat (e.g., strong hazardous currents).
Such tools can support the identification and setting of risk reduction objectives in local mitigation and resilience planning (e.g., identifying evacuation routes and safe centres, or critically exposed locations and assets) [2,[33][34][35]. These zones were developed using conservative attenuation-based ('bath-tub') models and are largely based on coarse elevation models (e.g., 90 m Shuttle Radar Topography Mission (SRTM) or greater), apart from Apia and Faleolo airport, which were the only locations with onshore Light Detection and Ranging (LiDAR) data at the time. The recent acquisition of 1 and 5 m gridded LiDAR topography and near-shore bathymetry for the whole of Samoa [36] now enables the development of more detailed hazard and risk assessments. For example, the work presented in [14] provides a characteristic set of 2009-type tsunami analogies (or hazard models), which can be applied to scenario-based exposure assessments (e.g., [37][38][39]), a key focus of this study.

Methods and Data
The components of the exposure model assessment approach used in this study are described in three sections. The tsunami inundation layers representing the 'hazard' component in the analysis are outlined in Section 3.1. The 'exposure' data layers capturing the distribution of buildings and roads in coastal areas are described in Section 3.2, with the exposure analysis framework described in Section 3.3.

Tsunami Hazard Layers
Numerical models for the 2009 SPT inundation in Upolu and Savaii [14] provide a characteristic set of maximum likely local event scenarios to quantify and compare the exposure characteristics at different sociopolitical levels. The models developed in [14] were simulated using BG-Flood (Block-adaptive on Graphics processing unit Flood model) [40]. BG-Flood was configured with a nested grid scheme to simulate tsunami initiation from the near-simultaneous earthquake sequence through to propagation and inundation. LiDAR topography and nearshore bathymetry with a spatial grid resolution of 5 m [36], provided the digital elevation model (DEM) used in simulating the inundation and runup.
Twenty-three scenarios of maximum runup and inundation were the output for the whole country at 50 m grid resolution, with thirty-four scenarios output at 10 m grid resolution for southeast Upolu. The individual scenarios represent the range of different proposed source configurations available in the published literature as synthesized by [14]. Each scenario differed in terms of slight variations to the source parameters such as rise time and rupture dimensions, with a uniform quadratic bottom roughness (cf = 0.0001) applied for the 50 m inundation grid. A spatially variable roughness formulation was applied to the higher-resolution 10 m inundation grid, whereby a high roughness length (z 0 = 0.01 m) was applied to all areas higher than 15 m below mean sea level, which included coral reefs as well as onshore built and vegetation features [14]. A low roughness length (z 0 = 0.0001 m) was used for all other areas lower than 15 m below mean sea level [14]. Since the 50 m model results were used to force the 10 m model (i.e., nested grid scheme), the tsunami waveform is quite similar in both grids. The main difference between the two model suites are due to roughness, as well as topographic representation which denotes the average elevation over each grid cell.
Model outputs of inundation extent and flow depth in digital raster format for both the 10 and 50 m resolution models showed good agreements between modelled-to-observed tsunami water level records (e.g., DART buoys and Apia Tide Gauge), as well as physical observations of runup and inundation (e.g., [14,16,17,41,42]). The maximum envelope for all 23 inundation scenarios at 50 m resolution provided the hazard layer for the whole country, with maximum flow depth providing the intensity parameter in the exposure analysis. For southeast Upolu, the maximum envelope and associated flow depths for all 34 inundation scenarios at 10 m resolution was used ( Figure 2). physical observations of runup and inundation (e.g., [14,16,17,41,42]). The maximum envelope for all 23 inundation scenarios at 50 m resolution provided the hazard layer for the whole country, with maximum flow depth providing the intensity parameter in the exposure analysis. For southeast Upolu, the maximum envelope and associated flow depths for all 34 inundation scenarios at 10 m resolution was used ( Figure 2).
Using the maximum envelope for both the 50 and 10 m resolution model suites was necessary for the purposes of this study to provide an overall, albeit conservative, representation of the 2009 SPT hazard given that no single scenario adequately produced consistent modelled-to-observed matches between the DART buoys, tide gauge, and runup) [14]. In addition, the tsunami evacuation zones developed in [20] were used as a benchmark for identifying and comparing built features located on land which could potentially be exposed to inundation associated with a range of potential tsunamigenic sources.  Using the maximum envelope for both the 50 and 10 m resolution model suites was necessary for the purposes of this study to provide an overall, albeit conservative, representation of the 2009 SPT hazard given that no single scenario adequately produced consistent modelled-to-observed matches between the DART buoys, tide gauge, and runup) [14]. In addition, the tsunami evacuation zones developed in [20] were used as a benchmark for identifying and comparing built features located on land which could potentially be exposed to inundation associated with a range of potential tsunamigenic sources.

Buildings, Roads and Sociopolitical Boundary Layers
The sum of building and road features used in the analysis are shown in Table 1. Buildings at risk to tsunami comprised of an amalgamated dataset based on: (1) the 2009 building outlines for the Apia area obtained from the Pacific Catastrophe and Risk Financing Initiative (PCRAFI) database; (2) Open Street Maps (OSM) building outlines for the rest of the country (current as of October 2020); and (3) manually digitized/updated buildings outlines based on 2020 Google satellite imagery. A composite of these datasets was required for the optimal coverage of the whole country, as no single dataset adequately captured the full distribution of buildings in Samoa. For example, the PCRAFI dataset covers the Apia urban area only, with the OSM data covering much of rural Samoa, but not Apia. Manual digitization and/or updating of coastal building outlines located within the maximum envelope of the 50 m resolution tsunami model was carried out using the QGIS 3.10.4 package using similar digitization techniques as presented in [19]. Physical attributes of building shape/outline and location were assigned. Importantly, the exposure analysis considers all building types without compartmentalization into specific typologies and/or use categories. This is due to the difficulties in assigning such attributes remotely without ground-truthing. A road centerline dataset incorporating the physical attributes of length and location was obtained from the OSM database (accessed in October 2020). These were visually checked for consistency with Google satellite imagery and then split into 10 m road lengths using the split-by-distance function in QGIS 3.10.4. Buildings and roads, respectively, were spatially joined to district and village boundary polygons obtained from the Samoa Ministry of Natural Resources and Environment (MNRE). This enabled individual features to be labelled with the districts and villages corresponding to their physical location.

Exposure Analysis
The natural hazards impact and loss modelling software, RiskScape [22] was configured to enumerate building and road exposure to tsunami flow depths for both the 10 and 50 m resolution model representations ( Figure 3). RiskScape has been applied in several tsunami exposure and impact analyses in Samoa and internationally (e.g., [6,19]). In this study, the RiskScape framework combines hazard intensity (e.g., flow depth), with building and road exposure information to determine the spatial distribution of tsunami hazard exposure. Flow depth increments of 0.5 m ranging from 0.0 to >3.0 m were defined. Buildings and roads intersecting inundated grid-cells were classed as exposed (i.e., grid cell > 0.0 m flow depth). Building and road exposure were enumerated at 0.5 m incremental flow depths at national, district and village levels. While higher resolution increments (e.g., 0.1 m) can particularly be useful for land-use or mitigation planning, the results are presented in 0.5 m increments which denote a minimum likely threshold, whereby evacuation potential begins to be severely compromised (e.g., [44,45]). Features/clusters exposed to >1.0 m flow depth are typically indicative of critical areas, whereby evacuation potential is mostly fully compromised (e.g., [44,45]). In addition, the sum of the buildings and road features exposed within the tsunami evacuation zones is presented for comparison for the 50 and 10 m scenarios for the 2009 SPT. presented in 0.5 m increments which denote a minimum likely threshold, whereby evacuation potential begins to be severely compromised (e.g., [44,45]). Features/clusters exposed to >1.0 m flow depth are typically indicative of critical areas, whereby evacuation potential is mostly fully compromised (e.g., [44,45]). In addition, the sum of the buildings and road features exposed within the tsunami evacuation zones is presented for comparison for the 50 and 10 m scenarios for the 2009 SPT.  [14] and evacuation zones [20] are combined with exposure layers (distribution of buildings and road features in districts and villages) to enumerate the sum of the exposed features. The RiskScape analysis framework is described by [22] and has been previously demonstrated in several inundation-related impacts and exposure analyses (e.g., [6,19,46]).

Results and Interpretations
Our results are presented in four sections: (1) national-or island-level overview of the total exposure of buildings and roads to the 2009 SPT scenario relative to the total number of buildings/roads within the tsunami evacuation zones; (2) district-level exposure based on the 50 m grid tsunami inundation; (3) village-level exposure based on the 10 m grid tsunami inundation; and (4) comparisons of district-and village-level exposure for 50 and 10 m inundation grids.

National Building and Road Exposure
A total of 12,040 buildings are located within the inundation extent of the 50 m grid tsunami model, of which approximately 74% are located on Upolu (which also includes Manono island in the Aiga o le Tai district) ( Table 2). In comparison, 25,483 buildings are located within the maximum extent of all tsunami evacuation zones; approximately 111% greater than the number of buildings exposed in the 2009 scenario. Savaii has a slightly greater proportion of buildings within the tsunami evacuation zones (51%) compared to Hazard layers depicting maximum extent and flow depths at 50 and 10 m grid resolution [14] and evacuation zones [20] are combined with exposure layers (distribution of buildings and road features in districts and villages) to enumerate the sum of the exposed features. The RiskScape analysis framework is described by [22] and has been previously demonstrated in several inundation-related impacts and exposure analyses (e.g., [6,19,46]).

Results and Interpretations
Our results are presented in four sections: (1) national-or island-level overview of the total exposure of buildings and roads to the 2009 SPT scenario relative to the total number of buildings/roads within the tsunami evacuation zones; (2) district-level exposure based on the 50 m grid tsunami inundation; (3) village-level exposure based on the 10 m grid tsunami inundation; and (4) comparisons of district-and village-level exposure for 50 and 10 m inundation grids.

National Building and Road Exposure
A total of 12,040 buildings are located within the inundation extent of the 50 m grid tsunami model, of which approximately 74% are located on Upolu (which also includes Manono island in the Aiga o le Tai district) ( Table 2). In comparison, 25,483 buildings are located within the maximum extent of all tsunami evacuation zones; approximately 111% greater than the number of buildings exposed in the 2009 scenario. Savaii has a slightly greater proportion of buildings within the tsunami evacuation zones (51%) compared to Upolu, and approximately 314% more buildings within evacuation zones compared to the total number exposed in the 2009 scenario.

District Building and Road Exposure
The distribution of buildings exposed to different flow depths within each district using the 50 m grid inundation model are shown in Figure 4a. The Aiga I le Tai and Falelatai/Samatau districts in southwest Upolu exhibit the greatest numbers of exposed buildings. These districts, along with Lefaga/Faleseela and Siumu in south central Upolu, and Falealili, Lotofaga, Lepa and Aleipata Itupa I Lalo in southeast Upolu, exhibit the highest proportionate exposure distributions relative to the total number of exposed buildings per district. That is, >80% of the total building exposure in each of these districts is exposed to flow depths of more than 0.5 m. On Savaii, districts located in the southeast have the greatest numbers of exposed buildings, with >80% of the total number of buildings exposed in Palauli East, Palauli West, Palauli i le Falefa, and Satupiatea subjected to flow depths >0.5 m.
Upolu, and approximately 314% more buildings within evacuation zones compared to the total number exposed in the 2009 scenario.

District Building and Road Exposure
The distribution of buildings exposed to different flow depths within each district using the 50 m grid inundation model are shown in Figure 4a. The Aiga I le Tai and Falelatai/Samatau districts in southwest Upolu exhibit the greatest numbers of exposed buildings. These districts, along with Lefaga/Faleseela and Siumu in south central Upolu, and Falealili, Lotofaga, Lepa and Aleipata Itupa I Lalo in southeast Upolu, exhibit the highest proportionate exposure distributions relative to the total number of exposed buildings per district. That is, >80% of the total building exposure in each of these districts is exposed to flow depths of more than 0.5 m. On Savaii, districts located in the southeast have the greatest numbers of exposed buildings, with >80% of the total number of buildings exposed in Palauli East, Palauli West, Palauli i le Falefa, and Satupiatea subjected to flow depths >0.5 m.
Districts on Upolu showing the longest lengths of exposed roads include Falealili in the southeast of the island, followed by Aana Alofi III and Aiga I le Tai in the southwest (Figure 4b). In several districts along Upolu's south coast, >80% of the exposed road lengths are at risk to flow depths >0.5 m, with >90% at risk in Lepa. On Savaii, districts located on the southeast of the island in general exhibit the longest lengths of exposed roads. Faasaleleaga-I has the longest length of exposed roads, with between 75 and 80% exposed to flow depths >0.5 m. (a)

Village Building and Road Exposure in Southeast Upolu
The variability in the total exposure distributions between the 10 and 50 m grid inundations models relative to the total number of features within evacuation zones on southeast Upolu are shown in Table 3. Compared with the total number of buildings exposed in the 10 m grid scenario, approximately 520 more are exposed in the 50 m scenario (35% increase) and 1380 (128% increase) located within the maximum evacuation zone limit. Similar trends are observed for the exposed roads, with a 32% increase in exposure observed in the 50 m grid scenario and 116% increase within the evacuation zones. Differences are observed between the 10 and 50 m grid scenarios for each village ( Figure 5), with notable variations in the proportionate distribution of features exposed to different flood depths. For example, some villages exhibit a greater number of exposed buildings in the 50 m scenario (e.g., Lalomanu), some exhibit less (e.g., Saleapaga), while in others the difference is negligible (e.g., Malaela) (Figure 5a). Similarly, villages located in the Aleipata Itupa I Lalo district exhibit greater numbers of buildings exposed to flow depth >1.0 m in the 50 m compared with the 10 m scenario. Overall, a total of 1138 buildings are exposed to flow depth >0.5 m in the 10 m grid scenario compared with 1427 exposed in the 50 m scenario-an approximate 25% difference. Comparable exposure differences between the 10 and 50 m grid scenarios are also observed for roads (Figure 5b), with an overall difference of approximately 26%. Districts on Upolu showing the longest lengths of exposed roads include Falealili in the southeast of the island, followed by Aana Alofi III and Aiga I le Tai in the southwest (Figure 4b). In several districts along Upolu's south coast, >80% of the exposed road lengths are at risk to flow depths >0.5 m, with >90% at risk in Lepa. On Savaii, districts located on the southeast of the island in general exhibit the longest lengths of exposed roads. Faasaleleaga-I has the longest length of exposed roads, with between 75 and 80% exposed to flow depths >0.5 m.

Village Building and Road Exposure in Southeast Upolu
The variability in the total exposure distributions between the 10 and 50 m grid inundations models relative to the total number of features within evacuation zones on southeast Upolu are shown in Table 3. Compared with the total number of buildings exposed in the 10 m grid scenario, approximately 520 more are exposed in the 50 m scenario (35% increase) and 1380 (128% increase) located within the maximum evacuation zone limit. Similar trends are observed for the exposed roads, with a 32% increase in exposure observed in the 50 m grid scenario and 116% increase within the evacuation zones. Differences are observed between the 10 and 50 m grid scenarios for each village ( Figure 5), with notable variations in the proportionate distribution of features exposed to different flood depths. For example, some villages exhibit a greater number of exposed buildings in the 50 m scenario (e.g., Lalomanu), some exhibit less (e.g., Saleapaga), while in others the difference is negligible (e.g., Malaela) (Figure 5a). Similarly, villages located in the Aleipata Itupa I Lalo district exhibit greater numbers of buildings exposed to flow depth >1.0 m in the 50 m compared with the 10 m scenario. Overall, a total of 1138 buildings are exposed to flow depth >0.5 m in the 10 m grid scenario compared with 1427 exposed in the 50 m scenario-an approximate 25% difference. Comparable exposure differences between the 10 and 50 m grid scenarios are also observed for roads (Figure 5b), with an overall difference of approximately 26%.

Discussion
The analysis performed in this work has demonstrated the practicality of using a combination of available information tools to perform a rapid multiscale tsunami exposure analysis. In Samoa, integrating modelled tsunami inundation from the BG-Flood tool with open building and road datasets with the RiskScape exposure model demonstrates the ability of risk tools to quantify exposure at multiple sociopolitical scales. However, observed variability in exposure distributions due to hazard model resolution poses several questions regarding the use and suitability of such information in planning contexts.

Exposure Vairability
The variability in building and road exposure distributions between the 10 and 50 m grid scenarios on southeast Upolu demonstrates the effect of model resolution on maximum inundation extent and subsequent risk exposure. Furthermore, they highlight the influence of model resolution on inundation flow depth and consequent distribution of proportionate exposure for each village. For this case, the 50 m grid scenario with inundation extent further inland and coarser representation of flow depth exposes up to 35% and 31% more buildings and roads, respectively, compared with the 10 m grid scenario ( Figure 6). Villages exhibiting minimal differences between the 50 and 10 m grid scenarios (e.g., Lepa and Saleapaga) are typically characterized by narrow coastal areas backed by steep cliffs. In contrast, villages with notable differences (e.g., Satitoa) are generally located in wider/flat low-lying coastal regions. This likely indicates the influence of terrain/topographical gradient on the model grid and subsequent representation of inundation extent, implying that the 50 m grid scenario is more likely to overestimate the exposure in low-lying, flat coastal terrain compared to the 10 m scenario (e.g., in parts of southern and southwestern Upolu).  Despite these variabilities, the general consistency in the relative exposure of villages in each district for both scenarios suggests that the 50 m representation does provide a suitable, while conservative, level of detail to facilitate a comparison of exposure at the island/district scale. That is, it does provide a rapid identification of representative high-risk clusters/areas where targeted mitigation interventions might be required. This can also facilitate assessments of the differences in exposure characteristics which might influence differential response (e.g., such as during an evacuation) or recovery among districts (e.g., longer-term land development/use plans).
However, the differences do pose challenges for tsunami preparedness and response planning at the village (or community) level. For example, in Satitoa, an additional 50 buildings are exposed in the 50 m grid scenario compared with the 10 m scenario. Such situations could translate into an additional 100 or more people being threatened and requiring evacuation, which will have direct implications on logistics planning and resourcing in evacuation centres. Similar challenges are apparent when comparing the total number of buildings and roads exposed in the 2009 local scenario with the total features located within tsunami evacuation zones (Figures 6 and 7). For example, in a warning situation to a local event, such as the 2009 tsunami, people will be required to evacuate outside of the maximum evacuation zone extent (i.e., yellow or orange zone in areas where the yellow zone is not present).

Applications in Tsunami Risk Management
While the variability in exposure distributions discussed above provide useful tools in understanding some of the uncertainty and subsequent challenges this might have on planning, these should be treated with caution and within the context of acceptable risk. For example, the overall difference of approximately 25% between the total number of buildings exposed in the 10 m grid compared with the 50 m grid scenario on southeast Nevertheless, with the increasing availability of finer resolution global DEMs (e.g., 30 m SRTM) as well as decreasing trends in the costs to acquire LiDAR data, the multiscale tsunami risk screening approach presented in this study can be adapted for other Pacific SIDS which are currently lacking whole-of-country LiDAR topography. Furthermore, the presented framework provides a potentially rapid method for assessing tsunami exposure in other countries to assist with preparedness and mitigation planning.

Applications in Tsunami Risk Management
While the variability in exposure distributions discussed above provide useful tools in understanding some of the uncertainty and subsequent challenges this might have on planning, these should be treated with caution and within the context of acceptable risk. For example, the overall difference of approximately 25% between the total number of buildings exposed in the 10 m grid compared with the 50 m grid scenario on southeast Upolu could be considered as an acceptable level of uncertainty. That is, the coarser 50 m scenario would provide a more conservative estimate suitable for application in preparedness and scenario-based evacuation planning.
Indeed, the general consistency in the total exposure trends between districts and villages for both scenarios corroborates this. High risk areas which are represented here by the total number of building and road features exposed in each district/village in the 50 m scenario are consistent with high risk areas identified in the 10 m scenario. However, the variability in proportionate exposure to different flow depths in each village suggests that the finer resolution is desired for more localized evaluations (e.g., at the village/community scales). For example, the total number of buildings and length of road in each village exposed to flow depth >0.5 m which could significantly impede/compromise evacuation ability are more accurately represented in the 10 m scenario.
While the choice of using a specific model resolution is influenced by other factors, such as computational costs relative to available resources, it should be informed and treated within the context of the desired application. The results in this study suggest that the 50 m grid scenario provides an acceptable level of detail required to identify and compare the characteristics of exposure risk hotspots around the country, including vulnerable locations where targeted interventions might be required to reduce vulnerability. On the other hand, the variability observed between the two scenarios suggests that the 50 m grid would not provide the level of detail required in localized, property-or feature-scale, assessments (e.g., insurance planning). In addition, limitations in the known attributes of the feature data (buildings and roads) in this study highlight the need for continual updating and ground validation of asset datasets in relation to planning requirements. For example, while the exposure of an individual building can be determined based on its outline/shape feature, its usage (e.g., commercial, residential) and criticality (e.g., health centres, schools) in tsunami preparedness is not known. Future work should consider the attribution of such parameters to infrastructure datasets including detailed ground-truthing to improve the level of detail and application of future exposure and impacts analysis.
Over the longer-term, the application of these types of information analysis tools in tsunami risk reduction and planning should be complemented with updated asset stocks to monitor and factor changing risk patterns in tsunami risk reduction plans (e.g., updating of building stock/features in line with 5-or 10-yearly census updates) (e.g., [19]). Longer-term planning should also incorporate future SLR and land-subsidence projections (e.g., [47][48][49][50][51]) to account for the effects these pose to the range of potential tsunami inundation and associated risk distributions (e.g., [52]).

Summary and Conclusions
This study examined the influence of grid model resolution on tsunami inundation arising from an event such as the 2009 SPT and subsequent variability in exposure distributions in Samoa, as well as the implications these pose to tsunami hazard and risk reduction planning. Overall, a 35% and 31% difference in the total exposure of buildings and roads, respectively, is observed between the 50 and 10 m grid inundation representations for southeast Upolu. The coarser model is more likely to overestimate exposure by up to three times or more in wide low-lying flat terrain at the village scale. In contrast, the differences in exposure between villages along narrow coastal stretches backed by steep cliffs is min-imal (i.e., up to approximately 200 m or less inland of the coastline). Nevertheless, the overall consistency in the relative distribution of exposure between districts in both the 10 and 50 m scenarios suggests that the observed variability is minor. That is, the 50 m scenario does provide the necessary level of detail for comparing the risk characteristics between districts and identifying regions where targeted risk mitigation/reduction interventions might be prioritized in national/district-level planning. The observed variability in exposure at the village scale, particularly regarding differential exposure of features to different flow depths, indicates the need for using fine resolution models in detailed evacuation and response plans at the community level.
These observations lead us to infer that while higher-resolution models are recommended where data/resources permit (e.g., LiDAR, financial resources), the absence of such datasets should not preclude the application of coarser hazard resolutions in risk assessments. This is particularly relevant in regions where the costs of producing high resolution models might exceed the available resources to produce them (i.e., within SIDS contexts). We suggest the suitability of model resolution should be treated within the scope and scale of the desired decision-support application. While the observations presented in this study are specific to the Samoan context, particularly in relation to the distribution of the exposed human environment, the results could be applied to similar risk exposure assessments for other inundation hazards, such as coastal or fluvial flooding in analogous SIDS contexts.