DreamBlooms/HansPub-15000-to-45000 huggingface dataset:
system: [SYSTEM]

user: 2. 资料整理与分析处理方法 2.1. 地震目录的选取 地震地热说的研究原理与方法，完全依赖于具有高精度震源深度资料的地震目录。为此，本文采用从公网上获取并整理的2套地震目录 [ 12 ]。其中包括： 图1. 19-3号伊朗地震柱分支的地理位置图 图2. 伊朗地震柱分支的地震与火山分布(据ANSS地震目录，1963~2019.5，M ≥ 4.0；火山资料据smithsonian institution [ 15 ] ) ANSS复合地震目录(ComCat)，来自北加利福利亚地震数据中心(Northern California Earthquake Data Center)，时段1963~2019.5，共计包括全球4级以上地震记录456，888条。网址包括https://earthquake.usgs.gov/earthquakes/search/ (2004年以前)和 https://earthquake.usgs.gov/earthquakes/search/#site-sectionnav (2004年及以后)。 EMSC地震目录来自欧洲-地中海地震中心(The European-Mediterranean Seismological Centre)，时段2004.10~2019.5，主要包含经度−40˚至110˚E地区，其中0~39˚E地区震级1级以上，其它地区2-3级以上，共计241，545个地震记录。网址https://www.emsc-csem.org/Earthquake/?filter=yes 2.2. 关于震源深度 关于地震震源深度的精度，可以从ANSS地震目录获得 [ 12 ]。ANSS目录包含全球4级以上地震，大量地震的深度测定误差均在几公里的量级，个别偏远地区深部地震的深度误差可达20多公里。以上2套地震目录中相同地震的震源深度基本吻合。我们所需要得到的并非是震源的绝对深度，而是要从理论和方法上将0~650 km深度的地震活动合理地离散开来，即相对震源深度。 2.3. 地震目录的适用性分析 由2个地震目录所采集的研究区地震集分档统计见表1。由表1可见： 1) 两个地震目录所描述的震源深度，4级以上的地震基本一致。4级以下地震的最大深度可达384 km。 2) ANSS目录的4级以上地震数量偏少，EMSC目录3级以下地震有所偏置，只能两套地震目录彼此参照。 3) 鉴于目前地震目录的科技水平，地震地热说当前最为关注的是未来壳内强震活动的时间和地点，暂不考虑能量关系，因此不对相关地震集的震级-频度关系做过分要求。 Table 1 Statistics of the iranian branch of seismic cone tectonic No.19-3 from earthquake catalogue 地震目录 震级下限Mmin 震级上限Mmax 频次 最大震源深度Hmax (km) ANSS目录 4 7.7 8612 183 5 1147 164 6 94 80 7 12 80 EMSC目录 1 7.8 33，632 384 2 32，757 384 3 18，037 384 4 3535 250 5 321 148 6 23 87 7 4 87 表1. 由不同地震目录所采集的伊朗地震柱分支地震统计表 2.4. 资料分析处理方法 目前，地震预测首要解决的问题是地震的成因到底是在地球的外部、地球的地表，还是地球的内部？地震最为本质的前兆信息到底是什么？为此，地震地热说依赖地震柱构造和地震子柱构造为平台，发现了大量与地震预测相关的构造信息。 地震地热说的分析方法，首先按照地震柱构造定义选取研究区的地震目录，然后采用地震与火山活动的平面分布、三维立体分布和震源深度时序分布分别研究地震柱构造的总体活动特征以及区内可分离的地震子柱构造的活动特征，总结以往事件的前兆信息并为后续事件的预测提供工作方法与经验。在原始数据同等精度的前提下，本方法的工作结果是唯一的，可重复的。
system: [SYSTEM]

user: 5. 讨论 上世纪70年代，“地震是可以预测、可以预防的”，感动了很多中国人；从上世纪末起，地震却变成了“不可知论”。但是，人们相信迟早还会回到“地震可以预测”的轨道。这就是人类认识自然的波浪式前进、螺旋式上升的规律。自然辩证法则，任何力量也阻挡不住。 本文采用地震地热说的理论与方法研究壳内强震和火山活动与地震柱构造的关系，并将地震柱构造与地质力学的构造体系相结合，产生了一种特殊的效果。在大陆漂移的大背景下，构造体系某些部位会应力集中，应变积累，造成地表构造的某些薄弱环节，当地震柱构造由于壳下地震活动积累的能量达到饱和状态时，就会冲破这样的薄弱环节爆发地震或者火山。图4显示，越是靠近构造体系应力容易集中的部位，强震与火山越是活跃，相应的地震子柱构造也越发强势。所以，壳内强震与火山喷发是地震柱构造释放柱体能量的两种不同方式，断层，褶皱，构造体系，以及火山锥，都只是地震柱构造释放能量的工具，因为大自然的运动总是以最省力的方式来进行。 地震柱构造的运动表现为地幔深部热能的节律性上升运动。壳下地震活动只是热能上升的一种表征，还应该包括热能的传导、辐射与对流。然而，后者目前是无法检测到的，因此抓住壳下地震活动这个表征性符号便成为当前唯一的可能。近几年的研究经验表明，地震柱构造内的壳下地震投影到地表的密集区与近地表地质构造的联合分析，以及壳下地震震源深度时序图的活动趋势，可以为未来强震或者火山活动的地域与时间域提供重要信息，给人以希冀 [ 5 ] [ 6 ] [ 17 ] [ 38 ]。 自从有了地震柱构造的概念(Seismic Cone Tectonic)和地震地热说(Seismo-Geothermics)的理论分析，采取地震活动的平面分布和三维立体分布图像，将强震的震源体和火山喷发的管道剥离为一个个倒立的圆锥体，通过地幔深部加热，以壳下地震活动的表征方式源源不断地为地表灾变提供能量，自下而上，逐层驱动，最后推动地表已有构造的薄弱环节突发强震和/或火山。如此，地震与火山的预测似乎是有希望的。
system: [SYSTEM]

user: 6. 结论 本文依据ANSS和EMSC等地震目录资料，采用地震地热说的原理和方法，对19号地中海地震柱的19-3号伊朗地震柱分支进行数据处理，做出了伊朗分支的总体地震、壳下地震、壳内强震和火山活动的平面分布图像、三维立体图像、壳下地震时序图，以及它们与伊朗山字型构造体系关系的理论解释，可以得出如下一些结论： 1) 伊朗地震柱分支构造，总体上满足地震柱构造的概念，可以按照地震地热说的基本理论和方法探索壳内强震与火山活动的预测途径。然而，由于伊朗地震柱分支构造与伊朗山字型构造体系的叠加，显得格外复杂，因此需要划分成若干地震子柱构造对于各个壳内强震与火山活动的地点做专门研究。 2) 地震子柱构造是能够独立控制自身影响区内的壳内强震或火山活动规律及其与壳下地震关系的构造实体，是地震柱构造内的最小地质构造单元。本文所划分的6个主要地震子柱构造中，除I04地震子柱构造的壳下地震深度较浅外，其余深度均在100 km以上，而且其立体图像基本满足地震柱构造的概念要求，因此对各个地震子柱构造的描述是基本可靠的。 3) 本文依据地震子柱构造的时序图对各相关构造的未来活动趋势进行了预测，并就I01和I04地震子柱构造未来的强震活动给出了壳下地震活动具有警示意义的趋势虚线，值得关注。 4) 本文将地震地热说的地震柱构造与地质力学的构造体系相结合的研究思路，是对地震地热说理论与方法的重大拓展。 本文研究区内4级以下的地震缺失较多，因而只是初步研究。如果拥有2级以上地震的完整目录，研究就可以更加深入，研究结果也更加具有说服力。然而，本文研究方法面临的最大困难在于具有高精度震源深度资料的地震目录的获取。本文所选用的两套地震目录对于研究区的监测能力都是不够理想的。世界上也有一些地方，尽管拥有长达数十年的地震目录资料，可是其壳内地震与壳下地震的分辨率偏低，同样不能适用于本文的研究方法。所以，本文研究方法的改善，依赖于全球24个地震柱构造范围内所有社会力量的共同努力。这取决于相关国家和地区的地震研究机构对于提高壳内地震与壳下地震分辨率重要性的理解与实践。
system: [SYSTEM]

user: 摘 要 随着锰矿区的开采，越来越多的土壤受到了锰污染，关于锰毒和植物抗锰机制的研究对治理锰污染土壤具有重要的意义。本文综述了高浓度锰对植物不同组织生长发育和生理生化过程的毒害作用，并重点阐述了植物响应高浓度锰的策略：通过植物根系的限制吸收、锰离子的外排、金属螯合作用、细胞内区域化等被动逃避策略和启动抗氧化系统，合成锰的吸收、转运有关的关键蛋白质和基因表达等主动耐受适应策略。最后，本文综合以上研究进展，并提出了一些具有参考价值的展望，通过对植物抗锰机理进行不断深入的研究，以实现利用植物修复技术治理重金属土壤污染，为提高农作物对重金属胁迫的抗性提供了理论基础。 关键词 :锰胁迫，区域化，螯合，抗氧化系统，锰转运蛋白 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 随着工农业发展，土壤重金属污染逐渐严重，且其污染具有隐蔽性、持久性和累积性，有害重金属通过土壤–作物系统进入食物链严重威胁着动物和人体的健康 [ 1 ]。土壤中锰含量超标主要源自于矿区开采，据统计全世界锰矿石开采量每年高达600万吨 [ 2 ]。再者，土壤中锰本是惰性金属，由于错误的施肥习惯，忽视有机肥的应用，有机质严重缺乏，导致土壤酸化，使本为惰性的锰变成游离态的锰，活性加大，增加了植物吸收锰的机率 [ 3 ]。推进锰污染治理是全面贯彻落实习近平生态文明思想，推进生态优先、绿色发展的重大实践，是打好打赢污染防治攻坚战的必然要求。因此，锰污染区域必须经过特殊的修复技术，并进行风险评价和连续检测，确保了有害重金属不会在植物内富集并通过食物链转移而危害其他生物健康，才可用于农业生产 [ 4 ]，其中，植物修复技术具有成本低、不破坏土壤组成、重金属可循环利用和美化环境等特点，而利用超富集植物对重金属污染土壤进行修复是目前植物修复的主要研究热点。 锰是植物生长发育的一种必需元素，在组成叶绿体结构，作为水光解的催化剂和调节氧化还原电位等方面发挥了重要的生理作用 [ 5 ]。然而，当土壤中锰含量过高，植物会发生锰中毒。在亚热带和热带酸性土壤地区，锰毒害是仅次于铝毒害影响植物生长最主要的限制因素。首先，过高含量的锰会快速诱导活性氧物质含量升高而造成氧化胁迫，会通过破坏类囊体结构和光合电子传递链而严重影响光合系统而降低光合效率，最终导致叶片枯萎和坏死 [ 6 ] [ 7 ]。其次，植物主根和侧根的生长发育也受到锰毒害的强烈抑制，由于根长期暴露在高浓度锰的土壤中，会出现褐化和断裂的症状 [ 8 ]。再次，在微观层面上，高浓度的锰还会使脂质、蛋白质和核酸等生命大分子物质发生降解，最终破坏细胞代谢导致细胞死亡 [ 9 ]。总之，锰毒害从各个层面上对植物的生长发育造成了如此严重的影响，应该引起广大研究者们的重视和关注。 面对高锰胁迫对植物的影响和危害，植物并非消极被动，而会通过一系列生存策略来响应高锰胁迫。植物抵御对重金属胁迫具有两条基本途径；一是重金属排斥性，即重金属被植物吸收后又被排出体外，或者重金属在植物体内的运输受到阻碍；另一途径是重金属积累，但可自身解毒，即重金属在植物体内以不具生物活性的解毒形式存在，如结合到细胞壁上、离子主动运输进入液泡、与有机酸或某些蛋白质的络合等 [ 10 ]。本文主要综述了植物以限制吸收，细胞内区域化和螯合高浓度的锰离子来实现逃避高锰胁迫的目的，通过过氧化酶清除锰胁迫造成的积累的活性氧物质，表达锰运输蛋白进行转移和外排体内高浓度的锰离子，从而增强对高锰胁迫的适应性和耐受性，并对今后植物响应锰胁迫的研究方向提出了展望。
system: [SYSTEM]

user: 2. 植物被动逃避高锰胁迫 2.1. 植物限制吸收 锰通常富含于土壤中，所以锰通过植物的根吸收进入，而对植物有效的锰形式有三种：直接以Mn 2+ 形式存在的水溶性锰和交换态锰，以及以高价锰氧化物存在的易被还原为有效锰的易还原态锰 [ 11 ]。植物在面对锰毒害时，可通过将锰氧化成非有效价态的方式来降低植物对锰的吸收和运输，从而减少锰毒害，更增加了锰耐性。水稻根能将二价锰离子氧化成非有效性的四价锰化合物，从而限制水稻吸收锰 [ 12 ]。淹水土壤中的根可通过形成明显可见的锰氧化物胶膜的方式来有效地阻止植物过量吸收二价锰离子 [ 13 ]。另外，植物还可以通过多种方式来阻止Mn 2+ 进入植物体内。赵中秋等发现菌根尤其是外生菌根能够显著提高植物的重金属耐受性，除此之外植物的根系还可以分泌改变根际的pH值和氧化还原状态的物质以及根际微生物以减少或阻止有害重金属进入植物 [ 14 ]。在高浓度锰处理下，氧化菌的活性增加使锰被氧化成非有效价态，结果减少了植物对锰吸收，增强了耐受性，因此含菌根植物如大豆的锰毒症状数天后缓解了 [ 15 ] [ 16 ]。 2.2. 植物外排 当植物体内的重金属达到一定量后，有些植物能够将多余的金属排出体外，进而降低毒害。白羽扁豆的锰只能在老叶中积累，不能移动；且根际土壤中含有较高浓度的锰，说明根能通过某种机制将锰排出体外 [ 17 ]。在篦齿眼子菜(Potamogeton pectinatus L)中，当锰累积到上限时，Mn 2+ 被氧化为Mn 4+ ，通过叶片衰老排出体外 [ 18 ]。角菱(Trapa natans)也可通过叶片表皮毛分泌锰-酚类螯合物来解毒 [ 19 ]。因此，不同的植物面对锰毒害，其外排方式具有特异性。 2.3. 区域化存储 区域化作用是指植物将有毒的金属离子转运和贮存在特定的组织器官或细胞的特定区域，降低代谢活性区域的金属离子浓度。这是一种提高植物耐金属能力的重要机制 [ 20 ]。通过对商陆和水蓼研究发现，叶片是锰的主要累积器官，而叶的边缘和表皮层是锰积累和解毒的主要组织 [ 17 ] [ 18 ]。在叶片中，细胞质可溶性部分锰含量最高，而液泡中锰含量占叶片细胞的95%，可推测液泡是叶片中过量锰的积累细胞器 [ 21 ]。锰可以在液泡中固存，并被转移到内质网和高尔基体中，P型ATP酶家族的AtECA 3和CDF锰转运体AtMTP 11可通过高尔基装置参与锰的解毒作用 [ 22 ]。在胞囊藻PCC 6803中，位于细胞质膜和类囊体中的转运体SynPAM 71可以介导锰从细胞质出口转移到细胞质周围和腔室，从而将锰以与外膜结合或者储存到专门的锰储存蛋白(如MncA)的方式储存在不干扰细胞质代谢过程的地方，提高植物对锰的耐受性，在一定范围内维持锰的稳态 [ 23 ]。根的细胞壁也是重金属的贮存库，因为细胞壁上有大量的阳离子配位，所以锰能大量结合在根的细胞壁上，只有当金属与细胞壁结合达到饱和时，多余的锰才会进入到细胞质中，这会导致锰向地上部分转移的能力减弱，这可能也是植物耐锰的一种策略 [ 17 ]。在锰富集植物中，锰主要积累在液泡中；在非富集植物中，锰主要积累在根的细胞壁，因此，不同锰富集能力的植物具有不同的存储区域 [ 24 ]。 2.4. 螯合作用 植物通过螯合作用形成的各种锰–螯合分子可以降低锰的毒性，对于植物中锰的吸收，转运，外排也有直接影响，是植物应对锰胁迫的重要解毒和耐受机制之一。植物根际可以分泌有机酸到土壤中，在体外螯合锰以降低根系周围锰浓度，减少根对锰的吸收。重金属会与有机酸络合在液泡中以降低其植物毒性，从而不影响细胞的正常代谢活动 [ 17 ]。在商陆中锰主要与草酸络合，研究发现草酸与锰结合最为稳定 [ 24 ] [ 25 ]。在高浓度的锰处理下，贵州风毛菊根部通过增加苹果酸脱氢酶介导苹果酸合成和渗出量，以提高对锰毒害的耐受性 [ 26 ]。豇豆中80%的根际锰离子都能与柠檬酸结合 [ 27 ]。此外，植物体内的螯合剂还通过螯合作用形成各种Mn–螯合物大多通过金属转运蛋白转移到叶片的表皮层、叶肉细胞的液泡和细胞壁中，从而减轻锰毒对植物代谢活跃部位的伤害。有些植物形成的Mn–螯合物会排出体外使细胞免受锰毒。例如角菱中锰与酚类化合物螯合后可以通过特化的结构毛状体分泌到体外，减少植物中锰的积累 [ 28 ]。 金属硫蛋白(MTs)和植物螯合肽(PCs)是目前研究较为深入的螯合剂。MTs是一类含有较多半胱氨酸残基、分子质量较小的金属离子结合蛋白 [ 29 ]。MTs中富含大量巯基有利于它与重金属的结合。美洲商陆根部胱氨酸含量较高，而胱氨酸中含有活性巯基，且根部总巯基含量较高，巯基有利于与重金属结合，因此推测美洲商陆的极强环境适应性可能与巯基有关 [ 22 ]。大米和西瓜中的重组MTs可以有效地清除巯基自由基，有利于减轻高锰胁迫带来的氧化损伤 [ 30 ]。 PCs是一类富含Cys的肽，在自然界广泛存在，由谷胱甘肽通过植物螯合肽合成酶而催化合成。PCs中巯基与金属离子螯合，通过Hmt1，Hmt2膜转运蛋白、H + /Cd 2+ 逆向转运蛋白运输到液泡，实现对金属离子的固定 [ 31 ]。吴惠芳等研究发现龙葵和小飞蓬中PCs和谷胱甘肽含量随锰浓度增加而显著升高(p < 0.05) [ 32 ]。另外，有机酸也是限制植物体内锰浓度的重要螯合剂之一。薛生国等通过FTIR法研究发现，酸模叶蓼根和叶分别在2 920 cm −1 和1 630 cm −1 处表现为分泌有机酸不断螯合植物体内的锰，这表明有机酸与锰螯合可能是锰超累积植物叶片锰累积和解毒的重要机制 [ 33 ]。
system: [SYSTEM]

user: 4. 结论和展望 综上所述，锰对植物生长发育多方面和从整体植株，组织，到细胞分子等不同水平都产生了不同程度的影响和危害，而且锰毒害的机理十分的复杂且相互联系。植物在面对锰毒害时首先选择逃避策略，经过长期的适应，也形成了独特的耐受高锰环境的策略。从宏观方面讲，根系限制吸收和外排；从微观方面讲，细胞内的区域化存储、螯合、抗氧化和相关的转运蛋白及基因(图1)。植物发展出了许多耐性机制，但这些机制都具有一定的局限性，目前还没发现一套存在于植物体内的普遍性的耐锰机制。根据目前锰对植物的毒害和植物耐锰机制的研究现状，以后的研究方向可以注重以下几个方面： 图1. 植物对高锰胁迫的响应策略 [ 11 ] 1) 细胞壁上耐锰代谢物质的研究。细胞壁是细胞的第一道屏障，因此物质进出细胞首先接触的就是细胞壁。其中肯定有一些耐锰的机制以及涉及到的一些物质如蛋白和基因等。植物在体内的存储区域除了液泡和细胞壁之外，是否可以储存在其他细胞器中，比如叶绿体，毕竟植物的光合作用需要锰的参与。 2) 锰和其他离子可以相互作用，影响植物对锰吸收，其中的具体原理是怎样的还未研究清楚；植物螯合肽能够缓解重金属的毒害作用，但它的高剂量表达需Cd的激活。如何利用它增强植物对锰的耐受和积累等问题目前还有待进一步研究。黎晓峰等报道了氮源影响锰的吸收情况 [ 61 ]，但具体影响锰毒的机制还不清楚。例如氮素是否影响锰在植物中的有效形态，如何影响锰在植物体的转运等等，都需要进一步研究。 3) 关于特异性重金属转运蛋白和相关调控基因的研究。目前，许多研究者已经发现了一些锰转运蛋白和基因。但这些转运蛋白大都可以转运多种金属离子 [ 18 ]。转运蛋白基因的拓扑结构、表达机制，转运蛋白的具体转运机制以及蛋白之间的协同作用可以作为接下来的研究重点。另外，动物对锰胁迫的响应机制与植物的有所不同，将动物体内有关的耐锰基因或重金属蛋白转入植物体内也是有重要的理论意义和潜在的实用价值。因此，加强对锰转运蛋白基因家族功能以及调控机制的研究，为培育耐受高锰胁迫的优良农作物提供了基因资源，奠定了理论基础。
system: [SYSTEM]

user: 摘 要 为考察抑郁在儿童期虐待与自伤行为之间的中介作用及其性别差异，本研究采用儿童期虐待问卷、流调中心抑郁量表和青少年自我伤害行为量表在线调查了312名本科生。结果发现：1) 儿童期虐待、抑郁和自伤行为呈显著正相关；2) 在总样本中，儿童期虐待、抑郁正向预测自伤行为，抑郁在二者间起部分中介作用；3) 区分群组，抑郁的中介效应在男、女生中均不显著，儿童期虐待对自伤行为的路径系数存在显著性别差异。研究结果提示：儿童期虐待、抑郁和自伤行为密切联系，三者间的相互作用机制存在性别差异。今后研究需就变量间的相互作用机制及性别差异做进一步验证，这对自伤行为的预防和干预具有重要指导意义。 关键词 :儿童期虐待，抑郁，自伤行为，性别差异，大学生 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 自伤行为是指不以自杀为目的的直接地、蓄意地伤害自己躯体的行为，包括割伤自己、抓伤自己、打自己等多种行式(Ross & Heath， 2002; 江光荣等，2011)。自伤行为不仅存在于精神疾病、人格障碍患者人群中，也广泛存在于普通人群中(Winchel & Stanley， 1991)。自伤行为和自杀行为有着相同的影响因素，自伤行为是自杀行为的风险因素。Hawton等人基于对11，583名自伤者长达20年的追踪，报告自伤行为者第1年的自杀风险性为0.7%，是普通人群的66倍，且自杀风险会随时间发展而升高(Hawton， Zahl， & Weatherall， 2003)。自伤行为因其普遍性和严重性而受到心理、精神卫生和公共卫生领域的广泛关注。青少年阶段是自伤行为的高发阶段，国内最近一项纳入23项研究且包括73，677名大学生的Meta研究报告我国大学生的自伤行为检出率高达16.6% (潘珍等，2016)，且国外调查表明青少年自伤发生率还存在逐年上升趋势(Fleming et al.， 2013)。因此，研究青少年自伤行为的影响因素及产生机制对预防和干预大学生自伤、自杀行为具有十分重要的意义。 在众多的自伤行为影响因素中，童年期创伤经历，特别是儿童期虐待，被认为是自伤行为的重要环境诱发因素。家庭环境对个人成长的影响是持续的、长期的，儿童期虐待经历的影响也是深刻的。许多心理学家认为，儿童期虐待经历对儿童的伤害并不会消除，甚至会在个体成年后以其他不良形式表现出来，如抑郁、攻击、物质滥用、犯罪、自伤等(Dube et al.， 2003)。实证研究证实儿童期虐待经历与自伤行为密切联系，儿童期虐待经历显著预测之后的自伤行为(Bakken & Gunter， 2012; Blasczykschiep & Jaworskaandryszewska， 2014; 苏静等，2015)。苏静等人的研究证实各种形式的童年期虐待经历，如：性、情感和躯体虐待，均会增加中学生自伤行为的发生(苏静等，2015)。本研究将利用结构方程模型统合各种形式儿童期虐待综合分析儿童期虐待与自伤行为之间的关系。据此，本研究提出第一个假设：儿童期虐待能够预测自伤行为。 以往研究多停留在儿童期虐待与自伤行为的关系上，较少有研究就它们之间的发生机制进行实证探索。受虐待儿童长期生存在被虐待环境下，常常会感觉到一种无能为力感和无价值感，长此形成的认知缺陷，会使他们具有抑郁易感性特点。大量研究证实儿童期虐待经历对青少年和成年阶段的抑郁有着重要影响(Zlotnick et al.， 1995; 于增艳等，2017)。另外，研究证实有自伤行为者会比无自伤行为者表现出更多的负性情绪管理缺陷(Gratz & Roemer， 2008; Zielinski et al.， 2017; 唐杰等，2014)和情绪问题(吴晓瑞等，2016；谢唯等，2017)。自伤行为的体验回避模型认为：自伤行为是外在情景事件引发了个体的厌恶情绪，个体为了逃脱或缓解不愉快的情绪体验而采取的情绪调节策略(Chapman， Gratz， & Brown， 2006)。由此，可以理解为早期儿童期虐待经历对个体自我伤害行为的影响过程可能是以个体情绪问题为中介产生的间接作用。本研究提出第二个研究假设：抑郁在儿童期虐待与自伤行为之间起部分中介作用。 一般认为，由于两性别之间情绪易感性、情绪调节策略、情绪表达等情绪心理机制上的差异，外在情景事件诱发的情绪反应和行为应对在男、女之间会是明显区别的。一直以来，心理学研究者们对情绪、行为问题的性别差异研究热情仍未消减。基于大量的抽样研究，抑郁元分析研究结果表明男、女在抑郁水平上的得分无显著差异(吴洪辉，廖友国，2013)，自伤行为的元分析结果表明女生的自伤发生率要高于男生(潘珍等，2016)。各类形式的儿童期虐待在男、女之间的流行情况也存在着较大差异，并对之产生着不同的伤害影响，如：于增艳等人关于儿童期虐待经历与抑郁关系的元分析发现，性别会影响儿童心理虐待与抑郁之间的关系，女性群体中二者关系的效应量会更大(于增艳等，2017；Fisher et al.， 2009)。因此，性别在本研究的三个变量之间也可能起到调节作用。本研究提出的第三个研究假设，即：儿童期虐待、抑郁情绪和自伤行为三者的作用模型存在性别差异。
system: [SYSTEM]

user: 2. 方法 2.1. 被试 本研究抽样方法为方便取样。在湖南省某高校的大学生心理健康教育老师在课堂上邀请学生参与本次调查研究，在班级QQ群中公布“问卷星”网络链接，学生在课后自愿填写，共收到提交问卷312份。网络问卷设置为每题必填后才能进入下一题填写，在数据初步筛查时将总分达极端分的4名被试数据剔除，最终得到308分有效问卷。其中，男生129人，女生179人。被试的年龄范围为17~24岁，平均年龄为18.43 ± 1.26岁。 2.2. 测量 1) 儿童期虐待问卷(Childhood Trauma Questionnaire， CTQ) 该量表由Bernstein等人编制，中文版由赵幸福等人翻译和修订，用于调查16岁前的受虐经历(赵幸福等，2004)。该量表共有28个条目，量表包括情感虐待、躯体虐待、性虐待、情感忽视和躯体忽视5个分量表，和3个效度条目。采用“从不”至“总是”5级计分。本量表在中国大学生群体中广泛应用，且信效度良好(何影等，2010；王佳慧，刘爱书，2014)。该量表的Cronbach’s α系数为0.75。 2) 青少年自我伤害行为量表(Adolescents Self-Harm Scale， ASHS) 本研究采用冯玉修订的“青少年自我伤害行为量表”，对大学生自伤行为进行调查(冯玉，2008)。该量表共有19项关于自我伤害情形的陈述，包含自伤频率和伤害程度两个分量表，分别要求被试根据项目陈述对其自伤行为的发生频率和伤害程度进行评估。本研究仅选频率分量表进行分析，总分为所有项目的得分和，该量表的Cronbach’s α系数为0.82。 3) 流调中心抑郁量表(The Center for Epidemiologic Studies-Depression Scale， CES-D) 此量表由Randloff编制，可用于抑郁症筛查和抑郁症状严重程度评定(Radloff， 1977; 汪向东等，1999)。此量表包括20个条目，对量表四个反向计分条目重新计分后，所有条目得分加和作为量表总得分，得分越高表明被试抑郁程度越高。本研究中该量表的Cronbach’s α系数0.94。 2.3. 统计分析 因本研究所有量表数据均由被试自我报告获得，所以测量中可能存在共同方法偏差。因此本研究首先采用Harman单因素检验法对数据进行分析，以排除共同方法偏差对研究的影响。随后，使用SPSS23.0、AMOS22.0软件对数据进行相关及中介模型建立等统计处理；采用偏差校正的Bootstrap法检验中介效应的显著性。
system: [SYSTEM]

user: 3. 结果 3.1. 共同方法偏差检验 本研究对三个研究变量(儿童期虐待、自伤行为和抑郁)的所有项目进行未旋转探索性因素分析，共析出特征根值大于1的因子共16个，第一个因子解释总方差变异量的19.26%，小于Podsakoff等提出的40%的判断标准(Podsakoff & Mackenzie， 2003)。其次，三个研究变量所有项目的单因素模型的验证性拟合结果较差(χ 2 /df = 3.816; RMSEA = 0.096; CFI = 0.387; TLI = 0.367; SRMR = 0.115)，即三个研究变量的所有项目无共同因子提取。两种检验方法均证明本研究共同方法偏差问题不明显，因此在模型检验中将不考虑共同方法偏差的影响。 3.2. 儿童期虐待、抑郁、自伤行为的相关分析 Spearmen相关分析发现：性别与抑郁、自伤行为相关不显著(p > 0.05)，性别与儿童期虐待的躯体忽视因子呈负相关(p < 0.01)。Pearson相关分发现：抑郁与儿童期虐待各因子均显著正相关(r范围：0.146~0.366)；除躯体忽视之外，自伤行为与其它各儿童期虐待因子也均显著正相关(r范围：0.148~0.371)；抑郁与自伤行为呈显著正相关(p < 0.01)。研究变量间的各相关系数见表1。 3.3. 抑郁在儿童期虐待与自伤行为关系间的中介作用 采用潜变量结构方程模型对假设模型进行验证，结果发现：各观测变量对相应各潜变量负荷值在0.32~0.91之间(p < 0.001)，因子对相应潜变量的负荷均达可接受水平，测量模型良好，潜变量间的路径系数见图1。模型拟合各适配度指标：χ 2 /df = 3.037，p < 0.001，RMSEA = 0.081，SRMR = 0.053，CFI = 0.940，TLI = 0.923，即表明假设模型成立。 设定抽样次数为1000，采用偏差校正的非参数百分位Bootstrap方法进行中介效应检验及置信区间估计。儿童期虐待对自伤行为的总效应值为0.221 (95% CI：0.082~0.524)，其中，儿童期虐待对自伤行为的直接效应值为0.167 (95% CI：0.022~0.508)，经抑郁的中介效应值为0.054 (95% CI：0.010~0.136)，95% CI不包括0，中介效应显著。抑郁在儿童期虐待与自伤行为之间起部分中介作用，中介效应占总效应的0.054/0.221 = 24.43%。 Table 1 Scores of research variables and relevant analysis result M ± SD 2 2-1 2-2 2-3 2-4 2-5 3 4 1性别 a −0.13 * 0.03 −0.12 * −0.01 −0.10 −0.17 ** 0.03 −0.01 2儿童期虐待 35.18 ± 8.54 0.77 ** 0.59 ** 0.38 ** 0.87 ** 0.72 ** 0.33 ** 0.27 ** 2-1情感虐待 7.12 ± 2.47 0.52 ** 0.28 ** 0.46 ** 0.39 ** 0.37 ** 0.19 ** 2-2躯体虐待 5.84 ± 1.72 0.20 ** 0.30 ** 0.21 ** 0.17 ** 0.37 ** 2-3性虐待 5.41 ± 1.38 0.13 * 0.11 * 0.15 * 0.15 ** 2-4情感忽视 9.51 ± 3.94 0.54 ** 0.24 ** 0.20 ** 2-5躯体忽视 7.30 ± 2.53 0.18 ** 0.10 3抑郁 15.52 ± 10.18 0.28 ** 4自伤行为–频率 2.28 ± 4.40 表1. 研究变量的得分情况及相关分析结果(N = 308) 注： a 性别为二分类变量(1 = 男，2 = 女)； * p < 0.05， ** p < 0.01。 图1. 抑郁对儿童期虐待和自伤行为的中介作用模型 3.4. 模型的性别差异 为检验上述中介模型在不同性别之间是否存在差异，本研究以性别为类别变量进行多群组分析，男、女潜变量结构方程模型的路径系数见图2。进一步检验路径系数的性别差异，结果发现：儿童期虐待与抑郁情绪、抑郁情绪与自伤行为的路径系数无显著性别差异(t分别为0.29和0.83，均p > 0.05)，儿童期虐待与自伤行为的路径系数在男女之间存在显著差异(t为3.47，p < 0.05)。男、女模型儿童期虐待对自伤行为的总效应、直接效应及抑郁情绪的中介效应结果见表2，根据效应值95%CI的范围可知：男生模型中，儿童期虐待对自伤行为的直接效应和总效应均不显著，抑郁情绪的中介效应也不显著；女生模型中，儿童期虐待对自伤行为的直接效应和总效应均显著，但抑郁情绪在儿童期虐待与自伤行为的中介效应不显著。 图2. 男/女生儿童期虐待、抑郁情绪与自伤行为的潜变量模型图 Table 2 Summary of effects of childhood abuse on self-injury behavior in male and female model 直接效应(95%CI) 中介效应(95%CI) 总效应(95%CI) 男生 0.108 (−0.171~0.388) 0.114 (−0.012~0.239) 0.222 (−0.034~0.479) 女生 0.181 (0.056~0.305) 0.024 (−0.017~0.064) 0.204 (0.088~0.321) 表2. 男、女生模型中儿童期虐待对自伤行为的效应汇总
system: [SYSTEM]

user: 4. 讨论 4.1. 抑郁情绪在儿童期虐待与自伤行为之间的中介作用 共同方法偏差检验结果表明本研究的共同方法偏差问题不明显，数据具有一定的可靠性。相关分析发现儿童期虐待、抑郁情绪和自伤行为三者之间相关密切，研究结果与之前文献报告的结果相同(苏静等，2015；Zlotnick et al.， 1995; 吴晓瑞等，2016)，即大学生在儿童期遭受的虐待频率越高，现体验到的抑郁水平也越高，同时也会报告更高水平的自伤行为。儿童期受虐待经历与青少年、成年阶段的自伤行为密切关联，以往研究多停留于变量相关研究，较少对影响二者关系的心理机制进行实证探索。在众多自伤行为发生的理论解释模型中，情绪体验回避模型认为自伤行为可以让个体回避不想要的内在体验，这种体验回避均与负性情绪相关联，包括个体悲伤的想法、记忆、躯体上的感觉或其它厌恶的内在体验(冯玉，2008；陈静等，2014)。童年期遭受虐待，触发个体厌恶情绪，不良环境导致他们不知道如何恰当地表达和调节情绪。为逃脱不愉快的情绪体验，该理论认为个体会采用自伤行为来应对情绪困扰，发展过程中不断强化，久之变成一种自动化的应对模式。这也就是说，情绪因素是儿童期虐待影响自伤行为的中介变量。冯玉在其研究中也证实情绪调节在言语虐待与自伤行为之间起部分中介作用(冯玉，2008)。 基于整个研究样本，本研究通过结构方程建模的方法对提出假设中介模型进行了检验，验证本研究所提假设中介模型成立。研究结果表明抑郁在儿童期虐待与自伤行为之间起部分中介效应，即：儿童期虐待可直接影响自伤行为；此外，儿童期虐待还可以通过抑郁的中介作用间接影响自伤行为。情绪调节困难与情绪症状密切关联，因此本研究结果与冯玉的研究结果相符(冯玉，2008)。比较中介效应和直接效应，中介效应对自伤行为发生解释的贡献量小于直接效应的贡献量，占总效应的24.43%。抑郁中介效应显著，提示干预抑郁情绪可以缓解儿童期虐待对青少年自伤行为的影响。研究结论对解释自伤行为的影响因素和形成机制具有一定的理论意义，其结论也对预防和干预自伤行为具有一定的实践意义。 4.2. 研究变量得分及相互作用机制的性别差异 基于以往元分析结果，本研究推测儿童期虐待、抑郁情绪和自伤行为三个研究变量的得分存在性别差异，但本研究相关分析仅发现儿童期虐待与性别相关显著，抑郁情绪和自伤行为与性别相关不显著。据儿童期虐待各因子与性别的相关结果可知，男生报告自己在儿童期比女生会遭受更高频率的虐待，在具体虐待形式中，仅在躯体虐待与躯体忽视两个因子上表现出性别差异，但其他虐待形式，即：情感虐待、情感忽视和性虐待的性别差异却并不显著。这可能是基于男、女不同的生理特征差异，社会对男、女的性别期待不同而导致。女生应表现柔弱，男生应表现坚强，加之“男孩穷养，女孩富养”等类似教养观念，中国社会普遍倾向认同对男生进行棍棒式教育，而对女生则应多采用说理式教育，这可能是男生在躯体上受虐待高于女生的一个原因。基于大量抽样研究，元分析结果表明儿童期虐待、抑郁情绪和自伤行为的性别差异并没有多大一致性。三个研究变量的相关性会受到抽样的局限，所以，本研究中研究变量与性别的相关结果也具有合理性。 本研究中，男、女在儿童虐待上得分显著差异，在自伤行为上的得分差异不显著，即性别对儿童期虐待与自伤行为的关系起显著调节作用，这提示相关心理因素在儿童期虐待与自伤行为之间的作用机制会存在性别差异。本研究通过中介模型分析和中介模型的性别群组差异检验证实抑郁情绪在儿童期虐待与自伤行为之间的作用机制确存在性别差异。在男生群组中，儿童期受虐待经历显著预测抑郁情绪，抑郁情绪显著预测自伤行为，中介作用分析发现抑郁情绪的中介作用并不显著。此结果提示，男生儿童期受虐经历对其抑郁情绪、抑郁情绪对其自伤行为的影响仅为简单单向作用，男生青少年、成年阶段的自伤行为与其儿童期受虐经历并无显著直接联系。女生群组中儿童期受虐经历显著预测抑郁情绪、自伤行为，但抑郁情绪不能显著预自伤行为，抑郁情绪的中介作用也不显著。本研究中，女生抑郁情绪对其自伤行为的预测并不显著，这与以往二者关系的研究结果相异(Gratz & Roemer， 2008; 吴晓瑞等，2106)。以往，抑郁症或抑郁情绪与自伤行为的关系研究多表明抑郁是自伤行为的一个独特危险因素(谢唯等，2017)，而本研究中抑郁情绪和自伤行为则似同为儿童期虐待的不良应对策略，两者之间的因素关系却不被证实。 4.3. 本研究对自伤行为干预的意义 基于整个研究样本，本研究证实：儿童期虐待可直接作用于自伤行为；另一方面，儿童期虐待可以通过抑郁的中介作用间接作用于自伤行为。显著直接效应提示保护未成年儿童，创造健康成长环境，对青少年、成年阶段的心理健康意义重大。模型中介效应显著，自伤行为的体验回避模型被支持，这提示在未来的心理健康教育中，提升大学生的情绪管理策略和能力，改善情绪，对大学生的自伤行为预防和干预具有积极作用。然而，区分性别时，抑郁情绪在儿童期虐待与自伤行为之间起中介作用的假设在男、女生群体中却均不被支持，三个研究变量的相互作用模型存在较大的性别差异，这与总群体的研究结果大为不同。一方面，今后研究需加大样本量进一步对性别差异结果进行检验。另外，情绪认知、调控等方面可能与自伤行为更为密切。所以，今后研究有必要加其它情绪心理因素来进一步验证自伤行为的回避理论模型。若自伤行为发生的情绪机制存在性别差异，则提示自伤行为的干预在男、女之间需要根据情绪机制的差异区别对待。
system: [SYSTEM]

user: 5. 结论 本研究结论：1) 抑郁在儿童期虐待与自伤行为间起中介作用的假设模型在本研究总样本中成立，但区分男、女群体时，假设中介模型均不成立，研究变量间的相互关系也存在性别差异；2) 男生群体中，儿童期虐待显著预测抑郁，抑郁显著预测自伤行为，儿童期虐待对自伤行为的预测不显著；3) 女生群体中，儿童期虐待显著预测抑郁和自伤行为，抑郁不能显著预测自伤行为。
system: [SYSTEM]

user: 摘 要 数据缺失是一个常见但难以处理的问题。文章简要介绍了数据缺失的几种机制，以及处理缺失数据的一般性方法，并对各种缺失数据的处理方法的特点及适用情况进行了比较。 关键词 :缺失值，缺失机制，填补方法 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 在目前众多的研究领域中，经常会出现无回答的调查，这就不可避免的会带来数据缺失。数据缺失不仅会给分析数据带来一定的困扰，也会给分析结果带来偏差。 目前国内的大多数心理学期刊并没有对缺失数据如何处理给出明确的要求，但杂志社对文章质量的要求越来越高，当下一些国际性的心理学期刊对缺失数据的处理要求必须给出详细的说明，另外随着国内心理学研究与国际逐渐接轨融合，对于缺失数据的处理要求也会越来越严格，因此，如何减少或者消除缺失数据带来的偏差变得越来越重要，对缺失数据的研究也越来越受到研究者的重视。
system: [SYSTEM]

user: 2. 数据缺失机制 在处理缺失数据之前，我们需要知道数据缺失的机制，不同的缺失机制对应不同的分析方法，了解这些信息后才能选择恰当的方法来处理缺失数据。 在缺失数据的处理中，数据缺失机制的作用被研究者们长期的忽略，直到1976年Robin博士才给出明确的说明，他把缺失机制划分为三类(Rubin， 1976)：随机缺失(Missing at Random， MAR)、完全随机缺失(Missing Completely at Radom， MCAR)、完全非随机缺失(Not Missing at Random， NMAR)。 我们用p表示概率分布，R代表缺失数据指标，Yobs和Ymiss分别表示数据集中实际观测到的部分以及缺失的部分，Φ是反映R与数据之间关系的一个参数。 1) 完全随机缺失。MCAR缺失要求变量Y上数据缺失的概率与其他观测变量无关，也与Y本身的值无关，用公式表示为p(R|φ) (Rubin， 1976)。这种缺失机制下数据缺失是完全随机的行为，其丢失的概率是未知的，换言之数据发生缺失与否与变量的取值没有任何关系。 2) 随机缺失。MAR缺失是指数据的缺失取决于数据集中其他变量，与自身取值没有关系，用公式表示为p(R|Yobs， φ) (Rubin， 1976)。 3) 完全非随机缺失。NMAR缺失是指数据是否缺失与数据集中其它变量的取值没有关系，只与缺失变量自身的取值有关，用公式表示为p(R|Yobs， Ymis， φ) (Rubin， 1976)。 不同的缺失值机制意味着需要采用不同的处理方法。NMAR与MAR和MCAR的情况不同，因此下面讨论的方法主要集中在MAR和MCAR两种条件下，NMAR的处理方法可参见Enders的文章(Enders， 2010)。 Enders曾以下面以一组数据对缺失数据的三种机制做了的介绍(Enders， 2010)。 例如在一次测试中，要求IQ达到88分以上才能参加随后的人格测验，这样IQ分数为78、84、84、85和87的数据便缺失了。这种数据缺失与人格变量自身无关，但却与IQ有关，称为随机缺失MAR。 完全随机缺失MCAR情况下数据缺失是随机的，不符合任何规律。换句话说，变量缺失值的出现完全是个随机事件，例如下表中IQ为78、84、96、112和134上数据的缺失，称为完全随机缺失MCAR。 完全非随机缺失NMAR数据缺失与其他变量无关，与自身表现得分相关。例如在表1中，公司新录用了14名员工，其中5名员工由于表现较差，在试用期内被辞退，年终表现评定中，被辞退的5名员工的表现分缺失了，这种情况下的数据缺失即为完全非随机缺失NMAR。 Table 1 Job performance evaluation for missing values of MCAR， MAR and NMA IQ Complete MCAR MAR NMAR 78 9 - - 9 84 13 13 - 13 84 10 - - 10 85 8 8 - - 87 7 7 - - 91 7 7 7 - 92 9 9 9 9 94 9 9 9 9 94 11 11 11 11 96 7 - 7 - 99 7 7 7 - 105 10 10 10 10 112 10 - 10 10 134 12 - 12 12 表1. MCAR、MAR和NMAR缺失值时的工作绩效评估
system: [SYSTEM]

user: 4. 缺失数据处理方法比较评价 处理缺失值数据的方法有很多，但各有各的特点，每种处理方法都有其需要满足的要求。有研究者认为(庞新生，2010)大部分缺失值处理方法应满足如下假定：一是多元正态分布，二是所有变量独立同分布，三是数据随机缺失。在缺失数据处理之前，需要对上述条件做出判断。 一般而言，在使用缺失数据技术时，最需要考虑的是缺失数据处理方法所适用的数据缺失类型。如果缺失数据不是随机缺失的，数据分析可能将导致偏差；如果数据类型不是MCAR机制，将不适用与多重插补MI的方法；以模型为基础的方法如完全信息极大似然估计又需要数据类型满足MAR机制。因此，在对缺失数据处理前我们需要对样本分布是MCAR、MAR或是NMAR机制有一定的了解。 选择何种缺失数据的处理方法主要取决于数据的性质和质量、用途以及缺失数据的机制。据此，对文章上述各种缺失数据处理方法进行比较分析，主要内容见下表2。 Table 2 Missing data processing method compariso 处理方法 前提 难易程度 适用范围 稳健性 偏差 列表删除 MCAR 易 缺失5%内 很差 大 成对删除 MCAR 易 几乎不用 很差 大 均值插补 MCAR 易 少用 很差 低估方差及抽样误差 回归插补 MAR 易 广泛 差 低估方差及抽样误差 hot-deck插补 MAR 较难 少用 差 方差估计较好但抽样误差不易控制 EM算法 MAR 难 广泛 很好 弱偏或无偏 多重插补 MCAR 难 广泛 好 弱偏 表2. 缺失数据处理方法比较 根据上表，EM算法和多重插补适用范围广泛，处理缺失数据时，具有比其他方法更好的稳健性，其偏差也很小，不失为很好的选择，但是对于这些以模型为基础的方法操作起来存在困难，需要有一定的数学和计算机基础，适合具有处理缺失数据的专长和工具的研究者。 均值插补和回归插补虽然操作相比多重插补简单，但是其稳健性较差，常常会低估方差及抽样误差，在实际应用中并不是特别广泛。 最后，如果缺失数据非常少，数据缺失在5%之内，我们也可以采用列表删除的方法，而此时也不会引入大的误差，操作起来也比较容易简单。
system: [SYSTEM]

user: 5. 小结 综上所述，没有哪一种方法是既简单适用性又广误差又小的，每种方法都有各自的利弊，有着各自的适用的条件。因此，我们在选择方法时，应该注意了解调查数据的背景、特征等，尽可能地挖掘样本总体的各种辅助信息，此外对待缺失值处理方法也应该持一种科学谨慎的态度，毕竟方法在不断发展，现在的方法都不会一直是最优最好的方法，选择最合适的方法处理自己的数据才是最好的方法。
system: [SYSTEM]

user: 摘 要 华西秋雨是我国华西地区特有的降水现象，东亚西风急流也是影响我国降水的重要大气环流之一。研究东亚西风急流和我国华西秋雨的关系，对预测我国未来华西秋雨有重要意义。本文使用相关分析、EOF-CCA分析等方法研究东亚西风急流与华西秋雨的关系，其中使用了华西秋雨指数和东亚西风急流的相关指数，得出了以下结论：华西秋雨变化的整体趋势为华西秋雨强度随着时间增长而减弱，存在华西秋雨异常偏强年和异常偏弱年；秋季东亚西风急流整体趋势为急流位置维持于39˚N，144.3˚E，急流强度维持于36.3 ms −1 左右；东亚西风急流和华西秋雨的7年滑动相关呈周期性变化，滑动周期为15年；秋季东亚西风急流强度区的西部与华西地区的东北部呈正相关，与华西地区的西部呈负相关。 关键词 :东亚西风急流，华西秋雨，相关分析，EOF-CCA分析 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 华西秋雨是秋季的一种特殊天气现象，其降水强度不大，但雨期较长。华西秋雨的累计降水量较大，容易发生洪涝灾害，损害人民的财产，对人民的生产生活造成影响。研究华西秋雨有重要意义，我国有许多学者在这方面做了很多工作，高由禧 [ 1 ] (1958)等早在1958年就讨论了华西秋雨的地区分布和开始结束日期，并指出华西秋雨的开始结束日期与西风急流有关。白虎志 [ 2 ] (2004)等使用了华西秋雨指数，分析了华西秋雨的气候区，研究了华西秋雨的成因。柳艳菊 [ 3 ] (2012)等分析了华西秋雨的多雨区和影响华西秋雨的大气环流等因素。王红军 [ 4 ] (2018)等分析了2017年影响华西秋雨的大气环流，分析了阻塞高亚、副热带高压、高空急流以及华西秋雨水汽来源等因素。本文在前人的基础上，使用了蒋竹将 [ 5 ] (2014)等定义的华西秋雨指数MARI，华西秋雨指数是包含雨量与雨日的综合指标，可以很好地反映华西秋雨的特征。 东亚西风急流是重要的大气环流之一，对我国气候有重大影响，是我国降水的重要影响因素之一。叶笃正 [ 6 ] (1958)等分析了6月和10月的大气环流突变，指出西风带向北向南移动时，大气环流形成典型的夏季环流型和冬季环流型。陶诗言 [ 7 ] (1958)等分析了梅雨与高空急流的关系，指出当大气环流发生季节性变化时，梅雨发生或结束。兰明才 [ 8 ] (2011)等利用降水资料和再分析资料，分析了急流对东北夏季降水的影响。方晓洁 [ 9 ] (2009)等分析了急流对夏季降水的影响，当急流异常移动时，会影响降水。李崇银 [ 10 ] (2004)等探讨了东亚西风急流和东亚夏季风活动的关系，指出东亚高空西风急流的位置一般会发生两次北跳，第一次北跳与南海夏季风爆发有关，第二次北跳与江淮梅雨的起始有关。况雪源 [ 11 ] (2006)等探讨了东亚西风急流位置对长江中下游夏季降水的影响，指出急流异常偏北时，长江流域不易降水；急流异常偏南时，长江流域易引发洪涝。杜银 [ 12 ] (2008)等探讨了东亚西风急流对梅雨的空间分布的影响，指出梅雨期当东亚高空西风急流主体呈显著的西太平洋分布型态时，有利于长江下游地区形成集中的强降水区域。本文通过相关分析、EOF-CCA分析等方法分析东亚西风急流和我国华西秋雨的关系，意在分析东亚西风急流对华西秋雨的影响成因，从而对我国华西秋雨的研究做出一定贡献。
system: [SYSTEM]

user: 2. 资料与方法 2.1. 资料选取 本文使用四川省气象局提供的站点降水资料，根据《华西秋雨监测业务规定(试行)》规定的华西秋雨监测站点，选取1971年~2016年全年无缺测的330个气象站点(图1为站点分布图)；及NCEP/NCAR提供的1971年~2016年200 hPa纬向风逐月再分析资料，分辨率为2.5 × 2.5。 图1. 华西秋雨监测站点分布 2.2. 方法简介 2.2.1. 相关分析 相关分析 [ 13 ] 是气象学上一种常用的分析方法。对于两个大气变量，通过计算二者的相关性，可以知道二者的紧密程度。对于任意的两个大气变量，其相关系数为： r k l = 1 n ∑ i = 1 n x z k i x z l i = 1 n ∑ i = 1 n ( x k i − x ¯ k s k ) ( x l i − x ¯ l s l ) (1) 通过计算相关系数，就可以知道两个大气变量之间的相关性。相关系数的值在−1~1之间，其绝对值越大，相关性越强。相关系数为负，说明二者呈负相关，即异增异减，当一者增加(减少)时，另一者减少(增加)；相关系数为正，说明二者呈正相关，即同增同减，当一者增加(减少)时，另一者增加(减少)。 2.2.2. 滑动相关分析 滑动相关分析方法是在相关分析的基础上再对相关系数进行一步处理的方法。选择一段时间作为滑动窗口，如7年、11年、21年，从开始时间开始计算滑动窗口内的滑动相关系数，滑动窗口滑动1年，计算一个滑动相关系数，到滑动窗口滑动到结束时间停止，这就得到所要分析的滑动相关系数时间序列。 相关分析方法计算所得的只是大气变量整体之间的相关性，所算得的只有一个值，难以进行进一步的分析。而通过滑动相关分析所得是一段滑动相关系数，反映了两个大气变量不同时期的相关性。通过分析滑动相关系数时间序列可以得到两个大气变量之间的关系。 2.2.3. EOF-CCA分析方法 EOF-CCA分析是将主分量分析(empirical orthogonal function，缩写EOF)和典型相关分析(canonical correlation analysis，缩写CCA)结合的一种分析方法。黄嘉佑 [ 14 ] 等将这种方法在2000年用于对夏季汛期的预报，本文也是借鉴这种方法用于对华西秋雨和东亚西风急流的分析。 主分量分析 [ 13 ] (EOF)可以提取出要素场的主分量，头几个主分量占原变量的方差绝大部分，所以主分量可以反映出气象要素的主要特征。典型相关分析 [ 13 ] (CCA)可以两个变量组合为两个新的典型变量。研究气象要素场就可以代替为研究典型变量，这样就可以分析出两个气象要素场之间的关系。 EOF-CCA分析方法的具体公式如下，设预报要素场为 ( q × n ) 的横资料阵，即为 Y = ( q × n ) ；预报因子为 ( p × n ) 的横资料阵，即为 X = ( p × n ) ，设 q < n ；分别对Y、X做EOF (主分量分析)有： Y = A y F y X = A x F x 其中A y 、F y 分别为Y的因子荷载和标准化主分量，A x 、F x 分别为 的因子荷载和标准化主分量。 对主分量进行CCA (典型相关分析)，设有新的线性组合： V = P F y U = Q F x 其中P、Q为线性组合系数，V、U为典型因子，其中V、U满足： ( 1 / n ) V V ' = I ( 1 / n ) U U ' = I 两个典型因子的相关系数为典型相关系数，其对角阵R表示为： R = ( 1 / n ) V U ′ = P [ ( 1 / n ) F y F x ' ] Q ′ 两个主分量的相关阵W表示为： W ( 1 / n ) F y F x ' = P ′ R Q 可得： W W ′ = P ′ R R ′ P = P ′ Λ P W ′ W = Q ′ R ′ R Q = Q ′ Λ Q 则P、Q分别为 W W ′ 和 W ′ W 的特征向量构成； D、C为典型因子的因子荷载，可用下式求出： D = A y P ′ C = A x Q ′ 至此，完成了Y、X的EOF-CCA模型构建；V、D分别为Y的模型因子和因子荷载；U、C分别为X的模型因子和因子荷载；以下三个因素将对EOF-CCA的效果产生影响： 1) 主分量的模态选取数越多，越能模拟原要素场；但模态选取数越多，运算量越大。 2) 在EOF-CCA模型中典型因子的模态数与主分量模态数相同；典型因子模态数越多，越能反应两个主分量的之间的典型相关性，但计算量加大，且后面典型相关性减小，不具有代表性。 3) 两个原始要素场选取的格点数越多，结果越准确；选取的时间越长，结果越准确；选取的时间段越接近，结果越准确。 2.2.4. 东亚西风急流相关指数和华西秋雨指数计算 华西秋雨和东亚西风急流的研究有许多，但二者始终不具有统一而公认的指标，本文选取了华西秋雨指数(MARI)作为研究华西秋雨的工具，选取了三个东亚西风急流相关的指数作为研究东亚西风急流的工具，东亚西风急流中心指数和东亚西风急流轴线指数为描述东亚西风急流位置的指标，而东亚西风急流强度指数为描述东亚西风急流强度的指标。 本文根据廖琦 [ 15 ] 等对东亚西风急流轴、急流强度的定义和魏林波 [ 16 ] 等对东亚西风急流中心的定义，对东亚西风急流相关的指数给出以下定义： 1) 东亚西风急流中心指数：0˚N~60˚N，40˚E~180˚E内200 hPa纬向风风速最大值的平均经度； 2) 东亚西风急流轴线指数：20˚N~60˚N，70˚E~150˚E内200 hPa纬向风风速最大值的平均纬度； 3) 东亚西风急流强度指数：32.5˚N~47.5˚N，110˚E~165˚E (图2)内的急流的平均强度。 图2. 东亚西风急流强度区(虚线方框)与我国的相对位置 本文采用蒋竹将 [ 5 ] 等改进的华西秋雨指数 M A R I = R a u t R y e a r × L (2) 华西秋雨指数的计算范围为每年的9~11月，降水量表示为R aut ，年降水量表示为R year ，其中L为降水日数。通过计算L来计算华西秋雨的雨期开始和结束，L连续出现5个降水雨日则华西秋雨开始，L连续出现5个非降水雨日则华西秋雨结束。R aut 则为L中的降水量，R year 则为每年的总降水量。