北京城市大气混合层与气溶胶垂直分布观测研究

2002年3～10月在北京大学利用微脉冲激光雷达(MPL)观测了气溶胶时空变化。提出一种反演混合层高度的方法,这种方法减小了仪器订正的误差,反演的混合层高度与探空测量结果有很好的一致性。利用该方法计算了观测期间晴天无云天气条件下的混合层高度,分析了混合层高度及其增长率的日变化、季节变化,初步研究了混合层高度和近地面气溶胶分布的相互关系,分析了表征大气扩散能力的通风系数的日变化。结果表明,利用MPL监测城市混合层是可行的和优越的。     

激光雷达和探空反演的混合层高度对比分析

本文利用2016年12月至2017年11月期间晴朗少云天气下的成都微脉冲激光雷达观测数据反演的混合层高度,与温江探空资料确定的混合层高度进行了对比和误差分析,结果表明:基于探空资料和激光雷达数据反演的混合层高度具有较好的一致性,两者相关系数达0.75,激光雷达反演的混合层高度略低于基于探空资料确定的值,在混合层高度为1000～2000m时,两种方法计算所得的值偏差幅度最小,约为20%;在1000m以内和2000m以上,偏差幅度略有增大,为26%;两种方法反演的混合层高度变化趋势较为一致,均呈现出12月、1月较低,4月、5月较高的特点;混合层高度具有明显的日变化特征:上午混合层高度迅速增高,午后增长速度减慢并发展到最大高度,日落后迅速降低;混合层内相对湿度的增加、残留层的存在是导致激光雷达反演混合层高度时产生较大误差的原因之一。     

结合激光雷达和探空资料研究青藏高原地区混合层高度特征

利用那曲地区的微脉冲激光雷达探测资料,采用梯度法获取了那曲地区白天混合层高度随时间的演变信息及混合层特征参数,结果表明混合层在上午发展缓慢,中午以后发展迅速,14:00（北京时）前后达到稳定;强烈的对流热泡活动导致混合层高度起伏较大,参数化反演得到的卷夹层厚度达到0.4~0.5 km,卷夹比在0.2左右。利用探空资料结合日最大位温资料,采用气块法得到了高原地区7个站点的每日最大混合层高度数据集。通过对由激光雷达和探空资料得到的最大混合层高度结果进行对比,发现二者有很好的一致性（相关系数0.85,均值偏差0.11 km,均方根误差0.30 km,并通过0.05显著性水平的t检验）。最大混合层高度在7个站点均有明显的逐日变化特征。从年均值看,最大混合层高度与海拔高度之间没有明显相关关系。从季节均值看,格尔木与都兰站表现出明显的春高冬低的分布特征,而其它五个站点则表现为春高夏低,表明高原地区的盆地地形和山地地形对混合层高度有显著而不同的影响。通过定义热力稳定度和确定特征气压层高度,利用热力稳定度与最大混合层高度之间良好的线性关系,获取了一种简便地获取最大混合层高度的统计方法。     

银川市大气混合层高度变化特征及其与空气污染的关系

基于银川高空站2008～2017年的L波段秒级数据和地面观测数据,利用干绝热法计算银川2008～2017年逐时大气混合层高度,分析其变化特征,同时利用银川6种污染物的质量浓度和AQI指数,分析大气混合层高度与空气污染物的关系。结果表明:银川市的大气混合层高度（MLH）大部时间在600 m以下,占比为68%;银川MLH具有明显的单峰型日变化特征,07:00（北京时间,下同）最低,16:00最高;各月MLH值在282～936 m,4～6月MLH值最高,12月MLH值最低;季节变化方面,春季最高,夏季次之,冬季最低;年变化方面,2012年MLH平均值最高为621 m,2015年最低为566 m;银川市6种污染物除O₃外,其余5种污染物的质量浓度与MLH值都为负相关,O₃与MLH值的相关性最好;四季中,冬季污染物浓度与MLH值的相关性最好,夏季最差,秋季好于春季;银川MLH值与AQI指数呈负相关。     

基于激光雷达和微波辐射计观测确定混合层高度方法的比较

利用苏州地区2010年1月4,7,16日和2月4日4天的激光雷达及微波辐射计观测资料,比较了不同遥感手段探测晴天大气混合层高度的差异,发现试验期间该地区的混合层高度在300~1500 m之间。利用梯度法、标准偏差法、小波法从激光雷达数据中提取混合层高度并进行了对比,结果表明三种方法都能较好地反演混合层高度并且一致性较好,三者差异主要存在于大气边界层的发展和消亡阶段;梯度法和小波法结果无明显差异,而标准偏差法结果稍高于其他方法。在此基础上,利用微波辐射计探测的大气温度,使用温度梯度法估算大气混合层高度,并与激光雷达探测结果进行了比较,结果表明,大多数情况下激光雷达探测结果高于微波辐射计观测结果;两种遥感手段有较好的相关性,相关系数为0.76。激光雷达同微波辐射计结果存在差异,尤其是在边界层的发展和消散阶段,这是由两种遥感手段探测原理不同造成的。     

大气红外探测器(AIRS)温、湿廓线反演产品及边界层高度在黄土高原的验证

利用黄土高原地区加强观测期探空资料和大气红外探测器(Atmospheric Infrared Sounder,AIRS)反演的温度、相对湿度廓线资料,评估了AIRS反演产品在黄土高原的适用性,以及利用AIRS温度廓线计算边界层高度的可行性。结果表明,各种边界层高度计算方法相关性显著,平均高度差异一般不超过200 m;Richardson数临界值的选取对确定边界层高度的影响不大。AIRS反演的大气温度和相对湿度均能很好地反映环境温湿的变化;温度平均偏差在±1 K以内,均方根误差一般不超过2 K;AIRS反演的地面气温误差相对较大,平均偏差和均方根误差分别为-1.68 K和3.32 K,并且会影响边界层高度的确定;AIRS相对湿度平均偏差在±10%以内,均方根误差不超过20%。利用Parcel法估计的边界层高度结果显示,根据AIRS反演的温度廓线确定的AIRS边界层高度低于利用探空观测资料确定的边界层高度,但能较好地再现边界层高度的变化,在没有探空数据时,可以用AIRS反演的温度廓线确定边界层高度。     

COSMIC掩星资料反演青藏高原大气廓线与探空观测的对比分析

利用青藏高原地区COSMIC掩星资料反演的大气湿廓线Wet Prf数据和8个站点的探空数据,分析了COSMIC反演大气廓线和可降水量与探空观测的偏差,并考查了偏差随高度的变化特征。结果显示:(1)COSMIC反演的温度、压强和水汽压廓线与探空观测具有很好的正相关;与探空观测相比,COSMIC的温度、压强和水汽压的偏差为-0.2℃、1.7 h Pa和0 h Pa,均方差为1.8℃、1.6 h Pa和0.4 h Pa;COSMIC反演大气廓线与探空观测的偏差基本上在大气低层较大,然后随高度增加而减小。(2)COSMIC反演的可降水量与探空观测正相关较好;COSMIC反演的可降水量低于探空观测,两者的偏差为-5.0 mm,均方差为5.7 mm;两者的负偏差在大气低层最明显。(3)探空观测在近地层的不稳定性和COSMIC反演方法中背景模式在青藏高原地区描述大气状态的能力有限,是造成COSMIC反演大气廓线和探空观测的偏差在近地层较大的主要原因;COSMIC观测的折射率偏小导致其反演的可降水量偏低。     

浙北地区一次冷空气过程中的灰霾天气过程分析

利用常规气象观测资料、中尺度自动站资料、探空站资料和空气污染资料,对2014年1月3-4日浙北地区的一次大范围灰霾天气过程进行了综合分析,结果表明:冷空气南下和中层高度的干燥暖舌影响,是这次灰霾天气过程的天气背景;冷空气受到700 h Pa干暖舌的抑制作用,侵入700 h Pa以上高度的速度缓慢,导致地面虽然有较大风速,但在中层形成的逆温结构,使得这次灰霾天气过程能够维持;发生霾时的混合层高度比轻雾高,大气混合层高度的变化对霾的发展变化有较好的指示作用,可为霾的预报提供参考依据。     

相干多普勒测风激光雷达反演厦门市风场及边界层高度的可靠性分析

利用2019年11月5日至12月13日厦门国家高空气象观测站的无线电探空仪数据和同期布设的多普勒测风激光雷达资料,对比分析二者的风速风向和通过梯度法反演的边界层高度。结果表明:(1)风速和风向的决定系数R~2分别达到0.91和0.98,一致性良好,450～1300 m高度范围内,对比效果最佳;(2)利用探空数据的温度廓线和测风激光雷达的雷达回波强度信号通过梯度法分别反演边界层高度,二者一致性很高,只是在边界层出现急流、云层和明显的污染等复杂情况时,探空数据反演的结果会明显大于测风激光雷达;(3)在出现降雨时,雷达探测高度明显降低,无法有效反演边界层高度。测风激光雷达不仅可以满足边界层内风场的精细化连续探测,也可以反演边界层高度。     

珠三角地区日最大混合层高度及其对 区域空气质量的影响

大气混合层高度是影响大气扩散的主要因子之一,其对空气质量评估与污染物的存储量及分布起着重要作用。利用2014年4月至2018年3月珠三角地区香港（沿海站点）和清远（内陆站点）气象探空数据,采用干绝热曲线法估算代表大气垂直方向上大气混合能力的最大混合层高度,探讨沿海与内陆地区混合层高度的差异性,并将最大混合层高度估算值与地面观测的污染物浓度进行相关性分析。结果表明：珠三角沿海与内陆地区的混合层变化具有典型局地特征,沿海日最大混合层高度普遍低于内陆,两地平均高度分别为982 m和1198 m。区域混合层高度的空间差异性由多方面原因造成,其中温度日较差起到关键作用。由于海洋水体的气温调节作用,沿海地区温度日较差较小,因此混合层发展相对较低。珠三角地区各污染物浓度与混合层高度的相关性有较大差异,其中以CO为代表的一次污染物与混合层高度间呈显著负相关,以O₃为代表的二次污染物与混合层高度间则呈显著正相关,而颗粒物作为多源性污染物（既有一次排放,又有二次生成）,其与混合层高度之间的相关性较弱。     

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.