State of the Art of Seismic Zonation Map in United States, Site Effect Model and Basin Effect Model
-
摘要: 2018版《美国地震危险性图》对地震动模型美国中部和东部地区地震动场地效应模型及西部地区沉积盆地放大等进行了更新,改善了对美国各地地震灾害的描述,并增强了对美国中部、东部地区与西部地区地震动差异的理解。本文分析并论述了2018版《美国地震危险性图》中美国中部和东部地区地震动场地效应模型及西部地区沉积盆地效应模型的主要特征,对我国新一代地震动区划图的编制、修订进行了思考。Abstract: The 2018 US National Seismic Hazard Model has updated the site effect models of ground motion in the central and eastern United States, and the amplification effect of sedimentary basins in the western United States, improving the description of seismic disasters across the United States and the understanding of ground motion differences between the central and eastern United States and the western United States regions. In this paper, the main characteristics of the site effect models of ground motion in the central and eastern United States and the sedimentary basin effect model in the western United States are analyzed and discussed, which triggers the deliberation of seismic zonation map of next generation in China.
-
引言
跨断层形变监测是一种观测物理意义明确、测量精度高、便于组织实施以及数据产品稳定可靠的地震前兆监测方法之一(薄万举等, 1994, 1998;李杰等,2010;张晶等,2011;张希等,2014;马伶俐等,2015;李腊月等,2016;孙启凯等,2017)。通过多期次、短距离的跨断层水准和基线测量,能够准确记录断层活动的位移量,监测活动断裂的运动特征,更好地反映断裂带及其所处区域的形变和运移特征(牛安福等,2013)。针对川滇地块的地壳形变、运动学特征及其形成机制开展的研究工作产出了大量成果,但是,目前关于这些问题也存在不同认识(乔学军等,2004;武艳强等,2011;洪敏等,2014;吴中海等,2015;白玉柱等,2016)。本文利用川滇地块南部附近的跨断层观测资料,对区域断裂带的现今形变变化进行分析,研究川滇地块南部的整体运动学特征,综合提取断层活动的异常信息,为区域地震危险性分析提供一定的参考。
1. 区域地质概况
川滇地块位于青藏高原的东缘,南北地震带的南段,大地构造上处于扬子板块、松潘-甘孜地块、兰坪-思茅地块以及冈底斯-腾冲地块等构造单元的拼合汇聚地带,在地质历史上经历了多次不同时期不同规模的构造运动与岩浆活动,构造属性复杂。该区域地震活动强烈而频繁,地震与构造运动之间存在复杂的因果关系,是研究地震活动与断层运动关系的热点地区(图 1)。
自新生代以来,印度板块与亚欧板块的相互碰撞拼贴,造成地壳强烈缩短增厚,青藏高原得以快速隆升(Tapponnier等,2001;许志琴等,2011)。伴随着两大板块的俯冲拼合,高原也发生着同构造挤压伸展,物质和能量向北西、北东、东、南东等多个方向转移释放。这种持续的俯冲碰撞使得两大板块拼合区域的喜马拉雅造山带、雅鲁藏布江缝合带以及高原内部地块与周缘区域构造活动强烈,并形成一系列大断裂带,如北西走向的鲜水河断裂、红河断裂及近南北走向的小江断裂、安宁河-则木河断裂带等。很多学者对此区域开展了研究,并取得了大量成果(Wang等,1998;Burchfiel等,2003;陈长云等,2016)。其中,处于川滇地块南部的活动断裂,是本文关注的重点。
2. 跨断层监测概况
川滇地块南部断裂所在的云南地区的跨断层监测场地从20世纪80年代开始观测,至今30余年,主要为流动观测场地,共有11处,分别是楚雄、峨山、建水、剑川、丽江、石屏、通海、下关、永胜、宜良和羊街,主要分布于金沙江-红河断裂带、小江断裂带、南华-楚雄-建水断裂带及其次级断裂(剑川断裂、洱海断裂、通海-峨山断裂、建水-石屏断裂、丽江断裂)上(表 1),多为综合观测场地。每个场地包括2条短水准和2条短基线测线,分别各有1条直交和斜交测线,通过这种组合方式来共同监测场地所跨断层的三维运动变化。从场地观测历史来看,该区域的跨断层观测基本是每月1期,个别年份每2个月观测1期,数据连续且完整。
表 1 研究区跨断层监测场地概况Table 1. Basic information of the cross-fault monitoring sites in the area场地 经度/°E 纬度/°N 所跨断层 断层产状 断层性质 楚雄 101.52 25.03 楚雄断裂 NE∠58° 右旋挤压 峨山 102.50 24.12 通海-峨山断裂 S∠60° 右旋挤压 建水 102.81 23.63 石屏-建水断裂 S∠70° 左旋挤压 剑川 99.93 26.63 剑川断裂 W∠75° 左旋拉张 丽江 100.25 26.88 丽江断裂 SW∠70° 左旋挤压 石屏 102.45 23.77 石屏-建水断裂 NE∠70° 右旋拉张 通海 102.72 24.00 通海-峨山断裂 NE∠68° 右旋挤压 下关 100.30 25.62 红河断裂 NE∠60° 右旋拉张 永胜 100.72 26.70 程海断裂 W∠35° 左旋挤压 宜良 103.03 25.02 马街-南羊街断裂 SE∠68° 左旋挤压 羊街 103.10 25.55 小江断裂西支 SE∠73° 左旋挤压 3. 数据计算及分析
将收集到的研究区30余年的跨断层观测数据进行预整理,对数据格式做统一标准化处理。首先,消除跨断层观测资料中明确因干扰引起的突跳变化,对观测资料进行连续性检验,剔除重复的数据。然后,对于日常观测中因各种原因造成缺测的数据采用线性插值法修补,得到周期统一的、可靠的形变数据。对数据的处理分析采用以下方法:
云南区域内的跨断层场地同时布设有两条同桩的基线和水准测线,可采用下面3个公式来计算断层活动的三维分量:
$$\Delta S=\frac{\Delta {{L}_{1}}\text{sin}{{\alpha }_{2}}-\Delta {{L}_{2}}\sin {{\alpha }_{1}}}{\sin ({{\alpha }_{2}}-{{\alpha }_{1}})}$$ (1) $$\Delta R=\frac{\Delta {{L}_{1}}-\Delta S\cos {{\alpha }_{1}}}{\text{sin}{{\alpha }_{1}}}$$ (2) $$\Delta H=\frac{1}{2}(\Delta {{H}_{1}}+\Delta {{H}_{2}})$$ (3) 其中,$\Delta H$为水准垂直变化量(以上盘相对于下盘的下降为正),$\Delta L$为基线伸缩量(以伸长为正),$\Delta S$为断层活动水平走滑量(两盘做顺扭运动为正),$\Delta R$为断层活动水平张压量,${{\alpha }_{1}}$和${{\alpha }_{2}}$分别代表由断层线方向逆时针转动至与两条测线方向重合或者平行时的角度。当$\Delta S$为正值时,表示右旋活动,$\Delta R$为正值时,表示张性活动,若均为负值则表示断层活动性质相反(薄万举等,1998)。
利用MATLAB对计算得出的数据绘制曲线图,单幅曲线图上方均标注有其代表的场地名及运动方式(图 2)。通过曲线图可以发现跨断层场地所观测的不同断裂表现出了跟地质背景相同或者不同的运动学特征,其中丽江、石屏和下关场地所跨断层的现今活动特征跟地质背景一致;楚雄、剑川和永胜场地所跨断层在走滑和拉张特征上都表现出跟背景特征不一致的现象;峨山、建水和通海场地所观测的断层在走滑性质上和地质背景相符,而在张压性质上相反。根据现今运动特征,可将场地分为3种类型:
(1)运动特征与地质背景一致。丽江场地所跨的丽江次级断裂,在地质背景上是一处表现为左旋走滑特征的活动断裂,兼具挤压特点,现今计算的三维活动量,依然表现出相同特征,其中断层的水平张压量在2012年后出现大幅异常变化。石屏场地所跨的石屏-建水断裂,在地质背景上是一处右旋拉张特征的活动断裂,现今右旋运动速率在2013年后虽然有所放缓,但是仍然与地质背景一致。下关场地所处的是一条表现为右旋走滑兼有拉张特征的深断裂,2003年水平走滑量曾出现右旋剪切转折,现今曲线重回向上趋势,继续保持右旋走滑特征。羊街和宜良场地所处的小江断裂带西支断裂,现今数据仍显示出以继承性的左旋走滑兼具挤压特征为主。
(2)运动特征与地质背景不一致。楚雄场地所跨断层为具有右旋挤压背景特征的楚雄断裂,在2008年后现今运动特征以强烈左旋拉张为主,但自2014年左旋和拉张的运动速率都出现大幅放缓的异常现象。剑川场地所处断裂为红河断裂带北段的剑川断裂,是一处左旋走滑兼具拉张特征的活动断裂,现今三维运动特征表现为右旋走滑兼具挤压特征,其中的垂直位移量自2015年至今出现了大幅抬升的异常变化;水平走滑量自2011年后也出现大幅抬升,右旋运动加速。永胜场地所监测的断裂是具有左旋挤压特征的程海深断裂,从计算所得的三维活动量曲线可以看出,该断裂现在运动学特征为右旋拉张,其垂直位移量曲线从2014年开始大幅异常抬升,水平走滑量自1994年开始稳步抬升。
(3)与地质背景相比,走滑特征一致,张压特征相反。峨山场地和通海场地共同监测曲江断裂的不同部位,该断裂在地质背景上具有右旋挤压特征,两者现今走滑特征与背景一致,但是张压量曲线表现出相反的拉张运动的异常特点。其中,峨山场地走滑分量在1998年后大幅抬升,之后趋于稳定的右旋运动状态;通海场地的垂直位移量在1992年后出现转折,断层表现为张性活动,走滑分量在2003年大幅下降,断层剪切活动速率减小,之后趋于稳定,依然显示右旋运动。建水场地所跨的石屏-建水断裂的建水段,断层性质为左旋挤压,从计算得出的水平走滑量参数来看,现今的走滑运动特征与地质背景一致,垂直位移量在1999年后发生趋势改变,断层活动由压性转变为相反的拉张特征。
4. 讨论与结论
在印度洋板块和亚欧板块的互相作用下,巨大的青藏高原造就了地球的“第三极”,这种作用力至今仍在持续。板块的碰撞拼合产生了强烈的构造活动,许多大规模断裂带应运而生(Xu等,2016),同时,也孕育了一系列强震。对此一些学者做过许多研究,并取得大量有价值的科研成果(Tapponnier等,2001;许志琴等,2011;陈长云等,2016)。尽管目前对青藏高原的陆内侧向运动形式、时间和逃逸量存在着争议,但是对位于其东缘的川滇地块发生的晚新生代南东向的逃逸和围绕东喜马拉雅构造结顺时针的旋转有着一致的认识(乔学军等,2004;武艳强等,2011;洪敏等,2014;吴中海等,2015)。
通过对位于川滇地块南部主要断裂上的跨断层场地三维形变量的分析,可以看出部分断裂表现出继承性特征,如现今下关和石屏场地的三维活动量表现出跟以往相同的持续平稳的活动速率。而另一部分断裂表现出了异常于背景活动特征的现象,如楚雄、永胜和剑川场地在走滑分量上都表现出了反向活动特点。研究区具有以水平运动为主、垂直运动为辅的背景特征,现今看来,依然如此。通过计算分析得到的宜良、羊街、楚雄以及丽江场地的形变结果来看,羊街、宜良场地所控制的小江断裂带分支为川滇地块东部边界断裂,现呈现左旋走滑兼具压性的特征;楚雄场地监控的楚雄-南华断裂位于川滇地块的南边部,表现出左旋拉张的运动特征;丽江场地所处的丽江-小金河断裂作为川滇地块南北部的分界断裂,现今以左旋挤压运动为主。可以看出,川滇地块南部区域整体上表现出以水平顺时针旋转为主的运动方式。根据GPS观测资料得到的该区域现今水平运动速度场也显示了相同的运动特点(乔学军等,2004;武艳强等,2011;洪敏等,2014;王伶俐等,2015)。造成这种运动方式的深层次动力学机制,应该源于印度洋板块俯冲至亚欧板块下,碰撞拼合后巨大高耸的青藏高原崛起,在超级重力作用下,物质开始向周缘,尤其是向东缘流动,在遇到横亘的北东向龙门山造山带时受阻,流动速率下降,被消化、吸收,再次转移(张家声等,2003;崔效锋等,2006)。其中,相当一部分能量和物质向南东运移,规模相对较小的川滇活动地块随之发生南东向旋转,同时起着调节周围各大板块和地块之间运动平衡的“交通环岛”的作用。另外,印度洋板块至今持续强烈地北向运动,跟扬子板块形成了大尺度上的右旋运动,也推动着川滇地块南东向旋转运移。
根据观测场地所跨断层的三维运动特征及计算所得曲线可以看出,在左旋的小江断裂带、右旋的红河断裂带和右旋的南华-楚雄-建水断裂带共同汇聚的峨山、通海和建水场地区域,原先的地壳形变表现出压性特征,现今由于断裂性质的异常变化,该区域场地集体表现出拉张的运动特点。滇西北区域的永胜、丽江和剑川场地所跨断层也表现出了一定的异常运动特征,剑川断裂和程海断裂呈水平反向运动,除此之外,包括丽江场地所跨的丽江断裂在内,3条断裂的垂向运动特征也存在大幅异常变化。云南地区中强震发生前的一定时段或震时,区域内跨断层场地观测到的水准和基线数据多数会出现较明显的变化,且部分场地形变量会发生速率的急变或转折,如文山5.3级、孟连西7.3级、姚安6.5级、丽江7.0级地震震前和震时,部分场地就表现出了加速继承性运动和转折性运动的特征(张兴华等,1996;付虹等,1997;李忠华等,1998;施顺英等,2005;李瑞莎等,2009;刘强等,2010;王永安等,2011)。同时,云南地区2015年至今存在长时间缺少4级以上中强地震的背景,不符合区域以往的地震活动特征,现今应力积累较高。自2016年云龙5.0级地震后至今,区域中强震活动已平静300天左右;滇西北地区自2001年永胜6.0级地震至今,已14.7年未发生6级以上地震;滇中南的石屏、建水等区域从2001年江川5.1级地震起算,5级地震平静期已长达15.5年。基于此,该区域在中短期尺度可能存在一定中强震的危险。
综上所述,得到以下几点结论:
(1)川滇地块南部跨断层场地的形变特征表现出继承和异常于地质背景特征的现象。① 丽江、石屏、下关、宜良和羊街场地所跨断层的现今活动特征与地质背景一致。丽江断裂和小江断裂表现出持续地左旋挤压运动特征;石屏-建水断裂石屏段表现出右旋拉张特征;红河断裂表现出跟地质背景相同的右旋拉张特征。② 楚雄、剑川和永胜场地所跨断层在走滑和张压特征上都出现与背景特征相背的现象。楚雄断裂、剑川断裂和程海断裂均表现出相反的运动特征。③ 峨山、建水和通海场地所监测的断层现今活动状态与地质背景相比,在走滑性质上相符,张压性质上相反。通海-峨山断裂现今运动特征,表现出与地质背景相同的右旋走滑和与地质背景相反的拉张特征。石屏-建水断裂建水段现今依然呈左旋走滑运动,但是张压性质发生反转,表现出张性运动的活动特点。
(2)川滇地块南部整体上表现出以水平顺时针旋转为主的运动方式。
(3)根据跨断层观测资料所得出的断层三维运动特征及曲线可知,滇西北及滇中的峨山、通海、石屏和建水等区域地壳形变特征存在异常,应关注该区域中短期尺度内可能存在的中强震危险。
致谢: 对编辑和审稿老师的辛苦审阅和宝贵建议,表示衷心感谢。 -
图 1 场地放大系数分布情况(Petersen等,2020)
Figure 1. Distribution of site magnification factors (Petersen et al., 2020)
图 2 不同模型、不同场地类别条件下地震动中值加权平均值分布情况(Rezaeian等, 2021)
Figure 2. Comparison of the 2018 GMM medians at original models and smoothed models, a range of site classes (Rezaeian et al., 2021)
图 3 NGA-West2的4个地震动模型在T=5 s时盆地效应系数与盆地深度的关系(Powers等,2021)
Figure 3. The relationship between basin effect coefficient and basin depth at T=5 s for the four ground motion models of NGA-West2 (Powers et al., 2021)
图 4 原始发布和 2018年美国地质勘探局修改后实施的地震动模型在不同盆地深度条件下对应的加速度反应谱中值对比结果(Powers等,2021)
Figure 4. The models as published and the 2018 USGS implementation of the median spectral acceleration at different basin depth conditions(Powers et al., 2021)
图 5 原始发布和 2018年美国地质勘探局修改后实施的地震动模型在不同盆地深度条件下对应的加速度反应谱中值对比结果(Powers等,2021)
Figure 5. The models as published and the 2018 USGS implementation of the median spectral acceleration at different basin depth conditions (Powers et al., 2021)
表 1 场地类别和盆地深度组合(Powers等,2021)
Table 1. Site class and basin-depth combinations for plots ( Powers et al., 2021)
场地类别-盆地 Vs30/(m·s−1) Z1.0/km(ASK14、BSSA14、CY14模型) Z2.5/km(CB14模型) BC类-默认盆地($ {V}_{\mathrm{s}30} $-基础) 760 (0.048,0.041,0.041) 0.607 BC类-深盆地(长滩) 760 (0.704,0.704,0.704) 3.830 BC类-浅盆地(旧金山) 760 (0.025,0.025,0.025) 0.850 DE类-默认盆地($ {V}_{\mathrm{s}30} $-基础) 185 (0.497,0.513,0.513) 3.060 DE类-深盆地(长滩) 185 (0.704,0.704,0.704) 3.830 DE类-浅盆地(旧金山) 185 (0.025,0.025,0.025) 0.850 -
中华人民共和国住房和城乡建设部, 中华人民共和国国家质量监督检验检疫总局, 2010. GB 50011—2010 建筑抗震设计规范. 北京: 中国建筑工业出版社, 18—20Ministry of Housing and Urban-Rural Development of the People’s Republic of China, General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China, 2010. GB 50011—2010 Code for seismic design of buildings. Beijing: China Architecture & Building Press, 18—20. (in Chinese) 周健, 李小军, 李亚琦等, 2021. 中美建筑抗震设计规范中工程场地类别的对比和换算关系. 地震学报, 43(4): 521—532 doi: 10.11939/jass.20200164Zhou J. , Li X. J. , Li Y. Q. , et al. , 2021. Comparative analysis and transformation relations between China and the US site classification systems in building seismic code provisions. Acta Seismologica Sinica, 43(4): 521—532. (in Chinese) doi: 10.11939/jass.20200164 Aagaard B. T. , Brocher T. M. , Dolenc D. , et al. , 2008. Ground-motion modeling of the 1906 San Francisco earthquake, part II: ground-motion estimates for the 1906 earthquake and scenario events. Bulletin of the Seismological Society of America, 98(2): 1012—1046. doi: 10.1785/0120060410 Abrahamson N. A. , Silva W. J. , Kamai R. , 2014. Summary of the ASK14 ground motion relation for active crustal regions. Earthquake Spectra, 30(3): 1025—1055. doi: 10.1193/070913EQS198M Algermissen S. T., Perkins D. M., 1976. A probabilistic estimate of maximum acceleration in rock in the contiguous United States. Reston: U. S. Geological Survey. Algermissen S. T., Perkins D. M., Thenhaus P. C., et al., 1990. Probabilistic earthquake acceleration and velocity maps for the United States and Puerto Rico. Reston: U. S. Geological Survey. Atkinson G. M. , Boore D. M. , 2008. Erratum to empirical ground-motion relations for subduction zone earthquakes and their application to Cascadia and other regions. Bulletin of the Seismological Society of America, 98(5): 2567—2569. doi: 10.1785/0120080108 Atkinson G. M. , Boore D. M. , 2011. Modifications to existing ground-motion prediction equations in light of new data. Bulletin of the Seismological Society of America, 101(3): 1121—1135. doi: 10.1785/0120100270 Boore D. M. , Stewart J. P. , Seyhan E. , et al. , 2014. NGA-West2 equations for predicting PGA, PGV, and 5% damped PSA for shallow crustal earthquakes. Earthquake Spectra, 30(3): 1057—1085. doi: 10.1193/070113EQS184M Campbell K. W. , Bozorgnia Y. , 2014. NGA-West2 ground motion model for the average horizontal components of PGA, PGV, and 5% damped linear acceleration response spectra. Earthquake Spectra, 30(3): 1087—1115. doi: 10.1193/062913EQS175M Chiou B. S. J. , Youngs R. R. , 2014. Update of the Chiou and Youngs NGA model for the average horizontal component of peak ground motion and response spectra. Earthquake Spectra, 30(3): 1117—1153. doi: 10.1193/072813EQS219M Frankel A. , Wirth E. , Marafi N. , et al. , 2018. Broadband synthetic seismograms for magnitude 9 earthquakes on the Cascadia megathrust based on 3 D simulations and stochastic synthetics, part 1: methodology and overall results. Bulletin of the Seismological Society of America, 108(5 A): 2347—2369. doi: 10.1785/0120180034 Frankel A. D., Mueller C. S., Barnhard T. P., et al., 1996. National seismic-hazard maps: documentation June 1996. Reston: U. S. Geological Survey. Frankel A. D. , Carver D. L. , Williams R. A. , 2002 a. Nonlinear and linear site response and basin effects in Seattle for the M 6.8 Nisqually, Washington, earthquake. Bulletin of the Seismological Society of America, 92(6): 2090—2109. doi: 10.1785/0120010254 Frankel A. D., Petersen M. D., Mueller C. S., et al., 2002 b. Documentation for the 2002 update of the national seismic hazard maps. Reston: U. S. Geological Survey. Goulet C. A., Bozorgnia Y., Abrahamson N., et al., 2018. Central and eastern North America ground-motion characterization-NGA-east final report. Berkeley: Pacific Earthquake Engineering Research Center. Graves R. , Jordan T. H. , Callaghan S. , et al. , 2011. CyberShake: a physics-based seismic hazard model for southern California. Pure and Applied Geophysics, 168(3—4): 367—381. doi: 10.1007/s00024-010-0161-6 Hartzell S. , Cranswick E. , Frankel A. , et al. , 1997. Variability of site response in the Los Angeles urban area. Bulletin of the Seismological Society of America, 87(6): 1377—1400. doi: 10.1785/BSSA0870061377 Hashash Y. M. A. , Ilhan O. , Harmon J. A. , et al. , 2020. Nonlinear site amplification model for ergodic seismic hazard analysis in central and eastern North America. Earthquake Spectra, 36(1): 69—86. doi: 10.1177/8755293019878193 Idriss I. M. , 2014. An NGA-West2 empirical model for estimating the horizontal spectral values generated by shallow crustal earthquakes. Earthquake Spectra, 30(3): 1155—1177. doi: 10.1193/070613EQS195M Moschetti M. P. , Hartzell S. , Ramírez-Guzmán L. , et al. , 2017. 3 D Ground-motion simulations of Mw7 earthquakes on the Salt Lake City segment of the Wasatch fault zone: variability of long-period (T≥1 s) ground motions and sensitivity to kinematic rupture parameters. Bulletin of the Seismological Society of America, 107(4): 1704—1723. Pacific Earthquake Engineering Research Center (PEER), 2015. NGA-east: adjustments to median ground-motion models for central and eastern North America. Berkeley: Pacific Earthquake Engineering Research Center. Petersen M. D., Frankel A. D., Harmsen S. C., et al., 2008. Documentation for the 2008 update of the United States national seismic hazard maps. Reston: U. S. Geological Survey. Petersen M. D., Moschetti M. P., Powers P. M., et al., 2014. Documentation for the 2014 update of the United States national seismic hazard maps. Reston: U. S. Geological Survey. Petersen M. D. , Shumway A. M. , Powers P. M. , et al. , 2020. The 2018 update of the US national seismic hazard model: overview of model and implications. Earthquake Spectra, 36(1): 5—41. doi: 10.1177/8755293019878199 Petersen M. D. , Shumway A. M. , Powers P. M. , et al. , 2021. The 2018 update of the US national seismic hazard model: where, why, and how much probabilistic ground motion maps changed. Earthquake Spectra, 37(2): 959—987. doi: 10.1177/8755293020988016 Powers P. M. , Rezaeian S. , Shumway A. M. , et al. , 2021. The 2018 update of the US national seismic hazard model: ground motion models in the western US. Earthquake Spectra, 37(4): 2315—2341. doi: 10.1177/87552930211011200 Rezaeian S. , Powers P. M. , Shumway A. M. , et al. , 2021. The 2018 update of the US national seismic hazard model: ground motion models in the central and eastern US. Earthquake Spectra, 37(S1): 1354—1390. Shumway A., Petersen M. D., Powers P. M., et al., 2018. Additional period and site class maps for the 2014 national seismic hazard model for the conterminous United States. Reston: U. S. Geological Survey. Stephenson W. J. , Frankel A. D. , Odum J. K. , et al. , 2006. Toward resolving an earthquake ground motion mystery in west Seattle, Washington State: shallow seismic focusing may cause anomalous chimney damage. Geophysical Research Letters, 33(6): L06316. Stewart J. P. , Parker G. A. , Atkinson G. M. , et al. , 2020. Ergodic site amplification model for central and eastern North America. Earthquake Spectra, 36(1): 42—68. doi: 10.1177/8755293019878185 -