史延飛 辛存林 梁浩東 劉海博

摘要：地震前后氣體地球化學信息的研究在地震監(jiān)測預報中具有較高的應用價值。為研究2021年5月22日青海瑪多MS7.4地震與CO時空變化之間的關系，在以瓦里關大氣觀測站地面數(shù)據(jù)驗證大氣紅外探測儀（AIRS）反演數(shù)據(jù)可靠性的基礎上，提取瑪多MS7.4地震前后AIRS反演的CO數(shù)據(jù)，通過滑動均值法、差值法對瑪多地震前后不同尺度的CO濃度數(shù)據(jù)進行處理和分析。結果表明：利用衛(wèi)星遙感數(shù)據(jù)提取CO地球化學信息是可靠的。時間上，瑪多MS7.4地震前兩個月CO濃度開始波動，出現(xiàn)峰值，地震發(fā)生后恢復平靜;空間上，震中位置的CO濃度在近地面變化尤其明顯，震中及附近區(qū)域的CO濃度從3月開始逐漸升高，由離散分布逐漸向震中和發(fā)震斷裂帶聚攏靠近，到4月底達到最大異常18.60×10-9，異常高值中心的連線與發(fā)震斷層江錯斷裂走向、地表破裂分布一致。排除背景值和季節(jié)變化的影響，推斷CO濃度異常變化是地震引起的，主要歸因于地下氣體釋放和巖石擠壓碰撞產氣，氣體逸散后在大氣圈中發(fā)生的一系列化學反應起次要作用。

關鍵詞：青海瑪多; MS7.4地震; AIRS傳感器; CO柱濃度; CO體積混合比

中圖分類號： P315????? 文獻標志碼：A?? 文章編號： 1000-0844（2024）03-0692-12

DOI：10.20000/j.1000-0844.20230105002

CO anomalies before and after the 2021 Maduo MS7.4 earthquake in Qinghai Province

SHI Yanfei1， XIN Cunlin1， LIANG Haodong1， LIU Haibo2

（1. School of Geography and Environmental Science， Northwest Normal University， Lanzhou 730070， Gansu， China;2. School of Earth Sciences， Zhejiang University， Hangzhou 310027， Zhejiang， China）

Abstract：?The investigation of gas geochemical information before and after an earthquake has high application value in earthquake monitoring and prediction. Thus， this paper investigates the relationship between the spatio-temporal variation of CO and the MS7.4 earthquake that occurred in Maduo County， Qinghai Province， on May 22， 2021. The accuracy of CO data obtained from the Atmospheric Infrared Sounder （AIRS） was compared with ground data from the Waliguan atmospheric observation station. The CO data retrieved by AIRS before and after the MS7.4 earthquake were extracted， from which the CO concentrations before and after the earthquake were processed and analyzed using the sliding mean method and difference method. Results indicate the reliability of extracting CO geochemical information from satellite remote sensing data. The CO concentration began to fluctuate and peak two months before the earthquake and then recovered after the earthquake. Furthermore， the change in CO concentration in the epicenter and its vicinity was particularly obvious near the surface. This gradually increased from March， gathered near the epicenter and the seismogenic fault， and eventually reached the maximum value of 18.60×10-9 by the end of April. The line connected to the centers of the abnormally high values is consistent with the strike of the seismogenic fault and the surface rupture distribution. Barring the influence of background value and seasonal change， it is inferred that the abnormal change in CO concentration is caused by the MS7.4 earthquake. The CO anomaly can be mainly attributed to the underground gas release and gas production caused by rock extrusion and collision， while the gas chemical reactions in the atmosphere play a secondary role.

Keywords：Maduo in Qinghai Province; MS7.4 earthquake; ARIS; CO column concentration; CO volume mixing ratio

0 引言

地球是一個開放的不平衡體系，液態(tài)地核中流體濃度大，含有H2、HF、H2S、CH4、CO、CO2、He、Rn及Ar等氣體，4 500 Ma來從地核向外不斷地排放［1-2］。隨著地震活動強烈期構造應力增加，地球深部流體沿斷裂帶在巖石圈內遷移、聚集或分散，引起水位、水溫、水化學離子組分以及深部氣體發(fā)生顯著變化［3-4］。構造活動產生的裂縫為氣體向地表擴散提供了良好通道［5］，促使氣體釋放到大氣圈形成局地氣體濃度升高。釋放的氣體還會通過物理作用（沉降作用、對流傳輸作用）和化學反應影響區(qū)域大氣成分［6-7］。氣體異常變化是地震前兆的一種表現(xiàn)，斷層氣地球化學特征異常形態(tài)對斷層特性有較好的反映［8］，科學地監(jiān)測地震前后氣體變化，在地震監(jiān)測預報中具有重要的應用價值。衛(wèi)星高光譜技術具有覆蓋范圍廣、觀測周期短、受地表影響少的優(yōu)勢［9-10］，彌補了地面觀測的不足，利用遙感監(jiān)測地震氣體異常成為研究的熱點之一。

目前，已有研究者利用中分辨率成像光譜儀（Moderate Resolution Imaging Spectroradiometer，MODIS）、對流層污染測量儀（Measurement of Pollution in the Troposphere，MOPITT）及大氣紅外探測儀（Atmospheric InfraRed Sounder，AIRS）等衛(wèi)星高光譜數(shù)據(jù)，提取地表溫度（Land Surface Temperature，LST）、射出長波輻射（Outgoing Longwave Radiation，OLR）、氣溶膠光學厚度（Aerosol Optical Depth，AOD）、地表潛熱通量（Surface Latent Heat Flux，SLHF）以及氣體地球化學等信息，探究其與地震的關系，發(fā)現(xiàn)地震前后往往會存在LST、OLR、AOD、SLHF與氣體地球化學等眾多信息中的一種或多種異常［11-17］。目前，觀測到地震前后出現(xiàn)異常的氣體主要有H2、He、Rn、Ar、CH4、CO2、CO、O3等。崔月菊等［18］指出與地震有關的氣體異常一般表現(xiàn)為地震前后幾天到幾個月局地氣體濃度升高，大量研究表明地震活動會伴隨CO逸散。姚清林等［19］利用MOPITT數(shù)據(jù)研究發(fā)現(xiàn)2000年4月在青藏高原出現(xiàn)大面積CO濃度升高，分布上與多處孤立升溫具有較高的一致性，推斷異?，F(xiàn)象可能與2000年6月6日甘肅景泰MS5.9地震和同年6月8日緬甸MS6.9地震有關。Singh等［20］利用MOPITT傳感器數(shù)據(jù)，指出2001年1月26日Gujarat7.6級地震前，近地面CO濃度顯著升高，是地震前地應力集中對震中區(qū)的水文地質邊界產生影響使得地下氣體排放造成的。Cui等［21］利用AIRS資料研究了2010—2012年世界范圍35次7級以上地震前后的氣體變化，顯示出12次地震前后出現(xiàn)CO濃度異常。劉海博等［22］利用AIRS數(shù)據(jù)通過差值法與異常指數(shù)法研究與2014年2月12日新疆于田MS7.3地震有關的CO、O3信息時發(fā)現(xiàn)，兩者均在地震后濃度升高并沿發(fā)震構造呈線性分布，氣體濃度在地震當月呈最低值，推斷是由于震區(qū)斷裂帶在彈性擠壓狀態(tài)下出現(xiàn)閉鎖造成的，指出氣體逸散與地震構造活動關系密切。但是由于不同研究區(qū)受地理因素、構造背景等影響，地震前后氣體異常排放信息的提取缺乏針對性篩選，缺少多角度的分析。

2021年5月22日青?，敹郙S7.4地震（以下簡稱瑪多地震）發(fā)生后，眾多學者針對瑪多地震的構造背景［23］、震源機制［24］、地震破裂過程［25］、地震深部環(huán)境［26］以及地下流體［27-29］等進行了相應研究。鐘駿等［27］分析震中附近地下流體觀測資料發(fā)現(xiàn)瑪多地震前震中距500 km范圍內存在水溫、水位和氣氡異常。另外瑪多震前存在動水位突降，震后采集的水化學樣品檢測結果顯示靠近震中地表破裂帶的泉水存在氫同位素異常［28-29］。然而目前對瑪多地震的研究，缺乏相關氣體變化的研討，故本文主要研究與瑪多地震相關的CO地球化學信息，在以瓦里關大氣觀測站地面數(shù)據(jù)驗證AIRS反演數(shù)據(jù)可靠的基礎上，結合地震區(qū)域地形特征，基于AIRS傳感器的衛(wèi)星監(jiān)測數(shù)據(jù)，通過差值法和滑動均值法分析了CO的時空變化特征及其與瑪多地震的關系，并探討了CO異常原因。

1 研究區(qū)概況

2021年5月22日2時4分，青海果洛藏族自治州瑪多縣發(fā)生MS7.4地震，震中位于黃河鄉(xiāng)（34.59°N，98.34°E），震源深度17 km，同震地表破裂達210 km，發(fā)震斷層為NW走向的昆侖山口-江錯斷裂，以左旋走滑運動為主［25，30］。江錯斷裂北段兼有逆沖性質，南段以走滑性質為主，余震與主震震源機制較為一致，多為走滑型，序列活動主要受NEE—SWW向近水平擠壓應力場控制［31］。江錯斷裂屬于從主干斷裂分支、深入巴顏喀拉塊體內部的次級斷裂［32］（圖1），其活動性相比于主干斷裂帶明顯更弱，震級與較小的斷層滑動速度之間存在一定矛盾［23］，而且其所屬的巴顏喀拉塊體周圍存在地應力顯著增加的地震空區(qū)，至今構造變形方式與過程彌散連續(xù)［26］，該區(qū)域存在較高的研究價值。瑪多地震震中遠離工業(yè)和人口稠密區(qū)，CO受人類因素影響小，而且CO的大氣環(huán)境本底濃度較穩(wěn)定，能夠很好地反映局部環(huán)境變化特征［33-34］，是研究地震活動與氣體異常關系的良好案例。

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

2.1 數(shù)據(jù)

本文所使用的CO遙感數(shù)據(jù)來源于美國國家航空航天局（National Aeronautics and Space Administration，NASA）提供的大氣紅外探測儀（Atmospheric InfraRed Sounder，AIRS）level-3的8天平均、月平均標準產品中的降軌數(shù)據(jù)，在戈達德地球科學數(shù)據(jù)和信息服務中心（Goddard Earth Sciences Date and Information Services Center，GESD ISC）下載（http：//disc.sci.gsfc.nasa.gov/AIRS/data-holdings）。AIRS是搭載于NASA 2002年發(fā)射的一顆太陽同步極軌衛(wèi)星AQUA/EOS上的高光譜分辨率傳感器，該衛(wèi)星從2002年9月開始向地面提供觀測數(shù)據(jù)，可實現(xiàn)對全球85%的地區(qū)每天兩次覆蓋。AIRS紅外探測器用于CO掃描的通道有36個，光譜為4.50～4.58 μm，空間分辨率為1°×1°，是迄今為止最為精確和穩(wěn)定的高光譜紅外探測儀［35］。其數(shù)據(jù)存儲格式為NASA標準的磁盤存儲格式HDF（Hierarchical Data Format）類型，利用MATLAB軟件進行提取。

本文所使用的CO地面數(shù)據(jù)來源于瓦里關全球大氣本底站（Waliguan，WLG;36.287°N，100.896°E），距離本次震中約300 km，在世界溫室氣體數(shù)據(jù)中心（World Data Centre for Greenhouse Gases，WDCGG）下載（http：//gaw.kishou.go.jp/）。瓦里關站是世界氣象組織/全球大氣觀測（WMO/GAW）32個大氣本底觀測站之一，位于歐亞大陸腹地青藏高原東北緣，遠離人類活動密集區(qū)，大氣清潔。該站于1994年開始運行，并對溫室氣體、反應性氣體、太陽輻射、氣溶膠等信息進行長期的觀測［36］。

2.2 方法

AIRS反演CO的基本方法是偏導數(shù)趨于零法（Variational Path-Delay，VPD）［37-38］。針對AIRS通道集，通過改變大氣物理狀態(tài)，反復迭代計算獲得觀測輻射和由大氣狀態(tài)計算的輻射之間的最小平方和F［式（1）］，F(xiàn)越小則變量越接近真值。隨后對方程式進行全微分，尋找T、q及CO的值，使得F最小，分別求各自的偏導數(shù)趨于0，即可得ΔF→0［式（2）、式（3）］。

∑i［I（Vi）observed-I（Vi，T，q，CO）computed］2?≡F（，T，q，CO）（1）

ΔF=FΤΔT+FqΔq+FCOΔCO+ε（2）

FXi→0 （3）

式中：T、q及CO均為變量，分別代表氣溫、降水量和一氧化碳含量;ε是噪聲。

對于每個變量Xi對輻射殘差貢獻為FXiΔXi。

為研究地震前后CO在垂直方向的分布特征，選取震中位置不同高度層CO體積混合比（Volume Mixing Ratio，VMR）數(shù)據(jù)，鑒于8天體積混合比數(shù)據(jù)的穩(wěn)定性相對較差，故運用滑動均值法進行處理，分析地震前后一年時間段中震中位置在不同高度層的CO體積混合比變化特征。具體計算公式為：

μCO（l，m）=13∑m+1n=m-1VCO（l，m） （4）

式中：VCO（l，m）為震中m時段（2020年7月—2021年7月）l層（600～100 hPa大氣）的CO體積混合比;μCO（l，m）為經滑動均值法處理后地震前后CO體積混合比。

為研究地震前后CO在水平尺度的變化特征，選取2015年及2017—2020年歷史同期8天CO濃度數(shù)據(jù)，以排除2016年10月17日青海玉樹MS6.2地震的影響，采用滑動均值法計算瑪多地震前后各研究時段內的8天CO濃度背景值，以剔除地形地貌、季節(jié)變化等固有非震因素的影響［21］，然后通過差值法計算各研究時段CO濃度變化值。其具體計算公式為：

A（i，j，m，t）=G（i，j，m，t）-Gbac（i，j，m，t）（5）

Gbac（i，j，m，t）=∑t-1n=t-5［VCO（i，j，m-1，n）+VCO（i，j，m，n）+VCO（i，j，m+1，n）］15（6）

式中：VCO（i，j，m，t）為研究區(qū)內（i，j）處t年m時段CO濃度值;Gbac（i，j，m，t）為研究區(qū)內（i，j）處t年m時段的CO背景值;G（i，j，m，t）為地震前后（i，j）處t年m時段的CO濃度值;A（i，j，m，t）為地震前后（i，j）處t年m時段的CO濃度變化量。

3 結果分析

3.1 地面觀測資料與AIRS反演數(shù)據(jù)的比較驗證

由于瑪多地區(qū)平均海拔高達4 500 m［39］，AIRS反演的月尺度體積混合比數(shù)據(jù)在600 hPa及以上高度大量缺失，所以選擇2014—2021年震中區(qū)域AIRS反演的400 hPa和500 hPa高度層月尺度CO體積混合比數(shù)據(jù)（AIRS-500 hpa、AIRS-400 hPa）和CO柱濃度數(shù)據(jù)（AIRS-Tot CO）與瓦里關地面觀測站數(shù)據(jù)進行可靠性驗證。由圖2可見，AIRS反演的CO濃度與瓦里關觀測站CO濃度基本一致，總體上呈現(xiàn)出春夏高、秋冬低的周期變化，與瓦里關連續(xù)觀測多年CO春高秋低的變化特征較為吻合［40］。但2015年、2016年和2019年的8、9月，以及2020年的9、10、11月CO濃度較高，究其原因可能與天氣和地震有關。調查發(fā)現(xiàn)2015年與2016年的秋季多陰雨天，致使CO光化學氧化消耗減弱，從而在大氣中得到了一定的積累，因此，推斷2015年CO濃度高值由2015年10月12日34.36°N、98.209°E處MS5.3地震和天氣共同引起，2016年則主要是天氣所致;2019年和2020年CO濃度高值分別與2019年8月9日37.69°N、101.59°E處MS4.9地震和2020年12月24日34.15°N、98.32°E處MS4.5淺源地震有關。

計算AIRS反演的CO各濃度數(shù)據(jù)與瓦里關觀測站CO濃度的均值、均值差及標準差之差（表1）。結果顯示AIRS-500 hPa與瓦里關地面觀測站的均值差之差為3.49×10-9，月均標準差之差為3.97×10-9，均小于AIRS-400 hPa與瓦里關地面觀測站的均值差與月均標準差之差，說明AIRS-500 hPa數(shù)據(jù)比AIRS-400 hPa數(shù)據(jù)的系統(tǒng)誤差小、精度高。相關性計算結果顯示AIRS-Tot CO、AIRS-500 hPa、AIRS-400 hPa數(shù)據(jù)與瓦里關地面觀測站數(shù)據(jù)的相關系數(shù)分別為0.62、0.56、0.44，其中AIRS-Tot CO、AIRS-500 hPa與瓦里關觀測站的相關系數(shù)較高，一致性較強，表明AIRS反演的柱濃度數(shù)據(jù)與500 hPa壓力層數(shù)據(jù)較為準確，利用衛(wèi)星遙感數(shù)據(jù)提取的CO氣體地球化學信息是可靠的。

3.2 CO濃度變化特征

瑪多地震位于青藏高原東北部，震中像元（34°～35°N，98°～99°E）內高程梯度較低，底層的CO體積混合比數(shù)據(jù)可較好地表征近地面CO濃度，因此選擇震中像元內不同壓力層體積混合比數(shù)據(jù)，分析CO垂直結構變化特征。運用滑動均值法處理震中不同壓力層的8天CO體積混合比數(shù)據(jù)，得到600～100 hPa等壓面上CO分布，如圖3所示，CO在400 hPa以上變化劇烈，300～100 hPa整體較為穩(wěn)定，這可能與CO的分布特點、研究區(qū)地形等有關。近地面600 hPa的CO從2月17日開始出現(xiàn)波動，3月5日出現(xiàn)小幅度下降后又開始上升，3月21日達到最大值169.8×10-9，隨后開始逐漸下降并于6月趨于正常變化。400 hPa、500 hPa的CO整體變化趨勢與600 hPa的CO變化較為一致，在3月5日—5月16日CO呈現(xiàn)較高水平，但無明顯峰值，其余時間段呈現(xiàn)正常狀態(tài)。300 hPa及以下由于CO很難傳輸?shù)酱颂帲⑶沂艿降匦?、傳感器敏感度的影響，幾乎無變化，間接證明400 hPa以上CO劇烈變化是近地面貢獻的結果。Liang等［41］研究2008年3月21日和2014年2月12日于田兩次MS7.3地震前后CO變化，結果顯示CO主要聚集在近地面，本文所得結果與其一致。而且圖4顯示，震中CO柱濃度與近地面500 hPa處CO體積混合比變化五角星代表震中，黑色線條表示斷層，藍色框為發(fā)震時段趨勢基本一致，并通過了0.05的顯著性檢驗。這表明CO的增加主要是近地面的貢獻。

為了分析瑪多地震前后CO排放在水平方向上的演化過程，基于地源性CO排放具有近地面濃度升高的特征，底層CO體積混合比數(shù)據(jù)更能反映近地面CO變化情況。由于研究區(qū)域內600 hPa及以上高度層8天尺度CO體積混合比數(shù)據(jù)監(jiān)測空白較多，故選擇500 hPa高度層識別地源排放信息。首先進行去背景場處理，選取2015年及2017—2020年歷史同期數(shù)據(jù)，以排除2016年10月17日青海玉樹MS6.2地震的影響，運用滑動均值法計算各研究時段的背景值，再通過差值法計算研究區(qū)內大氣底層的8天CO體積混合比變化。不同高度層的CO體積混合比變化特征（圖3）顯示，在2021年3月至6月初CO體積混合比在近地面變化明顯，故選取該時間段內震中及附近區(qū)域進行時空變化特征分析。通過普通克里金插值得到近地面CO濃度空間分布圖（圖4），結果顯示，近地面CO最大異常值出現(xiàn)在地震前的4月份，異常分布于震中附近及斷裂帶沿線。相比于歷史同期背景值，CO濃度在整個研究時段內較高。地震發(fā)生前，3月CO異常逐漸增大并由離散分布逐漸聚攏到震中附近，到4月初異常幅度變小并向東移動，隨后到4月14日異常開始突然增強，主要分布在震中北部和沿斷裂帶沿線區(qū)域，延續(xù)到震前一個月左右異常突然消失，整個研究區(qū)域的CO恢復到正常水平。4月末CO異常突然增大達到地震前后的最大程度，最大異常超出背景值18.60×10-9，異常分布于發(fā)震斷層的南部。隨后異常由震中南部逐漸擴散轉移到震中東南部及北部地區(qū)，并逐漸減小。直至地震發(fā)生前后一周左右［圖4（j）］，CO在震中的北部又出現(xiàn)明顯的異常，但異常程度相比震前較小。地震發(fā)生后，異常逐漸減小，并在震中附近形成了一個明顯的低值區(qū)域。

為了排除季節(jié)變化和其他因素對CO濃度變化的影響，選取震中像元2014-2022年間4月份8天尺度CO柱濃度、500 hPa處CO體積混合比數(shù)據(jù)進行分析，結果（圖5）顯示，CO柱濃度與CO體積混合比在2021年4月（黑色方框）高于歷史同期值，在4月30日左右達到最大值，而在其他非地震年的4月CO柱濃度與CO體積混合比濃度值較小。僅2015年4月14日和2016年4月16日兩處（箭頭位置）CO柱濃度與CO的體積混合比較其他值偏高，可能分別與2015年2月4日34.73°N、98.99°E處MS4.1地震和2016年3月21日35.74°N、99.20°E處MS4.1地震有關。這說明2021年4月CO濃度的異常變化不是由季節(jié)變化引起的，很可能與地震有關。

4 討論

4.1 CO異常與地震的關系

對瑪多地震前后CO濃度變化分析發(fā)現(xiàn)：CO異常最大值出現(xiàn)在震前的4月份;垂直方向主要表現(xiàn)為近地面濃度升高，水平方向上異常主要集中在震中及附近區(qū)域和斷裂帶沿線。地震前后CO異常在空間分布上呈現(xiàn)出分散-聚攏-分散的變化，異常程度由弱增強然后減弱。地震前兩個月震區(qū)周圍應力緩慢加載，到4月位于兩分段傾角變化處的震源位置易于形成應力閉鎖區(qū)，在周緣塊體的推擠下自西向東運動，震源處的應力沿斷層面向東、向西擴展，形成了一個條帶狀的應力集中區(qū)［42］，應力的空間分布與CO異常分布一致，氣體異常程度與瑪多地震烈度由西向東逐漸減小的分布［32］具有較好的對應性。地震前一個月左右，CO異常達到最大，整體分布于震中南部，空間上與斷裂帶不完全重合，但與地下流體異常分布相似［27］，可能由深部流體運動、地形特征和氣體運移等因素造成。Qi等［39］運用微波亮度溫度（Microwave Brightness Temperature，MWBT）數(shù)據(jù)研究發(fā)現(xiàn)，瑪多地震有關的MWBT異常在地震前出現(xiàn)，空間上分布于震中的南部及西南部，本文中CO異常分布特征與Qi等［39］發(fā)現(xiàn)的MWBT空間分布結果相似。發(fā)震時段，整個研究區(qū)CO異常程度較小，發(fā)震后在震中形成了低值區(qū)域，推測是因為瑪多地震的發(fā)生使得震源區(qū)積累的應力得到了充分的釋放［43］，CO也得以釋放，故而震后氣體的釋放量大大減小。Jing等［44］通過異常指數(shù)法處理與瑪多地震有關的MWBT和OLR數(shù)據(jù)，發(fā)現(xiàn)MWBT的異常變化從3月開始出現(xiàn)于巴顏喀拉塊體及震中西南部，4月異常逐漸增大，并在空間上分布于發(fā)震斷裂帶周圍;OLR的變化在震前的5月出現(xiàn)于震中的西南部，隨后不斷向震中移動并向北擴散，最后在發(fā)震前半個月消失。綜上所述，瑪多地震相關的異常信息多出現(xiàn)在震前，空間上分布于震中及周圍，并與斷裂帶的走向存在較高的一致性。

4.2 CO背景含量及來源分析

大氣中CO的來源主要有地殼活動、海洋釋放、化石燃料燃燒、化學反應（甲烷氧化、非甲烷碳氫化合物氧化和光化學反應）和植被自燃等［22，45-46］。首先，此次瑪多地震區(qū)域處于寒溫帶高原區(qū)，該地遠離人口和工業(yè)稠密區(qū)，海洋、森林火災的發(fā)生率低［33］，所以CO濃度不受海洋釋放、化石燃料燃燒與植被自燃的影響。其次，周凌晞［47］和Zhou等［48］在研究瓦里關觀測站CO來源時發(fā)現(xiàn)，來自東北方約100 km處青海西寧市和東北方約500 km處黃河峽谷工業(yè)及人口集中地帶的污染氣團對該地CO濃度存在影響，但這種污染源對該區(qū)域的CO柱濃度造成抬升的現(xiàn)象大多發(fā)生在冬季，故此次CO異常現(xiàn)象受大氣傳輸?shù)挠绊戄^小。其次，Granier等［45］指出CH4氧化作用所產生的CO占對流層中CO含量的28%，因而大氣CO柱濃度變化與CH4的氧化作用關系密切。而本研究選取震中區(qū)2020—2021年間CH4柱濃度數(shù)據(jù)進行相關性驗證的結果表明，這期間CO與CH4柱濃度相關系數(shù)為-0.132，在地震前后短時間內（2021年3—7月）兩者的相關系數(shù)為-0.149，不存在相關性，所以瑪多地震前后CO異常不受CH4氧化作用的影響。然后，非甲烷碳氫化合物主要包括異戊二烯、乙烷、丙烷、乙烯、丙烯等，異戊二烯的排放一般出現(xiàn)在熱帶區(qū)域，其余幾種主要產生于工業(yè)生產和生物質燃燒［45］，所以非甲烷碳氫化合物的氧化對此次CO異常也不產生影響。最后，瑪多地震發(fā)生時間為5月，一般來說OH自由基濃度在夏季最高，CO作為OH自由基最主要的消耗者［49］，通常表現(xiàn)為春季達到最大、夏季降至最小。然而該區(qū)域夏季大氣光化學過程較弱，秋冬季晴好天氣相對較多［50］，所以瑪多地震區(qū)及附近OH自由基對CO柱濃度的影響較小。經分析發(fā)現(xiàn)地震前后CO濃度受外界環(huán)境的干擾較小，種種跡象表明此次CO濃度異常主要與地震活動有關。

4.3 CO異常原因分析

瑪多地震前CO異常主要歸因于地下氣體釋放、巖石擠壓產氣和14N化學反應。首先，地球作為一個復雜體，在地震這個復雜過程中，往往表現(xiàn)出不盡相同的地球內部氣體釋放類型與強度［51］。構造演化和地震孕育過程中地球內部物質遷移、能量釋放及應力變化可能導致活動斷裂帶釋放氣體，進而出現(xiàn)氣體地球化學異常現(xiàn)象［52］。地質深部構造顯示，巴顏喀拉塊體的N向俯沖與瑪多地下的特殊鍥狀形態(tài)導致區(qū)域內深部地殼軟物質流擠壓、垂直上涌變形［23］。孕震過程中，地下流體中的水化學類型和水-巖作用程度發(fā)生顯著變化［28-29］，隨著構造應力的不斷積累，貫穿地表裂隙數(shù)量增多，地球內部碳質氣體（CO或CO2）的排放增強［44］，氣體在巖石圈和大氣圈的作用下導致了不同海拔氣體異常的強度不同［53］。在瑪多地震前表現(xiàn)為震中附近局部地區(qū)尤其是近地面的CO濃度升高，因此推斷由地下作用釋放的氣體為CO異常的主要貢獻者。水平方向上異常高值的空間分布與發(fā)震斷層走向一致，與地下流體東北部與西南部的異常具有較高的對應性［27］。表明活動斷裂帶為地球深部氣體的集中釋放地帶，地震前后的氣體異常現(xiàn)象與地震本身有密不可分的聯(lián)系，并由地震活動觸發(fā)［54］。其次，研究表明起源于地幔超基性巖礦物包裹體和地球深部的巖漿噴出地表所形成的火山巖玻璃包裹體中會含有部分CO［51］。該地區(qū)存在的流紋巖［55］在摩擦擠壓、碰撞過程中，會產生CO并通過地殼裂隙釋放到大氣中。最后，大氣圈中時刻在發(fā)生化學反應，瑪多地震發(fā)生前后存在地震電磁擾動［56-57］，導致電離層物質和電磁波輻射增強，引起等離子體（電子、離子）化學、物理參數(shù)發(fā)生變化［58-59］，促進14N的生成，14N經過一系列化學反應產生CO，反應機制為n+14N→14C+H、214C+O2→214CO［60］。

5 結論

通過分析2021年5月22日青?，敹郙S7.4地震前后震中附近的CO異常變化，得到以下結論：

（1） AIRS反演結果與瓦里關地面觀測站的CO月尺度濃度變化均表現(xiàn)出春夏高秋冬低的周期變化，反演的CO柱濃度與500 hPa壓力層體積混合比與地面觀測結果的相關性較高。表明利用AIRS反演的CO數(shù)據(jù)提取地球化學信息是可靠的。

（2） CO異常在地震前兩個月左右出現(xiàn)異常，震后逐漸恢復平靜，主要表現(xiàn)為近地面濃度升高?；诓钪捣ㄈケ尘皥鎏幚淼慕Y果顯示，震中及附近區(qū)域近地面CO在地震前兩個月濃度逐漸升高并由離散分布逐漸向震中聚斂靠近，到4月底沿斷裂帶NW向分布在震中南部，達到最大異常18.60×10-9，震后逐漸恢復平靜。氣體異常高值的空間分布與地震烈度分布對應性較好，和發(fā)震斷裂帶江錯斷裂走向具有一致性。

（3） 地震孕育過程中地下流體中的氣體以及巖石擠壓與破裂產生的氣體會通過地裂縫逸散至大氣圈，引起近地面CO濃度升高，此外14N在大氣圈發(fā)生化學反應也會產生CO。推斷瑪多地震前CO的異常為巖石圈和大氣圈綜合作用的結果，地下釋放為主要因素，大氣圈的作用為次要因素。

本文在結合地面監(jiān)測數(shù)據(jù)與衛(wèi)星監(jiān)測數(shù)據(jù)的同時考慮了研究區(qū)域的地形特征，表明地震前后震中附近CO出現(xiàn)異常的特征可用于地震前兆研究，豐富了地球化學信息在地震預報與監(jiān)測中的應用，對地震的監(jiān)測具有一定的意義。同時還有利于對江錯斷裂帶構造活動的研究，為巴顏喀拉塊體地震空區(qū)的進一步研究提供了參考數(shù)據(jù)。但是本文選擇的氣體參數(shù)較為單一，在今后的研究中需要結合多參量信息，以便更準確地監(jiān)測預報地震。

參考文獻（References）

［1］ 陳彧，徐瑞松，蔡睿，等.遙感技術在地震研究中的應用進展［J］.地球物理學進展，2008，23（4）：1273-1281.

CHEN Yu，XU Ruisong，CAI Rui，et al.Application progress of remote sensing technology in earthquake research［J］.Progress in Geophysics，2008，23（4）：1273-1281.

［2］ 杜樂天.地球排氣作用的重大意義及研究進展［J］.地質論評，2005，51（2）：174-180.

DU Letian.Significance of earth degassing and its research progress［J］.Geological Review，2005，51（2）：174-180.

［3］ 劉耀煒，任宏微，張磊，等.魯?shù)?.5級地震地下流體典型異常與前兆機理分析［J］.地震地質，2015，37（1）：307-318.

LIU Yaowei，REN Hongwei，ZHANG Lei，et al.Underground fluid anomalies and the precursor mechanisms of the Ludian MS6.5 earthquake［J］.Seismology and Geology，2015，37（1）：307-318.

［4］ 崔月菊.大地震前后CO、O3和CH4遙感地球化學異常特征［D］.北京：中國地質大學，2014.

CUI Yueju.Characteristics of remote sensing geochemical anomalies of CO，O3 and CH4 before and after the earthquake［D］.Beijing：China University of Geosciences，2014.

［5］ KING C Y，ZHANG W，ZHANG Z C.Earthquake-induced groundwater and gas changes［J］.Pure and Applied Geophysics，2006，163（4）：633-645.

［6］ 強祖基，孔令昌，王弋平，等.地球放氣、熱紅外異常與地震活動［J］.科學通報，1992，37（24）：2259-2262.

QIANG Zuji，KONG Lingchang，WANG Yiping，et al.Earth outgassing，thermal infrared anomaly and seismic activity［J］.Chinese Science Bulletin，1992，37（24）：2259-2262.

［7］ 崔月菊，李靜，王燕艷，等.遙感氣體探測技術在地震監(jiān)測中的應用［J］.地球科學進展，2015，30（2）：284-294.

CUI Yueju，LI Jing，WANG Yanyan，et al.Application of gas remote sensing technique to earthquake monitoring［J］.Advances in Earth Science，2015，30（2）：284-294.

［8］ 趙振燊.甘東南地震重點危險區(qū)主要活動斷裂帶斷層氣地球化學特征［D］.蘭州：中國地震局蘭州地震研究所，2012.

ZHAO Zhenshen.Geochemical characteristics of fault gas in main active fault zones in key danger areas of southeast Gansu earthquake［D］.Lanzhou：Lanzhou Institute of Seismology，CEA，2012.

［9］ 強祖基，杜樂天.地球排氣與森林火災和地震活動［J］.地學前緣，2001，8（2）：235-245.

QIANG Zuji，DU Letian.Earth degassing，forest fire and seismic activities［J］.Earth Science Frontiers，2001，8（2）：235-245.

［10］ 申旭輝，吳云，單新建.地震遙感應用趨勢與中國地震衛(wèi)星發(fā)展框架［J］.國際地震動態(tài)，2007，37（8）：38-45.

SHEN Xuhui，WU Yun，SHAN Xinjian.Remote sensing application in earthquake science and general proposal for earthquake satellite project in China［J］.Recent Developments in World Seismology，2007，37（8）：38-45.

［11］ 郭廣猛，曹云剛，龔建明.使用MODIS和MOPITT衛(wèi)星數(shù)據(jù)監(jiān)測震前異常［J］.地球科學進展，2006，21（7）：695-698，766.

GUO Guangmeng，CAO Yungang，GONG Jianming.Monitoring anomaly before earthquake with MODIS and MOPITT data［J］.Advances in Earth Science，2006，21（7）：695-698，766.

［12］ PULINETS S A，OUZOUNOV D，CIRAOLO L，et al.Thermal，atmospheric and ionospheric anomalies around the time of the Colima M7.8 earthquake of 21 January 2003［J］.Annales Geophysicae，2006，24（3）：835-849.

［13］ 崔麗華.汶川地震前的遙感信息異常及其機理研究［D］.唐山：河北理工大學，2009.

CUI Lihua.Study on the anomaly of remote sensing information before Wenchuan earthquake and its mechanism［D］.Tangshan：Hebei Polytechnic University，2009.

［14］ SINGH R P，MEHDI W，GAUTAM R，et al.Precursory signals using satellite and ground data associated with the Wenchuan earthquake of 12 May 2008［J］.International Journal of Remote Sensing，2010，31（13）：3341-3354.

［15］ 荊鳳，申旭輝，張鐵寶，等.與地震有關的活動斷裂帶紅外輻射變化特征［J］.國土資源遙感，2013，25（1）：56-60.

JING Feng，SHEN Xuhui，ZHANG Tiebao，et al.Variation characteristics in infrared radiation of active fault zone related to earthquakes［J］.Remote Sensing for Land & Resources，2013，25（1）：56-60.

［16］ 孫玉濤，崔月菊，劉永梅，等.蘇門答臘2004、2005年兩次大地震前后CO和O3遙感信息［J］.遙感信息，2014，29（2）：47-53.

SUN Yutao，CUI Yueju，LIU Yongmei，et al.Remote sensing anomalies of CO and O3 related to two giant Sumatra earthquakes occurred in 2004 and 2005［J］.Remote Sensing Information，2014，29（2）：47-53.

［17］ SHAH M，TARIQ M A，NAQVI N A.Atmospheric anomalies associated with MW>6.0 earthquakes in Pakistan and Iran during 2010—2017［J］.Journal of Atmospheric and Solar-Terrestrial Physics，2019，191：105056.

［18］ 崔月菊，杜建國，荊鳳，等.2008年汶川MS8.0地震前后川西含碳氣體衛(wèi)星高光譜特征［J］.地震學報，2016，38（3）：448-457，508.

CUI Yueju，DU Jianguo，JING Feng，et al.Mapping emission of carbon-bearing gases from the satellite hyperspectral data in western Sichuan before and after the 2008 Wenchuan MS8.0 earthquake［J］.Acta Seismologica Sinica，2016，38（3）：448-457，508.

［19］ 姚清林，強祖基，王弋平.青藏高原地震前CO的排放與衛(wèi)星熱紅外增溫異常［J］.地球科學進展，2005，20（5）：505-510.

YAO Qinglin，QIANG Zuji，WANG Yiping.CO release from the Tibetan Plateau before earthquakes and increasing temperature anomaly showing in thermal infrared images of satellite［J］.Advance in Earth Sciences，2005，20（5）：505-510.

［20］ SINGH R P，SENTHIL KUMAR J，ZLOTNICKI J，et al.Satellite detection of carbon monoxide emission prior to the Gujarat earthquake of 26 January 2001［J］.Applied Geochemistry，2010，25（4）：580-585.

［21］ CUI Y，DU J，ZHANG D，et al.Anomalies of total column CO and O3 associated with great earthquakes in recent years［J］.Natural Hazards and Earth System Sciences，2013，13（10）：2513-2519.

［22］ 劉海博，崔月菊，辛存林.探測與2014年新疆于田MS7.3地震相關的大氣CO和O3異常變化［J］.地震，2020，40（1）：99-111.

LIU Haibo，CUI Yueju，XIN Cunlin.Detecting anomalies of atmospheric total column CO and O3 related to the 2014 Yutian MS7.3 earthquake［J］.Earthquake，2020，40（1）：99-111.

［23］ 宋向輝，王帥軍，潘素珍，等.2021年瑪多MS7.4地震的深部構造背景［J］.地震地質，2021，43（4）：757-770.

SONG Xianghui，WANG Shuaijun，PAN Suzhen，et al.Deep seismotectonic environment of the 2021 Maduo MS7.4 earthquake［J］.Seismology and Geology，2021，43（4）：757-770.

［24］ 張志朋，李君，馮兵，等.2021年青?，敹郙S7.4地震序列精定位與震源機制研究［J］.地震工程學報，2022，44（1）：218-226.

ZHANG Zhipeng，LI Jun，F(xiàn)ENG Bing，et al.Precise location and focal mechanism solutions of the 2021 Maduo，Qinghai MS7.4 earthquake sequence［J］.China Earthquake Engineering Journal，2022，44（1）：218-226.

［25］ 李智敏，李文巧，李濤，等.2021年5月22日青?，敹郙S7.4地震的發(fā)震構造和地表破裂初步調查［J］.地震地質，2021，43（3）：722-737.

LI Zhimin，LI Wenqiao，LI Tao，et al.Seismogenic fault and coseismic surface deformation of the Maduo MS7.4 earthquake in Qinghai，China：a quick report［J］.Seismology and Geology，2021，43（3）：722-737.

［26］ 詹艷，梁明劍，孫翔宇，等.2021年5月22日青?，敹郙S7.4地震深部環(huán)境及發(fā)震構造模式［J］.地球物理學報，2021，64（7）：2232-2252.

ZHAN Yan，LIANG Mingjian，SUN Xiangyu，et al.Deep structure and seismogenic pattern of the 2021.5.22 Madoi （Qinghai） MS7.4 earthquake［J］.Chinese Journal of Geophysics，2021，64（7）：2232-2252.

［27］ 鐘駿，王博，周志華，等.2021年青?，敹郙S7.4地震前地下流體異常特征分析［J］.中國地震，2021，37（3）：574-585.

ZHONG Jun，WANG Bo，ZHOU Zhihua，et al.Analysis on anomaly characteristics of underground fluid before 2021 Maduo MS7.4 earthquake in Qinghai Province［J］.Earthquake Research in China，2021，37（3）：574-585.

［28］ 路暢，周曉成，李營，等.瑪多MS7.4地震地表破裂帶與東昆侖斷裂溫泉的水文地球化學特征［J］.地震地質，2021，43（5）：1101-1126..

LU Chang，ZHOU Xiaocheng，LI Ying，et al.Hydrogeochemical characteristics of ground-water in the surface rupture zone of Madoi MS7.4 earthquake and hot springs in the east Kunlun fault［J］.Seismology and Geology，2021，43（5）：1101-1126.

［29］ 蘇維剛，劉磊，孫璽皓.瑪多7.4級地震和門源6.9級地震前佐署地下流體異常特征分析［J］.地震工程學報，2022，44（3）：700-706，712.

SU Weigang，LIU Lei，SUN Xihao.Characteristics of underground fluid anomalies in Zuoshu station before Maduo MS7.4 and Menyuan MS6.9 earthquakes［J］.China Earthquake Engineering Journal，2022，44（3）：700-706，712.

［30］ 華俊，趙德政，單新建，等.2021年青?，敹郙W7.3地震InSAR的同震形變場、斷層滑動分布及其對周邊區(qū)域的應力擾動［J］.地震地質，2021，43（3）：677-691.

HUA Jun，ZHAO Dezheng，SHAN Xinjian，et al.Coseismic deformation field，slip distribution and coulomb stress disturbance of the 2021 MW7.3 Maduo earthquake using Sentinel-1 InSAR observations［J］.Seismology and Geology，2021，43（3）：677-691.

［31］ 趙韜，王瑩，馬冀，等.2021年青?，敹?.4級地震序列重定位和震源機制特征［J］.地震地質，2021，43（4）：790-805.

ZHAO Tao，WANG Ying，MA Ji，et al.Relocation and focal mechanism solutions of the 2021 Maduo，Qinghai MS7.4 earthquake sequence［J］.Seismology and Geology，2021，43（4）：790-805.

［32］ 蓋海龍，姚生海，楊麗萍，等.青?，敹唷?·22”MS7.4級地震的同震地表破裂特征、成因及意義[J].地質力學學報，2021，27（6）：899-912.

GAI Hailong，YAO Shenghai，YANG Liping，et al.Characteristics and causes of coseismic surface rupture triggered by the “5·22” MS7.4 earthquake in Maduo，Qinghai，and their significance[J].Journal of Geomechanics，2021，27（6）：899-912.

［33］ 周凌晞，溫玉璞，李金龍，等.瓦里關山大氣CO本底變化［J］.環(huán)境科學學報，2004，24（4）：637-642.

ZHOU Lingxi，WEN Yupu，LI Jinlong，et al.Background variation in atmospheric carbon monoxide at Mt. Waliguan，China［J］.Acta Scientiae Circumstantiae，2004，24（4）：637-642.

［34］ DANIEL J S，SOLOMON S.On the climate forcing of carbonmonoxide［J］.Journal of Geophysical Research （Atmospheres），1998，103（D11）：13249-13260.

［35］ 蔣耿明，金亞秋.中國風云二號SVISSR傳感器熱紅外波段與AIRS/Aqua的交叉輻射定標［J］.遙感技術與應用，2012，27（4）：504-510，495.

JIANG Gengming，JIN Yaqiu.Intercalibration of SVISRR/FY-2 thermal infrared channels with AIRS/Aqua channels［J］.Remote Sensing Technology and Application，2012，27（4）：504-510，495.

［36］ 王劍瓊，祁棟林，薛麗梅.瓦里關全球大氣本底站站址環(huán)境及部分觀測要素變化特征［J］.環(huán)境化學，2020，39（8）：2084-2092.

WANG Jianqiong，QI Donglin，XUE Limei.Characteristics of the station environment and some atmospheric compositions of China global atmosphere watch base line observatory at Mt. Waliguan［J］.Environmental Chemistry，2020，39（8）：2084-2092.

［37］ SUSSKIND J，BARNET C D，BLAISDELL J M.Retrieval of atmospheric and surface parameters from AIRS/AMSU/HSB data in the presence of clouds［J］.IEEE Transactions on Geoscience and Remote Sensing，2003，41（2）：390-409.

［38］ MCMILLAN W W.Daily global maps of carbon monoxide from NASA's atmospheric infrared sounder［J］.Geophysical Research Letters，2005，32（11）：L11801.

［39］ QI Y A，WU L X，DING Y F，et al.Extraction and discrimination of MBT anomalies possibly associated with the MW7.3 Maduo （Qinghai，China） earthquake on 21 May 2021［J］.Remote Sensing，2021，13（22）：4726.

［40］ 周凌晞，湯潔，ERNST M，等.中國西部本底大氣中CO的連續(xù)測量［J］.環(huán)境科學，2001，22（3）：1-5.

ZHOU Lingxi，TANG Jie，ERNST M，et al.Continuous measurement of baseline atmospheric carbon monoxide in western China［J］.Chinese Journal of Enviromental Science，2001，22（3）：1-5.

［41］ LIANG H D，XIN C L，LIU H B，et al.Remotesensing of atmospheric CO and O3 anomalies before and after two Yutian MS7.3 earthquakes［J］.Geofluids，2021，2021：1-14.

［42］ 吳立新，盧菁琛，毛文飛，等.基于斷層傾角分段變化的瑪多地震發(fā)震斷層構造應力場演化數(shù)值模擬分析［J］.地球物理學報，2022，65（10）：3844-3857.

WU Lixin，LU Jingchen，MAO Wenfei，et al.Sectional fault-inclination-change based numerical simulation of tectonic stress evolution on the seismogenic fault of Madoi earthquake［J］.Chinese Journal of Geophysics，2022，65（10）：3844-3857.

［43］ 祝愛玉，王永哲，李永華，等.基于InSAR地表形變約束的瑪多MS7.4地震孕育發(fā)生機理數(shù)值模擬研究［J］.地球物理學報，2021，64（12）：4548-4561.

ZHU Aiyu，WANG Yongzhe，LI Yonghua，et al.Numerical simulation on the mechanism of the Madoi，Qinghai MS7.4 earthquake constrained by InSAR deformation［J］.Chinese Journal of Geophysics，2021，64（12）：4548-4561.

［44］ JING F，ZHANG L，SINGH R P.Pronounced changes in thermal signals associated with the Madoi （China） M7.3 earthquake from passive microwave and infrared satellite data［J］.Remote Sensing，2022，14（11）：2539.

［45］ GRANIER C，PTRON G，MLLER J F，et al.The impact of natural and anthropogenic hydrocarbons on the tropospheric budget of carbon monoxide［J］.Atmospheric Environment，2000，34（29/30）：5255-5270.

［46］ SEILER W.The cycle of atmospheric CO［J］.Tellus，1974，26（1/2）：116-135.

［47］ 周凌晞.中國大陸地區(qū)主要溫室氣體本底特征研究［D］.北京：北京大學，2001.

ZHOU Lingxi.Study on background characteristics of main greenhouse gases in Chinese Mainland area［D］.Beijing：Peking University，2001.

［48］ ZHOU L X，TANG J E，WEN Y P，et al.The impact of local winds and long-range transport on the continuous carbon dioxide record at Mount Waliguan，China［J］.Tellus B：Chemical and Physical Meteorology，2003，55（2）：145.

［49］ NOVELLI P C，MASARIE K A，LANG P M.Distributions and recent changes of carbon monoxide in the lower troposphere［J］.Journal of Geophysical Research：Atmospheres，1998，103（D15）：19015-19033.

［50］ MA J Z，TANG J，ZHOU X J，et al.Estimates of the chemical budget for ozone at Waliguan observatory［J］.Journal of Atmospheric Chemistry，2002，41（1）：21-48.

［51］ 鄭樂平.溫室氣體CO2的另一源：地球內部［J］.環(huán)境科學研究，1998，11（2）：21-24.

ZHENG Leping.Another source of grenhouse gas CO2：the Earth's interior［J］.Research of Environmental Sciences，1998，11（2）：21-24.

［52］ ITALIANO F，BONFANTI P，DITTA M，et al.Helium and carbon isotopes in the dissolved gases of Friuli region （NE Italy）：geochemical evidence of CO2 production and degassing over a seismically active area［J］.Chemical Geology，2009，266（1/2）：76-85.

［53］ CUI Y，ZHENG C，JIANG L，et al.Variations of multiple gaseous emissions associated with the great Sumatra earthquakes in 2004 and 2005［J］.Chemical Geology，2023，618：121311.

［54］ CUI Y，OUZOUNOV D，HATZOPOULOS N，et al.Satellite observation of CH4 and CO anomalies associated with the Wenchuan MS8.0 and Lushan MS7.0 earthquakes in China［J］.Chemical Geology，2017，469：185-191.

［55］ 張雪亭.青海省大地構造格架研究［D］.北京：中國地質大學，2006.

ZHANG Xueting.Study on tectonic framework in Qinghai Province［D］.Beijing：China University of Geosciences，2006.

［56］ 李霞，馮麗麗，趙玉紅，等.2021年5月22日青?，敹郙S7.4地震前地磁加卸載響應比變化特征［J］.地球物理學進展，2022，37（3）：964-970.

LI Xia，F(xiàn)ENG Lili，ZHAO Yuhong，et al.Analysis on characteristics of geomagnetic load-unload response ratio method before Maduo MS7.4 earthquake［J］.Progress in Geophysics，2022，37（3）：964-970.

［57］ 辛建村，張晨蕾.漾濞縣MS6.4、瑪多縣MS7.4地震前地電場變化分析［J］.地震工程學報，2021，43（4）：818-825.

XIN Jiancun，ZHANG Chenlei.Geoelectric field changes before the Yangbi MS6.4 and Maduo MS7.4 earthquakes［J］.China Earthquake Engineering Journal，2021，43（4）：818-825.

［58］ 丁鑒海，申旭輝，潘威炎，等.地震電磁前兆研究進展［J］.電波科學學報，2006，21（5）：791-801.

DING Jianhai，SHEN Xuhui，PAN Weiyan，et al.Seismo-electromagnetism precursor research progress［J］.Chinese Journal of Radio Science，2006，21（5）：791-801.

［59］ 張學民，劉靜，錢家棟，等.西藏改則6.9級地震前的電離層電磁擾動［J］.地震，2008，28（3）：14-22.

ZHANG Xuemin，LIU Jing，QIAN Jiadong，et al.Ionospheric electromagnetic disturbance before Gaize earthquake with MS6.9，Tibet［J］.Earthquake，2008，28（3）：14-22.

［60］ PANDOW M，MACKAY C，WOLFGANG R.The reaction of atomic carbon with oxygen：significance for the natural radio-carbon cycle［J］.Journal of Inorganic and Nuclear Chemistry，1960，14（3/4）：153-158.

（本文編輯：張向紅）

基金項目：國家自然科學基金資助項目（41262001）

第一作者簡介：史延飛（1997-），女，碩士研究生，主要從事地質工程。E-mail：syf18693271702@163.com。

通信作者：辛存林（1967-），男，博士，教授，主要從事地質環(huán)境與地質礦產等方面研究。E-mail：xincunling@163.com。

史延飛，辛存林，梁浩東，等.2021年青?，敹郙S7.4地震前后CO異常變化研究［J］.地震工程學報，2024，46（3）：692-702.DOI：10.20000/j.1000-0844.20230105002

SHI Yanfei，XIN Cunlin，LIANG Haodong，et al.CO anomalies before and after the 2021 Maduo MS7.4 earthquake in Qinghai Province［J］.China Earthquake Engineering Journal，2024，46（3）：692-702.DOI：10.20000/j.1000-0844.20230105002