不同初始场及陆面方案对青藏高原中东部积雪消融过程的模拟研究

基于WRF(Weather Research and Forecasting)模式,本研究使用更为准确的气象站点及卫星遥感积雪资料替换初始场中积雪深度、雪水当量等积雪数据,对2014年2月17-27日青藏高原中东部一次积雪消融过程进行模拟研究,评估WRF模式中CLM(Community Land Model)、Noah-LSM(Noah land surface model)和Noah-MP(Noah-Multiparameterization Land Surface Model)3种陆面过程方案对该次积雪消融过程的模拟性能。结果表明：3种陆面过程方案均能较好地再现2 m气温、积雪深度和反照率的日变化趋势,但各试验模拟效果有一定差异。气象站点及卫星遥感积雪资料作为初始场时,CLM陆面过程方案模拟的2 m气温平均误差最小,为0.002℃;Noah-LSM陆面过程方案中2 m气温均方根误差（4.01℃）和平均绝对误差（3.30℃）最小,但昼夜温差较观测显著偏小;同时CLM陆面过程方案模拟的积雪深度均方根误差、平均误差和平均绝对误差均最小,分别为4.70 cm、-1.25 cm和2.75 ...     

太湖蒸发量对未来气候变化响应的初步模拟分析

与深水湖泊相比,太湖等浅水湖泊更容易发生富营养化和水资源危机,且对气候变化的响应更为敏锐。本文利用气候模式产品数据驱动CLM4-LISSS湖泊陆面过程模型,模拟分析未来(2010—2100年)RCP2.6、RCP4.5、以及RCP8.5不同温室气体排放情景下太湖蒸发量的变化特征及其影响因子。结果表明:(1)CLM4-LISSS模型湖表温度的观测值与模拟值的相关系数为0.94,均方根误差为0.85℃,准确的湖表气温模拟使得通量的结果也比较准确,潜热模拟与观测的相关系数在0.78,均方根误差为55.32 W·m~(-2);(2)2010—2100年,三种不同温室气体排放情景下太湖蒸发都呈现增加的趋势,但增量比例不同,RCP2.6,RCP4.5和RCP8.5情景下,蒸发量每10 a增加量分别为23.7 mm,29.2 mm和34.5 mm。蒸发量的增加速率随着辐射强迫的增加而增大,其变化主要受风速与水汽压差的乘积的影响。     

苏州-无锡-常州城市带热岛效应个例研究

应用WRF（weather research and forecasting）及其耦合的多层城市冠层模式BEP（building energy parameterization）,对2013年8月13日长江三角洲地区一次高温天气过程进行了模拟.此次过程盛行东南风,风向与苏州-无锡-常州城市带走向一致.模拟结果表明：苏州-无锡-常州城市带热岛效应明显,热岛强度向下游城市逐渐增加;在东南风作用下,三座城市的热岛连成一片,形成了一个更强大的热岛环流.夜晚,边界层逐渐趋于稳定,热岛环流减弱,有利于热岛温度向下游地区输送.热岛效应导致城市边界层高度明显上升.白天太湖产生强盛的湖风对其周边城市影响显著,来自太湖的冷气团导致无锡和常州边界层内热岛强度明显下降,抑制城市热岛向上发展,削弱了无锡与常州两城市热岛间的联系.白天太湖使得无锡和常州边界层高度明显下降.     

苏州—无锡—常州城市带热岛效应个例研究

应用WRF(weather research and forecasting)及其耦合的多层城市冠层模式BEP(building energy parameterization),对2013年8月13日长江三角洲地区一次高温天气过程进行了模拟。此次过程盛行东南风,风向与苏州—无锡—常州城市带走向一致。模拟结果表明:苏州—无锡—常州城市带热岛效应明显,热岛强度向下游城市逐渐增加;在东南风作用下,三座城市的热岛连成一片,形成了一个更强大的热岛环流。夜晚,边界层逐渐趋于稳定,热岛环流减弱,有利于热岛温度向下游地区输送。热岛效应导致城市边界层高度明显上升。白天太湖产生强盛的湖风对其周边城市影响显著,来自太湖的冷气团导致无锡和常州边界层内热岛强度明显下降,抑制城市热岛向上发展,削弱了无锡与常州两城市热岛间的联系。白天太湖使得无锡和常州边界层高度明显下降。     

基于三种WRF陆面过程方案的东南亚毁林增温响应研究

选取东南亚中南半岛地区作为模拟区域,基于高分辨率的遥感观测森林变化数据和WRF数值模式,设计毁林前后的两种情景对旱季气候进行模拟,评估NoahMP、CLM和Noah mosaic三种陆面过程方案对热带毁林增温响应的模拟能力。结果表明,CLM方案在模拟历史气温中有着更好的表现,Noah mosaic方案的结果存在明显低估。然而,对比毁林前后两种情景的模拟结果,本文发现,只有采用了"次格网"方式的Noah mosaic方案较好地模拟出毁林增温响应特征。在格网尺度采用"主导类型"计算方式的NoahMP方案没有合理地呈现出森林损失对区域气候的影响。理论上,CLM模式在计算中同时考虑格网内所有植被类型,然而本文发现CLM方案在主导类型不变的格网对森林损失比例不敏感,而且对毁林反馈的模拟结果与NoahMP方案的结果更接近。据此推测,在WRF模式耦合CLM方案的过程中,格网内参数处理方式可能产生了错误,实际采用的是"主导类型"方式。在模拟土地覆盖类型变化对气候的影响时,本文推荐使用Noah mosaic方案。同时,建议在未来版本的WRF模型中修正目前耦合的CLM方案关于次格网方法的处理方式,提供更合理的水热通量模拟。     

陆面过程参数化对太湖地区雷暴过程模拟的影响

以2010年8月发生在太湖地区的一次雷暴过程为例,利用WRF模式进行48 h的短期天气模拟,分析两个陆面参数化方案(Noah方案和RUC方案)对雷暴过程模拟的影响。对比模式结果与实况降水以及太湖地区两个站点的近地面要素表明:雷暴过程对陆面参数化方案的选取较为敏感,不同陆面参数化方案可影响雷暴发生的时间、地点及强度,两种方案的降水中心值差达40 mm以上,其中Noah方案所模拟的降水与实况更为接近。通过对两个方案模拟的物理量场的对比分析发现,RUC方案中对流发展滞后于Noah方案2 h;这表明陆面过程对雷暴等中尺度对流过程有显著的反馈作用,不同陆面参数化方案的使用影响雷暴发生的动力和热力作用,并改变雷暴的时空分布特征。陆面过程通过改变地面热通量输送影响边界层结构,使得水平和垂直方向上的风和温度等发生变化并产生辐合辐散,进而影响对流的启动时间和对流发展强度。由于对不同植被的参数化处理的差异,Noah方案对下垫面特征的描述能力优于RUC方案,尤其是对城市下垫面的处理,这也使得之后该方案模拟的雷暴发生时间更加接近于实况且雷暴过程更加强烈。     

太湖湖陆风背景下的苏州城市化对城市热岛特征的影响

利用南京大学区域边界层模式,选取苏州地区2006年8月12日晴天小风的天气作为背景天气条件进行算例的模拟。分析了该地区在湖陆风环流背景下的城市化进程对热岛特征的影响。结果表明:苏州地区的风场受大尺度系统、湖陆风环流和城市热岛环流系统的共同影响,且湖陆风环流的强度大于热岛环流。苏州地区的城市热岛影响高度可达400 m左右,且在背景风场的影响下,该地区的温度分布会体现出比较明显的下游效应特征,下游效应的距离可达10 km左右。这种下游效应在400 m高度以下比较明显,随高度的升高而逐渐减弱。20 a的城市化进程(1986—2006年)导致苏州城区地表感热通量增加了200 W/m²、潜热通量减少了200 W/m²、湍能增加了0.045 m²/s²、日平均热岛强度增加了0.4 ℃;热岛环流加强,并使白天混合层最大高度升高了400 m,但对夜间逆温层发展的影响较小。太湖湖陆风,使靠近太湖的苏州市郊的混合层、逆温层高度明显下降,但对距离较远的苏州市区的混合层、逆温层高度影响不大。且太湖的存在对于城市化进程导致的苏州地区热岛强度增加、地表湍流动能增加、感热通量增加、潜热通量减少的趋势有比较明显的缓解作用。     

昆明城市热岛效应的数值模拟研究

本文利用WRF（V3.9.1）模式中耦合Noah/SLUCM方案作为Control试验,研究了土地利用类型（Md04试验）、陆面过程（NoUCM试验）和湖泊（Nolake试验）对城市热岛强度及昆明城市气象要素水平、垂直的时空分布影响。主要结论如下:（1）四个试验城市热岛强度的平均日变化趋势相似,白天城市热岛强度较弱、夜间较强,在20时（北京时,下同）左右达到最大值。城市冠层（湖泊）对城市热岛有较明显的减（增）温,Control-NoUCM（Nolake）试验中,平均日最大差值为?0.79°C（+1.07°C）。（2）从能量平衡方程分析Control-Md04试验,感热（潜热）通量的差值为+46.18（?79.71）W m?2,潜热通量释放大于感热通量的绝对值。Control-NoUCM试验中,感热（潜热）通量的差值为?40.88（+29.60）W m?2;因NoUCM试验未考虑几何建筑物储热与遮挡,太阳辐射大部分被地表所吸收,导致感热通量偏大。（3）四种试验中,15（07）时边界层高度达到最大（小）值。NoUCM（Nolake）试验中城市边界层高度分别降低103 m（32 m）左右,而Md04试验中城市边界层高度增加102 m左右。（4）湖泊（滇池）对城市热岛环流影响的试验表明,湖泊上空垂直运动较弱,但水平方向湖陆风较大,这有利于向城市输送水汽,增加干空气湿度,使城市中空气的水汽含量增加,同时增大潜热能量释放,降低感热通量,减小了垂直温度梯度。     

大理苍山—洱海局地环流的数值模拟

利用耦合了湖泊模型的WRF_CLM模式模拟了秋季大理苍山—洱海地区的局地环流特征。结果表明:模式对近地面温度、风向、风速的模拟与观测基本一致,模拟结果能较好地再现该地区山谷风和湖陆风相互作用的局地环流特征。在秋季,大理苍山的谷风起止时间为08:00~17:00(北京时,下同),湖风起止时间为09:00~19:00。局地环流受高山地形及洱海湖面影响明显,山谷风形成早于湖陆风1 h,夜间山风、陆风强盛于白天谷风、湖风。白天苍山谷风与洱海湖风的叠加作用会驱动谷风到达2600 m的高度,而傍晚最先形成的苍山山风则会减弱洱海的湖风环流。夜间盆地南部在两侧山风、陆风的共同作用下,形成稳定而持续的气旋式环流。日出以后,对流边界层迅速发展,边界层高度逐渐增高。陆地17:00温度达到最高,边界层高度也达到峰值2000 m,之后逐渐降低。日落后形成稳定边界层,边界层高度在夜间基本保持在100 m。相对于陆地,湖面白天边界层高度低300 m,夜间边界层高度高100 m。     

LASG/IAP 大气环流谱模式对陆面过程的敏感性试验

将中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室最近发展的高分辨率全球大气环流谱模式SAMIL-R42L26, 分别与两个陆面模式--NCAR通用陆面模式CLM和简化的简单生物圈模式SSiB进行耦合.在比较了两个陆面过程模式, 指出CLM改进方面的基础上, 通过分析两个陆气耦合模式所模拟的陆气通量交换结果, 指出新版本的陆气耦合模式(SAMIL-R42L26与CLM耦合)对表面感热、温度、降水率、潜热通量和海平面气压场的模拟能力大大提高, 尤其对于夏季表面感热通量场, 使亚洲北部和东南部、格陵兰岛以及北美洲大部分地区的数值从100 W/m2 降低到接近60 W/m2, 与NCEP再分析资料一致.新版本的陆气耦合模式模拟陆地表面能量收支趋于平衡, 为下一步发展海-陆-气-冰耦合气候系统模式提供保障.采用CLM陆面模式, SAMIL-R42L26能较好地模拟亚洲季风区地表感热和潜热的季节演变趋势, 而采用SSiB陆面模式的结果, 则存在较大误差.文中的结果表明, 不同的陆面模式所模拟的大气下垫面(包括洋面)通量发生的变化, 通过陆气耦合过程产生的影响不仅仅是局地性的, 而且是全球范围的.     

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.     

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.     

IMPACT OF ATMOSPHERIC LOW-FREQUENCY OSCILLATION ON THE IMMIGRATION OF NILAPARVATA LUGENS(STL) IN HUNAN AND JIANGXI PROVINCES

Various features of the atmospheric environment affect the number of migratory insects, besides their initial population. However, little is known about the impact of atmospheric low-frequency oscillation(10 to 90 days) on insect migration. A case study was conducted to ascertain the influence of low-frequency atmospheric oscillation on the immigration of brown planthopper, Nilaparvata lugens(Stl), in Hunan and Jiangxi provinces. The results showed the following:(1) The number of immigrating N. lugens from April to June of 2007 through 2016 mainly exhibited a periodic oscillation of 10 to 20 days.(2) The 10-20 d low-frequency number of immigrating N. lugens was significantly correlated with a low-frequency wind field and a geopotential height field at 850 h Pa.(3) During the peak phase of immigration, southwest or south winds served as a driving force and carried N. lugens populations northward, and when in the back of the trough and the front of the ridge, the downward airflow created a favorable condition for N. lugens to land in the study area. In conclusion, the northward migration of N. lugens was influenced by a low-frequency atmospheric circulation based on the analysis of dynamics. This study was the first research connecting atmospheric low-frequency oscillation to insect migration.     

A COMPARATIVE ANALYSIS OF ATMOSPHERIC AND OCEANIC CONDITIONS BEFORE THE OCCURRENCE OF TWO TYPES OF EL NIO EVENTS

The atmospheric and oceanic conditions before the onset of EP El Ni?o and CP El Ni?o in nearly 30 years are compared and analyzed by using 850 hPa wind, 20℃ isotherm depth, sea surface temperature and the Wheeler and Hendon index. The results are as follows: In the western equatorial Pacific, the occurrence of the anomalously strong westerly winds of the EP El Ni?o is earlier than that of the CP El Ni?o. Its intensity is far stronger than that of the CP El Ni?o. Two months before the El Ni?o, the anomaly westerly winds of the EP El Ni?o have extended to the eastern Pacific region, while the westerly wind anomaly of the CP El Ni?o can only extend to the west of the dateline three months before the El Ni?o and later stay there. Unlike the EP El Ni?o, the CP El Ni?o is always associated with easterly wind anomaly in the eastern equatorial Pacific before its onset. The thermocline depth anomaly of the EP El Ni?o can significantly move eastward and deepen. In addition, we also find that the evolution of thermocline is ahead of the development of the sea surface temperature for the EP El Ni?o. The strong MJO activity of the EP El Ni?o in the western and central Pacific is earlier than that of the CP El Ni?o. Measured by the standard deviation of the zonal wind square, the intensity of MJO activity of the EP El Ni?o is significantly greater than that of the CP El Ni?o before the onset of El Ni?o.     

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