基于分位数映射法的黑河上游气候模式降水误差订正

区域气候模式降水弥补了高寒山区气象站点稀少的缺陷,是水文模拟的重要驱动变量。然而,高寒山区模式输出降水的总量和频率都存在较大不确定性。因此,改进了用于降水频率纠正的分位数映射法(Quantile Mapping,QM),对中尺度数值预报模式(Weather Research and Forecasting model,WRF)模拟的黑河上游日降水输出数据进行误差订正。选取第95分位和第98分位降水量为阈值,选择2004-2009年为建模时段,2010-2013年为验证时段,使用分段拟合的方法建立传递函数,侧重于对极端降水进行单独订正。研究结果表明:该方法不仅对降水空间分布有明显的改善,对极端降水也有很好的订正效果。订正前模式模拟日降水与台站之间的均方根误差为3.41 mm·d^-1,绝对偏差为115.67 mm·y^-1,订正后均方根误差减少为3.11 mm·d^-1,绝对偏差有明显改善,为60.3 mm·y^-1。订正后流域内年降水空间分布更加合理,年降水量也更接近于观测降水插值结果,其空间相关系数由0.74改善为0.94。春、夏季订正效果优于秋、冬季,其中夏季订正效果较为明显,订正前降水偏差百分比在-0.1~0.1以内的区域面积仅占流域总面积的28%,而订正后占比增加至66%。同时,该方法对极端降水有较好的订正效果,减小了日降水强度(SDII)和极强降水量(R99p)的模拟偏差,订正后的第95分位模拟降水与观测降水插值的相关系数由0.15提高到0.48。本研究为站点稀少的黑河上游提供了一种更有效的误差订正方案,有利于为寒区水文研究获取更精确的降水数据。     

分位数映射法在RegCM4中国气温模拟订正中的应用

将一种分位数映射法RQUANT,应用到一个区域气候模式（RegCM4）所模拟中国气温的误差订正中。从气候平均态、年际变率、极端气候及农业气候等多方面,评估了该方法对日平均气温、日最高气温和日最低气温模拟的订正效果。结果表明,该订正方法对模式模拟的日平均、日最高和最低气温气候平均态的订正效果都非常明显,中国大部分地区的订正结果与观测的偏差在±0.5℃之间。在降低极端气温指数和农业气候相关指数的模拟误差方面也有显著的效果,但对气温年际变率的订正效果有限。结合以往对降水订正的评估分析,该方法对模式模拟结果有较好的订正效果,可以应用于区域气候模式的气候变化模拟预估中,为气候变化及相关影响评估研究提供更适用和可靠的数据。     

QM和QDM方法对中国极端气候的高分辨率气候变化模拟的误差订正对比

采用分位数映射（Quantile Mapping, QM）和delta分位数映射（Quantile Delta Mapping, QDM）两种误差订正方法对区域气候模式RegCM4在中国区域内模拟的逐日气温和降水数据进行订正。模式数据是5种不同全球气候模式驱动下的区域模式气候变化模拟结果。计算订正前后的极端气候指数进行对比分析,包括日最高气温极大值（TXx）、日最低气温极小值（TNn）、连续干旱日数（CDD）和最大日降水量（RX1day）。结果表明,5组模拟结果和其集合平均（ensR）都显示气温指数的模拟效果高于降水指数,其中对TXx模拟最好,对CDD的模拟最差;经过订正后,针对不同模式的两种订正结果都能够有效地减小模式与观测的偏差并提高了空间相关系数,且两种方法的订正效果无明显差别。对RCP4.5情景下未来变化的分析中,QM在一定程度上改变了模式模拟的未来变化幅度和空间分布特征,QDM则能够有效地保留所有极端指数的气候变化信号。从全国平均来看,除CDD外,所有指数未来都呈现增加趋势,且QDM订正结果与订正前模式模拟的变化趋势更为接近。建议在气候变化模拟的误差订正中采用QDM方法。     

气候模拟日降水量的统计误差订正分析——以上海为例

对区域气候模式系统PRECIS在SRES A1B情景下模拟的上海日降水输出按季节进行了统计误差订正。该方法首先对降水日数进行比率订正,以消除模式产生的微小值降水。然后利用Γ分布拟合日降水量的累计概率分布,采用整体和分段拟合两种方法构建传递函数TF（Transfer Function）进行订正。选取1962年12月—1992年11月作为控制时段,构建TF并将其应用于验证时段（1992年12月—2002年11月）。该订正方案消除了模式产生的微小值降水,解决了模拟的小降水值偏多的问题,频率误差保持在1%以下,分段拟合订正相比整体拟合订正具有更强的对极端降水的订正能力;对冬、春季的订正效果比夏、秋季更显著。该方案不仅有效消除了平均值的漂移,而且显著订正了变率,同时提高了极端降水事件的再现能力,是一种相对完善的订正方案。      

CWRF模式极端降水模拟误差订正

基于1980—2015年6—8月CWRF模式（Climate-Weather Research and Forecasting model）14种方案的模拟结果和全国逐日降水观测资料,对比了Q-lin,Q-tri,RQ-lin,RQ-tri,SSP-lin和CDFt 6种误差订正方法对CWRF模式控制化方案（C1）模拟中国东部夏季日极端降水的订正效果,以CWRF模式14种方案日极端降水的模拟效果排名为基础,对比了模拟效果较好的4种方案集合、模拟较差的4种方案集合以及14种方案集合的订正效果,选出相对较好的订正方案进一步评估其成员集合后订正和成员分别订正后再集合的订正效果,结果表明:采用6种误差订正方法均可明显减少日极端降水模拟误差,其中RQ-lin方法订正效果最佳。CWRF模式对中国东部的极端降水指数均表现出较好的模拟能力,不同参数化集合方案得到14种方案成员先订正再集合与观测日极端降水平均值最为接近,研究结果对于改进模拟结果、提高其预测能力有重要应用价值。     

1951—2015年通辽市科尔沁区日气温特征分析

文章利用通辽国家基本气象站1951—2015年逐日最高气温、最低温度和日平均气温,计算了3个要素序列的算数平均值、中位数、众数及样本标准差等基本统计量,分析了每个气温等级出现的频次和不同气温段的气温日较差出现日数、日最高气温〉30℃日数和日最低气温≤-20℃日数等。结果表明：（1）该地极端最高气温达39.4℃,极端最低气温为-33.9℃,最高气温出现在7月12日,最低气温出现在1月17日。（2）全年26.1~28.0℃的最高气温出现频次和16.1~18.0℃的日最低气温出现频次最高。（3）一年中,日最高气温在0℃以上的天数占全年总日数的75%;日最低气温在0℃以下的日数占全年总日数的46%;日较差〉10℃的日数占全年总日数的69%。（4）1981—2015年平均气温比1951—1981年明显升高,表现为低温日数明显减少,而近10a不仅冬季低温日数减少而且夏季高温日数明显增多。     

基于Delta分位数映射法的青藏高原中东部IMERG卫星降水误差订正

可靠的降水资料对理解青藏高原水量平衡和水循环过程尤为重要。IMERG(Integrated Multi-satellitE Retrievals for Global Precipitation Measurement)是新一代卫星降水产品,具有更广的覆盖范围与更高的时空分辨率,但在高原复杂地形条件下仍然存在较大的不确定性。鉴于此,本研究应用Delta分位数映射法(Quantile Delta Mapping, QDM),对IMERG日降水数据进行偏差订正,使用2001-2010年的中国区域地面气象要素数据集(China Meteorological Forcing Dataset, CMFD)降水数据和IMERG日降水产品分季节建立传递函数,对2011-2014年的IMERG逐日降水进行订正。研究结果表明：(1)Delta分位数映射法能够有效订正IMERG的降水频率、数值和空间分布,对极端降水和负偏差较大区域的订正效果更为明显。订正后的IMERG降水概率分布更加接近观测概率分布,降水偏差也更符合正态分布,改进了对全年和季节降水空间分布的刻画,提高了月降水的精度。(2)订正后日降水量均...     

中国夏季极端气温与降水事件日数随平均气温变化的定量分析

极端高温、极端低温、极端降水事件日数对全球平均气温变化都有较强响应。本文利用全国1960-2005年549个站点逐日均一化温度资料以及559个站点逐日降水数据,定量分析了中国夏季极端气温与降水事件日数随气温变化的特征。结果表明:全球平均气温升高1℃,全国平均的极端高温、低温和降水日数的变化量分别为5.69,-5.3和0.69天;区域尺度上,全球平均气温升高1℃,东南沿岸和四川地区极端高温日数可增加8～10天,东北地区极端低温日数减少10天左右,西北地区极端降水日数可增加4～6天。基于以上结果,利用模式对未来全球平均气温的预估表明:2006-2099年我国平均的极端高温、低温、降水日数变化量在RCP8.5情景下最大,分别增加23天、减少22天和增加3天;区域尺度上,RCP8.5情景下2006-2099年我国东南以及西南地区极端高温日数可增加42天左右,我国北方大部分地区极端低温日数减少33天左右,而西北地区极端降水日数可增加16天左右。该研究结果表明一系列应对气候变化的措施势在必行。     

不同升温阈值下中国地区极端气候事件变化预估

本文基于耦合模式比较计划第五阶段(CMIP5)的18个全球气候模式的模拟结果,预估了全球平均气温在不同典型浓度路径(RCPs)下达到2℃、3℃和4℃阈值时,中国地区气温和降水的变化,并采用了具有稳定统计意义的27个极端气候指标定量评估了全球平均气温达到不同阈值时,中国地区极端气候事件的可能变化。结果表明,未来我国平均气温增幅将高于全球平均增暖,极端暖事件(如暖夜、暖昼、热带夜)明显增多,达到4℃阈值时,暖夜指数相比参考时段增加约49.9%。极端冷事件(如冷夜、冷昼、霜冻)减少。随全球气温升高,中国北方平均降水增多。在不同升温阈值下,中国地区降水的极端性都体现出增强的趋势,强降水事件发生频率(如中雨日数、大雨日数)和强度(如五日最大降水量、极端强降水量)都明显增加。随升温阈值的升高,这些变化幅度更大,在 RCP8.5 情景下全球升温 3℃和4℃时,中国平均五日最大降水分别增加 12.5mm和17.0mm。我国西南地区极端降水强度的增幅高于其他地区。     

概率调整法在气候模式模拟降水量订正中的应用

应用概率调整法订正区域气候模式系统PRECIS在SRES A1B情景下模拟的各季节全国日降水量。以第95百分位降水量为阈值,利用Γ分布分段拟合1962年12月—1972年11月的模拟值,构建传递函数,得到1991年12月—2001年11月的订正值。结果表明:全国平均日降水量空间分布的模拟改善明显,偏差百分率高于100%的格点比例从23.5%降低到1.0%;对各地区平均降水月循环的模拟结果改善,冷季降水较暖季更接近观测,提高拟合优度是改进订正方法的关键;多数地区连续干日数、连续5 d最大降水量及极端降水贡献率的空间强度、概率分布与空间相关性的订正效果显著。总体来说,该方法对模拟中国区域降水的平均态与极端降水均有明显改善,有助于气候评估工作的展开。     

CHARACTERISTICS AND TRENDS OF CLIMATIC EXTREMES IN CHINA DURING 1959–2014

The spatial and temporal variations of daily maximum temperature(Tmax), daily minimum temperature(Tmin), daily maximum precipitation(Pmax) and daily maximum wind speed(WSmax) were examined in China using Mann-Kendall test and linear regression method. The results indicated that for China as a whole, Tmax, Tmin and Pmax had significant increasing trends at rates of 0.15℃ per decade, 0.45℃ per decade and 0.58 mm per decade,respectively, while WSmax had decreased significantly at 1.18 m·s~(-1) per decade during 1959—2014. In all regions of China, Tmin increased and WSmax decreased significantly. Spatially, Tmax increased significantly at most of the stations in South China(SC), northwestern North China(NC), northeastern Northeast China(NEC), eastern Northwest China(NWC) and eastern Southwest China(SWC), and the increasing trends were significant in NC, SC, NWC and SWC on the regional average. Tmin increased significantly at most of the stations in China, with notable increase in NEC, northern and southeastern NC and northwestern and eastern NWC. Pmax showed no significant trend at most of the stations in China, and on the regional average it decreased significantly in NC but increased in SC, NWC and the mid-lower Yangtze River valley(YR). WSmax decreased significantly at the vast majority of stations in China, with remarkable decrease in northern NC, northern and central YR, central and southern SC and in parts of central NEC and western NWC. With global climate change and rapidly economic development, China has become more vulnerable to climatic extremes and meteorological disasters, so more strategies of mitigation and/or adaptation of climatic extremes,such as environmentally-friendly and low-cost energy production systems and the enhancement of engineering defense measures are necessary for government and social publics.     

ANALYSIS OF “BIMODAL PATTERN” STORM ACTIVITY CHARACTERISTICS OF THE BAY OF BENGAL

Storms that occur at the Bay of Bengal (BoB) are of a bimodal pattern, which is different from that of the other sea areas. By using the NCEP, SST and JTWC data, the causes of the bimodal pattern storm activity of the BoB are diagnosed and analyzed in this paper. The result shows that the seasonal variation of general atmosphere circulation in East Asia has a regulating and controlling impact on the BoB storm activity, and the “bimodal period” of the storm activity corresponds exactly to the seasonal conversion period of atmospheric circulation. The minor wind speed of shear spring and autumn contributed to the storm, which was a crucial factor for the generation and occurrence of the “bimodal pattern” storm activity in the BoB. The analysis on sea surface temperature (SST) shows that the SSTs of all the year around in the BoB area meet the conditions required for the generation of tropical cyclones (TCs). However, the SSTs in the central area of the bay are higher than that of the surrounding areas in spring and autumn, which facilitates the occurrence of a “two-peak” storm activity pattern. The genesis potential index (GPI) quantifies and reflects the environmental conditions for the generation of the BoB storms. For GPI, the intense low-level vortex disturbance in the troposphere and high-humidity atmosphere are the sufficient conditions for storms, while large maximum wind velocity of the ground vortex radius and small vertical wind shear are the necessary conditions of storms.     

LOW-FREQUENCY CIRCULATION AND EXTENDED-RANGE FORECAST IN CONNECTION WITH AUTUMN EXTREME HEAVY RAINFALL PROCESS IN HAINAN

Observed daily precipitation data from the National Meteorological Observatory in Hainan province and daily data from the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis-2 dataset from 1981 to 2014 are used to analyze the relationship between Hainan extreme heavy rainfall processes in autumn (referred to as EHRPs) and 10–30 d low-frequency circulation. Based on the key low-frequency signals and the NCEP Climate Forecast System Version 2 (CFSv2) model forecasting products, a dynamical-statistical method is established for the extended-range forecast of EHRPs. The results suggest that EHRPs have a close relationship with the 10–30 d low-frequency oscillation of 850 hPa zonal wind over Hainan Island and to its north, and that they basically occur during the trough phase of the low-frequency oscillation of zonal wind. The latitudinal propagation of the low-frequency wave train in the middle-high latitudes and the meridional propagation of the low-frequency wave train along the coast of East Asia contribute to the ‘north high (cold), south low (warm)’ pattern near Hainan Island, which results in the zonal wind over Hainan Island and to its north reaching its trough, consequently leading to EHRPs. Considering the link between low-frequency circulation and EHRPs, a low-frequency wave train index (LWTI) is defined and adopted to forecast EHRPs by using NCEP CFSv2 forecasting products. EHRPs are predicted to occur during peak phases of LWTI with value larger than 1 for three or more consecutive forecast days. Hindcast experiments for EHRPs in 2015–2016 indicate that EHRPs can be predicted 8–24 d in advance, with an average period of validity of 16.7 d.     

AN ATTRIBUTION ANALYSIS OF CHANGES IN POTENTIAL EVAPOTRANSPIRATION IN THE BEIJING–TIANJIN–HEBEI REGION UNDER CLIMATE CHANGE

Based on the measurements obtained at 64 national meteorological stations in the Beijing–Tianjin–Hebei (BTH) region between 1970 and 2013, the potential evapotranspiration (ET0) in this region was estimated using the Penman–Monteith equation and its sensitivity to maximum temperature (Tmax), minimum temperature (Tmin), wind speed (Vw), net radiation (Rn) and water vapor pressure (Pwv) was analyzed, respectively. The results are shown as follows. (1) The climatic elements in the BTH region underwent significant changes in the study period. Vw and Rn decreased significantly, whereas Tmin, Tmax and Pwv increased considerably. (2) In the BTH region, ET0 also exhibited a significant decreasing trend, and the sensitivity of ET0 to the climatic elements exhibited seasonal characteristics. Of all the climatic elements, ET0 was most sensitive to Pwv in the fall and winter and Rn in the spring and summer. On the annual scale, ET0 was most sensitive to Pwv, followed by Rn, Vw, Tmax and Tmin. In addition, the sensitivity coefficient of ET0 with respect to Pwv had a negative value for all the areas, indicating that increases in Pwv can prevent ET0 from increasing. (3) The sensitivity of ET0 to Tmin and Tmax was significantly lower than its sensitivity to other climatic elements. However, increases in temperature can lead to changes in Pwv and Rn. The temperature should be considered the key intrinsic climatic element that has caused the "evaporation paradox" phenomenon in the BTH region.     

Could the Recent Taal Volcano Eruption Trigger an El Ni?o and Lead to Eurasian Warming?

正The Taal Volcano in Luzon is one of the most active and dangerous volcanoes of the Philippines. A recent eruption occurred on 12 January 2020(Fig. 1a), and this volcano is still active with the occurrence of volcanic earthquakes. The eruption has become a deep concern worldwide, not only for its damage on local society, but also for potential hazardous consequences on the Earth's climate and environment.     

VARIABILITY OF INTER-ANNUAL RELATIONSHIP BETWEEN INDIAN OCEAN DIPOLE AND HUAXI REGION’S AUTUMN PRECIPITATION

The moving-window correlation analysis was applied to investigate the relationship between autumn Indian Ocean Dipole (IOD) events and the synchronous autumn precipitation in Huaxi region, based on the daily precipitation, sea surface temperature (SST) and atmospheric circulation data from 1960 to 2012. The correlation curves of IOD and the early modulation of Huaxi region’s autumn precipitation indicated a mutational site appeared in the 1970s. During 1960 to 1979, when the IOD was in positive phase in autumn, the circulations changed from a “W” shape to an ”M” shape at 500 hPa in Asia middle-high latitude region. Cold flux got into the Sichuan province with Northwest flow, the positive anomaly of the water vapor flux transported from Western Pacific to Huaxi region strengthened, caused precipitation increase in east Huaxi region. During 1980 to 1999, when the IOD in autumn was positive phase, the atmospheric circulation presented a “W” shape at 500 hPa, the positive anomaly of the water vapor flux transported from Bay of Bengal to Huaxi region strengthened, caused precipitation ascend in west Huaxi region. In summary, the Indian Ocean changed from cold phase to warm phase since the 1970s, caused the instability of the inter-annual relationship between the IOD and the autumn rainfall in Huaxi region.     

Revisiting the Concentration Observations and Source Apportionment of Atmospheric Ammonia

正While China’s Air Pollution Prevention and Control Action Plan on particulate matter since 2013 has reduced sulfate significantly, aerosol ammonium nitrate remains high in East China. As the high nitrate abundances are strongly linked with ammonia, reducing ammonia emissions is becoming increasingly important to improve the air quality of China. Although satellite data provide evidence of substantial increases in atmospheric ammonia concentrations over major agricultural regions, long-term surface observation of ammonia concentrations are sparse. In addition, there is still no consensus on     

全球贸易隐含碳变化及其影响分析

基于最新的GTAP8 (Global Trade Analysis Project）数据库,使用投入产出法,分析了2004年到2007年全球贸易变化下南北集团贸易隐含碳变化及对全球碳排放的影响。结果显示,随着发展中国家进出口规模扩张,全球贸易隐含碳流向的重心逐渐向发展中国家转移。2004年到2007年,发达国家高端设备制造业和服务业出口以及发展中国家资源、能源密集型行业及中低端制造业出口的趋势加强,该过程的生产转移导致全球碳排放增长4.15亿t,占研究时段全球贸易隐含碳增量的63%。未来发展中国家的出口隐含碳比重还将进一步提高。贸易变化带来的南北集团隐含碳流动变化对全球应对气候变化行动的影响日益突出,发达国家对此负有重要责任。     

COMPARATIVE ANALYSIS OF VARIOUS DATA OF AIR–SEA TEMPERATURE DIFFERENCE AND ITS VARIATION ACROSS SOUTH CHINA SEA IN THE PAST 35 YEARS

Using the International Comprehensive Ocean-Atmosphere Data Set(ICOADS) and ERA-Interim data, spatial distributions of air-sea temperature difference(ASTD) in the South China Sea(SCS) for the past 35 years are compared,and variations of spatial and temporal distributions of ASTD in this region are addressed using empirical orthogonal function decomposition and wavelet analysis methods. The results indicate that both ICOADS and ERA-Interim data can reflect actual distribution characteristics of ASTD in the SCS, but values of ASTD from the ERA-Interim data are smaller than those of the ICOADS data in the same region. In addition, the ASTD characteristics from the ERA-Interim data are not obvious inshore. A seesaw-type, north-south distribution of ASTD is dominant in the SCS; i.e., a positive peak in the south is associated with a negative peak in the north in November, and a negative peak in the south is accompanied by a positive peak in the north during April and May. Interannual ASTD variations in summer or autumn are decreasing. There is a seesaw-type distribution of ASTD between Beibu Bay and most of the SCS in summer, and the center of large values is in the Nansha Islands area in autumn. The ASTD in the SCS has a strong quasi-3a oscillation period in all seasons, and a quasi-11 a period in winter and spring. The ASTD is positively correlated with the Nio3.4 index in summer and autumn but negatively correlated in spring and winter.     

Retrieval of outgoing longwave radiation from COMS narrowband infrared imagery

Hourly outgoing longwave radiation(OLR) from the geostationary satellite Communication Oceanography Meteorological Satellite(COMS) has been retrieved since June 2010. The COMS OLR retrieval algorithms are based on regression analyses of radiative transfer simulations for spectral functions of COMS infrared channels. This study documents the accuracies of OLRs for future climate applications by making an intercomparison of four OLRs from one single-channel algorithm(OLR12.0using the 12.0 μm channel) and three multiple-channel algorithms(OLR10.8+12.0using the 10.8 and 12.0 μm channels; OLR6.7+10.8using the 6.7 and 10.8 μm channels; and OLR All using the 6.7, 10.8, and 12.0 μm channels). The COMS OLRs from these algorithms were validated with direct measurements of OLR from a broadband radiometer of the Clouds and Earth's Radiant Energy System(CERES) over the full COMS field of view [roughly(50°S–50°N, 70°–170°E)] during April 2011.Validation results show that the root-mean-square errors of COMS OLRs are 5–7 W m-2, which indicates good agreement with CERES OLR over the vast domain. OLR6.7+10.8and OLR All have much smaller errors(～ 6 W m-2) than OLR12.0and OLR10.8+12.0(～ 8 W m-2). Moreover, the small errors of OLR6.7+10.8and OLR All are systematic and can be readily reduced through additional mean bias correction and/or radiance calibration. These results indicate a noteworthy role of the6.7 μm water vapor absorption channel in improving the accuracy of the OLRs. The dependence of the accuracy of COMS OLRs on various surface, atmospheric, and observational conditions is also discussed.