欧洲扇区不同纬度电离层暴特征的统计分析
摘要：本文利用Madrigal数据库的TEC数据对2001—2010年间的156次单主相型磁暴事件,统计分析了欧洲扇区从赤道到极光带共5个纬度区域的电离层暴特征,结果表明：（1）电离层暴有明显的纬度分布特征,正负暴出现次数的比例随纬度的降低呈现明显的增加趋势,但夏季赤道地区趋势相反,正负暴比例比更高纬度的反而降低;（2）与主相相比,恢复相期间大部分纬度地区正暴数量减少,负暴数量增加,但赤道地区恢复相期间正暴数量反而增加;（3）中低纬地区电离层暴随磁暴MPO地方时分布特征明显,正暴所对应的MPO主要分布在白天,而MPO发生在夜间容易引起负暴;（4）电离层负暴主要发生在夜间,中、高纬地区负暴的开始时间存在‘时间禁区’,但不同纬度‘时间禁区’的地方时分布有一定差异,正暴分布则相对分散.第23太阳活动周期太阳风参数及地磁指数的统计分析
摘要：日冕物质抛射（Coronal Mass Ejection,简称CME）和共转相互作用区（Corotating Interaction Region,简称CIR）是造成日地空间行星际扰动和地磁扰动的两个主要原因,提供了地球磁暴的主要驱动力,进而显著影响地球空间环境.为深入研究太阳风活动及受其主导影响的地磁活动的时间分布特征,本文对大量太阳风参数及地磁活动指数的数据进行了详细分析.首先,采用由NASA OMNIWeb提供的太阳风参数及地磁活动指数的公开数据,通过自主编写matlab程序对第23太阳活动周期（1996-01-01—2008-12-31）的数据包括行星际磁场Bz分量、太阳风速度、太阳风质子密度、太阳风动压等重要太阳风参数及Dst指数、AE指数、Kp指数等主要的地磁指数进行统计分析,建立了包括269个CME事件和456个CIR事件列表的数据库.采用事例分析法和时间序列叠加法分别对两类太阳活动的四个重要太阳风参数（IMF Bz、太阳风速度、太阳风质子密度、太阳风动压）和三个主要地磁指数（Dst、AE、Kp）进行统计分析,并研究了其统计特征.其次,根据Dst指数最小值确定了第23太阳活动周期内的355个孤立地磁暴事件,并以Dst指数最小值为标准将这些磁暴进一步分类为145个弱磁暴、123个中等磁暴、70个强磁暴、12个剧烈磁暴和5个巨大磁暴.最后,采用时间序列叠加法对不同强度磁暴的太阳风参数和地磁指数进行统计分析.统计分析表明,对于CME事件,Nsw/Pdyn（Nsw表示太阳风质子密度,Pdyn表示太阳风动压）线性拟合斜率一般为正;对于CIR事件,Nsw/Pdyn线性拟合斜率一般为负,这可作为辨别CME和CIR事件的一种有效方法.从平均意义上讲,相较于CIR事件,CME事件有更大的南向IMF Bz分量、太阳风动压Pdyn、AE指数、Kp指数以及更小的Dstmin.一般情况下,CME事件有更大的可能性驱动极强地磁暴.总体而言,对于不同强度的地磁暴,Dst指数的变化呈现出一利用ERA-Interim资料对平流层Brewer-Dobson环流变化趋势的分析
摘要：选用每天12∶00UTC时次的逐日ERA-Interim再分析资料,根据transformed Eulerian-mean（TEM）方程通过积分剩余速度v^-＊,研究了1979—2011年间Brewer-Dobson（BD）环流的时空演变规律.并将其与downward control（DC）原理研究的结果进行比较,同时还探讨了平流层温度与BD环流之间的相互联系.结果表明,由TEM方程通过积分剩余速度v^-＊估算的BD环流与利用DC原理估算的环流相比较,在热带地区的形势更加明显.环流在热带对流层中上层上升至平流层中下层,最高可达1hPa等压面附近.然后在热带外向极向下运动,最后在中高纬度下沉回到对流层.BD环流的上升中心及质量通量均随季节的变化产生变动,环流在冬半球的形势显著地强于夏半球.在春季和秋季期间,环流呈现出南北两半球的对称形势.从全球尺度物质输送的角度来看,在过去的33a间平流层BD环流的长期变化趋势是减弱的,且在平流层中下层减弱是明显的.环流的减弱趋势与纬向平均温度的长期变化趋势相匹配.中国区域夏季再分析资料高空变量可信度的检验
摘要：利用全球探空资料（IGRA）对1989—2008年美国国家环境预报中心（NCEP）和大气研究中心（NCAR）再分析资料、NCEP和美国能源部（DOE）再分析资料、NCEP气候预测系统再分析资料（CFSR）、日本气象厅25年再分析资料（JRA-25）、欧洲数值预报中心再分析资料（ERA-Interim）和美国国家航空航天局（NASA）现代回顾性再分析资料（MERRA）的高空变量在中国地区对流层中高层的可信度进行了初步的检验.分析结果表明：再分析资料对中高层位势高度和温度的夏季平均气候态具有较好的再现能力,其EOF的时空变化特征与观测吻合也较好;再分析资料的绝对湿度值较观测结果要偏大,其中MERRA与观测最为接近.再分析资料不能很好地反映经向风的夏季平均气候态及年际变化特征,EOF的时空模态和观测偏离也较大.总体而言,NCEP/NCAR、NCEP/DOE及NCEP/CFSR对这些变量的再现能力较JRA-25、ERA-Interim和MERRA弱.城市化对石家庄站日气温变化的影响
摘要：利用1962—2011年逐日平均、最低和最高气温资料,对比分析了石家庄站和藁城站平均、最低和最高气温的概率分布特征.石家庄站是我国少有的自建站以来从未迁址的城市气象站,而附近的藁城站可近似看作乡村气象站.结果表明：城市化致使石家庄站1962—2011年平均、最低和最高气温的概率密度分布向高温方向偏移,其中对最低气温分布的影响尤其明显,对最高气温分布的影响很小;受城市化影响,石家庄站最低气温概率密度分布的高温部分增温比低温部分增温更加明显,最低气温分布形状更加扁平;相对于1962—1986年,1987—2011年石家庄站平均、最低和最高气温概率密度分布均向高温方向偏移,其中最低气温偏移最为明显,并导致1962—2011年整个分析时期最低气温分布出现非正态性;城市化对石家庄站气温分布的影响在冬、春季比夏、秋季更显著,最显著的城市化影响出现在冬季最低气温上;石家庄站基于最低气温的极端气温指数趋势受城市化影响严重,冷夜日数和暖夜日数的城市化影响分别为-1.13d/10a、1.48d/10a,但基于最高气温的冷昼和暖昼日数等极端气温指数变化趋势受城市化影响不明显.出现这种现象的主要原因是城市化对最低、最高气温分布的影响存在差异.三峡地区上地壳结构的远震虚震源反射地震成像
摘要：远震虚震源反射成像方法利用远震初至波在台网之下地表与地下界面间形成的反射波（PPdp震相）波形资料进行台网地区地震反射结构研究.此方法先用台网各台站的平均初至波形求取震源信号,再用该震源信号与各道地震记录作反褶积从而取得反射剖面.本文介绍了远震虚震源反射成像的基本原理和实现步骤,并以三峡地区的观测资料为例,得到三峡库首区上地壳的反射地震剖面.成像剖面中解释出的四川盆地和秭归盆地的底界面位置和形态与地表地质观测和大地构造背景吻合.为验证虚震源成像的能力,本文使用弹性波正演模拟合成地震数据,经过处理实际资料一样的步骤获得虚震源成像结果.正演模拟表明,叠加多个远震的反射地震剖面可以有效地改善成像的信噪比和连续性;虚震源成像需要选择特定的震源频率范围,以减轻高频噪音以及低频造成的成像问题.在研究上地壳结构时,建议在不适合主动源采集和缺乏低频信号的地区尝试远震虚震源方法.2013年8月香格里拉德钦—得荣Ms5.9地震序列震源机制与应力场特征
摘要：利用中国区域台网地震波形记录,采用CAP方法反演了香格里拉德钦（位于云南省）—得荣（属于四川省）2013年8月28日Ms5.1、8月31日Ms5.9地震及8次Ms〉4余震的震源双力偶断层面解和震源质心深度.结合震区地质构造、余震分布、烈度分布、动力学背景等资料,分析了此次地震序列的震源机制和应力场特征.反演结果表明,此次地震序列为节面倾角倾斜的正断层型地震,发震断层为NWW向活动构造带.序列中最大地震MS5.9和次大地震Ms5.1地震的破裂节面分别为走向299°、倾角53°、滑动角-73°;走向290°、倾角55°、滑动角-72°.震源区受到强烈的水平拉张力、垂直挤压力作用.Ms5.9地震后续余震T、P轴方位角随时间变化强烈,表明Ms5.9地震后震源区应力调整作用明显.震源区应力场反演结果显示,地震发生的构造带上最大主拉应力为NNE-SSW向,最大主压应力为NW-SE向,与GPS观测所反映的地表最大主应力分布方向基本一致,表明震源区的应力状态可能主要受到背景大尺度构造应力场的控制.此次地震序列填充了川滇地区震源机制及应力场的空间分布图像,1976年以来可靠的震源机制解资料表明香格里拉次级块体是川滇块体及周边区域显著的拉张作用区域.香格里拉次级块体和保山次级块体正断层地震的断层节面及震源应力轴分布的空间变化,与GPS观测反映的地表最大主拉应力分布较一致,其空间分布特征反映了在青藏高原物质挤出背景下,块体之间相互作用、地势差异等作用对构造活动的影响.应用综合震源机制解法推断鄂尔多斯块体周缘现今地壳应力场的初步结果
摘要：利用2007年8月1日至2013年7月21日发生在鄂尔多斯块体周缘的8499个地震的49844个P波初动符号资料,应用综合震源机制解法获得了鄂尔多斯块体周缘0.25°×0.25°的精细地壳应力场,所得应力场结果基本上覆盖了整个鄂尔多斯周缘地区.研究结果表明鄂尔多斯周缘地壳应力场具有以下特征：（1）在环绕鄂尔多斯周缘的银川—吉兰泰断陷带、河套断陷带、岱海断陷带、山西断陷带和渭河断陷带内,综合震源机制解结果以正断层型为主,且综合震源机制解节面走向大体与控制断陷带边界的主要断裂走向相一致,与鄂尔多斯周缘断陷带现今的拉张状态相一致.（2）在鄂尔多斯西南缘,综合震源机制解类型主要为逆冲、逆冲走滑和走滑型,反映了鄂尔多斯块体在西南缘受到青藏高原北东向挤压作用.鄂尔多斯西南缘的应力场的主压应力方向在远处为东向,源自于青藏高原向东北挤压作用,靠近鄂尔多斯块体表现为北东—南西向.（3）P轴方位在局部地区变化较大,但总体呈现规律性变化.P轴方位在鄂尔多斯块体西缘,从南向北,主压应力轴方位更加偏北;在其北缘,由西向东,主压应力轴方位更加偏东.在其南缘和东缘,主压应力轴方位变化不大,大体上平行于控制各断陷带主要断裂走向.P轴倾角在西南缘为近水平,在其周缘各盆地内P轴倾角近直立.（4）T轴方位总体表现为北西—南东向;在鄂尔多斯周缘各断陷带内,T轴走向大体与控制断陷带主要断裂走向以及断陷盆地走向相垂直.（5）鄂尔多斯块体在其西南角受到来自青藏高原的北东向挤压和其东北角深部物质上涌形成的北西—南东向拉张力联合作用,上述作用使得鄂尔多斯块体周缘地区除西南区为挤压区外,其余区域均为剪切拉张区,与先前研究认为鄂尔多斯周缘地区处于引张应力场作用相符合,较好地解释了环鄂尔多斯周缘的�基于速度-状态摩擦本构定律的三维DDA方法
摘要：原有三维非连续变形分析（DDA）方法采用常摩擦系数的Mohr-Coulomb定律作为切向破坏准则,然而当描述更大尺度构造块体的运动与变形时,常摩擦系数不再适用.速度-状态摩擦本构定律能够定量描述地震周期各阶段断层面剪应力变化,解释发震断层行为.本文将速度-状态摩擦定律与三维DDA方法相结合,首先推导了计算摩擦系数的实用公式,随后通过滑动-保持-滑动实验与速度步进实验算例对改进的三维DDA方法进行了验证.结果表明,应用速度-状态摩擦本构定律的三维DDA方法能够比较准确地模拟静摩擦的时间依赖性与动摩擦的速度依赖性,解决了将三维DDA方法在地学中应用的基本问题.利用三维高斯射线束成像进行地震定位
摘要：常规的地震定位方法通常需要拾取地震记录的初至,当初至不明显或被较高水平的噪声淹没时精度较低.本文采用基于三维高斯射线束的偏移成像方法对震源进行定位,较好地解决了该问题.通过三维高斯射线束对台站记录进行偏移归位,并将各台站成像结果的交点作为地震能量释放的中心位置;当各台站成像结果不能交于一点时,采用三维空间高斯滤波方法可实现震源位置的自动获取.提出的变网格计算方案极大地减少了计算量,显著地提高了成像精度和计算效率.利用首都圈地震台网数据,对涿鹿、滦县以及房山三个地震事件进行试算,结果表明：基于变网格三维高斯束偏移成像的地震定位方法自动化程度很高,而且具有较好的抗噪能力,特别适合处理低信噪比资料的地震定位问题.地幔柱与岩石圈相互作用过程的数值模拟
摘要：地幔柱的研究是地球科学研究的热点之一.本文主要集中研究地幔柱与岩石圈的相互作用过程.基于质量守恒方程、动量守恒方程和能量守恒方程,通过有限元数值方法可以计算得到地幔柱与岩石圈相互作用的温度场、速度场和有效黏度等的时空图.本文的流变本构模型主要基于非牛顿流体的有效黏度模型,通过数值模拟计算分析了地幔柱与岩石圈相互作用过程,着重讨论了地壳流变结构对此过程的影响.数值模拟结果显示,地幔柱与岩石圈的相互作用分为三个阶段：地幔柱上升期,时间持续到0.2 Ma,平均速度为2.75m·a^-1,地幔柱顶部地形开始向上隆起;地幔柱与岩石圈纵向作用期,时间从0.2 Ma到0.26 Ma,地幔柱上升的平均速度为0.83m·a^-1,地表地形隆升达到最大值;地幔柱与岩石圈横向作用期,0.26 Ma以后,岩石圈开始剪切变形,地幔柱水平运动速度为0.47m·a^-1,当剪切变形达到一定程度,岩石圈底部开始出现拆沉作用.当下地壳流变强度比较小时,上地壳的流变结构控制着地幔柱顶部地表地形隆起程度,流变强度越大,隆升高度越小;而下地壳的流变结构控制着地幔柱两侧地表地形的下沉幅度,下地壳流变强度越小,下沉幅度越大.最后,讨论了数值模拟对峨眉山大火成岩省地幔柱发展演化的应用.礼县—罗家堡断裂晚第四纪活动特征：兼论1654年礼县8级地震孕震机制
摘要：1654年7月21日甘肃省礼县8级地震发生在南北地震带的中北段,该地区的构造变形和构造活动与青藏高原向北东方向的扩展密切相关,复杂的构造几何特征主要受控于东昆仑断裂、西秦岭北缘断裂和一系列北东向断裂.礼县—罗家堡断裂为一条北东东向的左旋走滑活动断裂,错断了含有仰韶文化红色陶瓷片的一级阶地堆积物,阶地面上断层陡坎高约1.5m.沿断裂带发现冲沟的左旋位错量为3～10m,晚更新世黄土中残留的断层陡坎高4.5～8m.其中两条冲沟中发现裂点,高3.5～3.9m,距断层陡坎的距离为16m.礼县—盐关—罗家堡—天水镇一带发育大量的滑坡,长轴走向与礼县—罗家堡断裂一致,滑坡体后缘、侧壁陡峭,出露晚第四纪黄土,鲜有植被覆盖.礼县—罗家堡断裂为1654年8级地震的发震断裂.综合分析认为,受青藏高原向北东方向的扩展,被西秦岭北缘断裂、礼县—罗家堡断裂和岷县—宕昌断裂围限的礼县次级地块向南东滑动可能导致了1654年礼县8级地震的发生.断裂北、南两侧地壳电性结构的差异为强震的孕育提供了深部构造条件.用于区域重力场定量解释的多尺度刻痕分析方法
摘要：本文介绍一个把小波多尺度分析、表面刻痕分析以及位场频率域解释理论和反演方法结合起来的数据处理、反演解释和信息提取的方法系统.这一方法系统简称为区域重力场多尺度刻痕分析方法,应用于刻画地壳分层的三维密度结构、地壳变形带分布和构造单元分区.多尺度刻痕分析包含频率域重力场场源分层、重力场小波变换多尺度分解、场源分层深度及密度扰动反演、分层刻痕分析和构造边界定位四个子系统.文中扼要地介绍这四个子系统基本原理、方法技术及应用效果.从地球物理探测到大地构造学发现,是一个多学科综合研究的探索过程.要取得重大研究成果,必须研发和组合来自不同学科的多个新方法技术,使多学科综合研究有宽厚的理论支撑.本文介绍的四个子系统组合的理论支撑分别来自应用数学、地球物理学和信息科学.鄂尔多斯—中秦岭—四川东部的重力异常场与深部地壳结构
摘要：本文基于在跨越鄂尔多斯盆地、中秦岭造山带与渭河盆地、四川盆地东部长达1010km的综合地球物理探测剖面,实地采集的最新高精度重力位场数据,分析了不同构造单元的重力位场特征,构建了该剖面区域的地壳密度结构模型.进而通过分析地壳内部不同密度界面的横向差异特征及分区,确定了该剖面区域内的断裂分布,提出了鄂尔多斯盆地—中秦岭造山带—四川盆地的＂盆-山-盆＂型构造体系的新认识.西北太平洋岩石圈有效弹性厚度及其构造意义
摘要：本文引入滑动窗口导纳技术（MWAT）,计算西北太平洋岩石圈有效弹性厚度（Te）.首先,基于SIO V15.1海底地形模型,模拟研究了MWAT法计算Te的精度,表明当Te〈5km时,误差在±1km以内,当Te≥5km时,相对误差在10%以内.分别采用GEBCO、SIO V15.1和BAT_VGG海底地形模型,构建了西北太平洋Te,通过对获得的洋壳密度参数和实测导纳与模型导纳之差的均方根进行分析,结果表明,BAT_VGG模型更适用于Te计算.西北太平洋Te均值为13.2km,标准差为6.9km,以板块冷却模型为参考,主要分布在150℃～450℃等温线深度范围内.白垩纪和侏罗纪时期岩石圈Te分布在150℃～300℃等温线深度范围内,且未随海山加载时岩石圈年龄增大而增大,说明海山加载时岩石圈年龄不是影响其强度的唯一因素.南太平洋超级海隆活动,以及研究区域广泛存在的断裂带构造,都曾对本区域岩石圈演化产生过重要影响,可能是本地区岩石圈Te较小的构造原因.基于抛物线密度模型的频率域三维界面反演及其在川滇地区的应用
摘要：密度界面反演作为了解地球内部结构的一种重要方法,长期以来都是重力学研究的主要内容.本文结合抛物线密度模型及频率域算法的优点,将抛物线密度函数应用于Parker-Oldenburg算法,经过理论推导得到了抛物线密度模型的频率域公式,从而建立了基于抛物线密度模型的三维密度界面重力异常正反演的算法和流程.理论模型数据试验表明本方法快速、有效,适用于大多数浅部比深部增加更快的实际地壳密度.研究中还利用该方法对川滇地区重力异常进行了反演,获得了该区的莫霍面深度分布,并与接收函数研究结果进行对比分析,进一步验证了本文方法的正确性和有效性.华北地区上地幔及过渡带电性结构研究
摘要：采用远参考道和Robust技术,处理了华北地区14个地磁台站资料,得到了相干度超过0.8的地磁测深响应函数.并将其转换为大地电磁测深的响应函数,获取了10^5～10^7 s周期范围内的视电阻率和相位.应用ρ^＋理论对数据进行了一致性检验和反演,结果表明417km,850km深度附近可能存在电性间断面.同时采用基于一维最光滑模型的Occam反演方法得到了300～1000km范围的地幔电性结构,并与前人在其他地区的研究结果进行了对比.发现华北地区地幔过渡带的电导率在大兴安岭—太行山重力梯度带东西两侧表现不同,重力梯度带附近及西侧台站下方过渡带深度的电导率和北美的Tucson地区相当,而华北地区东部的电导率在地幔过渡带范围高出西侧约2～5倍,这很可能和太平洋板块的俯冲有关.两次近距离大震前成都台视电阻率重现性、相似性和各向异性变化
摘要：本文评价了四川汶川MS8.0、芦山MS7.0地震前后成都台地电观测环境,研究了该台视电阻率变化.结果为：（1）两次大震发生在该台以西的龙门山断裂带、震源机制和震源深度接近,是近距离大震,相应地,在两次地震前该台两个正交测道中的每一测道视电阻率变化均显示了中期异常及其变化过程的重现性和异常幅度的相似性;（2）在每次地震前,两个测道表现了异常变化形态、幅度和起始时间的差异性.其重现性、相似性证明这些异常与两次大震晚期孕育有关;差异性主要展示了与震源机制有直接联系的视电阻率各向异性变化,揭示了震前该台地下介质经历了强烈的电性各向异性变化的过程.
