中近磁尾等离子体片统计特性研究
摘要：本文使用Cluster-C1卫星的CIS仪器和FGM仪器测量得到的质子通量数据和计算的β数据,判断Cluster卫星在地球磁尾不同位置位于等离子体片内的概率.使用2001—2004年7—11月的Cluster-C1数据,分别在行星际磁场南向和北向时,得出X〈-10RE区域内卫星位于等离子体片的概率在Y-Dz平面的分布图（Dz是卫星到中性片的距离）.通过对比行星际磁场南向和北向时的卫星位于等离子体片的概率的分布图,我们发现等离子体片在行星际磁场南向时比在行星际磁场北向时要薄,并且这个效应在磁尾晨昏两侧比在午夜附近明显,同时我们还发现等离子体片在晨侧比在昏侧厚.基于COSMIC资料分析电离层F层不规则体结构
摘要：根据电离层不规则体的产生会导致周围电子浓度发生起伏变化的原理,利用2007年COSMIC掩星系统的TEC数据,通过平滑滤波得到TEC的扰动值ΔTEC的变化,利用其研究F层不规则体的时空变化特征.统计结果表明：扰动较大的掩星事件主要发生在磁纬±20°之间和高纬地区,春季和秋季带状分布较为明显,不同经度地区较强扰动的掩星事件的分布也有不同特征;较强ΔTEC的掩星事件主要发生在地方时午夜前和午夜后两个时段,发生的高度主要在250～400km范围内.这些结果与已知的F层不规则体的时空分布特征较为一致,说明利用TEC的扰动量来分析电离层F层不规则体结构是可行的.基于行星际/太阳风和地磁条件的紫外极光卵边界建模和预测
摘要：极光卵的尺度大小与太阳风-磁层-电离层能量耦合过程紧密相关,准确预测其大小对空间天气研究和预报具有非常重要的意义.本文基于模糊c均值聚类算法,从Polar卫星紫外极光图像中自动提取极光卵边界数据（～1215000个赤道向边界点和-3805000极向边界点）,统计分析其与太阳风等离子体、行星际磁场、地磁指数等之间的相关特性,并构建了以行星际、太阳风为模型参数（模型1）和以行星际、太阳风及地磁指数为模型参数（模型2）的2种极光卵边界多元回归模型.以模型预测的极光卵边界与实际极光卵边界之间的平均绝对误差作为模型评价标准,将本文预测模型与Carbary（2005）模型和Milan（2009）模型进行了对比.结果表明,模型2对极光卵极向、赤道向边界预测的平均绝对误差为1.55和1.66地磁纬度,优于Carbary和Milan模型（Carbary模型极向、赤道向边界的平均绝对误差为2.18和5.47地磁纬度,Milan模型极向、赤道向边界的平均绝对误差为1.71地磁纬度和1.90地磁纬度）.北京MST雷达探测中间层-低热层观测结果初步分析
摘要：北京MST雷达是子午工程建设的国内仅有的两部MST雷达之一,为研究其在中间层-低热层MLT区域的探测能力以及数据可靠性,本文应用北京MST雷达2012、2013两年高模式数据,从数据获取率、与廊坊流星雷达测风对比以及风场时空分布特征三个方面进行初步分析.结果是：（1）数据获取率日变化特征为：白天65-100km均可获取数据,数据获取率的高值区主要集中在70-80km,最大值可达80%;夜间主要集中在80-100km,数据获取率在30%及以下.表明该MST雷达白天可以探测到电离层D层和E层低层,夜间D层消失,只探测到E层低层.季节变化特征为：夏季白天可获取数据的时间和高度区间都比较大,春季次之,冬季最小.夏季白天以及日落后1h内可探测到120km.（2）对北京MST雷达与廊坊流星雷达2012年5月份、80-100km高度区间测量的水平风进行对比分析,二者测风结果在时空分布上有很好的一致性,表明MST雷达探测数据是可靠的.（3）2012年和2013年相应月份平均的纬向风、经向风时空分布特征有较高的一致性,并与HWM07模式结果也基本一致.上述初步分析结果表明,北京MST雷达对中间层-低热层60-120km高度区域已具备较强的探测能力,所得结果将可用于MLT过程揭示与驱动因子研究,并可与该高度上其他探测手段作综合研究.夜光云内小尺度重力波对冰晶粒径谱的影响规律研究
摘要：本文利用AIM卫星搭载的CIPS云图反照率和冰晶粒径数据,从中提取了2007/08南半球和2008年北半球共6489个小尺度重力波活动（波长5~150km范围）个例,对重力波区域与背景云层冰晶粒径谱进行对比分析,从而研究重力波对冰晶平均半径和谱宽的影响规律.结果表明,北半球重力波区域冰晶的平均半径和谱宽分别比背景云层小2.5nm和6.1nm,南半球则分别减小1.1nm和7.9nm.在随纬度的分布上,小于80°时,南北半球的平均半径扰动值均为负值,绝对值随纬度增大而减小,而大于80°时,负扰动转变为正扰动,且绝对值增加;谱宽扰动的绝对值也随着纬度增加而减小,但均为负值.在季节内随时间的分布上,南北半球重力波对冰晶平均半径和谱宽的扰动在始末阶段以负值为主,且绝对值较大,而在中期阶段正负值相当,且绝对值较小.这一特征与重力波引起冰晶粒径变化的振幅在纬度和时间上的分布趋势一致.重力波的波长均随纬度升高而减小,在季节的始末阶段较大,中期小,且南半球的平均波长和变化幅度都要明显大于北半球的,粒径扰动振幅随波长的变化率为南半球0.207nm·km-1,北半球的0.163nm·km-1.根据分析推断,重力波自身的扰动振幅应与其影响区域内的谱参数相对于背景云层的变化量有直接关系,且振幅越大,平均半径和谱宽的负扰动就越大.格点大气环流模式GAMIL2参数不确定性的量化分析与优化
摘要：物理过程参数化方案的不确定性是目前气候系统模式不确定性的重要来源之一.随着模式内在复杂度攀升,模拟场景多样化,参数化方案中基于先验的和人工的物理参数选取方法已经逐步成为限制模式模拟能力的瓶颈之一.为此,本文设计并提出了初选与寻优相结合的两步法参数优化方案.初选阶段用全因子采样方法对不确定参数空间进行初始敏感性分析,估计最优解所在区域;寻优步采用单纯型下山法,基于初选阶段确定的参数组合快速寻优.将两步法应用于中国科学院大气物理研究所（英文缩写：IAP）大气科学和地球流体力学数值模拟国家重点实验室（英文缩写：LASG）格点大气模式第2版：GAMIL2,选取其深对流方案和云量方案中的3个重要参数开展寻优,优化以综合减小模式降水、风场、温度、湿度、位势高度以及辐射通量的误差为目标.这些变量用GAMIL2标准版本标准化后形成单一的目标.结果显示,优化后的目标函数值比GAMIL2标准版本改进了7.5%.机理分析表明,调优后的参数优化了大气中的水汽凝结作用,进而减少模式的湿度偏差,改进云量的模拟效果;同时水汽凝结作用的变化通过大气内部动力和热力相互作用及响应影响温度、位势高度和风场的模拟.垂直重力梯度反演Moho面的频谱域公式及其应用
摘要：通过求解引力相等原则下的Fredholm积分方程,可以得到不规则单一密度界面（Moho面）的起伏.本文充分参考了前人的理论研究,推导出扰动垂直重力梯度确定Moho面深度的频谱域表达式,该式具有二次项迭代精度.运用此公式进行了全球Moho面的恢复计算,并将该结果与CRUST1.0模型和GEMMA Moho模型进行了对比和验证.地表负载及GPS测站分布对参考框架转换的影响分析
摘要：GPS数据处理通常采用Helmert七参数转换将瞬时站坐标转换到指定的框架下,但瞬时站坐标中尚未模型化的季节性地表负载会影响平移参数（地心运动）和尺度参数的估值,进而影响测站坐标;不均匀的测站分布会加剧这一影响.本文利用GRACE重力场系数仿真GPS地表负载的实验表明,基于网平移法采用实际的IGS站至少能够恢复90%的地心运动信号.地表负载及GPS实际测站的不均匀分布可以解释大约30%的GPS尺度的周年变化.相对于IGb08的所有框架站,目前采用91个全球均匀分布的核心站作为框架转换的基准是合理的.采用IGb08的所有框架站进行转换会导致U方向误差增加,特别是对框架站密集的欧洲区域（误差均值约为1mm）.因此框架转换时,应尽量选取均匀分布的测站,同时不估计尺度参数.东北地区背景噪声的Rayleigh和Love波相速度层析成像
摘要：本文利用中国数字地震台网位于东北地区的122个宽频地震台站的18个月记录的三分量连续地震噪声数据,采用互相关方法提取了Rayleigh和Love波经验格林函数,并利用时频自动分析技术获取了相应的相速度频散曲线.通过反演频散曲线,获得了Rayleigh和Love波周期为8-35s的二维相速度分布.结果表明,东北地区相速度的分布存在横向和垂向的不均匀性.短周期的相速度分布同地表地质构造密切相关,松辽盆地及山间沉积盆地呈现低速异常,而大兴安岭、小兴安岭及东部的一些山岭显示高速异常.随着周期的增加,位于中间的松辽盆地变为高低速相间,两侧的造山带呈现低速异常.这种异常的转变,可能是受构造活动或者莫霍面深度的影响.另外,在周期为20-35s频段内,Rayleigh和Love波同一周期的相速度在松辽盆地和位于吉林地区的郯庐断裂带表现不一致,表明可能存在径向各向异性.郯庐断裂带中南段及邻区Pn波速度结构与各向异性
摘要：郯庐断裂带是一条纵贯我国大陆东部NNE走向的巨型深断裂,其中南段及邻区（115°E—122°E,29°N—38°N）跨越了华北断块区、扬子断块区和华南褶皱系三大一级构造单元,由于其重要性和复杂性,长期以来一直是地学家们研究的热点.本文从国际地震中心（ISC）、中国地震台网及区域地震台网的地震观测报告中精心挑选出6381个Pn震相数据,用Pn波时间项层析成像法反演得到了郯庐断裂带中南段及邻区上地幔顶部Pn波速度结构和各向异性.结果显示,研究区上地幔顶部具有显著的横向非均匀性,相对于7.95km·s-1的平均速度而言,Pn波速度值在7.68-8.24km·s-1范围内变化.Pn波速度分布在郯庐断裂带中段和南段具有分段性：沿中段及周边存在一NE向低速异常带,低速可能是由于岩石圈的减薄和软流圈的高温物质沿郯庐带上涌导致;沿南段表现为一NNE向弱高波速异常带,作为高低速的边界带清晰地勾勒出了华北与扬子这两个不同块体,该边界在江苏域向华北地块NW方向凹进.Pn波速度各向异性的强弱与速度分布存在一定的相关性.总体上,如鲁西隆起及以南等低速区、茅山断裂附近的高低速过渡带,其速度各向异性较为强烈;而在具有高速异常的苏北盆地、合肥盆地等稳定区域下方其各向异性较弱.本文通过Pn波震相基本未能探测到郯庐断裂带中段的方位各向异性,推测是上地幔顶部被＂冻结＂下来的各向异性痕迹被软流圈热物质上涌这一强烈构造运动削弱所导致.南段具有与断裂伸展方向近乎平行的快波速方向.Pn波速度横向变化和强震活动存在一定关联.强震主要发生在Pn波低速异常区或高低速过渡带上.郯城8.5级地震震中位于中段和南段高低速过渡带,该区域也是速度横向变化最大的地方,最容易集中应力和产生应力差.利用接收函数反演青藏高原西部地壳S波速度结构
摘要：相对于宽阔的腹地,青藏高原西部南北向宽度仅约600km,却记录了印度和欧亚板块汇聚的深部过程及其响应.本文用22台宽频带流动地震台站在西缘构建了一条南北向探测剖面（～80°E,TW-80试验）.利用接收函数反演剖面下方S波速度结构,综合西部已有的宽频带探测结果,分析认为：印度板块向北俯冲可能已到达班公湖—怒江缝合带附近,俯冲过程中下地壳发生榴辉岩化;喀拉昆仑断裂带、班公湖—怒江缝合带、阿尔金断裂带均为切穿地壳的深断裂,莫霍面发生错断;喀拉昆仑断裂带和龙木错断裂带之间的中上地壳没有发现连续的S波低速体,说明可能缺乏解耦层,支持青藏高原西部地壳为整体缩短增厚模式.太行山东缘石家庄南部地壳结构及断裂活动性探测
摘要：采用深、浅地震反射和钻孔地质剖面相结合的方法,对太行山东缘石家庄南部的地壳结构和隐伏断裂的活动性进行了研究.深地震反射探测结果表明,该区地壳厚度33～38km,莫霍面从华北平原向太行山下倾伏.石家庄—晋县凹陷是受拆离断层控制的盆岭构造,太行山山前断裂为凹陷的西边界断裂,表现为上陡下缓的铲形断裂.石家庄—晋县凹陷中还发育北席断裂和栾城断裂,它们与太行山山前断裂一样受拆离断层的控制,未错断早更新世晚期以来沉积的地层不属于活动断裂.深地震反射剖面的中部还揭示了一个近垂直的穹窿状反射异常体,它可能起源于莫霍面,向上,穿过上、下地壳分界面,并延伸至上地壳.穹窿状反射异常体内部反射波视频率随深度增加而降低,在莫霍面附近的壳幔过渡带也出现明显的频率降低、界面扭曲和变形现象,推断它可能是上地幔岩浆上涌到地壳内部的侵入体.结合电磁测深结果可以发现,上地幔热物质的上涌和东、西向拉张可能是形成石家庄—晋县凹陷的动力学机制.探测结果为深入理解石家庄地区的深部地球动力学过程、华北克拉通破坏机制、深浅构造关系和地震构造提供了依据.青藏高原东北部多尺度重力场及其地球动力学意义
摘要：由于受到SN和EW双向挤压作用的综合影响,处于年轻高原与古老地块交接区的青藏高原东北部岩石圈强烈变形,构造活动十分活跃.为了整体性地了解青藏高原东北部重力场和深部动力学机制,本文基于EGM2008全球重力场模型数据,利用小波多尺度分析获得不同尺度的重力异常信息;同时反演了研究区的地壳厚度;通过构建穿越龙门山造山带和西秦岭造山带两条剖面的岩石圈密度结构模型,分析了地壳上地幔内不同介质的分布特征.研究结果表明,青藏高原东北部岩石圈显示十分复杂的塑性体的特征,其重力异常走向多以EW或SSE为主,反映了高原岩石圈物质向东运移的趋势;地壳厚度由西向东逐渐减薄,边缘造山带深部并没有发现＂山根＂痕迹,结合该地区低重力的特征,推测西秦岭—松潘构造结岩石圈发生过大规模地幔流底侵作用;地幔流上涌的动力可能来源于印度板块向欧亚大陆板块俯冲,激发了地幔流体侧向移动,在扬子地台和华北地台附近受到坚硬岩石圈的阻挡而被迫上移,并因此造成龙门山与西秦岭的隆升.2014年鲁甸6.5级地震成因、破裂特征和余震分布特征的库仑应力作用
摘要：根据2014年鲁甸MS6.5地震的区域构造和余震共轭分布特征,本文首先计算了1733年小江断裂带北段M7.75地震,1850年则木河断裂带M7.5地震和1974年马边MS7.1地震对鲁甸MS6.5地震震源机制解两个节面的黏弹性库仑应力作用,结果显示NNW向发生主破裂的包谷垴—小河断裂受到这3次地震,尤其是1850年M7.5地震明显的库仑应力作用,我们认为则木河断裂的高速左旋走滑运动以及7级以上强震的重复发生对于包谷垴—小河断裂的强震孕育和断裂演化方面具有促进作用;然后分析了鲁甸6.5级地震的共轭破裂与余震分布特征,并分别计算了两个共轭破裂面单独破裂对另一破裂面的库仑应力作用,结果显示NEE向破裂促进NNW向破裂的发生,而NNW向破裂后则阻碍了NEE向破裂的进一步发展,最终发展成以NNW向破裂为主的共轭破裂事件;最后计算了共轭破裂产生的库仑应力变化对余震的影响,认为位于NEE向破裂面西侧的余震集中分布主要是由于应力触发作用而形成.1989—1999大同地震序列的隐伏断层研究：库仑应力分析和余震JHD重定位
摘要：1989年到1999年,大同—阳高地区发生了一系列MS≥5的中强地震.本文基于前人对1989年三次MS≥5地震的震源机制反演的结果,通过建立不同断层模型,利用库仑应力方法,计算前震对于主震,以及前震和主震对于余震的库仑应力触发关系,提出了一种可能的破裂模型,即1989年前震沿北西西方向发生左旋破裂,之后主震和余震沿北北东方向发生右旋破裂.根据这种破裂模式计算得出,前震发生后,主震震源处的库仑应力增加了约2×10^5 Pa,余震震源处的库仑应力出现下降;主震发生后,余震处的库仑应力出现回升,最后余震处的库仑应力几乎没有变化.基于大同地震台网的近场观测数据,用JHD（Joint Hypocentral determination）定位方法,对1999年11月1日MS=5.6地震后一个月的余震进行重定位,得到一条走向118°,倾角85°的左旋走滑断层,余震的深度分布在5km至20km范围内,显示该断层是隐伏断层.另外提出对主震震中位置约10km的修正.本文对1989年三次MS≥5地震序列和1999年MS=5.6地震余震空间分布的研究揭示该地区存在两条活跃的共轭隐伏走滑断层（1989年主震的北北东方向和1999年地震的北西西方向）,并且推断已知的大王村断裂和团堡断裂是地下这两条共轭的隐伏走滑断层构造/地震活动在地表的响应.汤加—克马德克俯冲带现今非均匀应力场特征及其动力学意义
摘要：汤加—克马德克俯冲带是太平洋板块向澳大利亚板块俯冲碰撞的动力作用区,是全球俯冲带动力学研究的热点区域.本研究基于EHB地震目录,对汤加—克马德克俯冲带（18.5°S—28.5°S）区域进行平面拟合,得到该范围内俯冲带走向约为196°,倾角约为48°;利用该俯冲带研究区域内Global CMT目录,对不同位置、不同深度进行区域应力张量反演,得到汤加—克马德克俯冲带研究区内精细的应力图像.结果显示：（1）俯冲带浅部（60-300km）应力结构非均匀特征明显,主应力轴倾伏角变化多样,并且最大主压应力轴方位在24°S左右发生明显偏转,我们推测这可能与洋底构造路易斯维尔海链俯冲有关;（2）中部（300-500km）最大主压、主张应力轴由北向南逐渐发生偏转,这可能与由北向南流动的地幔流对俯冲板片产生推挤作用有关,并且这种推挤作用向南逐渐减弱;（3）深部（500-700km）最大主压应力轴沿俯冲方向分布;（4）本文的结果还发现了主俯冲带深部西侧＂偏移＂板片与主俯冲带应力结构不同,表明＂偏移＂板片与主俯冲带是分离的.数值模拟静态应力扰动下的断层失稳：结果分析兼与Dieterich模型和Coulomb模型的对比
摘要：基于大量岩石力学实验,Dieterich和其他研究者（Dieterich,1978;Ruina,1983）首先提出了描述岩石摩擦过程的速率-状态摩擦定律（R-S摩擦定律）.如今R-S摩擦定律已成为研究地震成核等地震演化机制的有效手段.Dieterich（1992,1994）最早提出了描述受静态剪应力扰动后断层失稳时间提前或推后的余震触发机制的解析模型.现在Dieterich模型已经成为解释余震随时间衰减规律的Omori定律等地震观测现象的有力工具.与之相对应,广泛使用的Coulomb应力失稳模型也可以给出断层受到静态剪应力扰动后,断层失稳时间的提前和推后量.Dieterich模型和Coulomb应力失稳模型基于不同的物理方法,所以在进行地震危险性评估时,二者均有各自的局限性.本文利用R-S摩擦定律控制的1-D弹簧-滑块模型,模拟计算了理论地震循环以及在不同静态剪切应力扰动下,失稳时间的提前和推后量的变化情况,然后将计算得到的时间提前和推后量分别与Dieterich模型和Coulomb模型的相应计算结果进行了定量化对比和差异性分析,并给出了相应的解释.数值模拟结果显示,对于R-S摩擦定律在参数不变的条件下,断层模型失稳时间的提前和推后量的大小强烈依赖于静态剪应力扰动的大小和作用时间,而且绝对值相同的正、负向静态剪应力扰动造成的失稳时间的提前和推后量的变化情况并不完全一致.在震后松弛/滑移阶段和闭锁阶段,时间提前和推后量是常数,且随静态剪应力扰动绝对值的增大而增大,两者的比值接近于1.0,这与Dieterich模型和Coulomb模型的结果是一致的,相应的差值小于两模型结果的10%.而在自加速阶段,模拟计算结果则存在与Dieterich模型和Coulomb模型结果不同的特征.首先,在自加速阶段模拟计算结果均偏离Coulomb模型,而且时间提前和推后量的比值小于1.0,相异于Coulomb模型的论断.不过当受到正向静态剪应�三维层状黏弹性半空间中球面SH、P和SV波源自由场
摘要：采用刚度矩阵方法结合Hankel积分变换,求解了层状黏弹性半空间中球面SH、P和SV波的自由波场.首先,在柱坐标系下建立层状黏弹性半空间的反轴对称（柱面SH波）和轴对称（柱面P-SV波）情况精确动力刚度矩阵.进而由Hankel变换将空间域内的球面波展开为波数域内柱面波的叠加,然后将球面波源所在层的上下端面固定,求得固定层内的动力响应和固定端面反力,将固端反力反向施加到层状黏弹性半空间上,采用直接刚度法求得固端反力的动力响应,叠加固定层内和固端反力动力响应,求得波数域内球面波源动力响应.最后由Hankel积分逆变换求得频率-空间域内球面波源自由场,时域结果由傅里叶逆变换求得.文中验证了方法的正确性,并以均匀半空间和基岩上单一土层中球面SH、P和SV波为例分别在频域和时域内进行了数值计算分析.研究表明基岩上单一土层中球面波自由场与均匀半空间情况有着本质差异;基岩上单一土层中球面波位移频谱峰值频率与场地固有频率相对应,基岩面的存在使得基岩上单一土层地表点的位移时程非常复杂,振动持续时间明显增长;阻尼的增大显著降低了动力响应的峰值,同时也显著减少了波在土层的往复次数.
