• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    利用連續(xù)重力觀測約束2014智利Iquique地震的震源機制解

    2016-07-29 02:04:43江穎李輝劉子維申重陽張曉彤
    地球物理學報 2016年3期
    關鍵詞:標量臺站震源

    江穎,李輝,劉子維,申重陽,張曉彤

    1 中國地震局地震研究所(中國地震局地震大地測量重點實驗室),武漢 430071 2 中國科學院測量與地球物理研究所大地測量與地球動力學國家重點實驗室,武漢 430077

    ?

    利用連續(xù)重力觀測約束2014智利Iquique地震的震源機制解

    江穎1,2,李輝1,劉子維1,申重陽1,張曉彤1

    1 中國地震局地震研究所(中國地震局地震大地測量重點實驗室),武漢4300712中國科學院測量與地球物理研究所大地測量與地球動力學國家重點實驗室,武漢430077

    摘要2014年4月1日,智利北部Iquique地區(qū)近海發(fā)生MW8.1地震,地震發(fā)生之后,國際上一些著名的地震科研機構和學者采用不同的數(shù)據(jù)和方法計算得到此次地震的震源機制解,但這些結果存在較大差異.地球長周期自由振蕩的振幅主要依賴于地震矩的大小及地震斷層的破裂方式,可以很好地約束地震震源機制.因此,本文根據(jù)2014年Iquique地震現(xiàn)有的6個不同震源機制解模擬計算了該地震激發(fā)的自由振蕩信號,并與全國連續(xù)重力臺網中16個彈簧重力儀的觀測結果進行比對,基于1.5~5.3 mHz的球型簡正模分析和約束了Iquique地震的震源機制解.研究發(fā)現(xiàn),基于美國地質調查局WPhase Moment Tensor Solution反演的震源機制解的自由振蕩模擬值與實際觀測符合最好,其相應的震級能較好反映Iquique地震釋放的總能量,而利用海嘯數(shù)據(jù)反演的標量地震矩偏小,聯(lián)合遠場和近場長周期觀測數(shù)據(jù)反演可顯著改善震源機制解.另外,還基于格爾木重力臺站的模擬值與觀測值定量分析了不同震源機制解參數(shù)對自由振蕩振幅的影響.結果表明地震的標量地震矩M0對自由振蕩振幅的影響最大,而斷層走向、傾角、滑動方向角和震源深度對自由振蕩的振幅影響相對較小.

    關鍵詞Iquique地震; 震源機制解; 地球自由振蕩; 重力觀測

    1引言

    智利西部沿海處于納斯卡板塊和南美板塊之間,由于納斯卡板塊以約8 cm·a-1的收斂速度向南美板塊俯沖(Angermann et al.,1999),該地區(qū)大地震頻發(fā).2014年4月1日23∶46 (UTC)在智利西北部Iquique近海地區(qū)發(fā)生了震級為MW8.1的地震,震中位為19.642°S,70.817°W(圖1,紅色五角星位置),震源深度25 km (USGS National Earthquake Information Center-NEIC: http:∥earthquake.usgs.gov/regional/neic/).地震發(fā)生后,Iquique地震的震源性質和發(fā)震構造吸引了許多學者的關注和研究(An et al.,2014;Lay et al.,2014;Ruiz et al.,2014;Schurr et al.,2014;Kato et al.,2014;Yagi et al.,2014;Bürgmann et al.,2014;Hayes et al.,2014;Liu et al.,2015).國際上一些著名的地震科研機構和學者采用不同的數(shù)據(jù)和方法計算得到Iquique地震的震源機制解,見表1.其中震源機制解1是GCMT(The Global Centroid-Moment-Tensor Project)通過計算單位矩張量元的格林函數(shù),將格林函數(shù)與震源時間函數(shù)做卷積,運用最小二乘求得最優(yōu)矩張量解.GCMT采用體波(截斷周期為50 s),長周期面波(截斷周期為50 s)和地幔波(截斷周期為200 s)計算了Iquique地震的震源機制解.震源機制解2、3是美國地質調查局(USGS)的結果,USGSCMT(USGS Centroid Moment Tensor Solution)采用體波和長周期面波反演震源機制,USGSWPHASE(USGS WPhase Moment Tensor Solution)采用W-phase(100~1000 s)的波形資料反演震源機制.對于中級地震,USGSBODY(USGS Body-Wave Moment Tensor Solution)采用遠場臺站記錄的體波數(shù)據(jù)反演震源機制,比如2013年的蘆山地震.單純用體波數(shù)據(jù)反演的標量地震矩偏小(Jiang et al.,2014a).通常情況下USGSBODY不提供大地震的震源機制解結果,如2011年日本Tohoku MW9.1地震等,包括此次Iquique MW8.1地震.Yagi 等(2014)利用48個臺站的遠震P波數(shù)據(jù),考慮了格林函數(shù)的不確定性反演了Iquique地震的震源機制解.Lay等(2014)利用遠程P波、SH波及海嘯數(shù)據(jù)進行有限斷層反演,得到的標量地震矩為1.7×1021N·m.An等(2014)利用3個深海評估和海嘯臺站的海嘯數(shù) 據(jù),反演得到了Iquique地震的滑動分布,并得到該地震的標量地震矩為1.22×1021N·m.從上述研究結果來看,不同的研究機構和學者給出的Iquique地震的震源機制解相差較大,特別是標量地震矩M0之間的差異顯著,最大值約為最小值的2倍,而這種差異不容忽視.

    表1 智利Iquique地震震源機制解一覽表

    注:文獻來源震源機制解1為Global-CMT Project的結果;震源機制解2和3分別為USGS提供2個不同的結果:USGSCMT、USGSWPHASE;4為Yagi et al.(2014)的結果;5為Lay et al.(2014)的結果;6為An et al.(2014)的結果.

    地震在瞬時釋放巨大能量、產生大面積地表破裂的同時,激發(fā)大量的地震波,并且導致地球整體的震動.地球的自由振蕩一旦被激發(fā),將以駐波形式穿過地球的內部介質.由于自由振蕩能量的大小與震源的破裂方式和破裂程度密切相關(Geller and Stein,1979;Stein et al.,2005),因此地震后的自由振蕩信號可用于約束地震震級和檢驗地震的震源機制解,它也是目前可對震源機制解進行總體檢驗和約束的非常有效的方法.利用自由振蕩信號約束MW9.0級以上大地震的震源機制解已經得到了不錯的效果,比如2004-12-26蘇門答臘大地震后,利用應變儀、長周期地震儀及超導重力儀記錄到的球型低頻模態(tài)0S0、0S2、0S3、0S4等,徑向模0S0、1S0等與理論模擬結果對比,發(fā)現(xiàn)相應的標量地震矩是GCMT結果的2.5~2.67倍,并且得出斷層的破裂長度、持續(xù)時間及破裂速度(Stein et al.,2005;Park et al.,2005;Lambotte et al.,2006).2011-03-11Tohoku大地震后,薛秀秀等(2012)利用0S0振型對Tohoku大地震的多個震源機制解進行分析和約束.對于中級地震來講低頻態(tài)的信號不易提取,Jiang等(2014)根據(jù)2013-04-20蘆山地震的4種不同震源機制解計算了蘆山地震頻段為2.3~5 mHz的自由振蕩信號,發(fā)現(xiàn)利用遠場體波反演的地震矩偏小.Hu等(2014)為了評估蘆山地震標量地震矩的正確性,對比了17個臺站和5個矩張量解的觀測值與模擬值,也得到了同樣的結論.目前,利用自由振蕩信號約束震源機制解所用到的觀測數(shù)據(jù)均來自超導重力儀(SG)和寬頻帶地震儀,當震級較大時,SG會發(fā)生限幅情況,因此在計算過程中必須去掉同震時的觀測段,對于MW8.1 Iquique地震來說,中國大陸的幾臺超導重力儀均出現(xiàn)了限幅情況,而寬頻帶地震儀主要用于記錄高頻信息.

    通過數(shù)字地震觀測網絡、中國大陸構造環(huán)境監(jiān)測網絡、中國地震背景場等國家重大科學工程項目,在中國大陸建成了一個集數(shù)據(jù)的自動匯集與數(shù)據(jù)庫式管理,儀器的遠程監(jiān)控和數(shù)據(jù)的圖形化分析處理為一體的連續(xù)重力觀測臺網,其平均空間分辨率約400 km.觀測儀器包括PET/gPhone、DZW、GS15、TRG-1和GWR超導等,除了SG以外,這些彈簧重力儀也可以觀測到微小的重力變化,分布廣泛且連續(xù)的彈簧重力儀觀測數(shù)據(jù)為檢測地球自由振蕩信號及利用自由振蕩信號約束震源機制解提供了很好的數(shù)據(jù)基礎.本文選取連續(xù)重力觀測臺網中分布廣泛的16個彈簧重力儀觀測數(shù)據(jù)(如圖1),利用自由振蕩信號觀測約束Iquique地震的震源機制解,并有效的評估哪一個震源機制解最符合地震實際釋放的能量,基于符合較好的重力站觀測資料,定量研究和分析影響自由振蕩振幅的因素.

    圖1 震中位置和16個連續(xù)重力觀測臺站的位置

    2數(shù)據(jù)和方法

    地球自由振蕩中的球型振蕩不僅有橫向分量,還有徑向分量,因此會引起地球的體積變化,從而導致重力發(fā)生變化.中國連續(xù)重力觀測臺網提供采樣率為1 Hz的重力觀測數(shù)據(jù)用于地球動力學研究.Iquique地震發(fā)生后,連續(xù)重力觀測臺網的重力儀清晰的記錄了重力變化.提取自由振蕩信號之前,首先對各臺站的觀測資料進行預處理,利用Tsoft(Vauterin,1998)重力潮汐數(shù)據(jù)預處理程序,去掉尖峰、突跳、間斷等.采用Eterna調和分析軟件,計算大氣導納因子,并在重力觀測中扣除大氣的影響(Wenzel,1996).

    r0處的點源地震激發(fā)臺站的簡正模位移可以利用雙力偶地震矩心矩張量理論和各項同性的地球模型模擬計算.在r處的位移u(r,t)為(Gilbert,1973):

    ×(M(r0):εn(r0)),

    (1)

    其中,Sn(r)為在某站點的標準化的n階模本征位移,ωn為本征頻率,Qn為質量因子,εn(r0)為應變張量,M(r0)為點源矩張量.二階6分量Mi(i=1,…,6),均勻地震矩張量M和εni組成εn(r0).則公式(1)可轉化為

    (2)

    (3)

    ai(r,t)為轉換函數(shù),Mi由標量地震矩M0、走向、傾角、滑動角組成(Aidetal.,1980).

    地球自由振蕩計算值的振幅主要取決于地震發(fā)生的位置和地震的震源機制解.模擬計算臺站自由振蕩需要的參數(shù)有:地震發(fā)生的時間、位置、走向、傾角、滑動角、標量地震矩、臺站的位置等.在參考的震源機制解中,地震的震源機制模型將震源等效為雙力偶點源,根據(jù)雙力偶地震矩心矩張量理論,由地球格林函數(shù)對矩張量加權求和,可計算矩張量源激發(fā)產生的地球自由振蕩位移(江穎等,2014b,2015).

    3結果分析

    3.1觀測值與模擬值比較

    將地震激發(fā)的地球自由振蕩信號的模擬值與相應的觀測值相比對,可評估震源機制解是否能正確反映地震的實際破裂情況.由于不同的地球模型對自由振蕩振幅的影響小于0.7%(薛秀秀等,2012),幾乎可以忽略.基于PREM地球模型(Dziewonski et al.,1981),結合震源機制解,將格林函數(shù)對矩張量加權求和,可以得到接收臺站的合成地震圖(Geller and Stein,1979).

    本文選取分布廣泛的16個彈簧重力儀觀測數(shù)據(jù)進行分析,臺站包括:北京(BJ)、福州(FZ)、高臺(GT)、格爾木(GEM)、海拉爾(HLE)、拉薩(LS)、蘭州(LZ)、牡丹江(MDJ)、憑祥(PX)、瓊中(QZ)、獅泉河(SQH)、炭山(TS)、烏加河(WJH)、烏什(WS)、西安(XA)、漳州(ZZ),見圖1.為了避免環(huán)境噪聲在簡正模振幅上的影響,選取16個臺站采樣間隔為1 s,時間長度為16h的重力觀測數(shù)據(jù),觀測數(shù)據(jù)和模擬數(shù)據(jù)在經過漢寧窗處理后進行離散傅里葉變換獲取線性振幅譜.連續(xù)觀測的彈簧重力儀會存在較大的零漂,選取的數(shù)據(jù)長度為16 h,零漂的影響較小,在數(shù)據(jù)處理過程中已經將該影響去掉.圖2顯示了16個重力臺站重力儀的自由振蕩觀測值和基于6個震源機制解計算的模擬值在1.5~5.3 mHz頻段的結果.在選取的頻段內,彈簧重力儀測得的自由振蕩信號有較高的信噪比.從圖2可以明顯看出基于震源機制解6的模擬值明顯偏小,而基于GCMT和USGS的震源機制解模擬計算的自由振蕩信號和觀測值符合的較好.

    3.2震源機制解整體評估

    為了更詳細的反映比較結果,基于16個重力儀臺站選用1.5~5.3 mHz頻段中所有的球型自由振蕩的模擬值和觀測值進行比較,利用殘余方差F(不符合度)和比例因子S(符合度)作為判斷標準:

    (4)

    (5)

    這里Oi和Si是觀測模的振幅和模擬值的振幅,n是所有振型的數(shù)量.當比例因子接近1的時候,其 震源機制解對應的M0是可靠的,可以描述震源的大小和強度.

    圖2 在1.5~5.3 mHz頻段內16個彈簧重力儀觀測數(shù)據(jù)與基于6個震源機制解算的模擬值

    另外,計算了基于各震源機制解結果的F和S的平均值及標準偏差,如表2.由表2可知,基于震源機制解3的自由振蕩模擬結果與實測結果最接近,其不符合度F的平均值為0.090,比例因子S為0.992,最接近1,F(xiàn)和S的標準偏差也最小,分別為0.055和0.065,最能反映地震釋放能量的大小.基于震源機制解1和2的模擬值與觀測值的殘余方差也較小,分別為0.098和0.091,顯示了較高的符合度.但是,基于震源機制解6的模擬值與觀測值相差較大,其不符合度F的平均值為0.214,比例因子S達到1.610.震源機制解中各個參數(shù)對地球自由振蕩信號振幅的影響需要進行進一步的計算與分析.

    表2 基于6個震源機制解的F與S的平均值和標準差

    3.3震源機制解影響定量分析

    分別計算震源機制解(標量地震矩、走向、傾角、滑動角和深度)對自由振蕩振幅的影響.基于USGSWPHASE的震源機制解計算的模擬值與觀測值最符合,由圖3的結果發(fā)現(xiàn),16個臺站中格爾木臺站觀測到的自由振蕩信號較接近模擬值,其對應的F為0.047,S為0.996.因此,基于格爾木臺站接收到的1.5~5.3 mHz頻段的自由振蕩信號平均振幅做定量分析.

    基于USGSWPHASE的震源機制解結果,分別給定標量地震矩、走向、傾角、滑動角和深度一個變化范圍,計算不同的震源機制解因素對自由振蕩振幅的影響.根據(jù)表1,標量地震矩的變化范圍為1.22×1021~2.35×1021N·m,保證走向、傾角、滑動角和深度不變,分別計算不同的標量地震矩對應的自由振蕩信號平均振幅,結果如圖4a,隨著標量地震矩的不斷增加,自由振蕩的平均振幅也在明顯增大,增量達到3.33×103nms-2/Hz;走向的變化范圍為240°~358°,保證M0、傾角、滑動角和深度不變,分別計算不同的走向對應的自由振蕩平均振幅,結果如圖4b,隨著走向增大,自由振蕩平均振幅也在慢慢增大,增加了2.25×103nms-2/Hz;傾角的變化范圍為12°~18°,保證M0、走向、滑動角和深度不變,分別計算不同的傾角對應的自由振蕩平均振幅,結果如圖4c,隨著傾角的不斷增加,振幅在慢慢減小,減小了3.08×103nms-2/Hz;滑動角變化范圍為74°~109°,保證M0、走向、傾角和深度不變,分別計算不同的滑動角對自由振蕩平均振幅的影響,結果如圖4d,隨著滑動角的不斷增加,振幅也在慢慢增大,增量達到2.13×103nms-2/Hz;深度變化范圍為20~33 km,保證M0、走向、傾角和滑動角不變,分別計算不同的深度對應的內核平動振蕩振幅的最大值,結果如圖4e,在20~33 km的范圍內,振幅有小范圍波動,差異僅為0.51×103nms-2/Hz.通過以上數(shù)值實驗說明M0對自由振蕩振幅的影響最大,而走向、傾角、滑動角和深度對內核平動振蕩振幅也有一定影響,但是影響相對較小.因此,將6個震源機制解中的M0均設置為相同的數(shù)值2.35×1021N·m,模擬計算格爾木臺站接收到的自由振蕩 信號,結果如圖5,當M0為2.35×1021N·m時,基于6個震源機制解得到的自由振蕩信號相差不大,且和觀測值均符合的較好.進一步證明了,M0是影響自由振蕩振幅的最大因素.An 等(2014)僅利用海嘯數(shù)據(jù)反演得到的震源機制解中M0偏小是造成模擬值和觀測值偏差較大的主要原因.

    圖4 震源機制解對自由振蕩振幅的影響

    圖5 格爾木臺重力觀測數(shù)據(jù)與基于改進的6個震源機制解解算的模擬值(其中,M0均為2.35×1021 N·m)

    4結論與討論

    通過比較分析Iquique地震激發(fā)的自由振蕩信號的實測值與模擬值,結果表明An 等(2014)等僅利用海嘯數(shù)據(jù)反演得到的震源機制解中M0偏小是造成模擬值和觀測值偏差較大的主要原因,說明對于本次的Iquique地震僅用海嘯數(shù)據(jù)做反演,有一定的局限性,這是由于海嘯數(shù)據(jù)對淺層地震的分辨率較高,而對于相對較深的逆沖型地震做數(shù)據(jù)分析的時候,并不一定能得到最佳解.利用海嘯數(shù)據(jù)反演震源機制主要是快速的得到破裂面如何破裂,在哪個方向上會有更大的破壞,而對于標量地震矩的估算會有一定偏差.另外,Lay 等(2014)從112個全球地震臺站選取200~1000 s的周期范圍的160個地面運動記錄的W-phase數(shù)據(jù)反演了Iquique地震,得到的標量地震矩為1.99×1021N·m,而另外的一組結合海嘯數(shù)據(jù)得到的震源機制解M0的為1.7×1021N·m,海嘯數(shù)據(jù)的加入造成了反演的標量地震矩結果偏小.Yagi等(2014)僅采用了遠震P波數(shù)據(jù),反演結果同樣偏小,可能是由于遠場體波記錄不到一些緩慢破裂的過程,且遠場體波容易受到其他震相的污染,特別是在主震后發(fā)生大量余震的情況下.對于大地震,USGS的地震檔案中一般不提供利用體波反演的震源機制解.USGSWPHASE選取長周期(周期為100~1000 s)方位角為5°~20°的近場數(shù)據(jù)和30°~90°的遠場數(shù)據(jù)聯(lián)合反演得到的震源機制解最佳.因此,利用長周期面波信號可以較好的約束Iquique地震的標量地震矩.

    本文基于地球自由振蕩分析了Iquique地震的震源機制解.結果表明利用彈簧重力儀記錄的1.5~5.3 mHz球型簡正模觀測可以有效地對Iquique地震的震源機制解進行總體評估.從研究結果來看,地震的斷層走向、傾角、滑動方向角和震源深度對自由振蕩的振幅有一定影響,但是影響較小,M0對自由振蕩振幅的影響最大.另外,基于USGSWPHASE反演得到的震源機制解的自由振蕩模擬值與觀測值符合最好,其反演的標量地震矩可以較好地反映Iquique地震釋放的能量,而僅利用海嘯數(shù)據(jù)反演的標量地震矩偏小,聯(lián)合近遠場觀測數(shù)據(jù)反演可顯著改善震源機制解.因此,利用自由振蕩的方法可以較好地判斷地震震源機制和有效約束地震的震級.

    另外,地球自由振蕩觀測只能對震源機制解進行總體綜合評估,而不能單獨對各個參數(shù)的影響單獨評估.如果能利用其他大地測量手段,如近場GPS觀測約束震源傾角、滑動方向角,結合地球自由振蕩的觀測和模擬可更細致地評估震源機制解.

    致謝感謝國家重力臺網中心提供的全國廣泛分布的連續(xù)重力觀測數(shù)據(jù).

    References

    AidK,Richards P G.1980.Quantitative Seismology : Theory and Methods.San Francisco: University Books.

    An C,Sepúlveda I,Liu P L F.2014.Tsunami source and its validation of the 2014 Iquique,Chile,earthquake.Geophysical Research Letters,41(11): 3988-3994.

    Angermann D,Klotz J,Reigber C.1999.Space-geodetic estimation of the Nazca-South America Euler vector.Earth and Planetary Science Letters,171(3): 329-334.

    Bürgmann R.2014.Earth science: Warning signs of the Iquique earthquake.Nature,512(7514): 258-259.

    Dziewonski A M,Anderson D L.1981.Preliminary reference Earth model.Phys.Earth Planet.Inter.,25(4): 297-356.

    Geller R J,Stein S.1979.Time-domain attenuation measurements for fundamental spheroidal modes (0S6to0S28) for the 1977 Indonesian earthquake.Bull.Seismol.Soc.Amer.,69(6):1671-1691.

    Gilbert F.1973.Derivation of source parameters from low-frequency spectra.Philosophical Transactions of the Royal Society A: Mathematical,Physical and Engineering Sciences,274(1239): 369-371.

    Hayes G P,Herman M W,Barnhart W D,et al.2014.Continuing megathrust earthquake potential in Chile after the 2014 Iquique earthquake.Nature,512(7514): 295-298.

    Hu X G,Jiang Y.2013.Evaluation of the seismic moment of the April 20,2013 Lushan earthquake.Earthq.Sci.,26(3-4): 169-177.

    Hu X G,Jiang Y,Sun H P.2014.Assessing the scalar moment of moderate earthquake and the effect of lateral heterogeneity on normal modes—An example from the 2013/04/20 Lushan earthquake,Sichuan,China.Physics of the Earth and Planetary Interiors,232: 61-71.

    Jiang Y,Hu X G,Liu C L,et al.2014.Constraining the focal mechanism of the Lushan earthquake with observations of the Earth′s free oscillations.Science China Earth Sciences,57(9): 2064-2070.

    Kato A,Nakagawa S.2014.Multiple slow-slip events during a foreshock sequence of the 2014 Iquique,Chile MW8.1 earthquake.Geophysical Research Letters,41(15): 5420-5427.

    Lambotte S,Rivera L,Hinderer J.2006.Rupture length and duration of the 2004 Aceh-Sumatra earthquake from the phases of the Earth′s gravest free oscillations.Geophysical Research Letters,33(3): L03307.

    Lay T,Yue H,Brodsky E E,et al.2014.The 1 April 2014 Iquique,Chile,MW8.1 earthquake rupture sequence.Geophysical Research Letters,41(11): 3818-3825.

    Liu C L,Zheng Y,Wang R J,et al.2015.Kinematic rupture process of the 2014 Chile MW8.1 earthquake constrained by strong-motion,GPS static offsets and teleseismic data.Geophysical Journal International,2: 1137-1145.

    Park J,Song T-R A,Tromp J,et al.2005.Earth′s free oscillations excited by the 26 December 2004 Sumatra-Andaman earthquake.Science,308(5725): 1139-1144.Ruiz S,Metois M,Fuenzalida A,et al.2014.Intense foreshocks and a slow slip event preceded the 2014 Iquique MW8.1 earthquake.Science,345(6201): 1165-1169.

    Schurr B,Asch G,Hainzl S,et al.2014.Gradual unlocking of plate boundary controlled initiation of the 2014 Iquique earthquake.Nature,512(7514): 299-302.

    Stein S,Okal E A.2005.Seismology: Speed and size of the Sumatra earthquake.Nature,308:1127-1139.

    Vauterin P.1998.Tsoft: Graphical and interactive software for the analysis of Earth tide data.∥ Proc.13th Int.Sympos.on Earth Tides,Brussels,Observatoire Royal de Belgique,Série Géophysique.481-486.

    Wenzel H G.1996.The nanogal software: Data processing package ETERNA 3.3.Bull.Inf.Marées Terrestres,124: 9425-9439.

    Yagi Y,Okuwaki R,Enescu B,et al.2014.Rupture process of the 2014 Iquique Chile Earthquake in relation with the foreshock activity.Geophysical Research Letters,41(12): 4201-4206.

    附中文參考文獻

    江穎,胡小剛,劉成利等.2014a.利用地球自由振蕩觀測約束蘆山地震的震源機制解.中國科學: 地球科學,44(12): 2689-2696.

    江穎,徐建橋,孫和平.2014b.深內部地球結構對內核平動振蕩本征周期的影響.地球物理學報,57(4): 1041-1048,doi: 10.6038/cjg20140403.

    江穎,徐建橋,孫和平.2015.地球內核平動振蕩的地震激發(fā).地球物理學報,58(1): 92-102,doi: 10.6038/cjg20150108.

    薛秀秀,胡小剛,郝曉光等.2012.利用重力觀測約束2011日本Tohoku大地震的震源機制.地球物理學報,55(9): 3006-3015,doi: 10.6038/j.issn.0001-5733.2012.09.019.

    (本文編輯汪海英)

    基金項目科研院所基本科研業(yè)務費(中國地震局地震研究所所長基金)(IS201526225),大地測量與地球動力學國家重點實驗室開放基金(SKLGED2015-1-2-EZ),國家自然科學基金(41404064,41374088)資助.

    作者簡介江穎,女,1986年生,現(xiàn)主要從事地球自由振蕩、重力數(shù)據(jù)分析及其應用研究.E-mail:jiangyingchen@126.com

    doi:10.6038/cjg20160316 中圖分類號P315

    收稿日期2015-04-08,2016-01-11收修定稿

    Constraining the focal mechanism of the Iquique earthquake with observations of the continuous gravity stations

    JIANG Ying1,2,LI Hui1,LIU Zi-Wei1,SHEN Chong-Yang1,ZHANG Xiao-Tong1

    1KeyLaboratoryofEarthquakeGeodesy,InstituteofSeismology,ChinaEarthquakeAdministration,Wuhan430071,China2StateKeyLaboratoryofGeodesyandEarth′sDynamics,InstituteofGeodesyandGeophysics,CAS,Wuhan430077,China

    AbstractOn 1 April 2014,a magnitude MW8.1 interplate thrust earthquake ruptured a densely instrumented region of Iquique seismic gap in northern Chile.After Earthquake,the focal mechanism solutions of Iquique earthquake were provided by well-respected international earthquake research institutions and scholars based on different data and methods,which were quite different.The amplitudes of the Earth′s free oscillations have a close relationship to earthquake focal mechanisms.Focal mechanisms of large earthquakes can be well analyzed and constrained with observations of long period free oscillations.We compared free oscillations observed by 16 spring gravimeters of continuous gravity stations with synthetic normal modes corresponding to six different focal mechanisms for the Iquique earthquake,and the focal mechanisms solutions of Iquique earthquake were analyzed and constrained by spherical normal modes in a 1.5 to 5.3 mHz frequency band.It′s found that the synthetic modes corresponding to the focal mechanism determined by the United States Geological Survey WPhase Moment Tensor Solution show agreement to the observed modes,suggesting that earthquake magnitudes predicted in this way can reflect the total energy released by the earthquake.The scalar seismic moment obtained by tsunami data inversion is significantly underestimated.Focal mechanism solutions can be improved by joint inversion of far- and near-field data.In addition,based on the synthetic and observed normal modes at Golmud station,quantitative analysis was carried out to research the influences of free oscillation′s amplitudes caused by different focal mechanism parameters.The results show that scalar seismic moment M0 has a major influence on the amplitudes of free oscillations,and the strike,dip,rake and depth of the hypocenter have minor influences.

    KeywordsIquique earthquake; Focal mechanism solutions; Earth′s free oscillations; Gravimeter observations

    江穎,李輝,劉子維等.2016.利用連續(xù)重力觀測約束2014智利Iquique地震的震源機制解.地球物理學報,59(3):941-951,doi:10.6038/cjg20160316.

    Jiang Y,Li H,Liu Z W,et al.2016.Constraining the focal mechanism of the Iquique earthquake with observations of the continuous gravity stations.Chinese J.Geophys.(in Chinese),59(3):941-951,doi:10.6038/cjg20160316.

    猜你喜歡
    標量臺站震源
    中國科學院野外臺站檔案工作回顧
    氣象基層臺站建設
    西藏科技(2021年12期)2022-01-17 08:46:38
    一種高效的橢圓曲線密碼標量乘算法及其實現(xiàn)
    一種靈活的橢圓曲線密碼并行化方法
    震源的高返利起步
    可控震源地震在張掖盆地南緣逆沖斷裂構造勘探中的應用
    華北地質(2015年3期)2015-12-04 06:13:25
    基層臺站綜合觀測業(yè)務管理之我見
    西藏科技(2015年6期)2015-09-26 12:12:13
    同步可控震源地震采集技術新進展
    單調Minkowski泛函與Henig真有效性的標量化
    MDOS平臺臺站級使用方法及技巧
    欧美精品啪啪一区二区三区| 国产精品偷伦视频观看了| av中文乱码字幕在线| 亚洲精品久久午夜乱码| 黄片大片在线免费观看| 80岁老熟妇乱子伦牲交| 热99国产精品久久久久久7| 中出人妻视频一区二区| 男女免费视频国产| 精品亚洲成国产av| 少妇猛男粗大的猛烈进出视频| 后天国语完整版免费观看| 日韩欧美三级三区| 欧美成人午夜精品| 精品第一国产精品| 中文欧美无线码| 十八禁网站免费在线| 身体一侧抽搐| 免费观看精品视频网站| 欧美色视频一区免费| a级片在线免费高清观看视频| 一二三四在线观看免费中文在| 欧美另类亚洲清纯唯美| 最近最新免费中文字幕在线| 在线永久观看黄色视频| 亚洲三区欧美一区| 精品少妇久久久久久888优播| 欧美成狂野欧美在线观看| 18禁国产床啪视频网站| 亚洲精品在线美女| 久久青草综合色| 99久久国产精品久久久| 久久人妻福利社区极品人妻图片| 亚洲全国av大片| 色播在线永久视频| 欧美黄色淫秽网站| 男女之事视频高清在线观看| 久久精品aⅴ一区二区三区四区| 97人妻天天添夜夜摸| 嫩草影视91久久| 久久人妻av系列| 亚洲色图 男人天堂 中文字幕| 国产午夜精品久久久久久| 韩国精品一区二区三区| 亚洲av片天天在线观看| 老司机福利观看| 18在线观看网站| 欧美日本中文国产一区发布| 99久久人妻综合| 免费日韩欧美在线观看| 国产男女内射视频| 桃红色精品国产亚洲av| 女人久久www免费人成看片| 免费看a级黄色片| 精品福利观看| 91老司机精品| 日本wwww免费看| 最近最新中文字幕大全免费视频| 热99久久久久精品小说推荐| 国产精品综合久久久久久久免费 | 真人做人爱边吃奶动态| 免费不卡黄色视频| 亚洲黑人精品在线| 亚洲av美国av| 国产野战对白在线观看| 精品福利观看| 一级毛片精品| 亚洲综合色网址| 青草久久国产| 久久国产精品大桥未久av| 在线观看免费视频网站a站| 国产无遮挡羞羞视频在线观看| 丝袜美足系列| 色尼玛亚洲综合影院| 如日韩欧美国产精品一区二区三区| 国产免费现黄频在线看| 精品视频人人做人人爽| 飞空精品影院首页| 欧美精品啪啪一区二区三区| 看免费av毛片| 国产成人欧美在线观看 | 成人影院久久| 黄色片一级片一级黄色片| 午夜福利在线观看吧| 国产成人精品在线电影| 黄色女人牲交| 免费看a级黄色片| 国产一区二区三区视频了| 成年人午夜在线观看视频| 日日夜夜操网爽| www.精华液| 免费在线观看亚洲国产| 18在线观看网站| 亚洲av第一区精品v没综合| 欧美激情极品国产一区二区三区| 久久狼人影院| 亚洲精品国产区一区二| 久久久久国产精品人妻aⅴ院 | 一进一出好大好爽视频| 久久国产乱子伦精品免费另类| 国产一区二区三区视频了| 国产乱人伦免费视频| 99久久99久久久精品蜜桃| 久久午夜综合久久蜜桃| 免费在线观看黄色视频的| 在线天堂中文资源库| 黑人巨大精品欧美一区二区mp4| 国产精品久久电影中文字幕 | 曰老女人黄片| 国产97色在线日韩免费| av视频免费观看在线观看| 亚洲午夜精品一区,二区,三区| 狂野欧美激情性xxxx| videosex国产| 亚洲av成人不卡在线观看播放网| 精品一区二区三区av网在线观看| 99香蕉大伊视频| 日本黄色日本黄色录像| 亚洲伊人色综图| 伦理电影免费视频| 久久香蕉国产精品| 国产亚洲av高清不卡| 中文字幕精品免费在线观看视频| 欧美日韩亚洲国产一区二区在线观看 | 国产成人精品无人区| 一级a爱片免费观看的视频| 久久热在线av| 国产av又大| 久久亚洲精品不卡| 午夜福利影视在线免费观看| 99热国产这里只有精品6| 精品国产乱码久久久久久男人| 午夜视频精品福利| 色播在线永久视频| 国产区一区二久久| 亚洲一码二码三码区别大吗| 国产91精品成人一区二区三区| 手机成人av网站| 最新美女视频免费是黄的| 久久亚洲真实| 国产免费av片在线观看野外av| 欧美日韩av久久| 成人18禁在线播放| 欧美精品一区二区免费开放| 久热这里只有精品99| 一级a爱视频在线免费观看| 亚洲第一欧美日韩一区二区三区| 亚洲精品中文字幕一二三四区| av欧美777| 国产无遮挡羞羞视频在线观看| 精品国产乱子伦一区二区三区| 免费不卡黄色视频| 国产精品 国内视频| 99精品欧美一区二区三区四区| 久久久久久人人人人人| av电影中文网址| 国产极品粉嫩免费观看在线| 精品国产一区二区三区久久久樱花| 欧美av亚洲av综合av国产av| 国产精品99久久99久久久不卡| 91老司机精品| 国产单亲对白刺激| 亚洲五月天丁香| 国产欧美亚洲国产| 村上凉子中文字幕在线| 免费日韩欧美在线观看| 久久国产精品人妻蜜桃| 极品人妻少妇av视频| 中亚洲国语对白在线视频| 美女高潮到喷水免费观看| 伊人久久大香线蕉亚洲五| 午夜福利影视在线免费观看| 久久影院123| 精品久久久久久电影网| 最近最新中文字幕大全免费视频| 国产精品免费一区二区三区在线 | 成人手机av| 俄罗斯特黄特色一大片| 热99久久久久精品小说推荐| 飞空精品影院首页| 久久精品亚洲精品国产色婷小说| 两个人看的免费小视频| 国产极品粉嫩免费观看在线| 69av精品久久久久久| 亚洲国产看品久久| 动漫黄色视频在线观看| 免费在线观看视频国产中文字幕亚洲| 丰满人妻熟妇乱又伦精品不卡| 日本wwww免费看| 超色免费av| 亚洲av电影在线进入| 18禁国产床啪视频网站| 国产野战对白在线观看| 妹子高潮喷水视频| 成人永久免费在线观看视频| 欧美大码av| 成年女人毛片免费观看观看9 | 捣出白浆h1v1| 悠悠久久av| 国产精品乱码一区二三区的特点 | 日韩人妻精品一区2区三区| 日韩精品免费视频一区二区三区| 妹子高潮喷水视频| 午夜福利影视在线免费观看| 精品第一国产精品| x7x7x7水蜜桃| 亚洲人成电影观看| 欧美日韩亚洲综合一区二区三区_| 国产一卡二卡三卡精品| 999精品在线视频| 国产精品一区二区精品视频观看| xxx96com| 亚洲一区高清亚洲精品| 国产精品偷伦视频观看了| 91av网站免费观看| 人妻一区二区av| 国产又爽黄色视频| 黄网站色视频无遮挡免费观看| 免费在线观看日本一区| 国产激情久久老熟女| 亚洲人成电影免费在线| videosex国产| 亚洲成av片中文字幕在线观看| av线在线观看网站| 亚洲成人免费电影在线观看| 亚洲欧美一区二区三区久久| 亚洲欧美日韩高清在线视频| 成年女人毛片免费观看观看9 | 色94色欧美一区二区| 他把我摸到了高潮在线观看| 男女午夜视频在线观看| 黄色成人免费大全| 欧美性长视频在线观看| 中文字幕av电影在线播放| 国产aⅴ精品一区二区三区波| 777米奇影视久久| 露出奶头的视频| 老汉色av国产亚洲站长工具| av视频免费观看在线观看| 丁香欧美五月| 18禁观看日本| 一区二区三区激情视频| 美女国产高潮福利片在线看| 色综合婷婷激情| 激情在线观看视频在线高清 | 动漫黄色视频在线观看| 岛国在线观看网站| 日日夜夜操网爽| 久久国产精品人妻蜜桃| 国产亚洲欧美精品永久| 国产区一区二久久| 成年人黄色毛片网站| 69av精品久久久久久| 丝袜美腿诱惑在线| 人妻 亚洲 视频| 午夜福利视频在线观看免费| 国产欧美日韩一区二区三区在线| 国产亚洲av高清不卡| 韩国av一区二区三区四区| 精品久久久久久久久久免费视频 | 女警被强在线播放| 亚洲av熟女| 两个人免费观看高清视频| 精品人妻1区二区| 不卡一级毛片| 亚洲欧美精品综合一区二区三区| 国产亚洲精品一区二区www | 午夜精品久久久久久毛片777| 久久久久久久久久久久大奶| 中文字幕人妻丝袜制服| 欧美日韩中文字幕国产精品一区二区三区 | 免费在线观看日本一区| 国产日韩一区二区三区精品不卡| 国内久久婷婷六月综合欲色啪| 久久草成人影院| 日本一区二区免费在线视频| 大香蕉久久成人网| 男女之事视频高清在线观看| 国产精品二区激情视频| 国产精品久久久人人做人人爽| 日韩 欧美 亚洲 中文字幕| 免费黄频网站在线观看国产| 婷婷精品国产亚洲av在线 | 亚洲欧美一区二区三区久久| 国内久久婷婷六月综合欲色啪| 女人被狂操c到高潮| av福利片在线| 欧美日韩成人在线一区二区| 免费看十八禁软件| 欧美日韩乱码在线| 五月开心婷婷网| 久久中文看片网| 人妻一区二区av| 日韩欧美三级三区| 在线永久观看黄色视频| 他把我摸到了高潮在线观看| 在线永久观看黄色视频| 日韩欧美一区二区三区在线观看 | 在线国产一区二区在线| 精品国产乱子伦一区二区三区| 搡老乐熟女国产| 成人18禁高潮啪啪吃奶动态图| 国产精品 欧美亚洲| 午夜久久久在线观看| 在线观看免费日韩欧美大片| 一个人免费在线观看的高清视频| 成年人免费黄色播放视频| 99精品欧美一区二区三区四区| 久久精品国产亚洲av香蕉五月 | 精品国产美女av久久久久小说| 纯流量卡能插随身wifi吗| 伊人久久大香线蕉亚洲五| 午夜免费成人在线视频| 色在线成人网| 国产又色又爽无遮挡免费看| 免费在线观看日本一区| 黄片播放在线免费| 国产精品 欧美亚洲| 久久午夜综合久久蜜桃| 国产成人av激情在线播放| 国产激情欧美一区二区| 精品国内亚洲2022精品成人 | 如日韩欧美国产精品一区二区三区| 国产色视频综合| 香蕉国产在线看| 久久香蕉激情| 亚洲国产精品一区二区三区在线| 大型av网站在线播放| 亚洲五月婷婷丁香| 精品熟女少妇八av免费久了| 波多野结衣一区麻豆| 成人手机av| 一边摸一边抽搐一进一小说 | 国产成人精品无人区| 一个人免费在线观看的高清视频| 国产伦人伦偷精品视频| 亚洲一码二码三码区别大吗| 人成视频在线观看免费观看| 夜夜爽天天搞| 嫁个100分男人电影在线观看| 后天国语完整版免费观看| 别揉我奶头~嗯~啊~动态视频| 很黄的视频免费| 中文字幕色久视频| 国产精品久久电影中文字幕 | 在线观看免费午夜福利视频| av片东京热男人的天堂| 欧美中文综合在线视频| 午夜成年电影在线免费观看| 久久午夜综合久久蜜桃| 亚洲精品av麻豆狂野| 美国免费a级毛片| 亚洲自偷自拍图片 自拍| 咕卡用的链子| 女人爽到高潮嗷嗷叫在线视频| 激情在线观看视频在线高清 | 久久精品亚洲精品国产色婷小说| 一本大道久久a久久精品| 好看av亚洲va欧美ⅴa在| 欧美精品亚洲一区二区| 妹子高潮喷水视频| 精品国产超薄肉色丝袜足j| 亚洲精品国产区一区二| 免费看a级黄色片| 国产av一区二区精品久久| 91麻豆精品激情在线观看国产 | 久久久精品国产亚洲av高清涩受| 国产精品秋霞免费鲁丝片| 热99re8久久精品国产| 69精品国产乱码久久久| 国产精华一区二区三区| 国产精品美女特级片免费视频播放器 | 午夜福利影视在线免费观看| 下体分泌物呈黄色| 亚洲精品中文字幕在线视频| 色在线成人网| 丰满的人妻完整版| av欧美777| 99精品欧美一区二区三区四区| 91成年电影在线观看| 少妇裸体淫交视频免费看高清 | 欧美黑人精品巨大| 80岁老熟妇乱子伦牲交| 日韩三级视频一区二区三区| 少妇被粗大的猛进出69影院| 99国产精品免费福利视频| 精品久久久久久久毛片微露脸| a级毛片在线看网站| av福利片在线| 欧美成人免费av一区二区三区 | 欧美乱妇无乱码| 操美女的视频在线观看| 久久这里只有精品19| 成人三级做爰电影| 夜夜躁狠狠躁天天躁| 国产日韩一区二区三区精品不卡| 色播在线永久视频| 国产精品1区2区在线观看. | 欧美乱码精品一区二区三区| 久久久久国产精品人妻aⅴ院 | 亚洲三区欧美一区| 韩国精品一区二区三区| 涩涩av久久男人的天堂| 欧美午夜高清在线| 久久精品国产清高在天天线| 午夜福利欧美成人| 91精品三级在线观看| 国产成人精品无人区| 欧美国产精品一级二级三级| 麻豆乱淫一区二区| 日韩欧美三级三区| 69精品国产乱码久久久| 热re99久久国产66热| av中文乱码字幕在线| а√天堂www在线а√下载 | 动漫黄色视频在线观看| 丰满的人妻完整版| 俄罗斯特黄特色一大片| 国产av又大| 狂野欧美激情性xxxx| xxxhd国产人妻xxx| 国产激情欧美一区二区| 国产在视频线精品| 亚洲欧美日韩另类电影网站| 18禁观看日本| 亚洲伊人色综图| 黄片小视频在线播放| 亚洲熟妇熟女久久| 国产成人精品无人区| 黄色片一级片一级黄色片| 黄色片一级片一级黄色片| 欧美中文综合在线视频| 看黄色毛片网站| 侵犯人妻中文字幕一二三四区| 悠悠久久av| 日韩欧美一区二区三区在线观看 | 亚洲av成人不卡在线观看播放网| 一区二区日韩欧美中文字幕| 欧美老熟妇乱子伦牲交| av视频免费观看在线观看| 久久中文字幕人妻熟女| 最新的欧美精品一区二区| e午夜精品久久久久久久| 国产99久久九九免费精品| 亚洲午夜精品一区,二区,三区| 自线自在国产av| 1024香蕉在线观看| 少妇裸体淫交视频免费看高清 | 亚洲国产精品sss在线观看 | 精品熟女少妇八av免费久了| 十分钟在线观看高清视频www| 免费av中文字幕在线| 亚洲av成人一区二区三| 999精品在线视频| 亚洲五月婷婷丁香| 99riav亚洲国产免费| 亚洲免费av在线视频| 最近最新免费中文字幕在线| 中文字幕人妻熟女乱码| 国产人伦9x9x在线观看| 午夜福利影视在线免费观看| 国产无遮挡羞羞视频在线观看| 免费黄频网站在线观看国产| 久久久久精品人妻al黑| 久久亚洲真实| 免费不卡黄色视频| 人人妻人人添人人爽欧美一区卜| 美女国产高潮福利片在线看| 成人国产一区最新在线观看| 1024香蕉在线观看| 男女免费视频国产| 一进一出抽搐动态| 亚洲成a人片在线一区二区| 欧美丝袜亚洲另类 | 午夜精品国产一区二区电影| 日韩 欧美 亚洲 中文字幕| 脱女人内裤的视频| 欧美亚洲日本最大视频资源| 免费观看人在逋| 黄片播放在线免费| 日韩成人在线观看一区二区三区| 久久精品91无色码中文字幕| 99国产精品99久久久久| 久久亚洲精品不卡| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲专区字幕在线| 人人妻人人爽人人添夜夜欢视频| 国产亚洲欧美98| 日韩中文字幕欧美一区二区| 久久中文字幕人妻熟女| 制服人妻中文乱码| 婷婷精品国产亚洲av在线 | 国产精品电影一区二区三区 | 日本撒尿小便嘘嘘汇集6| 水蜜桃什么品种好| 国产精品一区二区在线观看99| 免费av中文字幕在线| 乱人伦中国视频| 在线观看午夜福利视频| 91九色精品人成在线观看| 亚洲,欧美精品.| 亚洲成国产人片在线观看| 亚洲精品国产区一区二| 最近最新中文字幕大全免费视频| 欧美最黄视频在线播放免费 | 超碰成人久久| 成年人免费黄色播放视频| www.熟女人妻精品国产| 涩涩av久久男人的天堂| 如日韩欧美国产精品一区二区三区| 99国产精品一区二区三区| 成年人免费黄色播放视频| 国产激情久久老熟女| 久久国产精品男人的天堂亚洲| 老司机福利观看| 欧美黑人精品巨大| 国产高清国产精品国产三级| 99国产精品免费福利视频| 一级毛片精品| 黄色视频,在线免费观看| 91在线观看av| 欧美激情高清一区二区三区| 国产精品乱码一区二三区的特点 | 99久久人妻综合| 女警被强在线播放| 一进一出抽搐gif免费好疼 | 亚洲色图 男人天堂 中文字幕| 亚洲三区欧美一区| 欧美精品一区二区免费开放| 女人被躁到高潮嗷嗷叫费观| 亚洲国产欧美一区二区综合| 超碰97精品在线观看| 很黄的视频免费| 巨乳人妻的诱惑在线观看| 欧美大码av| 亚洲熟妇熟女久久| 9191精品国产免费久久| 狠狠狠狠99中文字幕| 老熟妇乱子伦视频在线观看| av网站免费在线观看视频| 久热爱精品视频在线9| 国产成人影院久久av| 热re99久久精品国产66热6| 十八禁人妻一区二区| 亚洲精品一卡2卡三卡4卡5卡| 1024视频免费在线观看| 久久香蕉精品热| 久久午夜亚洲精品久久| 国产亚洲一区二区精品| 母亲3免费完整高清在线观看| 亚洲专区字幕在线| 露出奶头的视频| 国精品久久久久久国模美| 在线永久观看黄色视频| 成熟少妇高潮喷水视频| 久久国产乱子伦精品免费另类| 午夜老司机福利片| 午夜视频精品福利| 久久九九热精品免费| 侵犯人妻中文字幕一二三四区| 久久99一区二区三区| 中文字幕精品免费在线观看视频| 十八禁网站免费在线| 丰满迷人的少妇在线观看| 国产免费av片在线观看野外av| av网站在线播放免费| 亚洲 国产 在线| 免费在线观看黄色视频的| 久久婷婷成人综合色麻豆| 男人的好看免费观看在线视频 | 一本一本久久a久久精品综合妖精| 国产精品欧美亚洲77777| 免费日韩欧美在线观看| 99精国产麻豆久久婷婷| 色综合婷婷激情| 久久香蕉国产精品| 国产99久久九九免费精品| 亚洲久久久国产精品| 黑人欧美特级aaaaaa片| 1024香蕉在线观看| 国产片内射在线| 黄色女人牲交| 黑人猛操日本美女一级片| 久热这里只有精品99| 校园春色视频在线观看| 国产真人三级小视频在线观看| 国产单亲对白刺激| 亚洲国产精品sss在线观看 | 露出奶头的视频| 午夜日韩欧美国产| 一个人免费在线观看的高清视频| 看片在线看免费视频| 19禁男女啪啪无遮挡网站| 亚洲色图综合在线观看| 国产欧美亚洲国产| 免费在线观看日本一区| 精品人妻熟女毛片av久久网站| 久久久久久久久免费视频了| 深夜精品福利| 高潮久久久久久久久久久不卡| 亚洲久久久国产精品| xxxhd国产人妻xxx| 19禁男女啪啪无遮挡网站| 少妇猛男粗大的猛烈进出视频| 午夜福利免费观看在线| 9191精品国产免费久久| 叶爱在线成人免费视频播放| 久久精品亚洲精品国产色婷小说| 狠狠婷婷综合久久久久久88av| 免费在线观看视频国产中文字幕亚洲| 黄频高清免费视频| 色老头精品视频在线观看| av网站免费在线观看视频| 少妇 在线观看| 黄片大片在线免费观看|