Variations in population vulnerability to tectonic and landslide-related tsunami hazards in Alaska

Natural Hazards - Effective tsunami risk reduction requires an understanding of how at-risk populations are specifically vulnerable to tsunami threats. Vulnerability assessments primarily have been...     

Correction to: Beat-the-wave evacuation mapping for tsunami hazards in Seaside,Oregon, USA

Natural Hazards - Due to a procedural error in construction of Figs.&;nbsp;8 and 9, listed minimum speeds to beat the tsunami wave in areas of Seaside seaward of Neawanna Creek are too high. The...     

Pacific Basin tsunami hazards associated with mass flows in the Aleutian arc of Alaska

We analyze mass-flow tsunami generation for selected areas within the Aleutian arc of Alaska using results from numerical simulation of hypothetical but plausible mass-flow sources such as submarine landslides and volcanic debris avalanches. The Aleutian arc consists of a chain of volcanic mountains, volcanic islands, and submarine canyons, surrounded by a low-relief continental shelf above about 1000–2000 m water depth. Parts of the arc are fragmented into a series of fault-bounded blocks, tens to hundreds of kilometers in length, and separated from one another by distinctive fault-controlled canyons that are roughly normal to the arc axis. The canyons are natural regions for the accumulation and conveyance of sediment derived from glacial and volcanic processes. The volcanic islands in the region include a number of historically active volcanoes and some possess geological evidence for large-scale sector collapse into the sea. Large scale mass-flow deposits have not been mapped on the seafloor south of the Aleutian Islands, in part because most of the area has never been examined at the resolution required to identify such features, and in part because of the complex nature of erosional and depositional processes. Extensive submarine landslide deposits and debris flows are known on the north side of the arc and are common in similar settings elsewhere and thus they likely exist on the trench slope south of the Aleutian Islands. Because the Aleutian arc is surrounded by deep, open ocean, mass flows of unconsolidated debris that originate either as submarine landslides or as volcanic debris avalanches entering the sea may be potential tsunami sources.To test this hypothesis we present a series of numerical simulations of submarine mass-flow initiated tsunamis from eight different source areas. We consider four submarine mass flows originating in submarine canyons and four flows that evolve from submarine landslides on the trench slope. The flows have lengths that range from 40 to 80 km, maximum thicknesses of 400–800 m, and maximum widths of 10–40 km. We also evaluate tsunami generation by volcanic debris avalanches associated with flank collapse, at four locations (Makushin, Cleveland, Seguam and Yunaska SW volcanoes), which represent large to moderate sized events in this region. We calculate tsunami sources using the numerical model TOPICS and simulate wave propagation across the Pacific using a spherical Boussinesq model, which is a modified version of the public domain code FUNWAVE. Our numerical simulations indicate that geologically plausible mass flows originating in the North Pacific near the Aleutian Islands can indeed generate large local tsunamis as well as large transoceanic tsunamis. These waves may be several meters in elevation at distal locations, such as Japan, Hawaii, and along the North and South American coastlines where they would constitute significant hazards.     

Microfossil analysis of sediments representing the 1964 earthquake, exposed at Girdwood Flats, Alaska, USA

Diatom, pollen, foraminifera and thecamoebian assemblages from an outcrop of peat and silt at Girdwood Flats, in the upper Turnagain Arm of the Cook Inlet, Alaska, record four phases of relative land and sea-level changes. The first phase is the development of freshwater swamp above high marsh sediments during relative land uplift, caused by strain accumulation along the locked portion of the Alaska-Aleutian subduction zone. In second phase, the top 2 cm of the peat, all microfossil groups record pre-seismic relative sea-level rise (relative land subsidence). The third phase is rapid land subsidence, 1.7 m, during the earthquake of March 1964 that initiated intertidal silt accumulation above the peat. The final phase is the colonisation of mudflat by salt marsh communities during post-seismic land uplift. The microfossil data compare favourably with sequences from Washington, Oregon and British Columbia that record late Holocene submergence events caused by earthquakes. The comparable changes in microfossil assemblages record the different phases of relative land and sea-level changes and the magnitude of land subsidence caused by each earthquake (expressed relative to the tidal range at the site). These results raise the question whether preseismic sea-level rise represents any kind of warning of large earthquakes.     

Kinematics, mechanics, and potential earthquake hazards for faults in Pottawatomie County, Kansas, USA

Many stable continental regions have subregions with poorly defined earthquake hazards. Analysis of minor structures (folds and faults) in these subregions can improve our understanding of the tectonics and earthquake hazards. Detailed structural mapping in Pottawatomie County has revealed a suite consisting of two uplifted blocks aligned along a northeast trend and surrounded by faults. The first uplift is located southwest of the second. The northwest and southeast sides of these uplifts are bounded by northeast-trending right-lateral faults. To the east, both uplifts are bounded by north-trending reverse faults, and the first uplift is bounded by a north-trending high-angle fault to the west. The structural suite occurs above a basement fault that is part of a series of north–northeast-trending faults that delineate the Humboldt Fault Zone of eastern Kansas, an integral part of the Midcontinent Rift System. The favored kinematic model is a contractional stepover (push-up) between echelon strike-slip faults. Mechanical modeling using the boundary element method supports the interpretation of the uplifts as contractional stepovers and indicates that an approximately east–northeast maximum compressive stress trajectory is responsible for the formation of the structural suite. This stress trajectory suggests potential activity during the Laramide Orogeny, which agrees with the age of kimberlite emplacement in adjacent Riley County. The current stress field in Kansas has a N85°W maximum compressive stress trajectory that could potentially produce earthquakes along the basement faults. Several epicenters of seismic events (<M2.0) are located within 10 km of the structural suite. One epicenter is coincident with the northwest boundary of the uplift. This structural suite, a contractional stepover between echelon northeast-trending right-lateral faults, is similar to that mapped in the New Madrid Seismic Zone, and both areas currently feature roughly east–west maximum compressive stress trajectory. Based on these similarities, the faults in Pottawatomie County have the potential for seismicity. The results demonstrate that mechanical analysis of minor structural features can improve our knowledge of local earthquake hazards.     

Regional differences in cognitive dissonance in evacuation behavior at the time of the 2011 Japan earthquake and tsunami

This paper constructs an evacuation decision-making model that takes cognitive dissonance into consideration. The purpose of this construction is to clarify the psychological mechanism for the evacuation behavior of residents during an emergency, based on Akerlof and Dickens (Akerlof and Dickens Am Econ Rev 72:307–319, 1982) "The economic consequences of cognitive dissonance". Specifically, we empirically explore people’s psychological mechanism (e.g., cognitive dissonance) for evacuation behavior when a tsunami disaster occurs. As a result, we show that the level of anxiety depends on the area where residents live and that the average anxiety of residents is mostly correlated with the level of damage of past disasters, and that it is affected also by the ages of residents. Since the level of anxiety largely affects an individual’s evacuation behavior, this result can indicate for what kinds of people intervention and assistance are required based on the level of anxiety. A high level of anxiety basically promotes evacuation. Since our results show that anxiety is increased by the experience of tsunamis, education having people virtually experience tsunamis may increase evacuation rates efficiently.

     

Asperity distribution of the 1964 Great Alaska earthquake and its relation to subsequent seismicity in the region

The 1964 Alaska earthquake was the second largest seismic events in the 20th century. The aim of this work is the use of surface deformation data to determine asperity and slip distributions on the fault plane of the Alaska earthquake: these distributions are calculated by a Monte Carlo method. To this aim, we decompose the fault plane in a large number of small square asperity units with a side of 25 km; this allows us to obtain plane surfaces with an irregular shape. In the first stage, each asperity unit is allowed to slip a constant amount or not to slip at all, providing the geometry of the dislocation surface that best reproduces the observed displacements. To this purpose, a large number of slip distributions have been tried by the use of the Monte Carlo method. The slip amplitude is the same for all the asperities and is equal to the average fault slip inferred from the seismic moment. In the second stage, we evaluate the slip distribution in the dislocation area determined by the Monte Carlo inversion: in this case, we allow unit cells to undergo different values of slip in order to refine the initial dislocation model. The results confirm the previous finding that the slip distribution of the great Alaska earthquake was essentially made of two dislocation areas with a higher slip, the Prince William Sound and the Kodiak asperities. Analysis of the post-1964 seismicity in the rupture region shows a strong correlation between the larger earthquakes (M_w≥6) and the distribution of locked asperities following the 1964 event, which can be considered as an independent test of the validity of the model. We do not find slip values higher than 25 m for any of the patches, and we determine two separate high-slip zones: one correspondent to the Prince William Sound asperity, and one (18 m slip) to the Kodiak asperity. The slip distribution connected with the 1964 shock appears to be consistent with the following seismicity in the region.     

Hazard assessment of the Tidal Inlet landslide and potential subsequent tsunami, Glacier Bay National Park, Alaska

An unstable rock slump, estimated at 5 to 10 × 10⁶ m³, lies perched above the northern shore of Tidal Inlet in Glacier Bay National Park, Alaska. This landslide mass has the potential to rapidly move into Tidal Inlet and generate large, long-period-impulse tsunami waves. Field and photographic examination revealed that the landslide moved between 1892 and 1919 after the retreat of the Little Ice Age glaciers from Tidal Inlet in 1890. Global positioning system measurements over a 2-year period show that the perched mass is presently moving at 3–4 cm annually indicating the landslide remains unstable. Numerical simulations of landslide-generated waves suggest that in the western arm of Glacier Bay, wave amplitudes would be greatest near the mouth of Tidal Inlet and slightly decrease with water depth according to Green’s law. As a function of time, wave amplitude would be greatest within approximately 40 min of the landslide entering water, with significant wave activity continuing for potentially several hours.     

Environmental hazards from population pressure in the Jinta Oasis, arid northwestern China

Oasis is not only the most concentrated area of human activities in arid region but also the largest area where artificial disturbances occur at the regional scale. The Jinta Oasis is a very typical artificial oasis in arid region of China and is set as one of the national land resource developing and representative areas in China. With the continuous increase in population and livestock number in Jinta Oasis, such severe problems of environmental degradation as serious land pollution and desertification, water environmental degradation, and vegetation degeneracy occur within the whole oasis.     

Modeling a tsunami from the Nicoya,Costa Rica,seismic gap and its potential impact in Puntarenas

Although subduction zones around the world are known to be the source of earthquakes and/or tsunamis, not all segments of these plate boundaries generate destructive earthquakes and catastrophic tsunamis. Costa Rica, in Central America, has subduction zones on both the Pacific and the Caribbean coasts and, even though large earthquakes (Mw = 7.4–7.8) occur in these convergent margins, they do not produce destructive tsunamis. The reason for this is that the seismogenic zones of the segments of the subduction zones that produce large earthquakes in Costa Rica are located beneath land (Nicoya peninsula, Osa peninsula and south of Limón) and not off shore as in most subduction zones around the world. To illustrate this particularity of Costa Rican subduction zones, we show in this work the case for the largest rupture area in Costa Rica (under the Nicoya peninsula), capable of producing Mw ～ 7.8 earthquakes, but the tsunamis it triggers are small and present little potential for damage even to the largest port city in Costa Rica.The Nicoya seismic gap, in NW Costa Rica, has passed its ～50-year interseismic period and therefore a large earthquake will have to occur there in the near future. The last large earthquake, in 1950 generated a tsunami which slightly affected the southwest coast of the Nicoya Peninsula. We present here a simulation to study the possible consequences that a tsunami generated by the next Nicoya earthquake could have for the city of Puntarenas. Puntarenas has a population of approximately eleven thousand people and is located on a 7.5 km long sand bar with a maximum height of 2 m above the mean sea level. This condition makes Puntarenas vulnerable to tsunamis.     

Problems in the application of the SSHAC probability method for assessing earthquake hazards at Swiss nuclear power plants

From 2000 to 2004 a large scale probabilistic seismic hazard analysis (PEGASOS) was created and performed as a research project, sponsored by the Swiss NPP utilities to improve the assessment methodology for seismic risk in support of the plant-specific seismic PRAs. The project followed the methodology of the SSHAC [Senior Seismic Hazard Analysis Committee (SSHAC), 1997. Recommendations for Probabilistic Seismic Hazard Analysis: Guidance on Uncertainty and Use of Experts. NU-REG/CR-6372] procedures at its most elaborate way—level 4. Before practical implementation was to occur, a detailed review was performed including validation tests and analysis of uncertainty propagation. This paper presents the main results of the review. The review revealed that current PSHA (Probabilistic Seismic Hazard Analysis) methodology as based on logic trees, in conjunction with the SSHAC procedures, potentially leads to a significant overestimation of the seismic hazard in areas with low seismic activity due to the inherent possibilities of unconstrained accumulation of uncertainties. The preliminary results of the project did not pass any of our logical geological–scientific benchmark tests used in our attempts to perform a validation process of the PEGASOS analysis method. Some of the problems encountered are of generic nature and shall be studied carefully before making the decision of whether or not the Swiss nuclear power industry will adopt the recommended use of SSHAC-procedures as a basis for the evaluation of the seismic hazard for individual nuclear power plant seismic PRA without correction.     

Great earthquakes,seismicity gaps and potential for earthquake disaster along the Himalaya plate boundary

Analysis of the space-time patterns of seismicity in the Himalaya plate boundary has established the existence of three seismic gaps:

• 1.(1) The “Kashmir gap” lying west of the 1905 Kangra earthquake;
• 2.(2) the “Central gap”, situated between the 1905 Kangra and the 1934 Bihar earthquakes;
• 3.(3) the “Assam gap” between the 1897 and 1950 Assam earthquakes.

This study has shown that the above great earthquakes were preceded as well as followed by long periods (⩾ 19 years) of decreased levels of seismic activity in the epicentral regions. Remarkable decrease in the seismicity following the year 1970 has been observed in the western half of the Central gap as well as in the Assam gap. Local seismic investigation in the Assam gap confirms this feature and the seismicity suggests the existence there of an asperity.The local seismic investigations in Garhwal Himalaya have shown that the small earthquakes are confined to the upper 6–8 km of the crust and may have strike-slip motions. These earthquakes occur in a region where teleseismically recorded events were few.     

Assessment of predominant frequencies using ambient vibration in the Kachchh region of western India: implications for earthquake hazards

The Kachchh region is the second most seismically active region in India after the Himalaya. One of the disastrous Indian earthquakes of the millennium was the Bhuj earthquake of January 26, 2001, which caused about 14,000 casualties and huge property damage. The main reason for such devastation is due to lack of earthquake awareness and poor construction practices. Hence, an increase in the knowledge and awareness, based on improved seismic hazard assessment, is required to mitigate damage due to an earthquake. Natural predominant ground frequencies have been investigated in the Kachchh region of western India using ambient vibrations. The horizontal-to-vertical spectral ratio technique has been applied to estimate the predominant frequency at 126 sites. The ambient vibration measurements were conducted for about 1 h at each site in the continuous mode recording at 100 samples/s. We have validated the estimated predominant frequency with earthquake data recorded at six broadband stations in the region. It has been observed that geological time period has a significant effect on predominant frequency of the ground. The estimated predominant frequencies vary from 0.24 to 2.25 Hz for the Quaternary, 0.41–2.34 Hz for the Tertiary, 0.32–4.91 Hz for the Cretaceous, and 0.39–8.0 Hz for the Jurassic/Mesozoic. In the Deccan trap, it varies from 1.30 to 3.80 Hz. We found distinct variation of predominant frequencies of sites associated with hard rock and soft soil. The predominant frequencies were related to the thickness of the sediments, which are deduced by other geophysical and geological methods in the region. Our results suggest that frequencies of the region reveals the site characteristics that can be considered for studying the seismic risks to evolve a plan for disaster risk mitigation for the region.     

Deformation and displacement fields of the Himalayan arc,seismicity and large earthquake hazards in its eastern segment and vicinity

The Himalayan arc is a type of plate margin similar to an island arc and is a world-famous region of tectonic and seismic activities, where a series of large earthquakes have occurred in historical time. In this paper, the vertical deformation and horizontal displacement fields of the Himalayan arc are theoretically derived from the viewpoint of the collision between the Indian and Eurasian plates. In the light of the observed data, the seismicity, earthquake focal mechanism, seismotectonic and geomorphological features of the arc and its vicinity are reasonably explained. The characteristics of seismicity and the possibility of earthquakes with magnitude above 8 occurring in this region in the future are studied.     

Simulation of strong ground motion for a potential Mw7.3 earthquake in Kopili fault zone,northeast India

Natural Hazards - The occurrence of major natural disasters in recent years has impacted large cities worldwide and boosted the need of assessing urban resilience. As a key factor of resilience,...     

Geophysical methods applied to fault characterization and earthquake potential assessment in the Lower Tagus Valley, Portugal

The study region is located in the Lower Tagus Valley, central Portugal, and includes a large portion of the densely populated area of Lisbon. It is characterized by a moderate seismicity with a diffuse pattern, with historical earthquakes causing many casualties, serious damage and economic losses. Occurrence of earthquakes in the area indicates the presence of seismogenic structures at depth that are deficiently known due to a thick Cenozoic sedimentary cover. The hidden character of many of the faults in the Lower Tagus Valley requires the use of indirect methodologies for their study. This paper focuses on the application of high-resolution seismic reflection method for the detection of near-surface faulting on two major tectonic structures that are hidden under the recent alluvial cover of the Tagus Valley, and that have been recognized on deep oil-industry seismic reflection profiles and/or inferred from the surface geology. These are a WNW–ESE-trending fault zone located within the Lower Tagus Cenozoic basin, across the Tagus River estuary (Porto Alto fault), and a NNE–SSW-trending reverse fault zone that borders the Cenozoic Basin at the W (Vila Franca de Xira–Lisbon fault). Vertical electrical soundings were also acquired over the seismic profiles and the refraction interpretation of the reflection data was carried out. According to the interpretation of the collected data, a complex fault pattern disrupts the near surface (first 400 m) at Porto Alto, affecting the Upper Neogene and (at least for one fault) the Quaternary, with a normal offset component. The consistency with the previous oil-industry profiles interpretation supports the location and geometry of this fault zone. Concerning the second structure, two major faults were detected north of Vila Franca de Xira, supporting the extension of the Vila Franca de Xira–Lisbon fault zone northwards. One of these faults presents a reverse geometry apparently displacing Holocene alluvium. Vertical offsets of the Holocene sediments detected in the studied geophysical data of Porto Alto and Vila Franca de Xira–Lisbon faults imply minimum slip rates of 0.15–0.30 mm/year, three times larger than previously inferred for active faults in the Lower Tagus Valley and maximum estimates of average return periods of 2000–5000 years for M 6.5–7 co-seismic ruptures.