首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 109 毫秒
1.
龙门山断裂带多参数深部结构成像与地震成因研究   总被引:6,自引:4,他引:2       下载免费PDF全文
通过反演大量的纵、横波地震数据,获得了沿龙门山断裂带及周边区域的深部三维精细结构,结合前人二维大地电磁探测研究成果,提出龙门山断裂带地壳形变与深部速度结构和导电率不均匀性有关,探讨了2008年汶川和2013年芦山地震的诱发和产生与流体侵入及地壳形变的密切关系.本研究发现,2008年汶川地震发生在高速度、高泊松比和低电导率的区域,2013年芦山地震则位于高速度、低泊松比和低电导率的发震层.在上地壳中,四川前陆盆地的低速、低泊松比和低阻异常与松潘一甘孜地块的高速、高泊松比和高阻异常形成了鲜明的对比.在龙门山断裂带下方的两个低速和低阻块体,将龙门山断裂带分成南、中和北三段.我们的研究认为,这两个异常体与来自松潘甘孜地块的下地壳和(或)上地幔的局部熔融或流体侵入到龙门山断裂带的脆弱区有关.基于对汶川和芦山地震的余震分布特征及震源区的地震波速度、泊松比及电阻率参数分析,揭示了龙门山断裂带深部剧烈的地壳形变与流体应力积累对2008年汶川和2013年芦山地震的触发及其地震破裂过程具有重要的控制作用.  相似文献   

2.
运用非连续变形分析法与三维有限元法相结合的方法,以GPS资料作为位移速率和震源机制的约束条件,通过数值模拟研究了青藏高原及其东侧邻区构造地块的运动、变形、相互作用及其与近30年来发生于该区的大地震之间的关系。研究中引入了以应力与摩擦强度的比值定义的断层“失稳危险度”,通过数值模拟计算得到了研究区地壳块体边界断层的失稳危险度分布。结果表明,失稳危险度高的地段与近期该区发生的MS≥7.0地震所在的位置基本一致,其中龙门山断裂带上包括汶川和芦山大地震的发震断层均为失稳危险度最高值地区。计算得到的应变率强度分布图显示,青藏高原东部边缘整条地带均为应变率强度的陡变带,特别是以龙门山断裂带上的陡变最为明显,其西侧应变率强度为东侧的近4倍,而且,这个带位于宽度相同、走向与龙门山断裂带走向相一致的高应变能密度带中,表明这两次大地震前,作为其发震断层的龙门山断裂带已积累了相当高的应变能,失稳危险度高,处于力学上的不稳定状态。模拟计算得到在上地壳层中,2001年昆仑山口西MS8.1地震引起汶川、芦山地震发震断层的库仑破裂应力增加约0.016 MPa,相当于龙门山断裂带约两年的应力积累,也就是说,使汶川、芦山地震发震断层的失稳破裂提前了约两年。 此外,关于2008年汶川MS8.0地震的模拟计算表明,汶川地震的发生也使包括芦山地震发震断层的龙门山断裂带西南段和东昆仑断裂带东南端的库仑破裂应力增大,应变能积累增强,这说明汶川MS8.0地震的发生对已处于失稳危险度较高状态的2013年芦山地震和2017年九寨沟地震发震断层的提前失稳破裂起到了促进作用。   相似文献   

3.
2008年5月12日我国四川省汶川地区发生了震惊世界的MS8.0地震.历史上,同类地震在大陆内部极为罕见.该地震深部构造背景的研究对理解其成因极为重要.本文利用中国地震局地质研究所地震动力学国家重点实验室在川西地区布设的大规模密集流动宽频带地震台阵记录的远震P波波形数据和接收函数非线性反演方法,得到了沿北纬31°线的19个台站下方120 km深度范围内的S波速度结构及台站下方地壳的平均泊松比.该观测剖面穿越了主震区,总长度约为420 km. 我们的结果揭示了川滇地块、松潘-甘孜地块和四川盆地三个不同地块构造差异.上述三个地块的地壳结构特征可以概括为:(1)四川盆地前陆壳幔界面向西侧倾斜并有较为明显的横向变形,地壳厚度存在46~52 km的横向变化,中下地壳S波速度存在横向变化,地壳平均泊松比值较高(0.28~0.31),但在龙门山断裂带附近,显示了坚硬地壳的特征,地壳平均泊松比仅为0.2;(2)松潘-甘孜地块地壳厚度由西侧靠近鲜水河断裂的60 km,向东减薄为52 km,在14~50 km深度范围内存在S波速度2.75~3.15 km/s的楔状低速区,其厚度由西侧的~30 km向东逐渐减薄为~15 km,相应区域的地壳平均泊松比高达0.29~0.31; (3)鲜水河断裂西侧,川滇地块地壳结构相对简单,地壳厚度为58 km,并在26 km深度存在约10 km厚度的高速层,地壳内平均泊松比约为0.25;(4)汶川大震区在12~23 km深度上具有近乎4.0 km/s的S波高速结构,而其下方的地壳为低速结构,地壳平均泊松比0.31~0.32,汶川大震的余震序列主要分布在高速介质区域内. 本文的结果表明松潘-甘孜地块的地壳相对软弱;而且并不存在四川盆地向西侧的俯冲.我们认为在青藏高原东向挤压的长期作用下,四川盆地强硬地壳的阻挡作用可导致松潘-甘孜地块内部蓄积很大的应变能量以及上、下地壳在壳内低速层顶部边界的解耦,在龙门山断裂带附近形成上地壳的铲形逆冲推覆.汶川大地震及其邻近区域所具有的坚硬上地壳和四川盆地的阻挡作用为低应变率下的高强度应力积累创造了必要条件,而松潘-甘孜地块长期变形积累的高应变能构成了孕育汶川大地震的动力来源.  相似文献   

4.
徐晶  邵志刚  刘静  季灵运 《地球物理学报》2017,60(10):4056-4068
本研究基于分层黏弹介质模型,考虑同震位错效应和震后黏滞松弛效应,分析巴颜喀拉地块东端1976年松潘地震序列、2008年汶川8.0级地震、2013年芦山7.0级地震和2017年九寨沟7.0级地震等多次大地震的可能存在的触发关系,计算大地震引起的周边各主要断裂的库仑应力变化.结果显示,1976年松潘地震序列各次地震间关系密切,存在明显的相继触发作用;综合考虑同震和震后效应,汶川8.0级地震对同属于龙门山断裂带的芦山7.0级地震有触发作用,且震后效应影响不可忽略;1976年地震序列,特别是1976年8月16日7.2级地震促进了2017年8月8日九寨沟7.0级地震的发生;汶川地震对九寨沟地震的影响研究中,采用不同的汶川地震同震位错模型,计算结果有差异.综合考虑多次大地震对周边断裂带的影响,龙门山断裂带南段、鲜水河断裂带中南段、平武—青川断裂北段、灌县—安县断裂北段、文县断裂的累积库仑应力增加显著,巴颜喀拉地块东端的东昆仑断裂带东段、迭部—白龙江断裂带西段以及金沙江断裂带库仑应力亦有所增加.综合考虑各重要断裂带已有的大地震危险性分析结果和库仑应力变化计算结果,龙门山断裂带南段、鲜水河断裂带中南段、东昆仑断裂带玛沁—玛曲段和金沙江断裂带的发震紧迫性有所增强,需引起关注.  相似文献   

5.
汶川8.0级地震前区域地壳运动与变形动态过程   总被引:13,自引:6,他引:7       下载免费PDF全文
利用GPS等观测资料研究了2008年汶川8.0级地震前的区域地壳运动与构造变形、应变积累,以及大区域地壳运动微动态变化过程,结合同震位移场分布等分析、讨论了汶川地震前近10年区域地壳变形的表象所反映的大震孕震最后阶段的物理过程.结果表明,发生汶川8.0级地震的龙门山断裂带1999~2004年相对稳定的华南地块无明显的相对运动,2004~2007年有一定程度的相对运动显示,以右旋活动为主,年速率为1.6mm/a.龙门山断裂带西北侧的巴颜喀拉地块东部为右旋剪切为主兼有推挤的大尺度缓慢变形状态,右旋扭动变形率为2.1±0.2(10-8/a),地壳总体缩短变形率为-0.7±0.1(10-8/a).由GPS基准站资料计算的基线时间序列显示,2005年以来大尺度北东向地壳缩短的相对运动明显增强,青藏块体相对华南地块的北东向运动明显增强.对汶川大地震前应变积累的特殊性等问题进行了初步讨论,分析认为,汶川大地震的发生是四川盆地西缘的龙门山断裂带受到其西侧巴颜喀拉地块推挤导致大尺度、长时期、缓慢的地壳应变积累的结果.在发震前不仅龙门山断裂带为显著闭锁状态,且孕震区域可能存在地壳弹性变形趋于极限后难以发生变形的相持阶段,大尺度地壳运动与动力环境的动态变化对汶川大地震的发生有促进影响.  相似文献   

6.
川西龙门山及邻区地壳上地幔远震P波层析成像   总被引:31,自引:13,他引:18       下载免费PDF全文
本文利用川西地震台阵记录到的远震P波走时数据和非线性层析成像算法,获得龙门山地区400 km深度范围内的三维P波速度结构.为了适应川西地区复杂的地质结构,本文的层析成像方法采用了快速行进三维走时计算算法和Tarantola非线性反演算法.我们的结果揭示了川滇地块、松潘-甘孜地块和四川盆地三个不同地块构造差异及该区深部动力学特征.本文的研究表明:1)研究区地壳上地幔P波速度结构具有较为明显的分区特征,松潘-甘孜地块和川滇地块岩石圈速度较低,四川盆地岩石圈速度较高,四川盆地的岩石圈厚度从南250 km向北逐渐减薄至100 km.松潘-甘孜地块上地幔存在地幔上涌的特征.2)川滇地块和四川盆地仅是垂直接触关系,而在龙门山地区四川盆地前缘存在减薄的现象,并伴随松潘-甘孜地块上地幔低速物质有侵入四川盆地岩石圈下方的特征,这显示了四川盆地与松潘-甘孜地块和川滇地块的动力学关系的差异.3)以映秀为界,龙门山断裂带被从松潘-甘孜侵入的低速异常分为南北两段:龙门山南段和龙门山北段,汶川大地震及其余震序列均分布在龙门山断裂带的北段.在青藏高原向东挤压和地幔上涌的双重作用下造成松潘-甘孜地块隆升,由于汶川处于龙门山北段的最南端,应力容易在此集中.这些因素可能是汶川MS8.0地震的基本动力学背景.本文的结果不支持四川盆地的俯冲及层间流动的动力学模型.  相似文献   

7.
汶川Ms8.0地震:地壳上地幔S波速度结构的初步研究   总被引:27,自引:11,他引:16       下载免费PDF全文
2008年5月12日我国四川省汶川地区发生了震惊世界的Ms8.0地震.历史上,同类地震在大陆内部极为罕见.该地震深部构造背景的研究对理解其成因极为重要.本文利用中国地震局地质研究所地震动力学国家重点实验室在川西地区布设的大规模密集流动宽频带地震台阵记录的远震P波波形数据和接收函数非线性反演方法,得到了沿北纬31°线的19个台站下方120 km深度范围内的S泼速度结构及台站下方地壳的平均泊松比.该观测剖面穿越了主震区,总长度约为420 km.我们的结果揭示了川滇地块、松潘-甘孜地块和四川盆地三个不同地块构造差异.上述三个地块的地壳结构特征口J以概括为:(1)四川盆地前陆壳幔界面向西侧倾斜井有较为明显的横向变肜,地壳厚度存在46~52 km的横向变化,中下地壳S波速度存在横向变化,地壳平均泊松比值较高(0.28~0.31),但在龙门山断裂带附近,显示了坚硬地壳的特征,地壳平均泊松比仪为0.2;(2)松潘-甘孜地块地壳厚度由西侧靠近鲜水河断裂的60 km,向东减薄为52 km,在14~50 km深度范围内存在S波速度2.75~3.15 km/s的楔状低速区,其厚度由西侧的~30 km向东逐渐减薄为~15 km,相应区域的地壳平均泊松比高达0.29~0.31;(3)鲜水河断裂西侧,川滇地块地壳结构相对简单,地壳厚度为58 km,并在26 km深度存在约10 km厚度的高速层,地壳内平均泊松比约为0.25;(4)汶川大震区在12~23 km深度上具有近乎4.0 km/s的S波高速结构,而其下方的地壳为低速结构,地壳平均泊松比0.31~0.32,汶川大震的余震序列主要分布在高速介质区域内.本文的结果表明松潘-甘孜地块的地壳相对软弱;而且并不仔在四川盆地向西侧的俯冲,我们认为在青藏高原东向挤压的长期作用下,四川盆地强硬地壳的阻挡作用可导致松潘-甘孜地块内部蓄积很大的应变能量以及上、下地壳存壳内低速层顶部边界的解耦,在龙门山断裂带附近形成上地壳的铲形逆冲推覆.汶川大地震及其邻近区域所具有的坚硬上地壳和四川盆地的阻挡作用为低应变率下的高慢度应力积累创造了必要条件,而松潘-甘孜地块长期变形积累的高应变能构成了孕育汶川大地震的动力来源.  相似文献   

8.
在青藏高原东北缘至鄂尔多斯地块沿玛沁—兰州—靖边剖面进行62个测点的大地电磁观测,采用Robust技术对观测数据进行了处理和张量阻抗分解.分析了视电阻率和阻抗相位曲线、二维偏离度、区域走向.采用RRI二维反演技术进行了资料的反演解释,二维剖面的电性结构显示:(1)玛沁断裂带、兰州深断裂带、马家滩—大水坑断裂带将剖面分为4个电性区块:巴颜喀拉地块、秦祁地块、边界带和鄂尔多斯地块.(2)区块1、2和4的地壳电性结构有类似特点:上地壳为高阻层,下地壳上部为低阻带,下地壳下部到上地幔电阻率随深度逐渐升高.区块3电性成层性差、结构复杂,是现今构造活动较强烈的地区.(3)玛沁断裂带、海原断裂带和罗山—云雾山断裂带为较陡立的超壳断裂带;西秦岭北缘断裂带为壳内断裂带.  相似文献   

9.
2008年汶川8.0级地震前横跨龙门山断裂带的震间形变   总被引:17,自引:4,他引:13       下载免费PDF全文
利用区域GPS和水准测量资料,结合地震构造背景的分析,本文研究2008年汶川8.0级地震前横跨龙门山断裂带地区的震间地壳形变,探讨引起这种形变的活动构造与动力学模式,并由此认识汶川地震的孕育与成因机制.主要结果表明:1997~2007年期间,自龙门山断裂带中段朝北西约230 km的地带内存在垂直于断裂的水平缩短变形、以及平行于断裂的水平右旋剪切变形,缩短率为1.3×10-8/a (即:0.013 mm/km/a),角变形速率为2.6×10-8/a;同一地带在1975~1997年期间还表现出垂直上隆变形,上隆速率在龙门山前山断裂与中央断裂之间仅0.6 mm/a,而至龙门山后山断裂及其以西达2~3 mm/a.这些反映了在汶川地震之前至少10~30余年,龙门山断裂带中段的前山与中央断裂业已闭锁、并伴有应变积累.造成这种形变的主要原因是:以壳内的低速层为“解耦”带,巴颜喀拉地块上地壳朝南东的水平运动在四川盆地西缘受到华南地块的阻挡、转换成龙门山断裂带中段的逆冲运动;由于该断裂段的震间闭锁,致使西侧的巴颜喀拉地块的上地壳发生横向缩短以及平行断裂的右旋剪切变形.然而,龙门山断裂带北段在1997~2007年期间除了有大约0.9 mm/a的右旋剪切变形外,横向的缩短变形极微弱,这可能与该断裂段西侧的岷江、虎牙、龙日坝等断裂带吸收了巴颜喀拉地块朝东水平运动的大部分有关.另外,汶川地震前,横跨龙门山断裂带中段与北段的地壳形变特征的差异,与汶川地震时能量释放的空间分布吻合.  相似文献   

10.
2008年5月12日汶川MW7.9地震发生在龙门山断裂带。龙门山断裂带及其邻域的地壳上地幔三维速度结构的研究对于理解汶川大地震的动力学背景具有重要的意义。2006年10月至2009年10月,在国家重大基础研究项目(973)的支持下,中国地震局地质研究所地震动力学国家重点实验室在川西地区(26°~32°N,100°~105°E)布设了由297台宽频带数字地震仪组成的流动观测台阵(简称川西台阵)。根据川西台阵记录的环境噪声和远震波形数据,利用噪声成像技术和接收函数方法,我们研究了川西地区(29°~32°N,100°~105°E)地壳上地幔100km深度范围内的三维S波速度结构。本文得到的结果为研究川西高原和四川盆地的地壳结构提供了新的高分辨率观测证据。我们的结果表明:1)观测台阵覆盖的川滇地块、松潘-甘孜地块和四川盆地的地壳上地幔S波速度结构具有显着差异,龙门山断裂和鲜水河断裂带,作为地块间的边界断裂带,对两侧地壳结构具有明显的控制作用。2)观测台阵覆盖区域的地壳厚度存在明显差异,川滇地块的地壳厚度为60~64km,松潘-甘孜地块的地壳厚度为52~56km,四川盆地前陆的地壳厚度为46~52km,沿龙门山断裂带松潘-甘孜地块和四川盆地形成镶嵌结构,汶川地震震中处南北两侧的壳幔边界存在约6km的断错。3)四川盆地前陆低速特征表明相应区域存在厚度8~10km的沉积盖层,松潘-甘孜地块和川滇地块的中下地壳具有大面积分布的S波低速区,松潘-甘孜地块地壳平均泊松比高达0.29~0.31,汶川地震余震绝大多数分布在低速区上方的高速介质区域内,而四川盆地的中下地壳呈现整体性的高速特征,以汶川地震的震中为界,龙门山断裂带北段和南段的S波速度结构显示了明显的速度分段特征,其北段的S波速度总体上高于南段。4)本文给出的研究区地壳三维S波速度结构表明,川西高原中下地壳较为软弱,而四川盆地中下地壳的强度应明显高于松潘-甘孜地块,意味着四川盆地坚硬中下地壳可以阻挡松潘-甘孜地块向东的逃逸;另一方面,川西高原和川滇地块的中下地壳虽然均存在大面积的S波低速区,但松潘-甘孜地块内的地壳速度结构相对来说较为复杂,并形成了高、低速相间的结构特征,表明在四川盆地的阻挡作用下,该地块形成了折皱变形的结构。5)与S波低速区相应,松潘-甘孜地块和川滇地块中下地壳应处于部分熔融的状态,这对该区域存在中下地壳通道流(Channelflow)的推断是一个支持;但是,松潘-甘孜地块内是否存在中下地壳通道流仍有待进一步的深入研究。6)接收函数方位各向异性的偏振分析表明,以汶川地震震中为界,龙门山断裂西南侧处于挤压状态,而其东北侧的主压应力方向与断层走向大体平行,推断先存应力场可能驱动了汶川地震逆冲破裂之后沿龙门山断裂向北东方向的走滑破裂。  相似文献   

11.
12.
The determination of the coefficient of viscosity of eruptive products gives useful elements to the knowledge of possible variations of composition and physical conditions of the magmas.  相似文献   

13.
14.
Geomagnetism and Aeronomy - The magnetic evolution of the active region of the Sun at the level of the photosphere has been studied. Magnetic synoptic maps of the 2007 Kitt Peak Observatory were...  相似文献   

15.
Until recently, the existing data prevented the geophysicists from accurately dating the Bysy-Yuryakh stratum, which outcrops in the middle reach of the Kotuy River, constraining the time of its formation to a wide interval from the end of the Late Cambrian to the beginning of the Silurian. The obtained paleomagnetic data unambiguously correlate the Bysy-Yuryakh stratum to the Nyaian regional stage and constrain its formation, at least a considerable part of it, by the Tremadocian. This result perfectly agrees with the data on the Bysy-Yuryakh conodonts studied in this work and yields a spectacular example of the successful application of paleomagnetic studies in solving important tasks of stratigraphy and, correspondingly, petroleum geology. Within the Bysy-Yuryakh stratum, we revealed a large normal-polarity interval corresponding to the long (>1 Ma) period when the geomagnetic reversals were absent. This result, in combination with the data for the Tremadocian and Middle–Upper Cambrian sequences of the other regions, indicates that (1) the rate of occurrence of the geomagnetic reversals on the eve of the Ordovician Moyero superchron of reversed polarity was at most one reversal per Ma; (2) the superchron does not switch on instantaneously but is preceded by a certain gradual change in the operation conditions of the dynamo mechanism which, inter alia, manifests itself by the reduction of the frequency of geomagnetic reversals with the approach of the superchron. This finding supports the views according to which a process preparing the establishment of the superchrons takes place at the core–mantle boundary.  相似文献   

16.
Summary The results of two different methods of the determination of the non-linearity of the gravimeter scale are compared; the economy is discussed as well.  相似文献   

17.
Atmospheric temperature and humidity fields as well as information on other meteorological parameters are nowadays retrieved from radiance measurements recorded by operational meteorological satellites. Up to now, the inversion procedures used only take into account crude information on the topography of the Earth’s surface. However, the applied radiative transfer codes have to consider the Earth’s surface as the lower boundary of the atmospheric model and, therefore, need a more precise mean elevation and a classification of the roughness of the Earth’s surface. The influence of the topography of the Earth surface on retrieved temperature profiles is studied by using a physico-statistical inversion method. An objective analysis is made of the more precise mean elevation and derivation of roughness parameters using a new high-resolution digital elevation model (DEM) with a resolution of 500 m×500 m. By means of a geomorphological process and a newly developed topography rejection test, areas with a high surface roughness are localized and singled out. The influence of topography on the retrieved temperature profiles is illustrated by case studies. Changes are found predominantly in areas with a high variation of topography. Using the new high-resolution DEM and the topography rejection test, the geographical position of the calculated temperature profiles tends to be shifted towards areas with a small vertical variation of topography. The mean elevation determined by the new elevation model better characterizes the area observed. Hence, the temperature profiles can be calculated down to lower atmospheric levels. Furthermore, a guess profile better describing the atmospheric situation is selected by the more precise elevation. In addition, the temperature profiles obtained near the coast are improved considerably by the more precise determination of the surface property ‘sea’ and ‘land,’ respectively. Integration of an independent physical information such as topography leads, on average, to a slight improvement of the results of the physico-statistical inversion procedure. In some cases, however, significant improvements have been achieved regarding the desired accuracy of temperature profiles of the order of 1 K. In future, the spatial resolution of new high-resolution sounding instrumentation on the next generation of operational meteorological satellites will be increased. To exploit the resolving power of this new instrumentation, the different variation of the topography of the Earth surface, especially in regions with a high variation of topography, can be taken into account more precisely by using a high-resolution DEM.  相似文献   

18.
Summary A simple model is proposed suitable for studying the effect of the ground layer of the atmosphere, polluted by aerosol, on the albedo. This model is founded on solving the equation of transfer of radiative energy. The numerical results are discussed, particular attention being paid to the analysis of errors due to neglecting the multiple reflection of solar radiation on the aerosol particles. A method which would also include the multiple reflection is proposed, and the conditions under which the presence of the aerosol is responsible for an increase or decrease of the solar radiation balance on the Earth's surface, are analysed.  相似文献   

19.
Резюме Дается вывод формул для вычислений длины хорды, длины и азимутов нормаляных сечений и длины геодезической линии эллипсоида. Применяются тригонометрические функции данных величин и постоянные принятого ргфгргнц-эллипсоида, что является целесообразным при использовании вычислителQjных машин, особенно автоматов. Рабочие формулы рекомендуются в форме (8–11).

Dedicated to Professor František Fiala on His 85th Birthday

Address: Veveři 95, Brno.  相似文献   

20.
Summary The maintenance of the axisymmetric component of the flow in the atmosphere is investigated by means of a steady-state, quasi-geostrophic formulation of the meteorological equations. It is shown that the meridional variations in the time-averaged axisymmetric variables can be expressed as the sum of three contributions, one being due to the eddy heat transport, another to the eddy momentum transport, and a third to the convective-radiative equilibrium temperature which enters the problem through the specification of a Newtonian form of diabatic heating. The contributions by the large scale eddies are evaluated through the use of observed values for the eddy heat and momentum transports.The contributions from each of the three forcing mechanisms to the temperature and zonal wind fields are invstigated individually and found to be of about equal importance. The sum of the three contributions are also presented for the temperature, the zonal wind, the stream function associated with the mean meridional circulation and the corresponding vertical motion. Although the results fail to reproduce the main observed features of the lower stratosphere, they are found to be in good agreement with observations in the middle latitude troposphere. At any pressure level, for example, the computed mean zonal wind has a jet-like profile and the axis of the jet is found to slope to the south with height, as observed in the atmosphere.Based in part on a thesis submitted by the first author as partial fulfillment of the requirements for the Ph.D. degree at the University of Michigan. — Publication No. 194 from the Department of Meteorology and Oceanography, The University of Michigan.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号