次网格地形坡度坡向参数化及其对区域气候模拟的影响

随着数值模式水平分辨率的提高，模式下垫面特征的描述更加仔细，针对东亚地区复杂的地形和植被特征，发展适合非均匀下垫面的地表通量参数化方案，对改进数值模式和刻画高原附近地区复杂地形动力和热力。效应的能力非常必要。本文通过计算次网格地形坡度、坡向参数及其对到达地面短波辐射通量的影响，提出次网格地形热力效应的参数化方案，改进数值模式中复杂地形区域地面热量平衡的计算。利用p—σ区域气候模式的数值试验结果表明，次网格地形热力效应参数化方案的引入，对东亚地区夏季气候尤其是降水的模拟有明显的改进。     

祁连山夏季地形云结构和云微物理过程的模拟研究(Ⅰ): 模式云物理方案和地形云结构

对三维非静力中尺度模式ARPS的云降水微物理方案进行了改进, 利用改进后的ARPS模式模拟了祁连山地区夏季的两个地形云个例, 通过对各自模拟结果的对比分析并结合实况资料研究了夏季祁连山地区地形云的发展状况、动力场特征、降水特征以及云微物理结构特征.研究结果表明, 地形云的发展受地形影响很大, 地形的抬升促进了云和降水的发展, 地形的作用也改变了地面降水特征, 使云的宏、 微观物理结构发生较大变化.     

中国西部陆面过程次网格地形参数化的改进对区域气温和降水模拟的影响研究

地表作为大气模块的下垫面,为大气模块提供边界条件,地形对于模式结果的准确性起到至关重要的作用。现有的陆面过程模式在陆面同一网格内的次网格单元采用相同的大气强迫量,没有考虑次网格地形对网格内大气强迫量的影响,这关系到模式对气象要素和陆气交换量的模拟水平。本文在陆面模式NOAH处理次网格单元的同时,将输入的大气强迫量根据其与地形高度的关系进行修订,提出新的次网格地形的参数化方案,并引入到WRF（Weather Research and Forecasting）模式中进行数值试验,通过3组数值模拟试验,与未改进的方案和细网格方案分析比较,探讨新参数化方案对WRF 模式模拟结果的影响。结果表明:地形越复杂区域,次网格地形的影响越大。本文引入的新陆面次网格地形方案对天山山脉和昆仑山脉以及青藏高原南部的地表气温的模拟有较大改善,模拟的地表气温在大范围区域内都更贴近细网格方案。虽然新陆面次网格地形方案和细网格试验都对温度的模拟结果都有改善,但新陆面次网格地形方案对降水的模拟改善甚微,而细网格试验对降水模拟却有改进,这是由于细网格试验在陆面和大气网格都进行了细化,而新陆面次网格地形方案只考虑了陆面次网格的影响。具体来说,新陆面次网格地形方案对温度的模拟结果改进是通过改变地表向上长波和地表感热实现的。而细网格试验由于同时细化了大气和陆面的空间网格,对降水和温模拟的改进是通过综合改变地表能量平衡实现的。     

祁连山夏季地形云结构和云微物理过程的模拟研究(I):模式云物理方案和地形云结构

对三维非静力中尺度模式ARPS的云降水微物理方案进行了改进,利用改进后的ARPS模式模拟了祁连山地区夏季的两个地形云个例,通过对各自模拟结果的对比分析并结合实况资料研究了夏季祁连山地区地形云的发展状况、动力场特征、降水特征以及云微物理结构特征。研究结果表明,地形云的发展受地形影响很大,地形的抬升促进了云和降水的发展,地形的作用也改变了地面降水特征,使云的宏、微观物理结构发生较大变化。     

地形湍流拖曳力参数化及在GRAPES中的应用

分辨率的限制使得不能被模式识别的地形称为次网格尺度地形,次网格尺度地形在热力和动力方面对实际大气有着不可忽略的作用,其效应只能通过参数化的形式回馈给模式。分辨率的提高使得与较小尺度地形相联系的地形湍流拖曳力凸显其重要性。数值模式中地形湍流拖曳力的参数化对完善模式物理过程和改善模式近地层预报效果具有积极意义,其方法包括有效粗糙度法和直接参数化法,而GRAPES模式中并未以任何方法考虑次网格尺度地形的影响。该文通过单柱模式比较了有效粗糙度法和直接参数化法的优劣, 发现后者在有些方面优于前者。最后,将应用于实际的一个直接参数化方案接入GRAPES中尺度模式中,进行个例模拟,并与NCEP再分析资料进行对比,结果表明:考虑地形湍流拖曳力方案对模式预报具有改进作用,尤其对局地低层风场具有积极影响。     

定量研究雅安地形坡向坡度对降水分布的影响

利用四川省雅安市30 m分辨率基础高程数据,提取栅格的坡向和坡度参数,将雅安307个区域自动站在2017年汛期(6—9月)共50次的降水天气个例,分为16次大尺度降水和34次中小尺度降水,使用对应时次的欧洲中心细网格0.25°×0.25°再分析风场资料,根据不同的站点地形高度将风场合成平均风场,和各站点地形的坡向和坡度计算出其动力抬升作用,同时使用当天日照和天文太阳辐射值来计算地形的热力抬升作用,与对应降水过程的降水分布进行多元线性回归,根据回归的标准系数的大小确定各自变量对降水分布的影响,得出以下结论:(1)中小尺度降水中,地形的热力抬升作用对降水分布的影响作用最大,其次是海拔高度,地形的动力抬升作用在三者中对降水分布的影响最小;(2)在大尺度降水中刚好相反,地形的动力抬升对降水的分布影响作用最大,海拔次之,热力抬升作用在三者中影响作用最小;(3)日降水量最大值的站点海拔高度基本位于1 000 m左右,与抬升凝结高度对应较好;(4)从长期的统计来看,地形的动力作用和地表的植被情况对降水分布的影响最大。在实际预报工作用,根据不同的降水类型,关注不同的动力和热力作用对于判断降水分布大值区的位置有较好的参考作用。     

祁连山夏季地形云和降水宏微观结构的数值模拟

利用ARPS(Advanced Regional Predictional System)中尺度数值模式,对2007年7月19日低涡天气背景下发生在祁连山区的一次比较典型的地形云降水过程中云和降水的宏微观结构特征进行了深入的模拟研究和分析。结果表明,ARPS模式能够较好地模拟出地面降水分布及其发展演变特征;祁连山北坡陡峭地形的抬升作用是祁连山云系降水的主要动力机制;祁连山地形作用下云和降水的微物理结构随云的不同发展阶段呈现出不同的特征。     

台风“桑美”的数值模拟和地形敏感性试验

用WRF模式对0608号台风“桑美”进行了数值模拟研究,较为成功地模拟出了台风路径和降水,但模拟的台风中心气压远高于实况。为研究“桑美”登陆期间地形的抬升作用对其降水及结构的影响,通过改变特定区域内的地形高度设计了一组敏感性试验。结果表明,台风登陆过程中地形抬升作用对台风降雨量有显著的增幅作用;台风中心位势涡度、气流垂直上升速度、水平水汽通量散度明显增大;地形抬升机制在台风登陆时刻达到最强。     

采用次网格地形方案对白鹤滩水电站坝区一次大风天气的数值模拟

基于中尺度气象模式,采用次网格地形方案模拟了2019年3月19~20日四川、云南交界处白鹤滩水电站的一次大风过程,对10 m风速、风向和2 m温度空间分布及日变化的模拟结果进行检验评估,并结合地形分析了坝区易产生大风的原因。结果表明:此次西南大风天气是由高低空一致的西南气流配合地面“东高西低”的环流形势共同造成的。午后热低压发展强盛,地面风速增大,坝区出现8级大风。采用次网格地形方案后,风速的平均绝对误差和均方根误差最高可分别减小17%和14%,同时该方案有效地缓解了模式对白天风速的低估和夜间风速的高估,在大风和小风时段均对风速有较好的模拟能力,从而能更好地刻画风速的日变化特征。综合来看,次网格地形方案能显著改善风场的模拟效果,其中Jiménez方案更适用于坝区大风的模拟,但次网格地形方案对温度模拟没有改善作用。白鹤滩水电站的大风受局地地形影响极大,南北向狭长河谷地形产生的狭管效应使气流增速显著,坝区主体高度区河谷风的放大系数超过3.0,使得白鹤滩水电站极易出现灾害性大风天气。      

贵州南部山区一次大暴雨的数值模拟

利用WRF模式,对2008年5月27日夜间贵州南部山区大暴雨天气进行了数值模拟。复杂地形条件下,数值模拟有一定难度。本次数值模拟能够反映此次大暴雨的降水区域、降水特征和降水强度。在此基础上,利用高分辨率的数值结果,对该暴雨过程的中尺度气旋及地形对降水的作用进行了分析,结果表明：中β尺度气旋是造成此次局地暴雨的直接系统;气旋在结构上表现出中低层强烈的辐合和正涡度柱的耦合,出现了伸展到对流层顶的深厚强上升运动;造成大暴雨的辐合线生成在贵州西部地形梯度较大地方,辐合线上发展起来的中尺度气旋在850 hPa和地面上表现为生成在一个开口向南的喇叭口地形上,由于喇叭口地形对气流的辐合和抬升作用,水汽辐合和上升运动在该地达到最强,无疑对降水的落区和强度有重要影响。     

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.     

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.     

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.     

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.