2015—2017年南极长城站地区地震活动特征

郭钊玥, 常利军

郭钊玥,常利军. 2024. 2015—2017年南极长城站地区地震活动特征. 地震学报,46(6):936−952. DOI: 10.11939/jass.20230055
引用本文: 郭钊玥,常利军. 2024. 2015—2017年南极长城站地区地震活动特征. 地震学报,46(6):936−952. DOI: 10.11939/jass.20230055
Guo Z Y,Chang L J. 2024. Seismicity of the Antarctic Great Wall Station area during 2015−2017. Acta Seismologica Sinica46(6):936−952. DOI: 10.11939/jass.20230055
Citation: Guo Z Y,Chang L J. 2024. Seismicity of the Antarctic Great Wall Station area during 2015−2017. Acta Seismologica Sinica46(6):936−952. DOI: 10.11939/jass.20230055

2015—2017年南极长城站地区地震活动特征

基金项目: 国家自然科学基金(41930218,42374079)和中国地震局地球物理研究所自主立项项目(JY2023Z01)联合资助
详细信息
    作者简介:

    郭钊玥,在读硕士研究生,主要从事南极地震活动性等方面的研究,e-mail:2975931293@qq.com

    通讯作者:

    常利军,博士,研究员,主要从事地球深部构造、地震各向异性和地球动力学方面的研究,e-mail:ljchang@cea-igp.ac.cn

  • 中图分类号: 315.5

Seismicity of the Antarctic Great Wall Station area during 2015−2017

  • 摘要:

    基于交错网格搜索的单台定位方法对长城站地震台2015—2017年期间记录的宽频带三分量波形数据进行了单台地震定位研究,定位了112个高信噪比的近震,获得了此期间南设得兰群岛及其周边地区的地震分布图。结果显示,地震主要沿NE向展布于南设得兰群岛及邻区,最大震级为ML4.8,最小震级为ML1.3,震源深度范围为2—134 km。相对于稳定的南极大陆内部,长城站地区的地震活动性较强,在靠近南设得兰海沟的弧前区域和布兰斯菲尔德海峡弧后区域发现了可能与俯冲和岩浆活动有关的地震事件,弧前中源深度的地震事件说明该区域可能存在俯冲作用,这直接给出了菲尼克斯板块的俯冲活动仍在继续的地震学证据。

    Abstract:

    The South Shetland Islands are located at the northern part of the Antarctic Peninsula in the West Antarctica. The region is subject to the direct or indirect influences from the Pacific Plate, the Antarctic Plate, the South American Plate, the Scotia Plate and the Phoenix Plate, which has resulted in significant tectonic activities, including subduction, rifting and volcanic magmatism. The global earthquake catalog shows that this region exhibits a higher level of seismicity in comparison to the stable Antarctic continent interior. The extreme climatic conditions of the Antarctic continent present significant challenges to seismic observation, and the level of seismological research conducted on the continent is limited. The distribution of seismic stations in the South Shetland Islands region is relatively sparse, which limits the seismic monitoring capacity of the region and consequently affects our understanding of the seismicity in the region. The Chinese Antarctic Great Wall Seismic Station (CCZ) is located on King George Island in the South Shetland Islands. Following the reconstruction and resumption of observation in 2010, the quality of seismic observation has been remarkably enhanced, thereby considerably improving the monitoring capability of small earthquakes in the Great Wall Station area and providing data support for a comprehensive understanding of the seismicity in the South Shetland Islands region.

    We used the single station location method of interlaced grid search to obtain earthquake catalog based on the broadband three-component waveform data recorded at CCZ during 2015−2017. Firstly, the P wave and S wave arrivals of the local seismic events are identified through visual analysis, and the events with clear Pg and Sg waves are selected for amplitude measurements. The arrival time difference, incidence angle and back azimuth of the events can be calculated. The interlaced grid search method initially establishes a set of grids for searching and location, and searching for three elements of arrival time difference, incidence angle, and back azimuth, respectively. Then subsequently spatially translates the grids to obtain the location results of searching different grids. Finally, it evaluates the location results by calculating the standard deviation of the results obtained by the same set of data in the search of different grids. This method effectively overcomes the limitations of single set grid searching and location by a spatial intersection of multiple sets of grids. It has the potential to enhance the spatial sampling rate of earthquake location while simultaneously ensure the reliability and precision of the outcomes through the integration of multiple grid systems in space.

    In this study, 1 452 local seismic events were identified, and 112 events with high signal-to-noise ratios were located, of which 87% of the events have an epicentral distance standard deviation and a depth standard deviation smaller than the grid division size of 1.0 km, indicating that the differences of grid models have less influence on the results in the location method of interlaced grid search, and the location results are stable and reliable. In this study, the earthquake catalog of CCZ and the events distribution in South Shetland Islands and its surrounding areas are obtained. During this period, the comparison between the CCZ catalog and the ANSS catalog in terms of magnitude, focal depth, and the number of located events suggests that CCZ effectively improves the seismic monitoring capability of the area, and that the use of direct P wave initial amplitude for calculating the angle of incidence better constrains the focal depth of local seismic events. The distribution of events shows that the events mainly spread along the NE direction in the South Shetland Islands region, with the maximum magnitude of ML4.8, the minimum magnitude of ML1.3, and the focal depth ranging from 2 to 134 km. The seismicity in the Antarctic Great Wall Station area is stronger than that of the stable interior of the Antarctic continent. The spreading of the Bransfield Strait was accompanied by extensive volcanism, and this study identifies clusters of earthquakes on Deception Island and its eastern flanks, as well as in the southern part of Greenwich Island, that may be related to the volcanic magmatism caused by subduction. Eleven events with focal depths exceeding 100 km, which have occurred beneath the lithosphere, were detected in the forearc region close to the South Shetland Trench. The locations of these events suggest that the subduction may have occurred in this region, which provides direct seismological evidence that the subduction of the Phoenix Plate is still continuing.

  • 地下水作为地壳结构中最活跃的组分之一,具有较好的流动性和难压缩性。当在含水层周围形成一个封闭条件较好的承压系统时,地下水的升降就能十分客观地反映出周边区域地壳应力应变的状态。同时自1988年至今所出版的 《中国震例》 报告(张肇成等,19881990ab19992000陈棋福等,2002ab20032008蒋海昆等,2014)中的统计结果显示,中强地震前普遍存在地下水位异常,因此,我们可以通过记录到的地下水位动态序列来研究周边含水岩体在地震发生之前的受力、变形和破坏过程(汪成民等,1988付虹等,2014)。

    正常观测得到的地下水位动态变化是在多种因素影响下产生的综合物理量,既包括水位的宏观动态变化,也包括微观动态变化,其中宏观动态变化主要是由含水层中储水量的改变引起的变化,微观动态变化则主要是由周边岩体的应力应变状态改变而引起的变化,因此,地下水位的异常变化形态有可能是周边环境干扰等引起的干扰变化,也可能是区域构造活动引起的前兆异常变化(贾化周等,1995孙小龙等,2013)。如何排除这些干扰因素的影响,识别出与构造活动有关的地下水微动态,进而判断地下水异常产生的机理,是目前地下流体资料分析应用的关键科学问题之一(车用太等,2011)。

    一般情况下,水位的宏观动态变化主要与观测井周边的地下水补给、径流和排泄状态密切相关,在日常研究中最为常见的因素包括地下水开采和大气降水等(车用太等,1993晏锐,2008杨明波等,2009胡小静等,2016),经过多年的实践,研究人员提出并尝试了不同的分析方法来排除降水对水位数据变化的影响,包括卷积滤波法(张昭栋等,1993)和组合水箱模型法(王旭升等,2010孙小龙等,2013)等,在实际应用中均取得了一定的实用效果;对于引起水位微观动态变化的构造活动,除了用各类地球物理场观测资料来描述地壳形变和构造活动之外,许多研究人员也尝试在不同区域,利用数字化观测的承压井水位动态资料,通过计算含水层的部分介质参数,来反演该区域内的构造活动和应力应变状态,从而分析探寻与地震有关的部分前兆信息,并积累了一定的经验和成果(李永善,1979孙小龙等,2011丁风和等,2015b)。

    滇南地区高大、开远、易门3井的井水位自2014年以来出现的趋势转折上升现象到底是环境干扰所致还是由构造活动所引起的前兆异常信息,针对此问题,笔者拟从降水影响定量排除、观测环境调查、井孔含水层和区域应力场状态以及井孔水体化学组分等方面进行分析,探讨目前水位趋势上升异常的发生机理。

    高大井、开远井和易门井在大的构造区域上,集中位于川滇菱形块体的东南角。从活动断裂带来看,高大井位于通海—峨山断裂带上,地处1970年通海MS7.8大震的极震区;开远井位于个旧断裂和开远—蒙自断裂附近,地处小江断裂与建水—石屏弧形断裂交叉、复合部位;易门井位于南华—楚雄—建水断裂带上(图1)。井孔柱状图结果显示(云南省地震局,2005),高大井目前实际井深约为201.4 m,井孔观测的含水层位于166—194 m处,岩性属震旦系微红组石英岩,地下水类型为基岩裂隙承压水,岩芯分析结果显示整个含水层处于受控制的强挤压断裂破碎带内;开远井深224 m,井孔观测的含水层有两段,分别位于96—100.36 m和169.38—224 m处,岩性属三叠系的灰岩,地下水类型为岩溶裂隙承压水,岩心结构显示50—90 m左右存在断层碎屑物,可能是断裂穿过井孔所致;易门井井深约为155 m,井孔观测的含水层分别位于96.31—102.47 m和122.64—159.58 m处,岩性属昆阳群板岩,地下水类型为基岩裂隙承压水。

    图  1  观测井孔周边地质构造图
    Figure  1.  Geological structure around the observation wells

    高大井和开远井的井水位自观测以来,呈现出雨季上升、旱季下降的年变形态,含水层补给类型属降雨渗透补给型,2014年以来高大井和开远井的井水位出现趋势转折上升现象,截至2017年8月,水位持续上升幅值已分别达到1.5 m和4.3 m,目前水位保持在高位状态(图2)。

    图  2  高大井 (a)、开远井 (b) 和易门井 (c) 月均水位和月累积降水量的时间进程曲线
    Figure  2.  The time history curve of monthly mean water levels and cumulative precipitations of Gaoda well (a),Kaiyuan well (b) and Yimen well (c)

    从数据分析的结果来看,引起水位趋势上升最常见的两类原因为:一是周边观测环境改变引起的干扰变化,包括本区域降水量逐年增多、观测井周边地下水开采量的减少等引起观测井含水层内水量的增多,从而引起观测井水位的上升;二是区域应力场改变构造活动导致地壳的岩体产生变形,使得含水层固体骨架之间的孔隙压增大,从而导致观测井水位出现持续上升现象。

    从水位与降雨量的同轴曲线(图2)可知,3口井的水位在大多数年份均表现为雨季上升、旱季下降的年变形态,表明井水位与降水之间的关系较为密切,但从图中也可以看出,水位变化在个别年份也会表现出与雨季、旱季不一致的变化形态,例如:高大井水位在2015,2016年间出现过旱季水位持续上升的现象,开远井水位在2014,2016年间出现过旱季水位持续上升的现象,而易门井水位则在2013年连续3个月的雨季期间持续下降,表明这3眼井近几年的水位变化除了受降水影响以外,还存在其它的影响因素。

    为了进一步分析降水与水位变化之间的相关性,利用月累积降雨量与月均水位变化量作线性回归,拟合出二者之间的线性关系,yax±b。其中,y为水位观测的月均变化量,x为水位观测的月累积降雨量,当水位月变幅值超出拟合线一定范围时,则认为水位变化与降水之间的相关性较小。从月累积降雨量与月均水位变化量的线性回归拟合可以看出,这3眼井的水位与降水呈现了较好的线性相关,与云南地区部分观测井水位对降雨响应的滞后时间长达4个月相比,这3口井的响应过程相对更快,但仍然存在数天的滞后现象,不过依然保持着雨季上升、旱季下降的年变形态。但高大井水位在2007—2008年间和2012—2016年间均出现了水位变化量与降雨影响幅值明显不吻合的现象;开远井水位在2012,2014—2016年间出现了水位变化量与降雨影响幅值不吻合的现象;易门井水位在2007,2012—2016年间出现了水位变化量与降雨影响幅值不吻合的现象(图3)。

    图  3  高大井 (a)、开远井 (b) 和易门井 (c) 的月均水位变化量与降雨量线性回归示意图
    Figure  3.  Linear regression between monthly mean water level difference and rainfall of Gaoda well (a),Kaiyuan well (b) and Yimen well (c)

    王旭升等(2010)针对水位变化对降水响应滞后的观测井孔,提出了降雨-水位动态的组合水箱模型,其可通过当地降水量资料来反演地下水位的动态变化过程。该模型主要通过伽玛分布密度函数所建立的单位脉冲响应函数来处理降水补给滞后的影响,并将地下水位观测得到的实际变化形态与反演计算得到的模拟变化动态进行对比,从而来定量地剔除降水量对水位动态变化的影响。孙小龙等(2013)利用该模型,对姚安井水位资料的降雨干扰进行了定量剔除,得到了比较理想的反演结果。

    笔者也尝试利用该模型,根据当地的月累积降雨量变化对高大、开远和易门3口井的水位变化进行反演,得到3口井水位的模拟值,并与去趋势之后的原始水位观测值进行对比(图4),从图4可以看出,易门井水位年动态受降水影响非常显著,而高大和开远井的水位年动态明显偏离模拟值,因此,仅从降水量补给角度看,易门井水位年动态主要受控于降水量,降水增多引起其趋势上升,而高大井和开远井的水位年动态除了受降水影响之外,在很大程度上还受其它因素的控制。

    图  4  高大井 (a)、开远井 (b) 和易门井 (c) 水位的观测值与模拟值的动态变化曲线
    Figure  4.  Dynamic changes of observation and simulation water levels in Gaoda well (a),Kaiyuan well (b) and Yimen well (c)

    自2009年以来,云南地区干旱严重,滇南地区作为旱区之一,新打机井几百眼,主要用于农业灌溉等,其中也存在部分深井。早在2013年水位处于下降过程时,付虹等(2014)便对开远井周边的抗旱打井抽水情况进行过调查,并对距离开远井最近的(约5 km)1口深井的抽水情况进行了定量估算,结果显示该深井抽水的影响范围远小于5 km。笔者于2017年再次对这3眼井周边地下水的开采情况进行了详细调查,调查结果显示:高大井和开远井所在坝区在2009年左右的干旱年份均未曾开挖抗旱井,已有的居民用井开挖深度均在10 m左右,与观测井含水层不在同一层位,用水量多年基本保持不变,因此,由高大井和开远井附近含水层内水量在2014年前后并未增多可知,这两口井附近的地下水开采并未减少。

    易门井周边地下水开采情况调查结果显示,在观测井东北方约3 km处为大椿树工业园区,园区内建有多家水泥厂和家具建材生产、销售基地,共打水井21口,井深大多在150 m左右,长期使用水泵进行抽水,直至2016年10月由当地政府将21口抽水井统一封存。对易门井周边大椿树工业园区内深井抽水情况进行具体分析可知:① 从工业园区内深井开挖的进程来看,在2002—2003年间,随着工业园区内开挖的深井明显增多,水位出现明显的趋势下降,2009—2010年间,随着深井开挖的另一个高峰,水位趋势下降的速率更为明显,至2012年县政府限制深井的开挖和抽水之后,水位出现短暂的回升,趋势下降过程有所减缓,2014年之后完全停止新井的开挖,水位出现了明显的转折回升过程,至2016年10月,工业园区内抽水深井被全部封存,水位保持持续回升的过程。整个工业园区内深井开挖和抽水的发展进程与易门井水位趋势下降后转折回升的过程在时间上比较吻合(图5);② 从抽水层与观测井含水层的对比来看,工业园区附近海拔为1 610 m,而井深大多在100—150 m左右,因此大量抽水层均位于1 460—1 510 m左右,观测井海拔为1 580 m,基础资料显示,该井含水层位于井孔下96.31—102.47 m,因此其含水层位于海拔1 483—1 478 m处,对比抽水层与观测井含水层,抽水层恰好位于观测井的含水层处。结合深井开挖进程和抽水层位进行分析认为,易门井水位变化动态受大椿树工业园区内深井抽水的影响较为明显。

    图  5  易门井水位时序图
    Figure  5.  Daily mean water levels of Yimen well

    从上述关于降水影响和地下水开采影响的分析结果可以看出,易门井井水位年动态主要受控于降水过程,整个趋势下降变化与大椿树工业园区的深井开挖及抽水关系非常密切;高大井和开远井井水位的年动态和趋势上升变化不完全受控于降水影响,与周边抗旱井地下水开采使用亦无直接联系。

    根据已有的研究结果,水位上升过程除了受降水和地下水开采的影响之外,很有可能还与周边的构造应力状态有关,因此,笔者进一步定量计算了3口井井孔含水层系统的垂直向应力场的动态变化过程,并收集了滇南地区的GPS观测结果,对异常井井孔周边的构造活动和应力场状态进行分析。

    根据前人的研究成果(李永善,1979孙小龙等,2011丁风和等,2015b),在理想的承压含水层模式下,井孔含水层垂直向的应力变化量Δσz与含水层部分介质参数和井水位变化量之间存在定量关系,并且能在一定程度上反映出该区域的构造活动特征,即

    $\Delta {\sigma _{{\textit z}}} {\text{=}} \frac{{ {\text{-}}2\beta \rho {\rm{g}}}}{{n\beta {\text{+}} {\displaystyle\frac{1{\text{-}} n} {E}}}}\Delta H{\text{,}}$

    (1)

    式中:Δσz为含水层垂直向应力变化量;ρg为水的重度,在研究过程中取9.8×103 N/m3;ΔH为剔除长趋势和年周期变化后,含水层应力变化引起的压力水头变化量,即井水位变化量;n为含水层的孔隙度;β为水的体积压缩系数;E为含水层固体骨架的弹性模量,E=1/αα为固体骨架的体积压缩系数。

    上述参数中,nβα即为含水层的部分介质参数,它们之间的关系表达式可通过气压系数和潮汐因子之间的关系推算得到(丁风和等,2015ac)。具体的推算和拟合过程如下:在仪器设备和观测环境正常的情况下,井水位的潮汐因子Bg和气压系数BP可分别表示为

    ${B_{\rm{g}}} {\text{=}} {\text{-}} \frac{{1 {\text{-}} n}}{{\rho {\rm g}\left[ {(1 {\text{-}} n)\alpha {\text{+}} n\beta } \right]}}{\text{,}}{B_{\rm P}} {\text{=}} \frac{{n\beta }}{{\alpha {\text{+}} n\beta }}{\text{,}}$

    (2)

    式中:潮汐因子Bg选择地球固体潮中振幅最大、信噪比最高的M2波作为计算对象,通过维尼迪柯夫(Venedikov)潮汐调和分析获取;气压系数BP主要利用水位数据和同井观测的气压数据经过滤波和差分,再通过线性回归模型进行拟合后获取;最后经过滑动推算得到含水层孔隙度n与水的体积压缩系数β之间的关系,以及固体骨架的体积压缩系数α

    从物理意义及异常性质判定来讲,当井孔含水层系统所受应力增强,即Δσz>0时,井水位上升,水位埋深变小,其变化量ΔH<0;当井孔含水层系统所受应力减弱,即Δσz<0时,井水位下降,水位埋深变大,其变化量ΔH>0。

    笔者利用文中3眼水位观测井的气压系数和M2波潮汐因子等参数,反演得到井孔含水层部分介质参数,在水平层状含水层模式下,定量地计算了2007年以来高大、易门和开远3口观测井的垂直向应力场的动态变化过程(图6)。计算结果显示,自2014年以来,高大井、开远井、易门井的含水层系统的应力状态均以压应力为主,表现为水位上升、水位埋深变小,表明该区域内近几年以持续的应力增强活动为主。从历史震例可以看出(表1),2007—2008年间水位观测井孔的含水层出现了应力增强之后,在云南和川滇交界地区先后发生了2008年会理MS6.1地震和2009年姚安MS6.0地震;2012年再次出现应力增强之后,2013年发生了四川芦山MS7.0地震;2014年以来在应力的持续增强过程中,云南地区先后发生了鲁甸MS6.5和景谷MS6.6等多次6级地震,2017年井孔垂直向应力仍保持增强状态。历史震例统计情况表明,井孔含水层应力增高,确实可以在一定程度上反映出区域范围内构造应力增强的过程,同时也有利于地震的发生。

    表  1  井孔应力增强与周边地震的对应情况统计
    Table  1.  Statistical analysis of the correspondence between the stress enhancement of observation wells and surrounding earthquakes
    序号 井孔名称 异常开始时间 异常结束时间 异常幅值/105 Pa 发震日期 震例
    年−月−日
    1 高大 2007−06−03 宁洱MS6.4
    开远 2007−03 2007−05 0.036 1
    易门 2007−05 2008−04 0.045 7
    2 高大 2008−08−30 会理MS6.1
    开远 2008−02 2009−01 0.055 0
    易门 2007−05 2008−04 0.045 7
    3 高大 2009−05 2009−12 0.001 96 2009−07−09 姚安MS6.0
    开远 2007−03 2009−01 0.055 0
    易门
    4 高大 2012−01 2012−12 0.012 9 2013−04−20 芦山MS7.0
    开远 2012−06 2012−12 0.005 4
    易门 2012−03 2012−09 0.003 0
    5 高大 2014−01 2017−05 (未结束) 0.087 4
    开远 2013−10 2017−01 0.125 0 2014−08−25
    2014−10−07
    鲁甸MS6.5
    景谷MS6.6
    易门 2014−03 2017−05 (未结束) 0.049 0
    注:截至结算结果日期2017年5月,高大井和易门井的高应力状态仍未结束。
    下载: 导出CSV 
    | 显示表格
    图  6  高大井 (a)、开远井 (b) 和易门井 (c) 含水层垂直向应力变化量Δσz的动态变化
    Figure  6.  Dynamic changes of vertical stress Δσz in aquifer of Gaoda well (a),Kaiyuan well (b) and Yimen well (c)

    陈立德和付虹(2003)的研究结果表明,前兆异常均为区域应力水平增强的产物。当区域应力水平出现改变,构造活动增强时,地下水位就可能出现趋势性的异常变化。根据同期云南省地震局形变测量中心的GPS解算结果,布设于滇中—滇南块体中自北向南的昆明—通海基线和通海—蒙自基线自2014年以来呈现出明显的压缩状态(图7),表明整个块体的南北向挤压应变积累较强,以压性活动特征为主,目前的异常井孔恰好位于两条基线观测区间内(图1)。根据已有的研究成果(刘耀炜等,2010付虹等,2014),如果岩层被挤压,则岩层的孔隙率变小,孔隙水压增大,水流由含水层向井孔流动,从而使水位上升。据此分析认为,高大井和开远井井水位自2014年以来的趋势性上升过程可能与整个块体南北向强挤压为主的构造活动存在一定的相关性。

    图  7  昆明—通海 (a) 和通海—蒙自 (b) 的GNSS基线观测时间序列
    Figure  7.  Time series of Kunming-Tonghai (a) and Tonghai-Mengzi (b) of GNSS baseline observation

    为了进一步判定目前观测井水体中的水-岩相互作用程度,并希望以此来判断观测井水体中是否携带来源于深部构造活动的物质成分,笔者于2016年12月对3眼观测井进行了取样,取样深度均为水面以下1 m,并委托云南省地震局防灾研究所实验室和中国地震局防灾科技学院流体实验室分别对水样的水化学成份和氢氧同位素进行了测试,化验结果分析如下。

    从水化学三线图(图8a)可见,3口井的地下水类型均属于重碳酸型水,反映了浅层、年龄较轻的水文循环特征,但相比于易门井,高大井和开远井的Na++K+的相对含量明显高于易门井,而Ca2,Mg2+的相对浓度则远低于易门井。根据离子之间的反应过程可知,当观测井中水体与含水层中的岩体之间发生深度接触时,水体中的Ca2+,Mg2+会发生一系列被吸附和交换的过程,从而释放出更多的Na+,使得含水层水体中Ca2+,Mg2+含量明显减少,而Na+含量则显著升高。因此,水化学Piper图的结果表明,高大井和开远井的含水层水体与岩体之间的相互作用比较充分,二者之间的平衡程度也较高;而易门井则是以典型的大气水与岩石之间的第一阶段反应为主。

    图  8  地下水水样测试结果分析
    (a) 地下水Piper图;(b) Na-K-Mg三角图;(c) 氢氧同位素对比图
    Figure  8.  Analysis of groundwater test results
    (a) The diagram of groundwater Piper;(b) The triangle diagram of Na-K-Mg; (c) Hydrogen and oxygen isotopes of observation wells

    从水样的Na-K-Mg三角图(图8b)可见,高大井的地下水比较接近“部分平衡水”状态,其水-岩相互作用程度相对较高,携带了部分深部构造活动的信息;易门井的地下水则位于Mg2+端元附近的“未成熟水区域”,表明该井的含水层水体与岩体之间的相互作用程度非常弱,二者之间仍在发生着相互溶解和交换的过程,所携带的深部信息较少;开远井虽然也处于“未平衡水状态”,但相比易门井,已有所偏离Mg2+端元,因此也存在一定的水-岩相互作用。

    从氢氧同位素对比图(图8c)来看,3口井均位于西南地区大气降水线右侧,表现为δ18O富集特征(正向漂移)。对于该现象,根据现有的研究结果,推测其可能的形成原因主要有两种:一是由于该区域正处在降水量小而蒸发作用较强的时段,这在干旱、半干旱地区表现的较为突出;二是由于地下水体与岩体之间长期以来发生了较为充分的相互作用,使得水体中的氢氧同位素之间发生了平衡交换。考虑到云南地区近几年雨水较为丰富,蒸发过程偏弱,因此分析认为氢氧同位素正向漂移现象很有可能是深部构造活动引起水岩反应的结果。

    通过对滇南地区高大井、开远井和易门井3眼水位观测井自2014年以来出现的水位趋势转折上升现象进行观测环境、大气降水、构造活动影响以及水化学组分等方面的分析,得到如下认识:

    1) 自2014年以来高大井、开远井和易门井井水位出现趋势上升,对水位变化与降水的相关性分析显示,高大井、开远井和易门井在短期内曾多次出现水位变化明显偏离降雨影响范围的上升过程,表现为与同期降水不完全相关的特征;利用降水-水位动态的组合水箱模型进行了降水影响定量排除,结果显示,易门井井水位年动态主要受控于降水量,而高大和开远井的水位年动态除了受降水影响之外,其它因素对其影响更大。

    2) 周边观测环境调查结果的相关性分析显示,易门井井水位的整体趋势下降与转折上升过程与易门大椿树工业园区的深井开挖抽水密切相关,高大和开远井井水位的年动态和趋势上升变化不受周边抗旱井及地下水开采使用等环境的影响。

    3) 在水平层状含水层模式下,3眼观测井的垂直向应力场的反演结果显示,自2014年以来,高大井、开远井和易门井的含水层系统的应力状态均表现为应力持续增强。历史震例表明,井孔含水层应力增高,确实可以在一定程度上反映出区域范围内构造应力积累的过程。GNSS基线观测结果同样表明,自2014年以来,异常井孔所处的整个滇中—滇南块体自北向南以明显的压性活动为主。

    4) 水化学组分和氢氧同位素分析结果表明,高大井和开远井地下水的水-岩相互作用程度较高,水体中存在一定的来源于深部构造活动反应的物质成分;易门井地下水的水-岩相互作用比较弱,基本以浅层的溶解作用为主。

    5) 结合上述几个方面的分析,笔者认为高大井、开远井的井水位自2014年以来的趋势转折上升现象不完全为降水影响所致,从水体的循环深度较深且携带有深部构造活动的物质成份以及GPS显示井孔所处区域地块的挤压现象等初步判定,这两眼井的水位趋势上升与区域构造活动有关,因此需作为异常继续观察;而易门井虽然曾经出现过的高应力状态也与地震有关,但因水位年动态与降水影响完全一致,趋势下降和转折上升与大椿树工业园区的深井开挖抽水关系也非常密切,因此认为其与深部的构造活动关系不大。

    本文在撰写过程中,得到了宁夏回族自治区地震局丁风和副研究员、中国地震局地壳应力研究所孙小龙博士和中国地震台网中心周志华博士的指导和帮助,作者在此表示衷心的感谢!

  • 图  1   研究区构造背景

    图中绿色三角形是中国南极长城站(CCZ)地震台的位置,地震数据来自美国国家地质调查局地震信息中心提供的1990—2017年的美国国家地震监测台网系统(Advanced National Seismic System,缩写为ANSS)综合地震目录(USGS,Earthquake Hazards Program,2017),震源机制解来自全球质心矩张量项目提供的1990—2017年的全球质心矩张量(Centroid-moment tensors,缩写为CMT)目录(Erik et al,2013),断层数据来自德克萨斯大学地球物理研究所的UTIG数据库(Institute of Geophysics,University of Texas,2015)和Dziak等(2010)。右上角插图为斯科舍板块及其周边的地震与台站分布,其中蓝色三角形是全球地震台网(GSN)的台站,蓝色方框代表主图的位置

    Figure  1.   Tectonic setting of studied area

    The green triangle is the location of Chinese Great Wall Seismic Station (CCZ),Antarctic. The seismic data is from the ANSS Comprehensive Earthquake Catalog between 1990 and 2017,provided by USGS National Earthquake Information Center (USGS,Earthquake Hazards Program,2017). The focal mechanism solution is from CMT catalog between 1990 and 2017 provided by the Global Centroid-Moment-Tensor Project (Erik et al,2013). The fault data is from the University of Texas Institute for Geophysics database (Institute of Geophysics,University of Texas,2015) and Dziak et al2010). The insert in the upper right is the distribution of earthquakes and seismic stations in and around the Scotia Plate,the blue triangles are Global Seismographic Network Stations,and the blue rectangle represents the location of the main figure

    图  2   中国南极长城站地震台的近震波形记录

    (a) 可定位地震事件的三分量波形记录;(b) 可定位地震事件的P波初动信号三分量记录;(c) 不可定位地震事件的三分量波形记录;(d) 不可定位地震事件的P波初动信号三分量记录

    Figure  2.   Near earthquake waveforms recorded at CCZ Station

    (a) Three-component waveform recordings of a locatable seismic event;(b) Three-component recordings of P wave initial motion signal of a locatable seismic event;(c) Three-component waveform recordings of an unlocatable seismic event; (d) Three-component recordings of P wave initial motion signal of an unlocatable seismic event

    图  3   交错网格的构建

    (a) 经度-纬度网格的平移过程;(b) 震中距-深度网格的平移过程黑色网格表示起始网格位置,红色、蓝色和绿色网格分别表示每一次移动后的网格位置,三角形表示台站。此示意图按小格尺寸的1/2为步长进行平移

    Figure  3.   The construction of interlaced gird

    (a) The translation process of the longitude-latitude grid;(b) The translation process of the epicentral distance-depth grid The black grid represents initial grid position,the red,blue and green grids represent the grid position after each move respectively. The triangle represents seismic station. This diagram is translated by 1/2 of the small lattice size per step

    图  4   长城站地区2015—2017年地震分布情况

    左上角插图给出了2015—2017年ANSS目录中在CCZ台500 km范围内的地震分布,蓝色方框为主图区域

    Figure  4.   Earthquake distribution in the Great Wall Station area from 2015 to 2017

    The insert at the upper left corner shows the distribution of earthquakes within 500 km from the Great Wall Seismic Station in the ANSS catalog from 2015 to 2017. The blue rectangle is the area of the main figure

    图  5   ANSS目录中2016−05−30 15:35:09.08发生的地震在CCZ台的记录

    Figure  5.   CCZ waveform record,which is in the ANSS catalog at 15:35:09.08 on May 30,2016

    图  6   2015—2017年长城站地区地震频次

    Figure  6.   Earthquake frequency in the Great Wall Station area during 2015−2017

    图  7   2015—2017年长城站地区可定位地震震级随时间变化的情况

    Figure  7.   The change of locatable earthquake magnitude with time in the Great Wall Station area during 2015−2017

    图  8   测线AA′ (a)和BB′ (b)的震源深度分布剖面图

    图a,b上部为AA′和BB’测线纵截面测线位置的地形高程变化,剖面上投影了AA′测线南北20 km范围内地震和BB′测线所有可定位地震,横轴表示地震到测线的垂线与CCZ台的距离,以台站为原点,向西(北)为负,向东(南)为正

    Figure  8.   Focal depths distribution along profile of survey lines AA′ (a) and BB′ (b)

    The upper part of the figure is the topographic elevation change along the longitudinal section of the survey lines. The earthquakes within 20 km from the north to the south of AA′ survey line and all locatable earthquakes of BB′ survey line are projected on the profile. The horizontal axis represents the distance between the vertical line of the earthquake to the survey line and the CCZ Station,with the station as the origin,negative to the west (north) and positive to the east (south)

    表  1   中国南极长城站地区速度模型

    Table  1   The velocity model of the Chinese Antarctic Great Wall Station area

    深度/km vP/(km·s−1 vS/(km·s−1
    0 4.10 2.20
    2 5.70 3.10
    9 6.65 3.59
    20 7.28 3.94
    30 8.05 4.43
    50 8.30 4.52
    下载: 导出CSV

    表  2   ANSS目录震级和本方法计算结果的对比

    Table  2   The magnitude comparison of ANSS catalog and calculation result of this method

    发震时刻(UTC) 南纬/° 西经/° 深度/km Mb ML (本文震级)
    年−月−日 时:分:秒
    2 016−05−30 15:35:09.080 61.65 54.70 10.0 4.6 4.9
    2 016−07−02 20:02:00.060 61.35 53.93 10.0 4.6 4.9
    2 016−07−16 17:30:40.090 60.36 57.56 16.5 4.3 4.7
    2 016−10−17 00:13:56.420 63.29 61.91 10.0 4.9 4.9
    下载: 导出CSV

    表  3   2015—2017年CCZ台单台定位结果

    Table  3   Single station location results of CCZ Station during 2015−2017

    发震时刻(UTC) 南纬/° 西经/° 深度/km ML 震中距标准差/km 深度标准差/km
    年−月−日 时:分:秒
    2 015−01−09 02:52:57.142 62.32 58.33 47.0 2.1 0.6 0.8
    2 015−01−11 18:55:35.171 62.65 59.49 46.6 3.2 0.3 0.5
    2 015−01−13 02:26:02.496 62.41 58.36 44.3 2.3 0.4 0.6
    2 015−01−17 15:08:20.536 62.69 59.48 44.9 3.0 0.4 0.5
    2 015−01−19 18:21:52.508 62.25 59.30 43.6 2.0 1.0 0.3
    2 015−01−26 23:24:17.910 62.67 59.56 36.9 3.4 0.3 0.6
    2 015−02−01 19:15:40.762 62.69 59.29 2.1 3.5 0.9 1.4
    2 015−02−01 21:56:40.560 62.67 59.54 15.7 3.6 0.4 2.2
    2 015−02−06 00:16:49.919 62.33 58.98 14.8 2.1 2.4 1.1
    2 015−02−22 05:07:40.639 62.33 59.64 45.0 3.6 0.3 0.5
    2 015−02−12 01:23:56.584 62.67 59.66 27.3 3.9 0.3 1.4
    2 015−02−14 01:52:26.810 62.35 59.29 42.7 1.7 0.7 0.4
    2 015−02−16 04:36:05.947 62.61 59.45 54.1 2.7 0.4 0.5
    2 015−02−16 09:37:06.332 62.67 59.40 47.1 2.9 0.4 0.5
    2 015−03−01 19:38:22.289 62.26 59.44 70.4 3.4 0.8 0.3
    2 015−03−08 21:29:07.408 62.64 59.40 48.4 4.5 0.4 0.5
    2 015−03−09 17:03:42.891 62.64 59.42 49.1 4.2 0.4 0.5
    2 015−03−10 12:17:11.868 62.69 59.53 31.5 3.1 0.3 1.1
    2 015−03−27 08:49:32.449 62.67 59.38 41.6 3.7 0.4 0.5
    2 015−03−27 11:50:30.010 62.30 59.50 77.1 1.9 0.6 0.3
    2 015−04−05 10:38:31.014 62.60 59.61 48.1 3.0 0.4 0.5
    2 015−04−07 15:01:46.761 62.69 59.41 43.9 3.6 0.4 0.6
    2 015−04−11 07:13:53.799 62.05 59.10 62.9 1.6 0.8 0.3
    2 015−04−14 17:24:41.493 62.35 59.41 49.9 2.8 0.6 0.3
    2 015−05−09 05:12:09.465 62.45 59.30 48.9 3.0 1.0 0.8
    2 015−05−31 11:24:33.141 62.68 59.55 21.4 3.8 0.4 1.5
    2 015−07−04 09:14:21.220 62.87 60.16 75.9 3.4 0.3 0.5
    2 015−07−05 04:30:40.211 62.83 59.92 91.6 3.6 0.4 0.5
    2 015−07−12 08:00:12.972 62.33 60.07 134.1 3.7 0.6 0.3
    2 015−07−12 23:11:11.567 62.51 59.98 70.2 3.0 0.4 0.6
    2 015−07−20 03:05:51.509 62.64 59.54 43.4 3.8 0.4 0.5
    2 015−07−20 05:59:45.273 62.64 59.43 48.4 4.8 0.4 0.5
    2 015−07−20 07:04:04.974 62.04 58.50 47.6 2.5 0.9 0.6
    2 015−07−20 11:23:14.377 62.68 59.32 48.1 3.5 0.4 0.5
    2 015−07−22 11:32:19.363 62.00 59.36 105.8 1.7 0.6 0.3
    2 015−07−24 22:09:27.906 62.42 59.51 58.6 2.7 0.6 0.4
    2 015−07−27 18:10:25.901 62.28 59.19 26.7 2.2 0.9 0.4
    2 015−07−31 18:27:09.665 62.90 60.08 98.9 4.0 0.4 0.5
    2 015−08−02 05:54:17.571 61.98 59.22 117.9 1.3 0.7 0.3
    2 015−08−02 17:53:21.228 62.12 59.65 76.8 4.5 0.5 0.3
    2 015−08−13 17 07:03.185 61.74 58.54 69.3 3.5 0.4 0.5
    2 015−08−23 18:01:25.544 62.65 59.51 43.6 3.3 0.4 0.5
    2 015−08−24 05:33:26.643 62.66 59.60 34.9 3.1 0.3 0.8
    2 015−08−26 07:21:16.635 61.47 58.73 121.4 3.5 0.9 0.6
    2 015−08−29 21:02:32.865 62.01 58.98 58.3 1.9 2.0 0.3
    2 015−09−04 01:00:05.695 61.72 58.99 100.9 2.7 0.6 0.4
    2 015−09−06 03:22:49.571 62.28 59.33 37.1 2.3 0.7 0.3
    2 015−09−09 15:43:45.491 61.86 58.37 78.6 3.6 0.5 0.3
    2 015−09−12 08:08:14.483 62.28 59.29 59.2 2.7 0.8 0.3
    2 015−09−15 21:15:24.730 62.32 58.89 14.0 1.8 1.0 0.7
    2 015−09−19 17:37:15.264 62.59 59.69 55.8 3.7 0.3 0.5
    2 015−09−26 21:19:39.555 61.99 58.86 59.0 2.8 0.7 0.3
    2 015−09−27 00:13:48.033 62.13 59.42 85.1 2.1 0.7 0.3
    2 015−10−01 07:15:24.196 62.63 59.59 67.7 4.0 0.4 0.6
    2 015−10−06 07:04:50.301 61.75 58.16 83.8 3.1 0.6 0.8
    2 015−10−11 00:24:17.187 62.23 59.15 56.2 2.2 1.8 0.3
    2 015−10−12 00:09:00.274 61.94 58.25 62.3 2.7 0.7 0.9
    2 015−10−30 09:02:58.971 61.93 59.40 117.7 2.6 0.6 0.3
    2 015−11−05 10:50:09.692 62.56 59.61 76.2 3.9 0.8 0.5
    2 015−12−17 20:10:02.746 62.67 59.86 69.4 3.3 0.3 0.5
    2 016−03−05 17:48:06.327 62.28 59.53 55.9 2.5 0.6 0.3
    2 016−04−10 21:07:44.857 62.57 59.55 65.9 3.8 0.4 0.5
    2 016−05−16 15:46:42.359 62.28 59.41 59.9 3.4 0.7 0.3
    2 016−06−07 09:39:09.115 61.77 58.37 51.6 3.6 0.3 0.5
    2 016−06−08 16:44:36.534 62.09 59.82 80.8 2.8 0.6 0.4
    2 016−06−19 00:34:16.557 62.47 59.57 62.1 2.8 0.4 0.3
    2 016−06−21 00:53:16.317 62.65 59.81 71.7 4.3 0.5 0.6
    2 016−07−08 06:59:53.855 62.56 59.53 61.6 2.7 0.4 0.4
    2 016−07−15 12:52:28.546 62.56 58.33 40.1 4.4 0.3 0.8
    2 016−08−11 22:37:06.427 62.27 59.30 57.6 2.2 0.9 0.3
    2 016−08−19 13:08:28.771 62.41 59.17 19.6 1.3 0.6 0.6
    2 016−08−20 04:48:49.164 62.31 59.61 61.1 2.1 0.6 0.4
    2 016−08−22 10:08:51.241 62.00 58.29 61.3 2.4 0.4 0.3
    2 016−08−26 16:11:27.861 62.27 59.76 50.9 1.9 0.5 0.5
    2 016−09−07 11:30:14.309 61.94 59.75 121.3 2.8 0.6 0.3
    2 016−09−16 18:31:24.635 62.51 59.60 43.6 2.9 0.4 0.5
    2 016−11−01 23:20:23.770 62.34 59.85 88.4 3.7 0.6 0.4
    2 016−11−03 09:46:21.038 62.67 59.84 84.5 4.0 0.3 0.4
    2 016−11−11 04:42:38.088 62.36 59.22 33.5 1.8 0.8 0.4
    2 016−11−11 20:06:48.364 62.40 59.61 58.8 2.8 0.5 0.3
    2 017−01−16 00:38:56.890 62.60 59.48 54.1 4.0 0.4 0.5
    2 017−01−16 01:18:34.637 62.62 59.63 44.6 4.0 0.4 0.6
    2 017−01−16 01:19:41.900 62.58 59.55 53.9 3.7 0.4 0.5
    2 017−01−24 16:03:04.799 62.13 59.37 47.6 1.8 0.8 0.3
    2 017−01−25 17:45:15.526 62.25 59.25 32.4 2.2 1.0 0.3
    2 017−02−15 04:07:15.207 62.13 59.35 62.1 3.1 0.7 0.3
    2 017−02−24 00:32:35.366 61.93 58.65 64.9 3.0 0.6 0.4
    2 017−03−01 04:29:05.009 62.54 59.80 47.9 3.5 0.4 0.5
    2 017−03−02 08:12:21.446 61.97 59.93 102.3 3.0 0.5 0.3
    2 017−03−06 17:54:49.262 62.02 58.41 35.4 2.4 0.4 0.5
    2 017−03−07 10:27:33.329 62.33 59.89 74.6 2.7 0.5 0.3
    2 017−03−09 16:07:42.708 62.64 59.73 70.6 2.6 0.5 0.6
    2 017−03−15 17:59:53.782 62.41 59.67 42.9 2.5 0.4 0.5
    2 017−03−31 10:05:30.533 61.78 59.78 132.6 4.1 0.5 0.3
    2 017−04−12 08:25:59.180 61.66 58.74 59.7 4.5 0.4 0.5
    2 017−05−24 09:54:32.085 62.18 58.81 26.4 1.6 1.5 0.3
    2 017−05−30 13:34:00.388 62.27 59.75 76.7 3.1 0.5 0.3
    2 017−05−31 20:40:08.639 61.98 58.65 84.6 3.6 0.6 0.3
    2 017−07−25 05:34:31.746 62.67 59.82 63.0 3.9 0.3 0.5
    2 017−08−09 11:22:15.439 63.00 60.83 77.4 3.8 0.4 0.5
    2 017−08−23 03:36:11.790 62.38 59.74 81.8 3.3 0.5 0.3
    2 017−08−24 10:41:36.910 62.69 59.72 58.1 3.2 0.3 0.5
    2 017−08−29 22:47:54.601 62.31 59.18 35.3 1.7 1.3 0.3
    2 017−09−15 16:14:20.470 62.65 59.99 62.9 3.8 0.4 0.5
    2 017−10−02 17:14:29.735 62.26 59.16 38.2 1.8 1.0 0.3
    2 017−10−03 05:40:27.856 62.17 59.15 38.5 1.9 1.1 0.3
    2 017−10−28 00:51:54.113 62.37 59.56 59.3 2.4 0.7 0.4
    2 017−11−12 04:58:43.418 62.10 59.14 65.7 2.9 0.9 0.3
    2 017−11−13 02:23:55.693 61.96 60.30 103.8 2.5 0.6 0.4
    2 017−11−15 13:24:29.718 62.26 59.33 54.3 3.2 0.8 0.3
    2 017−11−24 07:37:47.287 62.93 60.30 102.8 4.2 0.4 0.6
    2 017−12−28 23:20:23.074 62.26 59.32 53.9 3.3 0.9 0.3
    下载: 导出CSV
  • 彼得·鲍曼. 2002. 新地震观测实践手册[M]. 中国地震局监测预报司,译. 北京:地震出版社:584−585.

    Borman P. 2002. New Manual of Seismological Observatory Practice[M]. Monitoring and Forecasting Department of China Earthquake Administration,trans. Beijing:Seismological Press:584−585 (in Chinese).

    常利军,丁志峰,孙为国. 2012. 第27次南极科学考察度夏期间长城站地震观测[J]. 极地研究,24(1):95–100.

    Chang L J,Ding Z F,Sun W G. 2012. Seismological observation at the Great Wall Station during the 27th Chinese Antarctic expedition in the summer[J]. Chinese Journal of Polar Research,24(1):95–100 (in Chinese).

    常利军. 2013. 我国近期南极地震科考进展[J]. 科学,65(2):17–20.

    Chang L J. 2013. Advances in Chinese Antarctic earthquake expedition[J]. Science,65(2):17–20 (in Chinese).

    陈圣源,刘方兰,梁东红. 1997. 南极布兰斯菲尔德海域地球物理场与地质构造[J]. 海洋地质与第四纪地质,17(1):77–86.

    Chen S Y,Liu F L,Liang D H. 1997. Geophysical field and geologic structure of the Bransfield sea area in the Antarctic circle[J]. Marine Geology &Quaternary Geology,17(1):77–86 (in Chinese).

    何高文,付少英,邓希光,赵庆献,王力峰,王海峰,关晓春,刘胜旋,庞云天,陈玺,柴祎,罗贤虎,朱本铎,朱振华,李强,苏丕波,何赵,吴聪,黄宁,何发光,蓝明华. 2018. 南极南设得兰群岛海域地质新进展:“海洋六号”南极科考初步成果[J]. 地球学报,39(6):643–656.

    He G W,Fu S Y,Deng X G,Zhao Q X,Wang L F,Wang H F,Guan X C,Liu S X,Pang Y T,Chen X,Chai Y,Luo X H,Zhu B D,Zhu Z H,Li Q,Su P B,He Z,Wu C,Huang N,He F G,Lan M H. 2018. New geological progress in offshore area of South Shetland Islands,Antarctic:Preliminary results of R/V “Hai Yang Liu Hao” expedition[J]. Acta Geoscientica Sinica,39(6):643–656 (in Chinese).

    林云芳,贾国强. 1989. 中国南极长城地震台与地磁台的建立[J]. 地震学报,11(2):213–219.

    Lin Y F,Jia G Q. 1989. Establishment of China’s Antarctic seismological and geomagnetic observatories[J]. Acta Seismologica Sinica,11(2):213–219 (in Chinese).

    刘瑞丰,陈翔,沈道康,郑秀芬,杨辉,陈宏峰. 2014. 宽频带数字地震记录震相分析[M]. 北京:地震出版社:51−59.

    Liu R F,Chen X,Shen D K,Zheng X F,Yang H,Chen H F. 2014. Seismic Phase Analysis of Broadband Digital Seismic Records[M]. Beijing:Seismological Press:51−59 (in Chinese).

    刘小汉,郑祥身,鄂莫岚. 1991. 南极洲大地构造区划和冈瓦纳运动[J]. 南极研究,3(2):1–9.

    Liu X H,Zheng X S,E M L. 1991. Review on the geotectonic division and Gondwana movement,Antarctica[J]. Antarctic Research,3(2):1–9 (in Chinese).

    吕文正,吴水根. 1989. 布兰斯菲尔德海槽:新生代裂谷[J]. 南极研究,1(4):28–35.

    Lü W Z,Wu S G. 1989. The Bransfield trough:A Cenozoic rift[J]. Antarctic Research,1(4):28–35 (in Chinese).

    马龙,邢健. 2020. 南极布兰斯菲尔德海峡及邻区地壳结构反演及构造解析[J]. 海洋与湖沼,51(2):265–273.

    Ma L,Xing J. 2020. Structure inversion and its tectonic interpretation in Bransfield Strait and the adjacent area,Antarctic[J]. Oceanologia et Limnologia Sinica,51(2):265–273 (in Chinese).

    万永革. 2016. 地震学导论[M]. 北京:科学出版社:316−318.

    Wan Y G. 2016. Introduction to Seismology[M]. Beijing:Science Press:316−318 (in Chinese).

    王公念,陶军,吴宣志. 1997. 南极布兰斯菲尔德海峡海磁异常和深部地质[J]. 极地研究,9(2):145–151.

    Wang G N,Tao J,Wu X Z. 1997. Marine magnetic anomalies and deep geology of the Bransfield Strait,Antarctica[J]. Chinese Journal of Polar Research,9(2):145–151 (in Chinese).

    王改云,邓希光,刘金萍,杜民. 2013. 西南极南设得兰群岛中新生代地层特征[J]. 极地研究,25(2):161–166.

    Wang G Y,Deng X G,Liu J P,Du M. 2013. Stratigraphic characteristics of the Mesozoic-Cenozoic in the South Shetland Islands,Antarctica[J]. Chinese Journal of Polar Research,25(2):161–166 (in Chinese).

    吴水根,吕文正. 1988. 德雷克海峡的扩张历史及其影响[J]. 南极研究,1(2):1–7.

    Wu S G,Lü W Z. 1988. The history of spreading of Darke Passage and its effects[J]. Antarctic Research,1(2):1–7 (in Chinese).

    姚伯初,王光宇,陈邦彦,陈圣源. 1995. 南极布兰斯菲尔德海峡的地球物理场特征与构造发育史[J]. 南极研究(中文版),7(1):25–35.

    Yao B C,Wang G Y,Chen B Y,Chen S Y. 1995. The characteristics of geophysical field and tectonical evolution in the Bransfield Strait[J]. Antarctic Research (Chinese Edition),7(1):25–35 (in Chinese).

    Almendros J,Carmona E,Jiménez V,Díaz-Moreno A,Lorenzo F. 2018. Volcano-tectonic activity at deception island volcano following a seismic swarm in the Bransfield Rift (2014−2015)[J]. Geophys Res Lett,45(10):4788–4798.

    Almendros J,Wilcock W,Soule D,Teixidó T,Vizcaíno L,Ardanaz O,Granja-Bruña J L,Martín-Jiménez D,Yuan X,Heit B,Schmidt-Aursch M C,Geissler W,Dziak R,Carrión F,Ontiveros A,Abella R,Carmona E,Agüí-Fernández J F,Sánchez N,Serrano I,Davoli R,Krauss Z,Kidiwela M,Schmahl L. 2020. BRAVOSEIS:Geophysical investigation of rifting and volcanism in the Bransfield Strait,Antarctica[J]. J South Am Earth Sci,104:102834.

    An M J,Wiens D A,Zhao Y,Feng M,Nyblade A,Kanao M,Li Y S,Maggi A,Lévêque J J. 2015. Temperature,lithosphere-asthenosphere boundary,and heat flux beneath the Antarctic Plate inferred from seismic velocities[J]. J Geophys Res:Solid Earth,120(12):8720–8742.

    Barker P F. 1982. The Cenozoic subduction history of the Pacific margin of the Antarctic Peninsula:Ridge crest-trench interactions[J]. J Geol Soc,139(6):787–801.

    Craddock C. 1982. Antarctic Geoscience[M]. Madison:University of Wisconsin Press:3−14.

    Dziak R P,Park M,Lee W S,Matsumoto H,Bohnenstiehl D R,Haxel J H. 2010. Tectonomagmatic activity and ice dynamics in the Bransfield Strait back-arc basin,Antarctica[J]. J Geophys Res:Solid Earth,115(B1):B01102.

    Erik L,Ekström G,Nettles M. 2013. Global CMT Catalog Search[DB/OL]. [2023−05−15]. https://www.globalcmt.org/CMTsearch.html.

    González-Casado J M,Robles J L G,López-Martínez J. 2000. Bransfield Basin,Antarctic Peninsula:Not a normal backarc basin[J]. Geology,28(11):1043–1046.

    González-Ferrán O. 1983. The Larsen Rift:An active extension fracture in West Antarctica[M]//Antarctic Earth Science. Cambridge:Cambridge University Press:344−346.

    Institute for Geophysics,University of Texas. 2015. Institute for Geophysics[DB/OL]. [2023−05−15]. https://www-udc.ig.utexas.edu/external/plates/data.htm.

    Janik T,Środa P,Grad M,Guterch A. 2006. Moho depth along the Antarctic Peninsula and crustal structure across the landward projection of the Hero Fracture Zone[M]//Antarctica:Contributions to Global Earth Sciences. Berlin,Heidelberg:Springer-Verlag:229−236.

    Jiménez-Morales V,Almendros J,Carmona E. 2022. Long-Term evolution of the seismic activity preceding the 2015 seismic crisis at Deception Island Volcano,Antarctica (2008−2015)[J]. Surv Geophys,43(3):959–994.

    Kaminuma K. 1995. Seismicity around the Antarctic Peninsula[C]//Proceedings of the NIPR Symposium on Antarctic Geosciences. Tokyo:National Institute of Polar Research:35−42.

    Kaminuma K. 2006. Seismicity in the Antarctic and surrounding ocean[J]. J Ind Geophys Union,10(1):15–24.

    Lawver L A,Sloan B J,Barker D H N,Ghidella M,von Herzen R P,Keller R A,Klinkhammer G P,Chin C S. 1996. Distributed,active extension in Bransfield Basin,Antarctic Peninsula:Evidence from multibeam bathymetry[J]. GSA Today,6(11):1–6.

    Li W,Yuan X H,Heit B,Schmidt-Aursch M C,Almendros J,Geissler W H,Chen Y. 2021. Back-arc extension of the Central Bransfield Basin induced by ridge–trench collision:Implications from ambient noise tomography and stress field inversion[J]. Geophys Res Lett,48(21):e2021GL095032.

    Olivet J L,Bettucci L S,Castro-Artola O A,Castro H,Rodríguez M,Latorres E. 2021. A seismic swarm at the Bransfield Rift,Antarctica[J]. J South Am Earth Sci,111:103412.

    Parera-Portell J A,de Lis Mancilla F,Morales J,Almendros J,Jiménez-Morales V. 2021. Structure of the crust and upper mantle beneath the Bransfield Strait (Antarctica) using P receiver functions[J]. Tectonophysics,802:228744.

    Pelayo A M,Wiens D A. 1989. Seismotectonics and relative plate motions in the Scotia Sea region[J]. J Geophys Res:Solid Earth,94(B6):7293–7320.

    Poli P,Cabrera L,Flores M C,Báez J C,Ammirati J B,Vásquez J,Ruiz S. 2022. Volcanic origin of a long-lived swarm in the Central Bransfield Basin,Antarctica[J]. Geophys Res Lett,49(1):e2021GL095447.

    USGS,Earthquake Hazards Program. 2017. Advanced National Seismic System (ANSS) Comprehensive Catalog of Earthquake Events and Products:Various[DB/OL].[2023−05−15]. https://doi.org/10.5066/F7MS3QZH.

    Robertson Maurice S D,Wiens D A,Shore P J,Vcra E,DoRman L M. 2003. Seismicity and tectonics of the South Shetland Islands and Bransfield Strait from a regional broadband seismograph deployment[J]. J Geophys Res:Solid Earth,108(B10):2461.

  • 期刊类型引用(9)

    1. 李琼,付虹,胡小静,杜家云,杨黎. 云南地区氡异常与地震的统计特征分析. 地震学报. 2024(02): 340-352 . 本站查看
    2. 刘洁,翟宏光,张国强,朱琳,邱玉荣. 陕西凤翔井水位和区域应力场变化与地震活动性关系研究. 地震. 2024(03): 124-137 . 百度学术
    3. 刘伟,白细民,吕少杰,史浙明,齐之钰,何冠儒. 基于井水位气压效应计算含水层的水力参数. 地震地质. 2023(03): 652-667 . 百度学术
    4. 王秀娟,周璇. 江苏睢宁苏02、03井水位异常变化机理探讨. 四川地震. 2023(03): 31-37 . 百度学术
    5. 解晓静,张艺潇,张帆,李盛,郑在壮. 基于数字化水位的琼东北地区构造应力场时序特征分析. 地震科学进展. 2022(03): 110-119 . 百度学术
    6. 芮雪莲,杨耀,官致君,林圣杰. 四川会理川-31井水位变化与降雨量关系及异常识别. 华北地震科学. 2021(01): 78-83 . 百度学术
    7. 解晓静,苏荣托雅,张帆,孙三健. 海南琼海加积井水位水温同步上升后转平的异常成因初探. 华南地震. 2019(S1): 37-45 . 百度学术
    8. 王喜龙,杨振鹏,刘天龙. 辽宁高七井静水位破年变异常特征分析. 防灾减灾学报. 2019(04): 7-13 . 百度学术
    9. 王熠熙,孙小龙,邵永新,刘双庆,李悦,李赫,王喜龙,龚永俭. 天津王3井水位异常机理. 地震学报. 2019(06): 757-768+808 . 本站查看

    其他类型引用(0)

图(8)  /  表(3)
计量
  • 文章访问数:  219
  • HTML全文浏览量:  65
  • PDF下载量:  83
  • 被引次数: 9
出版历程
  • 收稿日期:  2023-05-14
  • 修回日期:  2023-06-01
  • 录用日期:  2023-06-06
  • 网络出版日期:  2023-10-17
  • 刊出日期:  2024-11-19

目录

/

返回文章
返回