地球科学进展  2018 , 33 (7): 762-774 https://doi.org/10.11867/j.issn.1001-8166.2018.07.0762

研究简报

相对湿度及其变化的年循环研究进展

念达, 邓琪敏, 付遵涛*

北京大学物理学院大气与海洋科学系,北京 100871

Research Progress of Relative Humidity and Its Changing Annual Cycle

Nian Da, Deng Qimin, Fu Zuntao*

Laboratory for Climate and Ocean-Atmosphere Studies, Department of Atmospheric and Oceanic Sciences, School of Physics, Peking University, Beijing 100871, China

中图分类号:  P467

文献标识码:  A

文章编号:  1001-8166(2018)07-0762-13

通讯作者:  *通信作者:付遵涛(1970-),男,黑龙江尚志人,教授,主要从事非线性大气动力学与气候变化研究.E-mail:fuzt@pku.edu.cn

收稿日期: 2018-01-11

修回日期:  2018-05-21

网络出版日期:  2018-07-20

版权声明:  2018 地球科学进展 编辑部 

基金资助:  *国家自然科学基金项目“变化的相对湿度年循环及其影响”(编号: 41675049)资助.

作者简介:

First author:Nian Da(1993-), female, Qujing City, Yunnan Province, Ph.D student. Research areas include climate changes.E-mail:danian@pku.edu.cn

作者简介:念达(1993-),女,云南曲靖人,博士研究生,主要从事气候变化研究.E-mail:danian@pku.edu.cn

展开

摘要

对相对湿度(RH)年循环进行有效诊断可以对季节内雾霾的监测和调控,流行性疾病爆发空间分布和防治以及季风降水量的实时预报和决策等提供重要帮助。通过总结RH及其年循环的研究现状,系统介绍了当前几类常用的信号提取方法,并评估了这些方法提取年循环信号的可行性。得出以下结论:针对RH年循环的时变特征和非对称三角波、类方波、状态瞬变等复杂结构,以谐波为基础的提取方法难以成功,高抗噪的非线性信号提取方法是解决这个问题的突破口。因此,现阶段对于RH年循环的研究迫切需要完善并提高RH观测数据的质量,精确提取并量化其年循环的时变特征和复杂结构,最后结合动力学过程和统计学进行分析,以期使我国对RH年循环(频率、相位、振幅等)变化的物理机制研究达到更高水平。

关键词: 相对湿度 ; 年循环 ; 时间序列分解

Abstract

Studies on the characteristics of relative humidity annual cycle change include the frequency, phase, and amplitude of the time series and their changes. The effective diagnosis of the relative humidity annual cycle can provide important help in the monitoring and regulation of seasonal haze, the spatial distribution of epidemic outbreaks and their prevention and control, and the real-time forecasting and decision-making of monsoon precipitation. Different from studies on the trend of the relative humidity, the diagnosis in the annual cycle is scarce. This paper summarized the research status of relative humidity and its annual cycle, introduced several current methods for extracting common signals, and evaluated the feasibility of these methods for extracting annual cycle signals. Due to the time-varying characteristics of the annual cycle of relative humidity and complex structures such as asymmetric triangular waves, square-like waves, and state transients, harmonic-based extraction methods are difficult to succeed. The nonlinear signal extraction method with high noise immunity will solve this problem. At present, the relative humidity annual cycle studies urgently need to be improved and improve the quality of relative humidity observation data, accurately extract and quantify the time-varying characteristics and complex structure of the annual cycle. Moreover, combining the dynamic process and statistical analysis, we also need to study the physical mechanism of the change of relative humidity annual cycle (frequency, phase, amplitude, etc.) in China.

Keywords: Relative humidity ; Annual cycle ; Time series decomposition.

0

PDF (2130KB) 元数据 多维度评价 相关文章 收藏文章

本文引用格式 导出 EndNote Ris Bibtex

念达, 邓琪敏, 付遵涛. 相对湿度及其变化的年循环研究进展[J]. 地球科学进展, 2018, 33(7): 762-774 https://doi.org/10.11867/j.issn.1001-8166.2018.07.0762

Nian Da, Deng Qimin, Fu Zuntao. Research Progress of Relative Humidity and Its Changing Annual Cycle[J]. Advances in Earth Science, 2018, 33(7): 762-774 https://doi.org/10.11867/j.issn.1001-8166.2018.07.0762

1 引 言

水汽作为大气中最强的温室气体,它的温室贡献是二氧化碳的数倍[1]。大气中水汽通过云的形成和演化吸收辐射而影响气候系统中其他变量的变化。大气中的水汽也是水循环的主要组成部分,因为对流层低层的水汽是所有天气系统降水的主要水源。用于量化大气中水汽常见的变量有比湿(q)和相对湿度(Relative Humidity,RH)。与比湿相比,RH被认为是气候反馈分析中更利于分析的状态变量[2]。因为在变暖的气候条件下,来自水汽和云的重要气候反馈很大程度上依赖于RH的变化[3]。RH作为一个衡量大气水汽含量的物理量,对温室效应、水循环、气候反馈、云微物理、疾病传播和人类健康等众多方面都有重要影响。就定义来说,RH是大气中现存水汽量与在一定温度下空气能够保持的水汽量的比值[4]。作为当前气候条件下大气中最强的温室气体,Schmidt等[5]用模式计算得出,水汽对温室效应的贡献将占到约50%,当大气中CO2的含量加倍时,这个贡献不变。大气水汽的变化与降水事件紧密相关,大量全球气候模式预测大气水汽含量会因为响应全球平均温度增高而增加,而增加的水汽又会因温室效应导致地面温度持续增加,从而形成加强的正反馈,这也是中高纬度模式模拟降水增加的主要原因[6,7,8]。RH在观测资料、再分析资料以及气候模式结果中存在较多差异,直接影响了大气可降水量的模拟准确度及趋势研究的可靠性[9]。很多关于RH的研究也有助于理解雾和霾的形成[10,11,12],雾一般出现在RH较高的条件下,当气象观测站无法提供雾是否形成的直接观测数据时(比如夜晚),可以通过4个可测量参数是否达到条件来进行诊断,而RH高于90%就是其中的条件之一[11]。沿海地区雾的形成可能由夜晚海陆风带来的平流导致,也可能由于近地面大气的辐射冷却导致[13],因而对于RH的研究促进了对区域天气、气候状况的理解。RH也和疫病的爆发和传播有关,比如19世纪北爱尔兰2次马铃薯疫病爆发时正好对应地出现了当地RH高值[14]。也有研究指出,气溶胶传播流感病毒依赖于RH和温度,RH高于80%时病毒传播会被完全阻塞[15]

2 RH及其年循环研究现状

由于RH在横跨多个领域都具有重要性,在全球变暖、大气温室效应增强的背景下,RH对于气候变化的响应以及对应的天气学、气候学过程研究也越来越受到人们的关注。虽然RH如此重要,但是,与对气温和降水的广泛研究相比,从RH的角度探究气候变化的研究还非常少。目前开展的工作主要集中在3个方面:RH的气候学与变化趋势,RH的变化机理,RH的年循环及其变化。下面分别总结这3个方面的研究现状和发展动态。

2.1 RH的气候学与变化趋势

关于RH气候学和变化趋势的研究工作首先由Peixoto等[16]在1996完成的。他们首次给出了RH的气候学特征,同时也发现RH的全球分布在对流层的各个高度均不是纬向均匀的,而是在所有的纬度都具有不同强度的中心。从全球角度来看,表面RH有相对小的空间和年际变化,海洋平均为75%~80%,陆地为70%~80%(高地与沙漠地区为30%~60%),由于日气温变化大,夜间的RH要比白天的高2%~15%。近年来,更多的研究认为RH只有很小量级的趋势,也有研究认为RH没有显著变化,它保持近乎常数[16,17,18,19,20,21,22,23,24,25]。第四次IPCC报告指出:在较大空间尺度和较大时间尺度上RH是保持不变的[26]。大量的全球模式也得到了RH在全球尺度和年平均尺度上为常数的结果[27],在研究温室效应的水汽反馈时也通常假定RH不变进行模拟[28,29] 。对另一个湿度参量——比湿的观测资料研究中发现:比湿随着温度的变化几乎是线性增长的。这个结果与RH没有长期变化趋势时用Clausius-Clapeyron方程(C-C方程)推导出的比湿和温度的关系十分类似,且得到的线性增长率与C-C方程导出的增长率量级差不多[30]。Zhao等[31]利用无线电探空数据集研究发现中国地区对流层比湿有增长趋势,且与中低层温度变化有较强正相关,增长率稍高于C-C方程得到的增长率,RH仅有较小变化以及微弱趋势。但是,Willett等[24]提出在一些局部地区,虽然比湿随着温度正相关变化,增长率却大大超过了C-C方程得到的增长率,说明RH可能不是常数。也有研究指出:在特定的地区,局地温度和比湿的关系在RH平均值较高的地区接近C-C方程导出的关系,此时RH变化不大,但对于云出现频率较低的地区,RH可能变化较大[32]。也有研究根据全球近地表比湿数据集HadISDH得出:大尺度的全球平均比湿、南北半球平均比湿的变化趋势非常小或几乎为0[33],那么RH在这些区域很有可能不保持为常数。在1976—2004年,全球整体的表面RH变化小,但是在海洋上下降趋势显著,为-0.11%/10a~-0.22%/10a,陆地上,RH增加大(0.5%/10a~2.0%/10a)主要发生在美国中部和东部、印度和中国西部[30]

虽然表面上RH的全球平均变化非常小,但是在区域尺度上它的变化可以非常大,不同区域间的差别也会非常大[34]。近几十年来,多个研究发现RH的变化并不是保持常数那么简单,特别是在区域尺度上RH变化的情况更为复杂。Simmons等[25]研究发现近几十年陆地区域上RH普遍减少,净干旱区域包括美国、南美和欧洲。另外关于北美、欧洲多个区域包括波兰、北爱尔兰、西班牙、瑞士等国家以及中国台湾等多个地区的研究也发现了RH年平均值的减小趋势[14,16,21,34~37]。第五次IPCC报告结合最近10~15年观测和再分析数据得出:中低纬度地区RH有下降趋势[38]。但与此相反的是,在保加利亚高山区的3个峰顶站中有2个观测站的RH序列出现增加趋势,1个观测站出现减小趋势[39]。在伊朗北部和南部沿海区域,年平均RH分别以1.03%/10a和0.28%/10a的速率增长[40]。在海洋上,Pfahl等[41]通过观测记录得出RH和温度在热带外呈现正相关,Pearson相关系数在北大西洋最大,为0.6,在热带内海洋呈现负相关。从RH变化趋势的研究看来,目前关于大气湿度对全球变暖的响应过程及其机理仍然不甚清楚。在美国,夜间的湿度增长趋势要大于白天的趋势,比湿的增加与温度的上升趋势一致。RH的趋势要弱于比湿的趋势,且统计上是不显著的。但是有证据表明RH确实呈现出上升趋势,特别是在冬季和春季。不论从区域还是全球的角度来说,化石燃料的消耗带来的人为水汽都太弱,不足以解释观测中湿度的趋势[17]。在伊朗,利用参数与非参数趋势检验均得到年平均RH的变化趋势为:在北部和南部沿海分别以1.03%/10a与0.28%/10a增加[40]。在中国,大部分地区表面RH呈现出下降趋势,主要原因是比湿的增加小于饱和比湿的增加。但是从2003年开始,在中国的大部分地区饱和比湿显著增加,但是比湿却突然下降,从而导致RH有更大程度的下降[42]。中国地区对流层中低层年平均或季节平均的RH在大部分区域没有显现出明显的变化趋势[4]。高RH发生频次的趋势随高度和季节变化,在近地面层,东南部的春季和北方的夏季大多数观测站的RH都呈现下降趋势;在对流层低层,大多数观测站的RH在夏秋冬三季均呈显出上升趋势[43]。在西班牙,利用1961—2011年的月平均表面RH和比湿序列研究发现,RH呈现出大的下降(最大的下降发生在春季和夏季),而与此相反的是,在此期间比湿整体上没有变化[36]。从这里可以看出,区域尺度上RH(或比湿)的变化行为差别非常大,强烈依赖于当地的环流条件和水汽输送情况。

另外,仅根据RH年平均值的变化趋势不足以描述RH随着气候变化的响应特征。尤其是目前区域性天气、气候预报以及气候内变率如ENSO等方面的研究越来越需要RH在天气尺度和季节性尺度的变化信息[44]。目前对RH季节性的研究较为零散,多为分季节性讨论RH的趋势变化,地区不同结果不同。例如,根据西班牙1920—2011年RH序列的研究表明,春季、夏季RH减少趋势最为显著[36]。在中国地区的研究中,一般将中国划分为6个区域[4,42],水汽在夏季增多、冬季减少。这些研究都说明了RH变化的复杂性,而单独用RH趋势的变化无法全面说明和理解RH的变化,这也说明目前研究尚不清楚RH对气候变化的响应机理是什么。

2.2 RH的变化机理

RH的变化主要受气温和水汽输送2个方面因素的影响。研究发现在长时间尺度上,气温与湿度的变化主要由C-C方程决定,即饱和水汽压随着气温呈指数增长而维持RH近似为常数,而且这种关系在气温变化大的时候更为显著[41]。这也是大气环流模式中得到的水汽反馈与恒定RH条件下的结果接近的原因,RH的恒定不变是基于全球平均的长时间尺度的结果[3]。但是,事实上,RH的变化具有非常明显的空间结构和不同时间尺度的变率,其对应的变化机理差别也会非常大。在短的天气时间尺度,RH的变化可以影响雾的形成或者表面蒸发结构的变化。在日平均时间尺度上,气温与RH的变化关系随着空间位置的不同而不同,从而偏离C-C方程的限制。具体为在热带内区,气温和RH呈现强的负相关,而在中纬度地区常见气温与RH呈现强正相关[41]。在临近海洋的大陆架上,尽管RH的月平均值几乎维持在75%的常数水平,季节内却具有相当大的日变率。夏天比湿是冬天的2倍,高气温是导致夏季月平均RH低的原因。冬季具有最大的RH变率,从50%变化到100%。比湿的主导变率在天气时间尺度上常与9月和11月的过境锋面有关。RH不仅依赖于气象条件的高频变率,也依赖于海洋条件的低频变率。特别是海表温度,它不仅控制着表面的热通量,也控制着海洋环流动力学[11]。由此可知,RH分布是由动力场而不是局地气温强烈地控制的,RH的分布变化与大气环流的变化密切相关。因此,在季风区季风动力学可以通过RH的变化体现,季风开始和季风结束时RH具有显著的变化,这些变化与大尺度大气与海洋的驱动机制有关[45]。中国的气候变化深受季风的影响,湿度的季节循环首要受控于东亚季风系统,气温为中国北部和西部的主导因素,而与季风环流相关的水汽平流是南部地区的主导因素[19]。青藏高原地区的近地面RH呈现出显著的下降趋势,特别是夏季和秋季。这与高原地区气温的显著增加相对应,更与海洋水汽输送不足有关[46]

2.3 RH年循环及其变化

虽然已经开展了一些对于RH的研究,但是,大气湿度如何响应全球变暖的过程依然不清楚[36],表面RH如何响应季风系统的演化和变率也不清楚[46]。困扰这些问题的最根本原因是对大气水汽在不同尺度的行为和大气水汽如何被变暖过程影响的认识不足,因此对于RH的研究依然非常迫切[36]。特别是在RH的气候学与变化趋势及其机理这2个方面的研究中,忽略了RH周期性变化的贡献。虽然RH的周期性在已有的研究中有所提及,例如,在临近海洋的大陆架上气候态的月平均RH具有显著的年循环特征,高RH值出现在夏季[11];由AIRS得到的相对湿度气候学研究表明,南极洲RH具有可重复的年循环[47]。1838—2008年,瑞士Armagh Observatory的年平均和季节平均RH仅自1880年代开始显现出微弱的趋势,但是却存在显著的年代际或多年代际的变率,其中有些变率循环出现。变率具有2个准周期,分别是23.4~25.5年和36~51年[14],这些周期变化与大气环流[3]和温盐环流变化[14]及太阳磁循环有关[14]。上述各种尺度的周期性的结构必然会影响RH的气候学与变化趋势及其机理方面的研究结果。但是,关于RH年循环的详细研究目前还较缺乏。

因为气候变化不仅与气候变量的平均值有关,也与对应的季节循环的特征变化有关[48,49,50,51,52]。特别是对于地表气温来说,气温年循环时常占据了气温方差的90%左右(即地表气温的主要变率发生在年时间尺度上)。而研究发现气温的年循环并不是恒定不变的,气温变化更多地体现在年循环的振幅与相位的长期变化上[48,52~55]。从观测研究中可以发现北半球热带外地区气温年循环的振幅呈现出显著的衰减,而相位则出现了显著的提前[48,49]。地表气温年循环的振幅与相位的变化主要是由太阳辐射的年循环、局地条件、大气环流和二氧化碳的排放引起的[48,51,56~60]。研究发现:地表气温年循环的振幅变化可以影响地表气温的趋势估计[52],多种方法都表明近几十年欧洲部分地区春季开始提前了[61],且在中国北方地区该提前有40%~60%是由年循环本身的变化导致的[62]。因此,Thomson[48]建议在气候研究中利用去掉月平均等方式得到的气候距平需要重新考虑,毕竟年循环的振幅和相位都不是恒定不变的,特别是相位的变化本身具有年代际尺度的涨落[63,64,65,66]。那么,RH的年循环是否存在?如果存在,它的年循环是否也和地表气温的年循环一样是变化的?其变化对于RH长期趋势变化的贡献怎么样?这些问题都需要通过研究来回答。

年循环的提取方法问题在地表气温年循环的研究中已经有了比较深入的研究,考虑到中高纬度的地表气温年循环基本上具有正弦函数形式,用到的方法主要有逐年的正弦函数拟合[51,54,58~61]、复解调方法[48,61,67]、经验模态分解(Empirical Mode Decomposition,EMD)相关的方法[52,53,62,68,69]、基于回归的序列分解方法[50,70]等。因此,在下一节将总结部分目前正在使用的年循环提取方法,同时介绍一些新的信号提取方法,并且分析它们对于RH年循环提取的可行性。

3 年循环提取方法

RH年循环是周期性或准周期性信号,可以通过较多信号提取方法来获得。由于气温也是包含显著年循环的大气变量,在气温年循环的周期信号提取中成功应用的方法具有非常好的借鉴价值,并可以尝试将这些方法运用到RH年循环的提取。

3.1 气候态年平均法

气候态是对气候的一种描述和特定研究,也用作对一个区域内气候变量的特征值的定量描述[71]。气候态平均作为计算气候态的一种方法,是通过对一段时期内的天气状况做平均而得到的,比如用20年的月气温数据计算气候态,就是计算每个月20年的平均值,最后得到1年(12个月)的气候态结果[52,72]。一般来说,多年平均气候态是长时间的平均,30年以上的时间序列才可能计算气候平均态。传统意义上,RH的年循环被认为就是气候态年循环(Climatology Annual Cycle,CAC)[19],这与气温等其他气候变量类似。CAC是提取年循环运用最为广泛和便捷的方式。在周期信号占主导的气温年循环的提取中,解释方差可达80%~90%,对一个区域内长时期的气温演变状况总结的信息十分有效。在目前的RH年循环研究中,CAC仍然是主要研究手段,通过该方法能对RH的分布、变化情况有基本和概括性的了解。同时,也有相关研究通过对RH的气候态研究指出,引起RH周期性变化的原因除了太阳辐射的年周期变化外,可能还有当地季风、水汽循环等因素[19]

但CAC方法是基于假定气候变量年循环不变,因而,其忽略了年循环逐年之间的变化,且提取的信号中仍然残留大量高频变化[65,66]。根据目前的观测资料,气温、RH等气候变量的变化并不是每年都相同的,尤其是RH的时间序列表现出复杂的结构,忽略掉年循环的逐年变化并不利于我们分析RH年循环变化的原因以及对全球变暖的气候响应。且CAC提取的年循环逐年拼接后,在拼接端点处常会导致三角波和锯齿波的产生,引入了其他误差。特别是RH相对气温来说季节性周期信号较弱,得到的结果更易被误差干扰。另外,CAC在进行一年中某个月(或某天)的累计平均时,保留了在这个月(这天)多年数据的误差的平均。由于RH的时间序列振荡较大,误差较大,CAC方法不能很好地去除误差的影响,只能大致给出年循环的信息。

3.2 滑动平均

有研究采用滑动平均(moving average)的方法得出年循环,解决了CAC的年与年相连的端点问题。滑动平均过程是常用的数据滤波手段之一。对一序列a1,a2,…,an,取其中连续的q个数at-q,at-q+1,…,at-1,at为一个子集,则可得到此子集的平均值 z~t:

z~t=1qt-qtai,i=t-q,t-q+1,,t,(1)

式中:t为所取子集中最后一个数的位置,固定子集大小为q,在原数列上每次向前滑动一个值取子集并计算平均值,最终可得到该数列的滑动平均数列{ z~t}[73]。这样取出的滑动平均序列缺少第1至第q处的值,可以通过 z~i=1i1iai,i=1,2,,q-1得到,但会引入一定偏差,在运用分析中可以去除开头段来消除这种端点效应。

这种方法简单、快捷,在研究中也常被采用,比如Bonsal等[74]在气温年循环的研究中采用31天滑动平均的方法对原始数据进行滤波,将高频噪声过滤,并得到连续的年循环变化,以此讨论年循环相位的改变。

但由于原始数据振荡和高频变化对滑动平均的结果影响较大,在数据包含较大噪声时,滑动平均的结果也会受到较大影响,不能有效去除年循环[64]。在Qian等[62]的研究中,若使用日资料研究气温的变化,滑动平均的窗口较小时(5天滑动平均),得到的结果中包含大量高频波动;当窗口较大时(31天滑动平均),得到的结果较为平滑,但仍然存在较多高频变化引起的结构。在使用月平均资料时,由于高频变化较少,滑动平均的结果明显更好。对于RH时间序列来说,日资料包含大量高频涨落,用滑动平均来提取年循环的效果较差。对于月资料则可得到更为平滑的结果,但由于月资料每年只包含12个点,滑动平均可能会滤去某些月尺度的年内循环。

直接从数据的平均值得到年循环的方法虽然快捷方便,但误差较大,且必须结合其他分析方法才能进一步得到年循环的振幅、相位等特征信息。若要提取数据中的谐波成分,采用时频分析方法,则可以更为直接地得到年循环的有效信息。

3.3 傅里叶合成逐年拟合法

RH年循环包含谐波信号,采用时频分析方法也是提取年循环的思路之一,在RH年循环的研究中可以引入傅里叶合成逐年拟合法。傅里叶变换是使用最广泛的一种时频分析方法,也十分简单有效。Stine等[51]使用简单的傅里叶逐年拟合方法对气温变化的年循环进行了提取,并得到了气温变化的年循环信息。

他们采用月平均气温资料,用傅里叶变换计算每年循环中的正弦成分:

Yx=212t=0.511.5e2πit/12x(t+t0),t=0.5,,11.5,(2)

式中:x(t+t0)是去平均值的月平均温度的值,t分别对应12个月的时间。式中分子上的2是因为正频率和负频率都被考虑在内。则相位为ϕx=tan-1[Im(Yx)/Re(Yx)],振幅为Ax=|Yx|,Re为取实部算符,Im为取虚部算符,Yx是一个复数,Re(Yx)是Yx的实部,Im(Yx)是Yx的虚部。

这种方法得到的气温年循环的正弦成分对月温度资料的年内变化的解释率在陆地上达到96%,在海洋上达到90%[51]。有效地得到了气温逐年变化的年循环变化,从而得到较为可靠的年循环的振幅、相位随时间变化信息。但此方法对每一年的数据都进行拟合,在年与年连接处可能会引入锯齿波,存在过度拟合的情况。

3.4 复解调法

由于逐年拟合法是通过傅里叶变换对每年数据单独进行拟合,缺乏对于数据整体的分析,复解调法(Complex Demodulation)建立在傅里叶变换基础上,有效地弥补了这一点,且更为灵活[75]。RH年循环的提取也可以采用这种方法。

复解调法假设一个包含扰动的时间序列的周期成分为xt=Rtcos2π(f0t+ϕt) ,其中{Rt}是一个缓慢变化的振幅,{ϕt}是缓慢变化的相位。复解调的目的即为提取{Rt}和{ϕt}的近似序列。假设数据记录为xt,则相应的复解调为yt=xte-2πif0t,其中f0为中心频率,在年循环提取中,若用月资料则为1/12,所以Rt=|yt|,ϕt=tan-1[Im(yt)/Re(yt)][75]

这个方法在大气中应用较多,比如Thomson[48]采用复解调方法提取了气温的年循环,并分析其相位的变化,得出相位的平均变化和大气CO2浓度的对数变化是一致的。这个方法建立在傅里叶变换基础之上,傅里叶变换已经假定序列包含谐波成分以及满足叠加性原理。但是,在真实序列中,无论是气温还是RH序列,都包含非线性的变化。RH中不仅包含谐波,还可能存在方波和三角波(见下节详细讨论),而傅里叶变换对于非线性成分的提取的效果不佳,因而采用复解调法时要注意其对噪声中非线性成分的处理问题。

3.5 经验模态分解法/集合经验模态分解法

由于大气系统的复杂性,RH序列中存在非线性成分,提取年循环的方法也可以采用EMD和集合经验模枋分解法(Ensemble Empirical Mode Decomposition,EEMD)。EMD方法是Huang等[76,77]1998年提出的信号提取方法,其本质是从数据的特征时间尺度辨别其中的本征振荡模。EMD方法从数据x(t)中提取多个本征模函数(IMFs)ci,

x(t)=i=1nci+rn,ci(t)=ai(t)cos[ωi(t)dt],(3)

式中:rn为残余,要么是一个平均趋势,要么为一个常数。

EMD方法是一个强有力的离散分解方法,不同于以傅里叶分解为基础的时频分析方法,它不需要预设基函数,不需要提前给定中心频率,是自适应信号的时频处理方法。但EMD方法存在严重的模态混淆问题[63]:因为EMD方法的分解没有理论支撑,仅以极值分布为基础进行数据分解,所以噪声对结果影响很大,且分解出的模态不一定具有物理意义。也就是说,如果原始数据是白噪声,EMD方法也会分解得出IMFs[78]

为解决这个问题,Wu等[79]提出EEMD。在EMD方法基础上,EEMD方法通过在数据中加入有限噪声来去除混合模态问题并保留有物理意义的分解结果。步骤是:①在原始数据中加入白噪声序列;②用EMD方法分解加入白噪声的序列得到IMFs;③重复步骤①和②,但每次加入的白噪声序列不同;④计算IMFs的集合平均,并作为最终结果。方法的设计者认为白噪声集合会在时空集合平均中互相抵消,因此剩下的为所需的信号。

EEMD方法在气温年循环的提取中运用较多[52,62,63],提取出的气温调制年循环MAC中频率、相位、振幅均随时间变化。比如Qian等[62]用EEMD方法分析了1955—2003年中国观测的气温数据,得出了北京春天每10年约提前2.98天的结果,其中由于变化的年循环导致的提前是每10年1.85天。

但EEMD方法对于强噪声背景下的信号提取能力相对较弱,噪声仍然会在很大程度上影响提取的IMFs[80]。对于RH这种包含较强噪声的信号来说,EEMD提取结果的可靠性仍需进一步研究与分析。

3.6 非线性模态分解法

非线性模态分解法(Nonlinear Mode Decomposition, NMD)是抗噪声能力较强的信号提取方法之一。对于包含噪声、非线性成分的RH时间序列,NMD方法也可以运用到RH年循环的提取方面。NMD方法是Iatsenko等[80]2015年提出的,结合多种目前正在发展的时频分析技术、代用数据检验方法与谐波辨识(harmonic identification)方法,极大地提高了抗噪声性能,且能够自适应信号,使得信息分析更加完整。

NMD方法在信号分解时,假定信号包含多个谐波,且由于非线性作用会使得信号成分更加复杂。定义由相同行为引起的所有振荡成分之和,即全部的非线性模态(Nonlinear Modes,NM)为:

c(t)=A(t)ν[ϕ(t)]=A(t)hahcos[(t)+ϕh],       (4)

式中:ν[ϕ(t)]=ν[ϕ(t)+2π]为具有周期相位的周期函数,可以展开成傅里叶序列。信号s(t)就表示为:

s(t)=ici+η(t),(5)

式中:ci为第i个NM,η(t)为噪声。

NMD方法主要分为4个步骤:①基于时频分析技术(Time Frequency Representition,TFR)提取信号中的基础谐波成分;②基于基础谐波性质,找出其可能存在的所有候选谐波;③通过代用数据检验法从候选谐波中取出真的谐波,组成一组NM;④从信号中减去该组NM,重复前面3个步骤,直到达到停止条件。由于采用代用数据检验方法,可以更有效地去除噪声和没有物理意义的分解结果。

NMD方法可以在噪声较强的背景下提取信号,如分解人类血液流动信号、从人类EEG信号中去除人为信号等[75],且相对来说提取到的信号更稳健。Deng等[66]采用NMD方法提取出气温序列变化的年循环,较好地量化了气温距平统计结构,避免了CAC方法带来的年循环结构残留问题,对研究气温距平的统计性质提供有力基础。RH序列中包含非线性成分以及较多噪声,NMD方法对于RH的年循环提取来说是或许是一个较好的选择。

3.7 STL方法

另一种抗噪声能力较强的方法是STL( Seasonal trend decomposition procedure based on loess)[81]。STL方法中包含一系列局部加权平滑(locally weighted regression,or loess[82])应用,过程简洁,便于分析,且计算速率较快,可以分解含有缺测值的时间序列。对于RH的观测和再分析资料,常常面对处理缺测值和噪声的问题,STL的这些特性满足对于RH时间序列分析的需要,也是提取RH年循环方法中的良好选择之一。

将时间序列Yt 分解为3个成分:趋势成分Tt,季节成分St和残余成分Rt

Yt=Tt+St+Rt(6)

STL方法由内循环和外循环2个递归过程构成,每次内循环,季节成分和趋势成分更新一次,每次外循环通过内循环进行鲁棒性权重的计算构成,循环主要由一系列平滑操作完成,其中除了局部加权平滑(locally weighted regression,Loess[82])之外,都是相同的平滑操作。STL方法中有较多参数,并不是完全自适应的,通过参数的合适选取,对信号的季节和趋势成分的提取较为有效。

通过STL得到的年循环也随着时间变化,但其中季节成分包含的相位、振幅等信息还需要结合其他手段进行分析。因为能够将非线性时间序列中的谐波有效提取出来,有研究采用STL方法去除农业气象数据中的季节成分来研究剩下部分的性质[83]。也有研究采用STL方法分析2000—2014年的海温和表面风数据,STL方法得出的结果比标准线性回归方法更真实地得到了表面风和海温的趋势[84]。对于提取RH年循环来说,STL方法也是一个较好的选择。

4 待解决的重要问题

RH作为重要的气象变量,它的年内变化直接或者间接影响着多个大气活动的进程,如降水、云雾形成等。水汽作为大气中最强的温室气体,表征大气含水量的RH年际变化和整体趋势对气候效应也有重要的影响。在当前全球变暖的环境下,RH的响应和变化机理研究已经成为一个关注焦点。同时,我国南方的水汽变化是季风形成的一个关键环节,研究RH年循环包含了对RH的季节尺度、年际尺度的变化规律的探讨,对于提高目前的季风预报能力有巨大潜力。

未来对于RH年循环研究应该注重以下几个方面:

(1)完善RH数据资料并提高资料的可信度:从模式模拟、再分析数据集以及卫星观测3种渠道得到的全球湿度气候态大不相同,这为通过气候态研究RH年循环带来很多困难。Jiang等[85]研究得出,在对流层高层大气,通过CMIP5得到的结果比卫星的大气红外探测数据(AIRS)产品和微波临边探测器(MLS)的观测结果湿润2倍,但在对流层中层,CMIP5模式只比AIRS数据集和MLS数据集湿润10%。Chen等[86]发现,由于采用的数据和同化方式不同,欧洲中心的ERA-40数据集和美国国家环境预测中心的NCEP再分析数据集得到气候态在热带湿度的年内变率不同。因而RH的气候态会随着采用不同模型、不同数据等原因而不同,这大大降低了利用气候态表征年循环的可信度。

(2)RH年循环随时间变化的特征:气候态的年循环是固定的,不能体现出年循环随着时间变化的特点。Du等[32]采用AIRS数据,对自由对流层RH数据进行了多元线性回归和经验模分析,发现平均年循环的振幅在垂直分布上不同,而且RH随时间变化的年循环包含多个谐波,存在12个月、6个月、4个月、3个月的周期变化,其中12个月和6个月谐波的解释方差占了较大比例。除了年内变率,也有研究发现RH的年际变化。除了对近地面RH变化的研究,也有人研究对流层高层RH的变化,发现对流层高层大气也存在气候态年循环,且随着地区的分布不同[14]。除了气候态外,RH的均值也存在不同时间尺度的波动[87]。还有一些研究没有直接研究RH的年循环,而是通过每个季节的逐年变化讨论了其变化的季节性[34]

(3)RH年循环复杂结构的类型和原因:利用中国194个观测站1960—2013年的日均RH数据集,可以直接观察到RH时间序列中包含的复杂结构。在图1a中,月平均RH序列的年循环信号十分明显,其中还包含明显的半年周期信号。同时由序列的形状可以发现在单个年循环中,RH增加到达极大值时较快,但减弱到极小值时花费的时间较长,体现出RH年循环结构的非对称性特点。从图1b中可以看出RH序列存在类方波结构且年与年的连接处基本上是快速转换,这种结构一般出现在水汽容易饱和的地区。从图1c中可以看出RH序列中存在非对称的三角波结构;对于长时间RH的演变,时间序列的振幅随着时间在变化,存在显著的年代际变率,体现其非线性的特点。

(4)RH年循环有效提取方法的选择:从RH时间序列变化的复杂性可以看出现有研究提供的方法可能无法实现对RH年循环的有效提取,还需要更多的信号提取手段。研究发现:若利用多线性回归把日平均RH分解为确定性部分和随机性部分(残差项)如下:

ri(t)=ri0+k=1nai,kcos(ωkt+bi,k)+βiE(t)+αit+εi(t),(7)

式中:等式右端的前4项为确定性部分(分别为平均项、谐波变化项、ENSO作用项、线性趋势项),最后一项为残差项(随机部分)。对于全球平均的RH,方差贡献最大的是包含12个月和6个月周期的谐波项与残差项[32]。这也说明RH的年循环占据了RH方差的很大部分,不同于地表气温的是RH的高频扰动也占据了RH方差的一定部分。造成这一情形的原因是RH周期变化部分的形状并不总是正弦函数,这在前人的研究中已经涉及。例如,在瑞士Armagh Observatory年平均RH序列中发现的周期并不是正弦形式的,而是先陡峭上升到一个维持约7年的高RH水平,接着下降到一个持续30~40年的低RH水平,同时,北风和RH间存在小但是显著的反相关关系[14]。这也从上述残差序列的自相关函数的特征上得到体现,它要经过15天后指数才能衰减到0,即只有在15天或更长的时间尺度,残差序列才能很好地用白噪声刻画。这说明上述模型不能很好地捕捉到瞬变局域化的扰动(如东传的对流耦合波、热带外风暴等具有时间尺度显著长于一天)的时间演化[32]。除了这个模型,还可以采用自适应的EMD/EEMD方法进行信号提取,但对于方波信号以及强噪声背景,这种方法的提取结果仍不理想[80]。因此,RH年循环的提取方法问题值得深入研究。上节总结的提取方法都非常成熟,但是,这些方法应用的前提大都是年循环具有正弦函数形式。而RH年循环的形式很多时候都不是正弦函数,特别是由于受到季风环流等的影响,高RH与低RH之间的转换都是在非常短的时间完成的,RH年循环的形状更接近于方波或三角波(图2)。这使得基于三角函数的拟合或者EMD相关的方法在提取RH年循环时受到很多限制,需要遴选或提出更适合于RH年循环自身特点的提取方法,这是研究RH年循环的基础,也是难点。

图1   中国3个观测站月平均RH时间序列(a)观测站50434 1997年3月至2003年3月(6年)的月数据;(b)观测站59501 1977年10月至1985年10月(8年)的月数据;(c)观测站51573 1968年2月至1993年5月(25年)的月数据

Fig.1   Segments of monthly averaged relative humidity records of three stations in China(a) Records from 1997/3 to 2003/3 at station 50434; (b) Records from 1977/10 to 1985/10 at station 59501;(c) Records from 1968/2 to 1993/5 at station 51573

(5)RH相对于年循环的距平序列性质研究问题:RH的距平序列(相对于变化的年循环)性质研究问题。基于气候平均态定义的RH距平序列的研究[88,89,90,91],最早可以追朔到1994年Vattay等[88]针对地面湿度日平均资料的研究。他们发现RH距平时间序列的谱密度呈现1/fa的特征,具有标度行为,从这个特征以及从云分形的空间结构可以把RH的变化看成是自组织临界性的表现。但随后Garcia等[89]通过研究西班牙城市萨拉戈萨1951—1993年的RH日平均资料的谱密度,认为并不能得到1/fa的特征,不能看成自组织临界现象。Chen等[90]发现:采用Dtrended Fluctuation Analysis方法分析中国日均RH距平序列时,发现RH距平序列有明显不同于其他气象变量的统计性质。由于气象变量变化的非线性和非平稳性,对应气象变量的年循环的结构都是变化的,而这种年循环的结构变化在经典气象变量距平定义(假定恒定的年循环)中很难被去除掉[65,68,69]。如何量化经典的年循环和变化的年循环对于RH距平的统计性质影响[65,68,69]也是需要认真研究的内容。

图2   RH在传统定义下的年循环形式(a)正弦函数;(b)类方波;(c) 状态迅速转换

Fig.2   Annual cycle of relative humidity based on CAC(a) Sine function;(b) Rectangle-type wave; (c) States of abrupt transition

5 结 语

对于RH年循环的研究,迫切需要提高数据资料有效性,认识到RH年循环是随时间变化的,并且具有复杂结构,放弃普遍使用的以谐波为基础的信号提取手段,寻找新的方法。RH年循环所反映的季节尺度、年际尺度变化规律,以及距平序列性质,将可能进一步提高我们对RH年循环变化内因和外因的确定和分析,有助于我们对气候变化的理解,同时使我国季风预报、疾病监测以及污染调控能力达到一个更高的水平。

The authors have declared that no competing interests exist.


参考文献

[1] Trenberth K E, Jones P D, Ambenje P, et al.

Observations: Surface and atmospheric climate change

[M]//Solomon S, Qin D, Manning M, et al, eds. Climate Change 2007: The Physical Science Basis: Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, 2007: 235-336.

[本文引用: 1]     

[2] Held I M, Shell K M.

Using relative humidity as a state variable in climate feedback analysis

[J]. Journal of Climate, 2012, 25(8): 2 578-2 582.

DOI      URL      [本文引用: 1]      摘要

An approach to climate change feedback analysis is described in which tropospheric relative humidity replaces specific humidity as the state variable that, along with the temperature structure, surface albedos, and clouds, controls the magnitude of the response of global mean surface temperature to a radiative forcing. Despite being simply a regrouping of terms in the feedback analysis, this alternative perspective has the benefit of removing most of the pervasive cancellation between water and lapse-rate feedbacks seen in models. As a consequence, the individual feedbacks have less scatter than in the traditional formulation. The role of cloud feedbacks in controlling climate sensitivity is also reflected more clearly in the new formulation.
[3] Sherwood S C, Ingram W, Tsushima Y, et al.

Relative humidity changes in a warmer climate

[J]. Journal of Geophysical Research, 2010, 115(D9): 1 063.

DOI      URL      [本文引用: 3]      摘要

[1] Key climate feedback due to water vapor and clouds rest largely on how relative humidity 090407 changes in a warmer climate, yet this has not been extensively analyzed in models. General circulation models (GCMs) from the CMIP3 archive and several higher-resolution atmospheric GCMs examined here generally predict a characteristic pattern of 090407 trend with global temperature that has been reported previously in individual models, including increase around the tropopause, decrease in the tropical upper troposphere, and decrease in midlatitudes. This pattern is very similar to that previously reported for cloud cover in the same GCMs, confirming the role of 090407 in controlling changes in simulated cloud. Comparing different models, the trend in each part of the troposphere is approximately proportional to the upward and/or poleward gradient of 090407 in the present climate. While this suggests that the changes simply reflect a shift of the 090407 pattern upward with the tropopause and poleward with the zonal jets, the drying trend in the subtropics is roughly 3 times too large to be attributable to shifts of subtropical features, and the subtropical 090407 minima deepen in most models. 090407 trends are correlated with horizontal model resolution, especially outside the tropics, where they show signs of convergence and latitudinal gradients become close to available observations for GCM resolutions near T85 and higher. We argue that much of the systematic change in 090407 can be explained by the local specific humidity having been set (by condensation) in remote regions with different temperature changes, hence the gradients and trends each depend on a model's ability to resolve moisture transport. Finally, subtropical drying trends predicted from the warming alone fall well short of those observed in recent decades. While this discrepancy supports previous reports of GCMs underestimating Hadley cell expansion, our results imply that shifts alone are not a sufficient interpretation of changes.
[4] Xie B, Zhang Q, Ying Y.

Trends in precipitable water and relative humidity in China: 1979-2005

[J]. Journal of Applied Meteorology and Climatology, 2011, 50(10): 1 985-1 994.

DOI      URL      [本文引用: 3]      摘要

ABSTRACT Annual and seasonal trends of precipitable water (PW) and relative humidity (RH) at 850, 700, and 500 hPa are studied using the data from 106 radiosonde stations over China during the period 1979-2005. Analysis shows evidence of an increase in PW associated with the slight warming observed in the lower to midtroposphere over China. The northern part of China shows a significant upward trend of PW in summer, and drying of the atmosphere in winter is found in most regions over China. Annual and seasonal trends in RH at the 850-, 700-, and 500-hPa levels show no significant trends in most regions in China except for Xinjiang, which shows an upward trend, and central China, where there was a downward trend in RH at 500 hPa. It is found that changes in PW are coincident with the warming of the surface and the lower to midtroposphere. The RH in the lower to midtroposphere in most regions over China has remained steady during the most recent 30 years, as might be expected given the increasing of PW and the warming above the surface. The long-term trend of precipitation over China may be linked to the trends of PW and RH at the lower level and midlevel.
[5] Schmidt G A, Ruedy R A, Miller R L, et al.

Attribution of the present-day total greenhouse effect

[J]. Journal of Geophysical Research, 2010, 115(D20): 20 106

DOI      URL      [本文引用: 1]     

[6] Cubasch U, Meehl G A, Boer G J,et al.

Projections of future climate change

[M]//Houghton J T, Ding Y, Griggs D J, et al, eds. Climate Change 2001: The Scientific Basis: Contribution of Working Group I to the Third Assessment Report of the Intergovernmental Panel, 2001: 526-582.

[本文引用: 1]     

[7] Dai A.Global precipitation and thunderstorm frequencies.

Part I: Seasonal and interannual variations

[J]. Journal of Climate, 2001, 14(6): 1 092-1 111.

DOI      URL      [本文引用: 1]      摘要

Abstract Present and past weather reports from 651565000 stations around the globe and from the Comprehensive Ocean–Atmosphere Data Set from 1975 to 1997 were analyzed for the frequency of occurrence for and the percentage of the days with various types of precipitation (drizzle, nondrizzle, showery, nonshowery, and snow) and thunderstorms. In this paper, the mean geographical, seasonal, and interannual variations in the frequencies are documented. Drizzles occur most frequently (655%–15% of the time) over mid- and high-latitude oceans. Nonshowery precipitation is the preferred form of precipitation over the storm-track regions at northern mid- and high latitudes in boreal winter and over the Antarctic Circumpolar Current in all seasons. Showery precipitation occurs 655%–20% of the time over the oceans, as compared with < 10% over land areas except in boreal summer over Northern Hemisphere land areas, where showery precipitation and thunderstorms occur in over 20% of the days. Inferred mean precipitation inte...
[8] Dai A.Global precipitation and thunderstorm frequencies.

Part II: Diurnal variations

[J]. Journal of Climate, 2001, 14(6): 1 112-1 128.

DOI      URL      [本文引用: 1]      摘要

Three-hourly present weather reports from 15 000 stations around the globe and from the Comprehensive Ocean tmosphere Data Set from 1975 to 1997 were analyzed for diurnal variations in the frequency of occurrence for various types of precipitation (drizzle, nondrizzle, showery, nonshowery, and snow) and thunderstorms. Significant diurnal variations with amplitudes exceeding 20% of the daily mean are found over much of the globe, especially over land areas and during summer. Drizzle and nonshowery precipitation occur most frequently in the morning around 0600 local solar time (LST) over most land areas and from midnight to 0400 LST over many oceanic areas. Showery precipitation and thunderstorms occur much more frequently in the late afternoon than other times over most land areas in all seasons, with a diurnal amplitude exceeding 50% of the daily mean frequencies. Over the North Pacific, the North Atlantic, and many other oceanic areas adjacent to continents, showery precipitation is most frequent in the morning around 0600 LST, which is out of phase with land areas. Over the tropical and southern oceans, showery precipitation tends to peak from midnight to 0400 LST. Maritime thunderstorms occur most frequently around midnight. It is suggested that the diurnal variations in atmospheric relative humidity contribute to the morning maximum in the frequency of occurrence for drizzle and nonshowery precipitation, especially over land areas. Solar heating on the ground produces a late-afternoon maximum of convective available potential energy in the atmosphere that favors late-afternoon moist convection and showery precipitation over land areas during summer. This strong continental diurnal cycle induces a diurnal cycle of opposite phase in low-level convergence over large nearby oceanic areas that favors a morning maximum of maritime showery precipitation. Larger low-level convergence induced by pressure tides and higher relative humidity at night than at other times may contribute to the nighttime maximum of maritime showery and nonshowery precipitation over remote oceans far away from continents.
[9] Zhao T, Wang J, Dai A.

Evaluation of atmospheric precipitable water from reanalysis products using homogenized radiosonde observations over China

[J]. Journal of Geophysical Research: Atmospheres, 2015, 120(20): 10 703-10 727.

[本文引用: 1]     

[10] Vautard R, Yiou P, Van Oldenborgh G J.

Decline of fog, mist and haze in Europe over the past 30 years

[J]. Nature Geoscience, 2009, 2(2): 115-119.

DOI      URL      [本文引用: 1]      摘要

Surface solar radiation has undergone decadal variations since the middle of the twentieth century, producing global `dimming' and `brightening' effects. These variations presumably result from changes in aerosol burden and clouds, but the detailed processes involved have yet to be determined. Over Europe, the marked solar radiation increase since the 1980s is thought to have contributed to the observed large continental warming, but this contribution has not been quantified. Here we analyse multidecadal data of horizontal visibility, and find that the frequency of low-visibility conditions such as fog, mist and haze has declined in Europe over the past 30 years, for all seasons and all visibility ranges between distances of 0 and 8km. This decline is spatially and temporally correlated with trends in sulphur dioxide emissions, suggesting a significant contribution of air-quality improvements. Statistically linking local visibility changes with temperature variations, we estimate that the reduction in low-visibility conditions could have contributed on average to about 10-20% of Europe's recent daytime warming and to about 50% of eastern European warming. Large improvements in air quality and visibility already achieved in Europe over the past decades may mean that future reductions in the frequency of low-visibility events will be limited, possibly leading to less rapid regional warming.
[11] Virmani J I, Weisberg R H.

Relative humidity over the West Florida continental shelf

[J]. Monthly Weather Review, 2005, 133(6): 1 671-1 686.

DOI      URL      [本文引用: 4]      摘要

Observed relative humidity variations on the coastal ocean of the West Florida Continental Shelf (WFS) are examined over the 5-yr period 199809“2003. Despite considerable daily variability within seasons, the monthly mean values are nearly constant at about 75%. Summertime specific humidity is twice that during winter, so high air temperatures are responsible for the low summer monthly mean relative humidities. Winter has the greatest relative humidity variability; values range from less than 50% to over 100% as extratropical fronts move over the WFS. Saturation (and fog) occurs as warm moist air passes over colder water. Two different sensors, mounted on multiple moorings, were used to make these observations. Monthly mean values from the Rotronics MP-100F are higher than the Hygrometrix 1020SHT. In addition to sensor differences, a contributing cause to this offset appears to be the locations chosen for sensor deployment. NCEP reanalysis climatology over the WFS and land-based coastal data both show an annual cycle in monthly mean relative humidity, with higher values in summer, suggesting that the reanalysis field is influenced by land. Air09“sea fluxes over the WFS are sensitive to small spatial variability in the coastal ocean and atmosphere. The large grid spacing of the NCEP reanalysis does not capture this variability. The lack of coastal ocean data for assimilation biases the NCEP reanalysis fields toward land-based measurements. Increased spatial coverage via evolving Coastal Ocean Observing Systems should remedy this problem by providing required information for describing and understanding the complicated ocean09“atmosphere interactions that occur on continental shelves.
[12] Gultepe I.

Fog and boundary Layer Clouds: Fog Visibility and Forecasting

[M].Berlin: Birkhäuser Basel, 2007.

[本文引用: 1]     

[13] Pilié R J, Mack E J, Rogers C W, et al.

The formation of marine fog and the development of fog-stratus systems along the California coast

[J]. Journal of Applied Meteorology, 1979, 18(10): 1 275-1 286.

DOI      URL      [本文引用: 1]     

[14] Butler C J, García-Suárez A M.

Relative humidity at Armagh observatory, 1838-2008

[J]. International Journal of Climatology, 2012, 32(5): 657-668.

DOI      URL      [本文引用: 7]      摘要

We present a calibrated atmospheric humidity series for Armagh Observatory, Northern Ireland. Beginning in 1838 and continuing to the present day, this series is believed to be the longest presented so far. Daily wet-bulb and dry-bulb temperature readings are corrected for the time of reading and for known instrumental errors to determine mean monthly, seasonal and annual relative humidity (RH) at 09:00 GMT from 1838 to 2008 and at 21:00 GMT from 1844 to 1964.The mean seasonal and annual time series of RH at this site show only marginal evidence for a trend in RH since the 1880s, however, there is significant decadal and multi-decadal variability, some of which appears to be cyclic. Wavelet analysis indicates the presence of two quasi-periodic components in RH, one ranging from 23.4 to 25.5 years and the other from 36 to 51 years. We discuss potential physical origins of these cycles such as changes in the atmospheric and thermohaline circulation and the solar magnetic cycle. There is a small but significant anti-correlation between RH and northerly winds.Two conspicuous peaks in RH during the 19th century coincide with major infestations of potato blight in Ireland. Copyright 2011 Royal Meteorological Society
[15] Lowen A C, Mubareka S, Steel J, et al.

Influenza virus transmission is dependent on relative humidity and temperature

[J]. PLoS Pathog, 2007, 3(10): e151.

DOI      URL      PMID      [本文引用: 1]      摘要

Using the guinea pig as a model host, we show that aerosol spread of influenza virus is dependent upon both ambient relative humidity and temperature. Twenty experiments performed at relative humidities from 20% to 80% and 5 °C, 20 °C, or 30 °C indicated that both cold and dry conditions favor transmission. The relationship between transmission via aerosols and relative humidity at 20 °C is similar to that previously reported for the stability of influenza viruses (except at high relative humidity, 80%), implying that the effects of humidity act largely at the level of the virus particle. For infected guinea pigs housed at 5 °C, the duration of peak shedding was approximately 40 h longer than that of animals housed at 20 °C; this increased shedding likely accounts for the enhanced transmission seen at 5 °C. To investigate the mechanism permitting prolonged viral growth, expression levels in the upper respiratory tract of several innate immune mediators were determined. Innate responses proved to be comparable between animals housed at 5 °C and 20 °C, suggesting that cold temperature (5 °C) does not impair the innate immune response in this system. Although the seasonal epidemiology of influenza is well characterized, the underlying reasons for predominant wintertime spread are not clear. We provide direct, experimental evidence to support the role of weather conditions in the dynamics of influenza and thereby address a long-standing question fundamental to the understanding of influenza epidemiology and evolution. In temperate regions influenza epidemics recur with marked seasonality: in the northern hemisphere the influenza season spans November to March, while in the southern hemisphere epidemics last from May until September. Although seasonality is one of the most familiar features of influenza, it is also one of the least understood. Indoor crowding during cold weather, seasonal fluctuations in host immune responses, and environmental factors, including relative humidity, temperature, and UV radiation have all been suggested to account for this phenomenon, but none of these hypotheses has been tested directly. Using the guinea pig model, we have evaluated the effects of temperature and relative humidity on influenza virus spread. By housing infected and na07ve guinea pigs together in an environmental chamber, we carried out transmission experiments under conditions of controlled temperature and humidity. We found that low relative humidities of 20%–35% were most favorable, while transmission was completely blocked at a high relative humidity of 80%. Furthermore, when guinea pigs were kept at 5 °C, transmission occurred with greater frequency than at 20 °C, while at 30 °C, no transmission was detected. Our data implicate low relative humidities produced by indoor heating and cold temperatures as features of winter that favor influenza virus spread.
[16] Peixoto J P, Oort A H.

The climatology of relative humidity in the atmosphere

[J]. Journal of Climate, 1996, 9(12): 3 443-3 463.

DOI      URL      [本文引用: 3]      摘要

Abstract The present paper deals with the analysis of the time-average relative humidity fields in the atmosphere. Twice-daily estimates of relative humidity are used. After some theoretical considerations on the relations between relative humidity, other moisture parameters, and temperature, a critical analysis of the various sources of data is made considering their possible limitations. Various methods of computing relative humidity are formulated and discussed. The global distribution of relative humidity at various levels shows that it is not zonally uniform with centers of various intensities at all latitudes. The global maps show maxima in the equatorial zone and minima in the dry subtropical belts around 30°N and 30°S. The land–sea contrast and variations related to the orographic relief are also apparent. The general pattern of relative humidity is similar at all levels but its magnitude decreases with altitude. The seasonal analyses show a similar pattern as the annual analyses but are slightly ...
[17] Gaffen D J, Ross R J.

Climatology and trends of US surface humidity and temperature

[J]. Journal of Climate, 1999, 12(3): 811-828.

DOI      URL      [本文引用: 2]     

[18] New M, Hulme M, Jones P.Representing twentieth-century space-time climate variability.

Part II: Development of 1901-96 monthly grids of terrestrial surface climate

[J]. Journal of Climate, 2000, 13(13): 2 217-2 238.

DOI      URL      [本文引用: 1]     

[19] Wang J X L, Gaffen D J.

Late-twentieth-century climatology and trends of surface humidity and temperature in China

[J]. Journal of Climate, 2001, 14(13): 2 833-2 845.

DOI      URL      [本文引用: 4]      摘要

Climatological surface temperature and humidity variables for China are presented based on 6-hourly data from 196 stations for the period of 1961-90. Seasonal and annual means for daytime, nighttime, and the full day are shown. The seasonal cycle of moisture is primarily controlled by the east Asia monsoon system, with dominant factors of temperature change in northern and western China and of moisture advection associated with monsoon circulations in the southeast.Trends during 1951-94 are estimated for each station and for four regions of the country, with attention paid to the effects of changes in instrumentation, observing time, and station locations. The data show evidence of increases in both temperature and atmospheric moisture content. Temperature and specific humidity trends are larger at nighttime than daytime and larger in winter than summer. Moisture increases are observed over most of China. The increases are several percent per decade for specific humidity, and several tenths of a degree per decade for temperature and dewpoint. Increasing trends in summertime temperature and humidity contribute to upward trends in apparent temperature, a measure of human comfort.
[20] Worley S J, Woodruff S D, Reynolds R W, et al.

ICOADS release 2.1 data and products

[J]. International Journal of Climatology, 2005, 25(7): 823-842.

DOI      URL      [本文引用: 1]      摘要

The International Comprehensive Ocean–Atmosphere Data Set (ICOADS), release 2.1 (1784–2002), is the largest available set of in situ marine observations. Observations from ships include instrument measurements and visual estimates, and data from moored and drifting buoys are exclusively instrumental. The ICOADS collection is constructed from many diverse data sources, and made inhomogeneous by the changes in observing systems and recording practices used throughout the period of record, which is over two centuries. Nevertheless, it is a key reference data set that documents the long-term environmental state, provides input to a variety of critical climate and other research applications, and serves as a basis for many associated products and analyses. The observational database is augmented with higher level ICOADS data products. The observed data are synthesized to products by computing statistical summaries, on a monthly basis, for samples within 2° latitude × 2° longitude and 1° × 1° boxes beginning in 1800 and 1960 respectively. For each resolution the summaries are computed using two different data mixtures and quality control criteria. This partially controls and contrasts the effects of changing observing systems and accounts for periods with greater climate variability. The ICOADS observations and products are freely distributed worldwide. The standard ICOADS release is supplemented in several ways; additional summaries are produced using experimental quality control, additional observations are made available in advance of their formal blending into a release, and metadata that define recent ships' physical characteristics and instruments are available. Copyright 08 2005 Royal Meteorological Society
[21] Van Wijngaarden W A, Vincent L A.

Examination of discontinuities in hourly surface relative humidity in Canada during 1953-2003

[J].Journal of Geophysical Research, 2005, 110(D22): 3 093-3 109.

DOI      URL      [本文引用: 2]      摘要

[1] Hourly values of relative humidity recorded at 75 stations across Canada were examined. Data were checked for possible discontinuities arising because of changes in procedures and instruments. It was found that the replacement of the psychrometer by the dewcel has produced a decreasing step in relative humidity at a number of stations. The historical records were closely examined to retrieve the dewcel installation date, and a procedure based on regression models was applied to determine if it corresponds to a significant step. Results show that there are more stations experiencing a dewcel step in the winter than in the summer. Examination of the trends also reveals that the step often accentuates the decreasing trends originally observed during winter and spring. However, significant steps taken into account, it appears that the relative humidity still decreased by several percent in the spring during 1953 2003 in western Canada. It seems that the southern and coastal stations are not as much affected by this change of instruments.
[22] Vincent L A, van Wijngaarden W A, Hopkinson R.

Surface temperature and humidity trends in Canada for 1953-2005

[J]. Journal of Climate, 2007, 20(20): 5 100-5 113.

DOI      URL      [本文引用: 1]      摘要

Annual and seasonal trends in temperature, dewpoint, relative humidity, and specific humidity are presented for the period 195309“2005. The analysis uses hourly observations from 75 climatological stations across Canada. Data were examined for discontinuities due to changes in instruments and observing practice. It was found that the main discontinuity corresponds to the replacement of the psychrometer by the dewcel in the early 1970s, which created an artificial negative step in relative humidity and dewpoint at many locations. After accounting for these discontinuities, the results of trend analysis show evidence of an increase in air moisture content associated with the warming observed in the country. During winter and spring, the significant warming in the western and southern regions is accompanied by an increase in dewpoint and specific humidity and by a decrease in relative humidity; in summer, warming is observed in the southeast and it is associated with significant positive trends in dewpoint and specific humidity. Although there is no strong evidence of a greater nighttime warming in Canada over 195309“2005, the nighttime dewpoint and specific humidity trends are slightly larger than the daytime trends, especially during the spring and summer.
[23] Willett K M, Jones P D, Gillett N P, et al.

Recent changes in surface humidity: Development of the HadCRUH dataset

[J]. Journal of Climate, 2008, 21(20): 5 364-5 383.

DOI      URL      [本文引用: 1]      摘要

Water vapor constitutes the most significant greenhouse gas, is a key driver of many atmospheric processes, and hence, is fundamental to understanding the climate system. It is a major factor in human “heat stress,” whereby increasing humidity reduces the ability to stay cool. Until now no truly global homogenized surface humidity dataset has existed with which to assess recent changes. The Met Office Hadley Centre and Climatic Research Unit Global Surface Humidity dataset (HadCRUH), described herein, provides a homogenized quality controlled near-global 5° by 5° gridded monthly mean anomaly dataset in surface specific and relative humidity from 1973 to 2003. It consists of land and marine data, and is geographically quasi-complete over the region 60°N–40°S. Between 1973 and 2003 surface specific humidity has increased significantly over the globe, tropics, and Northern Hemisphere. Global trends are 0.11 and 0.07 g kg611 (10 yr)611 for land and marine components, respectively. Trends are consistently larger in the tropics and in the Northern Hemisphere during summer, as expected: warmer regions exhibit larger increases in specific humidity for a given temperature change under conditions of constant relative humidity, based on the Clausius–Clapeyron equation. Relative humidity trends are not significant when averaged over the landmass of the globe, tropics, and Northern Hemisphere, although some seasonal changes are significant. A strong positive bias is apparent in marine humidity data prior to 1982, likely owing to a known change in reporting practice for dewpoint temperature at this time. Consequently, trends in both specific and relative humidity are likely underestimated over the oceans.
[24] Willett K M, Jones P D, Thorne P W, et al.

A comparison of large scale changes in surface humidity over land in observations and CMIP3 general circulation models

[J]. Environmental Research Letters, 2010, 5(2): 025210.

DOI      URL      [本文引用: 2]     

[25] Simmons A J, Willett K M, Jones P D, et al.

Low-frequency variations in surface atmospheric humidity, temperature, and precipitation: Inferences from reanalyses and monthly gridded observational data sets

[J]. Journal of Geophysical Research, 2010, 115(D1): 1-21.

[本文引用: 2]     

[26] Randall D A, Wood R A, Bony S, et al.Climate models and their evaluation[M]//Climate Change 2007: The Physical Science basis. Contribution of Working Group I to the Fourth Assessment Report of the IPCC (FAR). Cambridge: Cambridge University Press, 2007: 589-662.

[本文引用: 1]     

[27] Manabe S, Wetherald R T.

Thermal equilibrium of the atmosphere with a given distribution of relative humidity

[J]. Journal of the Atmospheric Sciences, 1967, 24(3): 241-259.

DOI      URL      [本文引用: 1]      摘要

Radiative convective equilibrium of the atmosphere with a given distribution of relative humidity is computed as the asymptotic state of an initial value problem.The results show that it takes almost twice as long to reach the state of radiative convective equilibrium for the atmosphere with a given distribution of relative humidity than for the atmosphere with a given distribution of absolute humidity.Also, the surface equilibrium temperature of the former is almost twice as sensitive to change of various factors such as solar constant, COcontent, Ocontent, and cloudiness, than that of the latter, due to the adjustment of water vapor content to the temperature variation of the atmosphere.According to our estimate, a doubling of the COcontent in the atmosphere has the effect of raising the temperature of the atmosphere (whose relative humidity is fixed) by about 2C. Our model does not have the extreme sensitivity of atmospheric temperature to changes of COcontent which was adduced by M ller.
[28] Manabe S, Wetherald R T.

The effects of doubling the CO2 concentration on the climate of a general circulation model

[J]. Journal of the Atmospheric Sciences, 1975, 32(1): 3-15.

DOI      URL      [本文引用: 1]     

[29] Held I M, Soden B J.

Water vapor feedback and global warming 1

[J]. Annual Review of Energy and the Environment, 2000, 25(1): 441-475.

DOI      URL      [本文引用: 1]     

[30] Dai A.

Recent climatology, variability, and trends in global surface humidity

[J]. Journal of Climate, 2006, 19(15): 3 589-3 606.

DOI      URL      [本文引用: 2]     

[31] Zhao T, Dai A, Wang J.

Trends in tropospheric humidity from 1970 to 2008 over China from a homogenized radiosonde dataset

[J]. Journal of Climate, 2012, 25(13): 4 549-4 567.

DOI      URL      [本文引用: 1]     

[32] Du J, Cooper F, Fueglistaler S.

Statistical analysis of global variations of atmospheric relative humidity as observed by AIRS

[J]. Journal of Geophysical Research, 2012, 117(D12): 13 300.

DOI      URL      [本文引用: 4]      摘要

[1] Atmospheric water vapor plays a key role in the climate. Numerical model calculations suggest that global mean relative humidity (RH) stays approximately constant in global warming scenarios. Here, we analyze the September 2002 to April 2011 daily mean free tropospheric relative humidity (RH) data from the Atmospheric Infrared Sounder (AIRS), version 5 level 3 data at a spatial resolution of 100°/100° longitude/latitude. We perform a multiple linear regression analysis with annual harmonics, trend and an El-Ni01±o090009Southern Oscillation (ENSO) index. For the mean annual cycle, we find strong compensation of local variability upon global averaging, leaving an amplitude for the global mean of about 4% RH (peak to peak) at 300 hPa, with the minimum in February. The amplitude decreases to about 2.5%RH at 500 hPa (with minimum shifted to boreal summer/early fall) and less than 2%RH further below. For ENSO, the local changes compensate strongly upon global averaging. Computational limitations restrict detailed analysis of the residual to the zonal mean residual, which we interpret with a stochastic model that takes the latitudinal covariance into account. We find that the relation between zonal mean residual RH variations (understood as the consequence of 090004weather090005 and subseasonal variability) and their global mean is equivalent to that of about 9 independent random time series with appropriate variances. The residual contributes more to the variance of the global mean than the harmonics and ENSO on all free tropospheric levels except 400 hPa. Our results apply to the AIRS version 5 data as reported, and possible problems in that data are discussed.
[33] Willett K M, Williams Jr C N, Dunn R J H, et al.

HadISDH: An updateable land surface specific humidity product for climate monitoring

[J]. Climate of the Past, 2013, 9(2): 657.

DOI      URL      [本文引用: 1]      摘要

HadISDH is a near-global land surface specific humidity monitoring product providing monthly means from 1973 onwards over large-scale grids. Presented herein to 2012, annual updates are anticipated. HadISDH is an update to the land component of HadCRUH, utilising the global high-resolution land surface station product HadISD as a basis. HadISD, in turn, uses an updated version of NOAA's Integrated Surface Database. Intensive automated quality control has been undertaken at the individual observation level, as part of HadISD processing. The data have been subsequently run through the pairwise homogenisation algorithm developed for NCDC's US Historical Climatology Network monthly temperature product. For the first time, uncertainty estimates are provided at the grid-box spatial scale and monthly timescale. <br><br> HadISDH is in good agreement with existing land surface humidity products in periods of overlap, and with both land air and sea surface temperature estimates. Widespread moistening is shown over the 1973???2012 period. The largest moistening signals are over the tropics with drying over the subtropics, supporting other evidence of an intensified hydrological cycle over recent years. Moistening is detectable with high (95%) confidence over large-scale averages for the globe, Northern Hemisphere and tropics, with trends of 0.089 (0.080 to 0.098) g kg???1 per decade, 0.086 (0.075 to 0.097) g kg???1 per decade and 0.133 (0.119 to 0.148) g kg???1 per decade, respectively. These changes are outside the uncertainty range for the large-scale average which is dominated by the spatial coverage component; station and grid-box sampling uncertainty is essentially negligible on large scales. A very small moistening (0.013 (???0.005 to 0.031) g kg???1 per decade) is found in the Southern Hemisphere, but it is not significantly different from zero and uncertainty is large. When globally averaged, 1998 is the moistest year since monitoring began in 1973, closely followed by 2010, two strong El Ni??o years. The period in between is relatively flat, concurring with previous findings of decreasing relative humidity over land.
[34] Fatichi S, Molnar P, Mastrotheodoros T, et al.

Diurnal and seasonal changes in near-surface humidity in a complex orography

[J]. Journal of Geophysical Research, 2015, 120(6): 2 358-2 374.

DOI      URL      [本文引用: 3]      摘要

Abstract Changes in near-surface humidity are evaluated in relation to changes in other meteorological variables such as air temperature, precipitation, and sunshine duration for a 31 ear period (1981 2011) using hourly time series recorded in Switzerland. Trends in meteorological variables are analyzed at the seasonal and subdaily scale, and changes in water vapor are tested for Clausius-Clapeyron scaling. Results show a marked seasonality of climatic changes with a trend toward warmer, clearer, and drier near-surface atmosphere strengthening from January to June. During this period considerable negative trends in relative humidity are detected. An abrupt shift in climatic trends occurs during the month of July, after which warming trends are considerably smaller and relative humidity remains constant or increases. An evaluation of reanalysis products demonstrates strong consistency with station observations and when combined with teleconnection indices supports the hypothesis that shifts in the general circulation patterns rather than local feedbacks are the principal drivers of observed seasonality of climate change. However, local feedbacks can play a role in enhancing summer humidity and convective activity (lightning and evening precipitation). Subdaily changes are significant with more pronounced daytime trends of temperature and humidity during spring and the opposite during summer. The strong seasonality and diurnal variability of changes in near-surface meteorology over the last three decades in Switzerland warns against climate change impact assessments at local scale that consider uniform changes in air temperature and assume Clausius-Clapeyron scaling for near-surface humidity.
[35] Wypych A.

Twentieth century variability of surface humidity as the climate change indicator in Kraków (Southern Poland)

[J]. Theoretical and Applied Climatology, 2010, 101(3/4): 475-482.

DOI      URL     

[36] Vicente-Serrano S M, Azorin-Molina C, Sanchez-Lorenzo A, et al.

Temporal evolution of surface humidity in Spain: Recent trends and possible physical mechanisms

[J].Climate Dynamics, 2014, 42(9/10): 2 655-2 674.

DOI      URL      [本文引用: 4]      摘要

We analyzed the evolution of surface relative humidity (RH) and specific humidity (q) in Spain, based on complete records available from the State Meteorological Agency of Spain. The surface RH records used span the period 1920 2011, but because of spatial and temporal constraints in the dataset we used a subset of the data, covering the period 1961 2011. The subset contained 50 monthly series of RH, which were created through a process of quality control, reconstruction and homogenization. The data shows that there was a large decrease in RH over mainland Spain from 1961 to 2011, which was greatest in spring and summer. In contrast, there was no overall change in the specific humidity in this period, except in spring, when an increase was observed. The decrease in RH affected the entire country, but the changes in specific humidity were less homogeneous. For specific humidity there was a general increase in the northern and eastern parts of Spain, whereas negative trends dominated in the central and southern areas, mainly during the summer months. The results suggest that an increase in the water holding capacity of the atmosphere as a consequence of warming during recent decades has not been accompanied by an increase in the surface water vapor content, probably because the supply of water vapor from the main terrestrial and oceanic areas has been constrained. We discuss the implications of these findings for evapotranspiration processes, precipitation and water management in Spain.
[37] Shiu C J, Liu S C, Chen J P.

Diurnally asymmetric trends of temperature, humidity, and precipitation in Taiwan

[J]. Journal of Climate, 2009, 22(21): 5 635-5 649.

DOI      URL      [本文引用: 1]      摘要

In this work, 45 years (1961 2005) of hourly meteorological data in Taiwan, including temperature, humidity, and precipitation, have been analyzed with emphasis on their diurnal asymmetries. A long-term decreasing trend for relative humidity (RH) is found, and the trend is significantly greater in the nighttime than in the daytime, apparently resulting from a greater warming at night. The warming at night in three large urban centers is large enough to impact the average temperature trend in Taiwan significantly between 1910 and 2005. There is a decrease in the diurnal temperature range (DTR) that is largest in major urban areas, and it becomes smaller but does not disappear in smaller cities and offshore islands. The nighttime reduction in RH is likely the main cause of a significant reduction of fog events over Taiwan. The smaller but consistent reductions in DTR and RH in the three off-coast islands suggests that, in addition to local land use changes, a regional-scale process such as the indirect effect of anthropogenic aerosols may also contribute to these trends. A reduction in light precipitation (10 mm h-1) are found over Taiwan and the offshore islands. The changes in precipitation are similar to the changes of other areas in Asia, but they are different from those of the United States, Europe, and the tropical oceans. The latter do not show any reduction in light precipitation.
[38] Hartmann D L, Klein Tank A M G, Rusticucci M, et al. Climate Change 2013: the Physical Science Basis: Working Group I Contribution to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change[M]. Cambridge: Cambridge University Press, 2013.

[本文引用: 1]     

[39] Nojarov P.

Variations in precipitation amounts, atmosphere circulation, and relative humidity in high mountainous parts of Bulgaria for the period 1947-2008

[J]. Theoretical and Applied Climatology, 2012, 107(1/2): 175-187.

DOI      URL      [本文引用: 1]     

[40] Talaee P H, Sabziparvar A A, Tabari H.

Observed changes in relative humidity and dew point temperature in coastal regions of Iran

[J]. Theoretical and Applied Climatology, 2012, 110(3): 385-393.

DOI      URL      [本文引用: 2]      摘要

Abstract) time series at ten coastal weather stations in Iran during 1966–2005. The serial structure of the data was considered, and the significant serial correlations were eliminated using the trend-free pre-whitening method. The results showed that annual RH increased by 1.03 and 0.2865%/decade at the northern and southern coastal regions of the country, respectively, while annual increased by 0.29 and 0.15°C per decade at the northern and southern regions, respectively. The significant trends were frequent in the series, but they were observed only at 2 out of the 50 RH series. The results showed that the difference between the results of the parametric and nonparametric tests was small, although the parametric test detected larger significant trends in the RH and time series. Furthermore, the differences between the results of the trend tests were not related to the normality of the statistical distribution.
[41] Pfahl S, Niedermann N.

Daily covariations in near-surface relative humidity and temperature over the ocean

[J]. Journal of Geophysical Research, 2011, 116(D19): 1 441-1 458.

DOI      [本文引用: 3]      摘要

[1] Changes in atmospheric relative humidity in concert with temperature changes in a future climate may have large consequences for the water vapor feedback, the hydrological cycle, and its interaction with weather systems. This study contributes to the basic understanding of the relationship between temperature and humidity by investigating the processes leading to synoptic-scale covariations of the two variables close to the ocean surface. Daily data from in situ observations between 1000°S and 5000°N and global ERA-Interim reanalyses are used. Correlations between temperature and both specific and relative humidity are calculated. The results from the two data sets appear to be greatly consistent. They show strong anticorrelation between temperature and relative humidity (RH) in the inner tropics with minimum correlation coefficients below 0908080.8. In midlatitudes, there are large areas where the correlation coefficient of temperature and RH is positive and greater than 0.6. The anticorrelation in the tropics is found to be related to convective precipitation, which, on the one hand, leads to local temperature decrease due to vertical mixing and reduction of solar radiation by clouds. On the other hand, rainfall is associated with an increase in boundary layer RH. Over the midlatitude ocean, daily temperature variations are mainly controlled by meridional transport, as shown with the help of backward trajectories. Moreover, advection of cold air typically goes along with vertical moisture transport, either due to large-scale subsidence or turbulent mixing, causing a reduction of near-surface RH. All together, this dynamical effect induces the positive temperature-RH correlation.
[42] Song Y, Liu Y, Ding Y.

A study of surface humidity changes in China during the recent 50 years

[J]. Acta Meteorologica Sinica, 2012, 26(5): 541-553.

DOI      URL      [本文引用: 2]      摘要

Abstract) shows basically an exponential growth, according to the Clausius-Clapeyron equation. The nationwide average moistening rate in winter is obviously less than the annual average rate and the summer rate. There are some regional differences in trends in different regions of China. For example, in central and eastern parts of China exhibits a reducing trend in summer, consistent with the weakening trend of temperature in these areas; (2) except parts of South China and Jianghuai Region in eastern China, unanimously increasing trends of annual and winter specific humidity () are found in most of China, especially in western China. In summer, except parts of Northeast China, Northwest China, and some areas over the Qinghai-Tibetan Plateau, the decreasing and drying trends are significant in most of China, which is not consistent with the global mean situation; (3) the surface relative humidity (RH) in most of China shows a reducing trend. One of the major reasons for the reduction of RH is that the increasing rates of are smaller than those of . Nonetheless, upward trends of RH in central and eastern China mainly due to the cooling temperature and rising in these regions are observed in summer, leading to more precipitation. From about 2003 or so, has remarkably increased while has sharply decreased in most parts of China; therefore, RH has reduced to a great extent. This may be closely related to the persistent growth of drought areas in China in the recent 10 years.
[43] Mao R, Gong D Y, Zhao T, et al.

Trends in the frequency of high relative humidity over China: 1979-2012

[J]. Journal of Climate, 2015, 28(24): 9 816-9 837.

DOI      URL      [本文引用: 1]      摘要

High relative humidity (HRH) is defined as a relative humidity of at least 80%, which is often associated with the occurrence of cloud layers. Thus, the frequency of HRH and its changes in the troposphere may be related to the occurrence frequency of cloud layers and their changes. In this study, trends in the frequency of HRH (defined as days with relative humidity 09‰0680%) over China from the surface to the midtroposphere (09‰06400 hPa) from 1979 to 2012 were analyzed using a homogenized humidity dataset for spring (March-May), summer (June-August), autumn (September-November), and winter (December-February). The results for the ground level indicate decreasing trends at most stations in southeastern China in spring and in northern China in summer. In the lower troposphere (850 and 700 hPa), most stations over China exhibit positive trends in summer, autumn, and winter. For the midtroposphere (500-400 hPa), increasing trends dominate over China in spring, summer, and autumn. Finally, six reanalysis datasets, the NCEP-NCAR, NCEP-DOE, CFSR, ERA-Interim, MERRA, and JRA-55 datasets, were compared with the observed increasing trends in HRH frequency in the low-to-middle troposphere. Similar increasing trends in HRH frequency in the reanalysis datasets and the homogenized humidity data are observed in certain seasons and for certain regions. These results are consistent with the increasing low-to-middle cloud amounts in recent decades.
[44] McCarthy M P, Toumi R.

Observed interannual variability of tropical troposphere relative humidity

[J]. Journal of Climate, 2004, 17(16): 3 181-3 191.

DOI      URL      [本文引用: 1]      摘要

Relative humidity fields from the High-Resolution Infrared Radiation Sounder (HIRS) flown on NOAA series satellites since 1979 have been used to study the seasonal aspects of the interannual variability of relative humidity in the tropical troposphere. The El Ni o Southern Oscillation (ENSO) is the only statistically identifiable physical mechanism of such variability. Boreal winter (December February) relative humidity variations during an ENSO event follow patterns of anomalous convection and large-scale upper-level circulation. During El Ni o (La Ni a) regions of large negative (positive) relative humidity anomalies exist at subtropical latitudes over the Pacific Ocean. These are not always balanced by increases (decreases) in humidity near the equator. NCEP NCAR reanalysis temperatures are used to separate observed changes in relative humidity into contributions from tropospheric temperature versus the contribution from changes in water vapor content. The authors find that at subtropical latitudes variations in temperature contribute between 50% and 70% of the observed change in relative humidity. It is also shown that large relative humidity anomalies exist over the equatorial Indian, Atlantic, and far east Pacific Oceans during the summer season (June August) following an ENSO event. Ocean atmosphere dynamics coupled with the seasonal cycle of relative humidity explain the existence of the long-lasting effects of ENSO in the atmosphere. The authors argue that observed linear trends in regional and tropical mean relative humidity are unlikely to be due solely to ENSO or a simple intensification of the hydrological cycle.<HR ALIGN="center" WIDTH="30%">
[45] Broman D, Rajagopalan B, Hopson T.

Spatiotemporal variability and predictability of relative humidity over West African Monsoon region

[J]. Journal of Climate, 2014, 27(14): 5 346-5 363.

DOI      URL      [本文引用: 1]     

[46] You Q, Min J, Lin H, et al.

Observed climatology and trend in relative humidity in the central and eastern Tibetan Plateau

[J]. Journal of Geophysical Research, 2015, 120(9): 3 610-3 621.

DOI      URL      [本文引用: 2]      摘要

Abstract Monthly surface relative humidity (RH) data for 71 stations in the Tibetan Plateau (TP) provided by the National Meteorological Information Center/China Meteorological Administration are compared with corresponding grid points from the National Center for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR hereafter) reanalysis. Mean climatologies, interannual variabilities, and trends calculated by the Mann-Kendal method are analyzed during 1961–2013. The annual regional long-term mean surface RH is 55.3%, with a clear maximum in summer (66.4%) and minimum in winter (44.9%). Compared with observations, NCEP/NCAR overestimates RH in all seasons, especially in spring (18.2%) and winter (17.8%). Mean annual regional surface RH has decreased by 610.23%65decade611 and even more rapidly in summer (610.60%65decade611) and autumn (610.39%65decade611). The reduction of surface RH is also captured by the NCEP/NCAR reanalysis at the surface, 400, 500, and 60065hPa. A particularly sharp reduction of RH since the mid-1990s is evident in both data sets, in line with rapid warming over the plateau. This suggests that moisture supply to the plateau from the Arabian Sea and the Bay of Bengal is limited and that variability and trends of surface RH over the TP are not uniquely driven by the Clausius-Clapeyron relationship.
[47] Gettelman A, Walden V P, Miloshevich L M, et al.

Relative humidity over Antarctica from radiosondes, satellites, and a general circulation model

[J]. Journal of Geophysical Research, 2006, 111(D9): 1 435-1 453.

DOI      URL      [本文引用: 1]      摘要

[1] Radiosonde measurements are used to validate measurements of relative humidity (RH) over Antarctica from the Atmospheric Infrared Sounder (AIRS) satellite instrument. Radiosonde observations are corrected for most known biases but still have a solar heating dry bias of up to 8% relative to other instruments. AIRS reproduces the observations of temperature and relative humidity with good fidelity. There is a 20% moist bias to the data in the upper troposphere relative to radiosonde measurements, but it is within the standard deviation of the measurements. Probability distribution functions of RH from radiosondes and AIRS are similar, suggesting that variability over Antarctica is well reproduced by the satellite. AIRS data are also compared to simulations from the Community Atmosphere Model version 3 (CAM3) and are found to be significantly moister than the model, although the model does not allow supersaturation with respect to ice or liquid water. A climatology from AIRS indicates that it has a repeatable annual cycle over Antarctica. Supersaturation with respect to ice is very common over the continent, particularly in winter, where it might occur almost half the time in the troposphere. This may affect the quantity and isotopic composition of ice over Antarctica.
[48] Thomson D J.

The seasons, global temperature, and precession

[J]. Science, 1995, 268(5 207): 59.

DOI      URL      PMID      [本文引用: 7]      摘要

Analysis of instrumental temperature records beginning in 1659 shows that in much of the world the dominant frequency of the seasons is one cycle per anomalistic year (the time from perihelion to perihelion, 365.25964 days), not one cycle per tropical year (the time from equinox to equinox, 365.24220 days), and that the timing of the annual temperature cycle is controlled by perihelion. The assumption that the seasons were timed by the equinoxes has caused many statistical analyses of climate data to be badly biased. Coherence between changes in the amplitude of the annual cycle and those in the average temperature show that between 1854 and 1922 there were small temperature variations, probably of solar origin. Since 1922, the phase of the Northern Hemisphere coherence between these quantities switched from 0 to 180 and implies that solar variability cannot be the sole cause of the increasing temperature over the last century. About 1940, the phase patterns of the previous 300 years began to change and now appear to be changing at an unprecedented rate. The average change in phase is now coherent with the logarithm of atmospheric CO$_2$ concentration.
[49] Thomson D J.

Dependence of global temperatures on atmospheric CO2 and solar irradiance

[J]. Proceedings of the National Academy of Sciences, 1997, 94(16): 8 370-8 377.

DOI      URL      PMID      [本文引用: 2]      摘要

Changes in global average temperatures and of the seasonal cycle are strongly coupled to the concentration of atmospheric CO2. I estimate transfer functions from changes in atmospheric CO2and from changes in solar irradiance to hemispheric temperatures that have been corrected for the effects of precession. They show that changes from CO2over the last century are about three times larger than those from changes in solar irradiance. The increase in global average temperature during the last century is at least 20 times the SD of the residual temperature series left when the effects of CO2and changes in solar irradiance are subtracted.
[50] Barbosa S M.

Changing seasonality in Europe’s air temperature

[J]. The European Physical Journal-Special Topics, 2009, 174(1): 81-89.

DOI      URL      [本文引用: 2]      摘要

Climate change is expected to involve not only changes in the mean of climate parameters, but also in the characteristics of the corresponding seasonal cycle. However, the discrimination from an observational record of long-term changes in the mean and low-frequency variations in the seasonal pattern is a challenging task, requiring the application of specific statistical methods. In this work, a time series decomposition method based on autoregression is applied in order to obtain a flexible description of seasonal variability from European temperature records. The method is based on the dynamic linear model representation for an autoregressive process and is particularly useful for isolating time-varying cycles in climate time series, allowing to retrieve fluctuations in the amplitude and phase of the periodic components and to assess their statistical significance. This approach is utilised in the analysis of long time series of daily mean temperature from the ECA (European Climate Assessment) project. Seasonality in Europe’s air temperature is characterised by an annual cycle with a stable phase but considerable inter-annual and inter-decadal variability. In particular, the annual amplitude was highest in the 1940’s and exhibits a distinct minimum around 1975, coincident with the climatic regime shift of the mid-1970’s.
[51] Stine A R, Huybers P, Fung I Y.

Changes in the phase of the annual cycle of surface temperature

[J]. Nature, 2009, 457(7 228): 435-440.

DOI      URL      PMID      [本文引用: 5]      摘要

Abstract The annual cycle in the Earth's surface temperature is extremely large-comparable in magnitude to the glacial-interglacial cycles over most of the planet. Trends in the phase and the amplitude of the annual cycle have been observed, but the causes and significance of these changes remain poorly understood-in part because we lack an understanding of the natural variability. Here we show that the phase of the annual cycle of surface temperature over extratropical land shifted towards earlier seasons by 1.7 days between 1954 and 2007; this change is highly anomalous with respect to earlier variations, which we interpret as being indicative of the natural range. Significant changes in the amplitude of the annual cycle are also observed between 1954 and 2007. These shifts in the annual cycles appear to be related, in part, to changes in the northern annular mode of climate variability, although the land phase shift is significantly larger than that predicted by trends in the northern annular mode alone. Few of the climate models presented by the Intergovernmental Panel on Climate Change reproduce the observed decrease in amplitude and none reproduce the shift towards earlier seasons.
[52] Qian C, Fu C, Wu Z.

Changes in the amplitude of the temperature annual cycle in China and their implication for climate change research

[J]. Journal of Climate, 2011, 24(20): 5 292-5 302.

DOI      URL      [本文引用: 6]     

[53] Qian C, Zhang X.

Human influences on changes in the temperature seasonality in mid-to high-latitude land areas

[J]. Journal of Climate, 2015, 28(15): 5 908-5 921.

DOI      URL      [本文引用: 1]      摘要

ABSTRACT The annual cycle is the largest variability for many climate variables outside the tropics. Whether human activities have affected the annual cycle at the regional scale is unclear. In this study, long-term changes in the amplitude of surface air temperature annual cycle in the observations are compared with those simulated by the climate models participating in phase 5 of the Coupled Model Intercomparison Project (CMIP5). Different spatial domains ranging from hemispheric to subcontinental scales in mid- to high-latitude land areas for the period 1950–2005 are considered. Both the optimal fingerprinting and a nonoptimal detection and attribution technique are used. The results show that the space–time pattern of model-simulated responses to the combined effect of anthropogenic and natural forcings is consistent with the observed changes. In particular, models capture not only the decrease in the temperature seasonality in the northern high latitudes and East Asia, but also the increase in the Mediterranean region. A human influence on the weakening in the temperature seasonality in the Northern Hemisphere is detected, particularly in the high latitudes (50–70N) where the influence of the anthropogenic forcing can be separated from that of the natural forcing.
[54] Eliseev A V, Mokhov I I.

Amplitude-phase characteristics of the annual cycle of surface air temperature in the Northern Hemisphere

[J]. Advances in Atmospheric Sciences, 2003, 20(1): 1-16.

DOI      URL      [本文引用: 1]     

[55] Dwyer J G, Biasutti M, Sobel A H.

Projected changes in the seasonal cycle of surface temperature

[J]. Journal of Climate, 2012, 25(18): 6 359-6 374.

DOI      URL      [本文引用: 1]      摘要

When forced with increasing greenhouse gases, global climate models project a delay in the phase and a reduction in the amplitude of the seasonal cycle of surface temperature, expressed as later minimum and maximum annual temperatures and greater warming in winter than in summer. Most of the global mean changes come from the high latitudes, especially over the ocean. All 24 Coupled Model Intercomparison Project phase 3 models agree on these changes and, over the twenty-first century, average a phase delay of 5 days and an amplitude decrease of 5% for the global mean ocean surface temperature. Evidence is provided that the changes are mainly driven by sea ice loss: as sea ice melts during the twenty-first century, the previously unexposed open ocean increases the effective heat capacity of the surface layer, slowing and damping the temperature response. From the tropics to the midlatitudes, changes in phase and amplitude are smaller and less spatially uniform than near the poles but are still prevalent in the models. These regions experience a small phase delay but an amplitude increase of the surface temperature cycle, a combination that is inconsistent with changes to the effective heat capacity of the system. The authors propose that changes in this region are controlled by changes in surface heat fluxes.
[56] Mann M E, Park J.

Greenhouse warming and changes in the seasonal cycle of temperature: Model versus observations

[J]. Geophysical Research Letters, 1996, 23(10): 1 111-1 114.

DOI      URL      [本文引用: 1]      摘要

Thomson [1995] argues that an enhanced green-house effect may be altering the seasonal cycle in temperature. We compare trends in the amplitude and phase of the seasonal cycle in observational temperature data in the northern hemisphere with the response of two general circulation models to increased COconcentrations. Sizeable amplitude decreases are observed in both models and observations. Significant phase delays (ie, later seasonal transitions) are found in the simulations, opposite to the phase advances isolated in the observations. The retreat of winter sea ice in high-latitude regions appears to explain the models' response to COincrease. Much of the variability in the observational data is not predicted by the models.
[57] Tesouro M, Gimeno L, Nieto R, et al.

Interannual variability of the annual cycle of temperature over Northern Africa

[J]. Studia Geophysica et Geodaetica, 2005, 49(1): 141-151.

DOI      URL      摘要

In this study, the imprints of two major atmospheric variability modes - ENSO and NAO - on the annual cycle of temperature over Northern Africa, a region sensitive to both modes, are investigated. Results from adjusting the annual cycle from daily data on a high resolution grid, indicate that both NAO and ENSO are able to influence significantly the amplitude and phase of the seasonal cycle and, consequently, that interannual trends found in amplitude and phase can be not exclusively due to greenhouse gases effects.
[58] Thomson D J.

Climate change: Shifts in season

[J]. Nature, 2009, 457(7 228): 391-392.

DOI      URL      PMID      [本文引用: 1]      摘要

It's cold in winter and hot in summer. But the latest analysis illustrates the need to put observational data at the forefront of attempts to achieve a more detailed understanding of the annual temperature cycle.
[59] Stine A R, Huybers P.

Changes in the seasonal cycle of temperature and atmospheric circulation

[J]. Journal of Climate, 2012, 25(21): 7 362-7 380.

DOI      URL      摘要

ABSTRACT The vast majority of variability in the instrumental surface temperature record is at annual frequencies. Systematic changes in the yearly Fourier component of surface temperature have been observed since the midtwentieth century, including a shift toward earlier seasonal transitions over land. Here it is shown that the variability in the amplitude and phase of the annual cycle of surface temperature in the northern extratropics is related to Northern Hemisphere atmospheric circulation as represented by the northern annular mode (NAM) and the Pacific orth America mode (PNA). The phase of the seasonal cycle is most strongly influenced by changes in spring atmospheric circulation, whereas amplitude is most strongly influenced by winter circulation. A statistical model is developed based on the NAM and PNA values in these seasons and it successfully predicts the interdecadal trends in the seasonal cycle using parameters diagnosed only at in-terannual time scales. In particular, 70% of the observed amplitude trends and 68% of the observed phase trends are predicted over land, and the residual trends are consistent with internal variability. The strong relationship between atmospheric circulation and the structure of the seasonal cycle indicates that physical explanations for changes in atmospheric circulation also extend to explaining changes in the structure of the seasonal cycle.
[60] McKinnon K A, Stine A R, Huybers P.

The spatial structure of the annual cycle in surface temperature: Amplitude, phase, and Lagrangian history

[J]. Journal of Climate, 2013, 26(20): 7 852-7 862.

DOI      URL      [本文引用: 1]     

[61] Paluš M, Novotná D, Tichavský P.

Shifts of seasons at the European mid-latitudes: Natural fluctuations correlated with the North Atlantic Oscillation

[J]. Geophysical Research Letters, 2005, 32(12): 161-179

[本文引用: 3]     

[62] Qian C, Fu C, Wu Z, et al.

The role of changes in the annual cycle in earlier onset of climatic spring in northern China

[J]. Advances in Atmospheric Sciences, 2011, 28(2): 284-296.

DOI      URL      [本文引用: 5]      摘要

Climatic changes in the onset of spring in northern China associated with changes in the annual cycle and with a recent warming trend were quantified using a recently developed adaptive data analysis tool, the Ensemble Empirical Mode Decomposition. The study was based on a homogenized daily surface air temperature (SAT) dataset for the period 1955–2003. The annual cycle here is referred to as a refined modulated annual cycle (MAC). The results show that spring at Beijing has arrived significantly earlier by about 2.98 d (10 yr) 611 , of which about 1.85 d (10 yr) 611 is due to changes in the annual cycle and 1.13 d (10 yr) 611 due to the long-term warming trend. Variations in the MAC component explain about 92.5% of the total variance in the Beijing daily SAT series and could cause as much as a 20-day shift in the onset of spring from one year to another. The onset of spring has been advancing all over northern China, but more significant in the east than in the west part of the region. These differences are somehow unexplainable by the zonal pattern of the warming trend over the whole region, but can be explained by opposite changes in the spring phase of the MAC, i.e. advancing in the east while delaying in the west. In the east of northern China, the change in the spring phase of MAC explains 40%–60% of the spring onset trend and is attributable to a weakening Asian winter monsoon. The average sea level pressure in Siberia (55°–80°N, 50°–110°E), an index of the strength of the winter monsoon, could serve as a potential short-term predictor for the onset of spring in the east of northern China.
[63] Wu Z, Schneider E K, Kirtman B P, et al.

The modulated annual cycle: An alternative reference frame for climate anomalies

[J]. Climate Dynamics, 2008, 31(7/8): 823-841.

DOI      URL      [本文引用: 3]      摘要

In climate science, an anomaly is the deviation of a quantity from its annual cycle. There are many ways to define annual cycle. Traditionally, this annual cycle is taken to be an exact repeat of itself year after year. This stationary annual cycle may not reflect well the intrinsic nonlinearity of the climate system, especially under external forcing. In this paper, we re-examine the reference frame for anomalies by re-examining the annual cycle. We propose an alternative reference frame for climate anomalies, the modulated annual cycle (MAC) that allows the annual cycle to change from year to year, for defining anomalies. In order for this alternative reference frame to be useful, we need to be able to define the instantaneous annual cycle: we therefore also introduce a new method to extract the MAC from climatic data. In the presence of a MAC, modulated in both amplitude and frequency, we can then define an alternative version of an anomaly, this time with respect to the instantaneous MAC rather than a permanent and unchanging AC. Based on this alternative definition of anomalies, we re-examine some familiar physical processes: in particular SST re-emergence and ENSO phase locking to the annual cycle. We find that the re-emergence mechanism may be alternatively interpreted as an explanation of the change of the annual cycle instead of an explanation of the interannual to interdecadal persistence of SST anomalies. We also find that the ENSO phase locking can largely be attributed to the residual annual cycle (the difference of the MAC and the corresponding traditional annual cycle) contained in the traditional anomaly, and, therefore, can be alternatively interpreted as a part of the annual cycle phase locked to the annual cycle itself. In addition to the examples of reinterpretation of physics of well known climate phenomena, we also present an example of the implications of using a MAC against which to define anomalies. We show that using MAC as a reference framework for anomaly can bypass the difficulty brought by concepts such as ecadal variability of summer (or winter) climate for understanding the low-frequency variability of the climate system. The concept of an amplitude and frequency modulated annual cycle, a method to extract it, and its implications for the interpretation of physical processes, all may contribute potentially to a more consistent and fruitful way of examining past and future climate variability and change.
[64] Qian C, Wu Z, Fu C, et al.

On multi-timescale variability of temperature in China in modulated annual cycle reference frame

[J]. Advances in Atmospheric Sciences, 2010, 27(5): 1 169-1 182.

DOI      URL      [本文引用: 2]     

[65] Qian C, Wu Z, Fu C, et al.

On changing El Niño: A view from time-varying annual cycle, interannual variability, and mean state

[J]. Journal of Climate, 2011, 24(24): 6 486-6 500.

DOI      URL      [本文引用: 4]     

[66] Deng Q, Nian D, Fu Z.

The impact of inter-annual variability of annual cycle on long-term persistence of surface air temperature in long historical records

[J]. Climate Dynamics, 2018, 50: 1 091-1 100.

DOI      URL      [本文引用: 3]     

[67] Thompson R.

Complex demodulation and the estimation of the changing continentality of Europe's climate

[J]. International Journal of Climatology, 1995, 15(2): 175-185.

DOI      URL      [本文引用: 1]      摘要

Six long records of European mean monthly temperature have been analysed for changes in continentality of climate using the time series technique of complex demodulation. Variations in amplitude and phase of the annual temperature cycle as estimated by the complex demodulation, in combination with low-pass filtering, are found to be coherent across much of Europe, from distinctly maritime to distinctly continental locations. The most significant departure of continentality found during the last 460 years is the maritime period of the 1920s when cool summers accompanied mild winters along with a retardation of the seasons throughout Europe.
[68] Vecchio A, Carbone V.

Amplitude-frequency fluctuations of the seasonal cycle, temperature anomalies, and long-range persistence of climate records

[J]. Physical Review E, 2010, 82(6): 066101.

DOI      URL      [本文引用: 3]     

[69] Capparelli V, Vecchio A, Carbone V.

Long-range persistence of temperature records induced by long-term climatic phenomena

[J]. Physical Review E, 2011, 83(4): 046103.

DOI      URL      PMID      [本文引用: 3]      摘要

Abstract The naming game (NG) describes the agreement dynamics of a population of agents that interact locally in a pairwise fashion, and in recent years statistical physics tools and techniques have greatly contributed to shed light on its rich phenomenology. Here we investigate in details the role played by the way in which the two agents update their states after an interaction. We show that slightly modifying the NG rules in terms of which agent performs the update in given circumstances (i.e., after a success) can either alter dramatically the overall dynamics or leave it qualitatively unchanged. We understand analytically the first case by casting the model in the broader framework of a generalized NG. As for the second case, on the other hand, we note that the modified rule reproducing the main features of the usual NG corresponds in fact to a simplification of it consisting in the elimination of feedback between the agents. This allows us to introduce and study a very natural broadcasting scheme on networks that can be potentially relevant for different applications, such as the design and implementation of autonomous sensor networks, as pointed out in the recent literature.
[70] Grieser J, Trömel S, Schönwiese C D.

Statistical time series decomposition into significant components and application to European temperature

[J]. Theoretical and Applied Climatology, 2002, 71(3): 171-183.

DOI      URL      [本文引用: 1]     

[71] Climate Program Office>About CPO>Glossary[EB/OL].[2017-05-11]..

URL      [本文引用: 1]     

[72] Zhao Tianbao, Qian Cheng.

Comparison between traditional anomaly and variability of temperature in China in modulated annual cycle reference frame

[J]. Climatic and Environmental Research,2010, 15(1): 34-44.

[本文引用: 1]     

[赵天保, 钱诚.

传统距平与变年循环参照系下的中国气温变率比较

[J]. 气候与环境研究, 2010, 15(1): 34-44.]

[本文引用: 1]     

[73] Box G E P, Jenkins G M, Reinsel G C, et al.

Time Series Analysis: Forecasting and Control

[M]. New Jersey: John Wiley & Sons, 2015.

[本文引用: 1]     

[74] Bonsal B R, Prowse T D.

Trends and variability in spring and autumn 0 C-isotherm dates over Canada

[J]. Climatic Change, 2003, 57(3): 341-358.

DOI      URL      [本文引用: 1]     

[75] Bingham C, Godfrey M, Tukey J.

Modern techniques of power spectrum estimation

[J]. IEEE Transactions on Audio and Electroacoustics, 1967, 15(2): 56-66.

DOI      URL      [本文引用: 3]      摘要

are the complex Fourier coefficients). Also discussed are raw and modified Fourier periodograms, bandwidth versus stability aspects, and aims and computational approaches to complex demodulation. Appendixes include a glossary, a review of complex demodulation without fast Fourier transform, and a short explanation of the fast Fourier transform.
[76] Huang N E, Shen Z, Long S R, et al.

The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis

[C]//Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences. The Royal Society, 1998, 454(1 971): 903-995.

[本文引用: 1]     

[77] Huang N E, Shen Z, Long S R.

A new view of nonlinear water waves: The Hilbert Spectrum 1

[J]. Annual Review of Fluid Mechanics, 1999, 31(1): 417-457.

DOI      URL      [本文引用: 1]      摘要

We survey the newly developed Hilbert spectral analysis method and its applications to Stokes waves, nonlinear wave evolution processes, the spectral form of the random wave field, and turbulence. Our emphasis is on the inadequacy of presently available methods in nonlinear and nonstationary data analysis. Hilbert spectral analysis is here proposed as an alternative. This new method provides not only a more precise definition of particular events in time-frequency space than wavelet analysis, but also more physically meaningful interpretations of the underlying dynamic processes.
[78] Wu Z, Huang N E.

Statistical significance test of intrinsic mode functions

[M]//Huang N E, Shen S S P, eds. Hilbert-Huang Transform and Its Applications Edited. Singapore: World Scientific, 2005: 149-169.

[本文引用: 1]     

[79] Wu Z, Huang N E.

Ensemble empirical mode decomposition: A noise-assisted data analysis method

[J]. Advances in Adaptive Data Analysis, 2009, 1(1): 1-41.

DOI      URL      [本文引用: 1]     

[80] Iatsenko D, McClintock P V E, Stefanovska A.

Nonlinear mode decomposition: A noise-robust, adaptive decomposition method

[J]. Physical Review E, 2015, 92(3): 032916.

DOI      URL      PMID      [本文引用: 3]      摘要

The signals emanating from complex systems are usually composed of a mixture of different oscillations which, for a reliable analysis, should be separated from each other and from the inevitable background of noise. Here we introduce an adaptive decomposition tool—nonlinear mode decomposition (NMD)—which decomposesa given signal into a set of physically meaningful oscillations for any wave form, simultaneously removing the noise. NMD is based on the powerful combination of time-frequency analysis techniques—which, together with the adaptive choice of their parameters, make it extremely noise robust—and surrogate data tests used to identify interdependent oscillations and to distinguish deterministic from random activity. We illustrate the application of NMD to both simulated and real signals and demonstrate its qualitative and quantitative superiority over otherapproaches, such as (ensemble) empirical mode decomposition, Karhunen-Loeve expansion, and independent component analysis. We point out that NMD is likely to be applicable and useful in many different areas of research, such as geophysics, finance, and the life sciences. The necessary MATLAB codes for running NMD arefreely available for download.
[81] Cleveland R B, Cleveland W S, Terpenning I.

STL: A seasonal-trend decomposition procedure based on loess

[J]. Journal of Official Statistics, 1990, 6(1): 3.

URL      [本文引用: 1]      摘要

ABSTRACT Abstract: STL is a filtering procedure for decomposing a time series into trend , seasonal , and remainder components. STL has a simple design that consists of a sequence of applications of the loess smoother; the simplicity allows analysis of the properties of the procedure and ...
[82] Cleveland W S, Devlin S J.

Locally weighted regression: An approach to regression analysis by local fitting

[J]. Journal of the American Statistical Association, 1988, 83(403): 596-610.

DOI      URL      [本文引用: 2]      摘要

Locally weighted regression, or loess, is a way of estimating a regression surface through a multivariate smoothing procedure, fitting a function of the independent variables locally and in a moving fashion analogous to how a moving average is computed for a time series. With local fitting we can estimate a much wider class of regression surfaces than with the usual classes of parametric functions, such as polynomials. The goal of this article is to show, through applications, how loess can be used for three purposes: data exploration, diagnostic checking of parametric models, and providing a nonparametric regression surface. Along the way, the following methodology is introduced: (a) a multivariate smoothing procedure that is an extension of univariate locally weighted regression; (b) statistical procedures that are analogous to those used in the least-squares fitting of parametric functions; (c) several graphical methods that are useful tools for understanding loess estimates and checking the assumptions on which the estimation procedure is based; and (d) the M plot, an adaptation of Mallows''s Cp procedure, which provides a graphical portrayal of the trade-off between variance and bias, and which can be used to choose the amount of smoothing.
[83] Krzyszczak J, Baranowski P, Zubik M, et al.

Temporal scale influence on multifractal properties of agro-meteorological time series

[J]. Agricultural and Forest Meteorology, 2017, 239: 223-235.

DOI      URL      [本文引用: 1]      摘要

Abstract Abstract Scale issues become very important when applying weather time series. We address problems associated with transferring meteorological data across time scales by comparing multifractal properties of hourly and daily meteorological time series. The multifractal detrended fluctuation approach revealed that temporal aggregation of agro-meteorological time series can impact on their multifractal properties. The most apparent evidence of changing the time scale on multifractal properties was found for precipitation. It was the least noticeable for the wind speed time series. The change from hourly to daily time scale had an effect on the long-range correlations and the broadness of the probability density function. The contribution of these two components to series multifractality was smaller than before data aggregation. Our results confirm the loss of unique multifractal features at daily time scale as compared to hourly time series.
[84] Eladawy A, Nadaoka K, Negm A, et al.

Characterization of the northern Red Sea’s oceanic features with remote sensing data and outputs from a global circulation model

[J]. Oceanologia, 2017, 59(3): 213-237.

DOI      URL      [本文引用: 1]     

[85] Jiang J H, Su H, Zhai C, et al.

Evaluation of cloud and water vapor simulations in CMIP5 climate models using NASA “A-Train” satellite observations

[J]. Journal of Geophysical Research, 2012, 117(D14):110 871.

[本文引用: 1]     

[86] Chen J, Del Genio A D, Carlson B E, et al.

The spatiotemporal structure of twentieth-century climate variations in observations and reanalyses. Part I: Long-term trend

[J]. Journal of Climate, 2008, 21(11): 2 611-2 633.

DOI      URL      [本文引用: 1]     

[87] Ryoo J M, Waugh D W, Gettelman A.

Variability of subtropical upper tropospheric humidity

[J]. Atmospheric Chemistry and Physics, 2008, 8(1): 1 041-1 067.

DOI      URL      [本文引用: 1]      摘要

Analysis of Atmospheric Infrared Sounder (AIRS) measurements for five northern winters shows significant longitudinal variations in subtropical upper tropospheric relative humidity (RH), not only in the climatological mean values but also in the local distributions and temporal variability. The largest climatological mean values in the northern subtropics occur over the eastern Pacific and Atlantic oceans, where there is also large day-to-day variability. In contrast, there are smaller mean values, and smaller variability that occurs at lower frequency, over the Indian and western Pacific oceans. These differences in the distribution and variability of subtropical RH are related to differences in the key transport processes in the different sectors. The large variability and intermittent high and low RH over the Eastern Pacific and Atlantic oceans, and to a smaller extent over the Indian ocean, are due to intrusions of high potential vorticity air into the subtropics. Intrusions seldom occur over the eastern Indian and western Pacific oceans, and here the subtropical RH is more closely linked to the location and strength of subtropical anticyclones. In this region there are eastward propagating features in the subtropical RH that are out of phase with the tropical RH, and are caused by modulation of the subtropical anticyclones by the Madden-Julian Oscillation.
[88] Vattay G, Harnos A.

Scaling behavior in daily air humidity fluctuations

[J]. Physical Review Letters, 1994, 73(5): 768.

DOI      URL      [本文引用: 2]      摘要

We show that the daily average air humidity fluctuations exhibit nontrivial 1/ f behavior which is different from the spectral properties of other meteorological quantities. This feature and the fractal spatial structure found in clouds make is plausible to regard air humidity fluctuations as a manifestation of self-organized criticality. We give arguments why the dynamics in air humidity can be similar to those in sandpile models of self-organized criticality.
[89] Garcia J M, Gimenez L M, Pacheco A F.

On fluctuations of the mean daily relative humidity and self-organized critical phenomena

[J]. Journal of Geophysical Research, 1997, 102(D8): 9 487-9 491.

DOI      URL      [本文引用: 2]      摘要

We analyze the spectral density of the fluctuations in the relative humidity of Zaragoza, Spain in the period 19510900091993. A similar study has already been carried out by Vattay and Harnos [1994] using data from Hungary. These authors claim that this spectrum shows a l/0408 behavior. Our results indicate a different trend, which seems to show that relative humidity at ground level is not a self-organized critical phenomenon.
[90] Chen X, Lin G X, Fu Z T.

Long-range correlations in daily relative humidity fluctuations: A new index to characterize the climate regions over China

[J]. Geophysical Research Letters, 2007, 34(7): L07804.DOI:10.1029/2006GL027755.

[本文引用: 2]     

[91] Lin G X, Chen X, Fu Z T.

Temporal-spatial diversities of long-range correlation for relative humidity over China

[J]. Physica A, 2007, 383: 585-594.

DOI      URL      [本文引用: 1]      摘要

Long-range correlations of daily relative humidity anomaly records from 191 weather stations over China during 1951鈥2000 are analyzed by means of fluctuation analysis (FA) and detrended fluctuation analysis (DFA). The information about trends in the relative humidity records can be obtained by comparing the FA curve with DFA curves. The daily relative humidity fluctuations are found to be power-law correlated and their average scaling exponent is higher than that of the temperature fluctuations, indicating that the relative humidity fluctuations take different statistical behavior from other meteorological quantities and there exists a stronger persistence in the relative humidity fluctuations. Furthermore, it is also found that these power-law scaling properties vary from station to station and show both spatial and temporal diversities, which may be explained by a proposed mechanism.

/