西北干旱和半干旱地区地表蒸发计算方法评估应用

针对我国西北地区干旱和半干旱的气候特征,以及相对简单的地表覆盖特点,论证了互补相关模型计算西北地区蒸发量的可行性。通过中国科学院山东禹城试验站农田观测资料,分析了植被多样性分布对互补相关模型的影响。得出植被的多样性以及植被物候变化的非同步变化是造成互补相关模型中关键参数——大尺度平流参数季节变化的重要原因。并利用遥感反演参数和气象台站观测数据计算了覆盖类型相对单一的西北地区蒸发。结果表明,蒸发与地表覆盖类型的空间分布一致,与气象台站20 cm蒸发皿观测数据呈现空间互补特征,其季节变化与区域降水、温度和植被物候变化相吻合。       

非均匀陆面条件下区域蒸散量计算的遥感模型

非均匀陆面条件下的区域蒸散计算是一个复杂的问题。文中首先在利用遥感资料求取地表特征参数 (如植被覆盖度、地表反照率等 )的基础上 ,建立了裸露地表条件下的裸土蒸发和全植被覆盖条件下植被蒸腾计算模型 ,然后结合植被覆盖度 (植被的垂直投影面积与单位面积之比 )给出非均匀陆面条件下的区域蒸散计算方法。实测资料验算表明该模型具有较高的计算精度。文章最后利用该模型对中国北方地区的蒸散量进行了计算 ,并对该研究区蒸散的特点进行了分析     

区域气候模式(RegCM2)与水文模式耦合的数值试验

建立不均匀的地表径流算法，修改RegCM2的径流方案，设计了一个适合与气候模式RegCM2耦合、能模拟水文站流量的汇流模式，模拟了1998年6、7、8月降水的空间分布，分析了该径流方案对降水、地表热量通量、地表径流、土壤湿度产生的影响。结果表明：（1）该方案在模拟1998年长江流域降水空间分布上有一定的合理性，在一定程度上改善了降水量的模拟，其影响大致是总降水量的10％；（2）地表径流方案改变了地面向大气输送的热量通量，这种作用随时间发生变化，这种变化与地表水分的再分配有关；（3）本方案计算的土壤渗透率较强，在暴雨初期，产生径流较少，而在暴雨后期土壤湿度增大，产生的地表径流较大，这一点更符合洪水形成的特点；（4）水文一气候耦合模式模拟了两个长江水文站的流量，模拟值基本反映了实测值的变化趋势，也表明耦合模式基本能反映1998年夏季长江流域大暴雨期间的地表水文过程。     

旱作春小麦农田蒸散与能量平衡

利用蒸渗计等仪器设备的观测资料，分析了旱作春小麦农田的蒸散耗水状况。得在降水基本正常的年份，农田水分收支基本平衡，农田休闲期蒸发耗水约占期间降水的７０％。旱作小麦的实际蒸散耗水峰值期与小麦生理需水峰值期并不一定吻合。与有灌溉条件的春小麦相比，平均日蒸散量偏小。各生育期的平均叶面积系数与平均日蒸散强度存在较好的线性关系。小麦主要生育期的地表能量平衡分析表明：潜热耗能占主导地位，乳熟－黄熟生育时段以感热、潜热耗能为主的出现几率各为５０％。由于旱作春小麦株冠不能完全覆盖棵间裸地，地表向下的热量传输耗能也占有一定的比重。     

基于遥感资料的陆面水循环模拟及检验

利用遥感资料和10年逐日气象资料,采用逐步插值方法(SIA),修改VIC-3L(Three-Layer Variable Infiltrat1on Capacity hydrological model)模式中多种覆盖类型蒸发的计算方法,计算植被对地表水循环的贡献.模拟地处南水北调工程水源区的褒河(陕南汉江支流)流域地表水的空间分布,将应用遥感资料与<中国资源环境数据库>资料的模拟结果进行比较,并用实测资料进行检验,结果表明:(1)有效的气候要素空间插值方法,是陆面水循环模拟的基础;(2)4 km分辨率DEM(Digital Elevation Model data)得到褒河流域的河网,与实际状况相吻合;(3)利用VIC-3L模拟褒河流域的陆面水文过程中,水分分量的分布合理,物理成因清晰,量值大小准确;(4)基于遥感资料模拟的地表水平衡分量的空间分布更为合理,尤其是改善了河流流量峰值的模拟;(5)遥感植被类型、叶面积指数较好地反映了实际地表覆盖状况和叶面积指数的季节变化情况.     

植被归一化植被指数对气候因子的空间敏感性

基于全球土地利用类型和覆盖度,利用生长季多年平均（1982～2015年）归一化植被指数（Normalized Difference Vegetation Index,NDVI）和气候平均态（气温、降水量）数据,讨论了全球植被格局与气候因子之间的关系,建立了两者之间的多元回归模型,并分析了植被对气温和降水气候态敏感性的特征。植被与气候因子在气候梯度上存在明显的对应关系,回归模型可较好拟合气候态NDVI的全球分布格局,拟合与观测NDVI的相关系数达0.90。其中,常绿阔叶林、混交林、常绿针叶林、落叶阔叶林、农田和木本稀树草原空间分布的拟合能力较好（r>0.8）。不同土地覆盖类型的NDVI对气温、降水气候态的空间敏感性特征不同。整体而言,植被对气温和降水的敏感性呈现反相关关系（r=-0.6）。不同土地覆盖类型对气温表现出正/负敏感性,寒带灌木对气温的敏感性最强,而农作物、草原、裸地对气温负敏感性较大;植被对降水的敏感性均表现出正敏感性,其中落叶针叶林、草原和稀树草原对降水的空间敏感性较强。     

辽宁省植被覆盖度时空变化特征及其对气候变化的响应

基于MODIS-NDVI遥感数据集计算辽宁省2001—2019年植被覆盖度,并结合MODIS土地覆盖产品和辽宁省61个气象观测站气温、降水资料,重点探讨辽宁省5种主要植被类型的植被覆盖度时空变化特征及其对气候变化的响应。结果表明:(1)辽宁省多年平均植被覆盖度为0.48,且呈现"东高西低"的空间分布特征。近19 a来,辽宁省绝大部分地区植被覆盖度呈上升趋势,整体上每10 a增加0.036;主要植被类型作物、草原、落叶阔叶林、多树草原和稀树草原的植被覆盖度均呈显著上升趋势,其中草原增加速率最大,作物增加速率最小。(2)作物的植被覆盖度在暖温带半干旱区与降水存在正相关,与气温存在负相关,而在暖温带半湿润区则与降水和气温均呈现正相关;草原的植被覆盖度对降水响应比气温强烈,而落叶阔叶林、多树草原和稀树草原对气温较为敏感。(3)主要植被类型的植被覆盖度对气温、降水的时滞响应不同。作物和草原在生长季内对上一月降水有时滞响应,而落叶阔叶林、多树草原和稀树草原在生长季末期对上一月气温、降水有时滞响应。     

西藏地区植被覆盖特征与气象因子的相关分析

利用NOAA-AVHRR数据,以归一化植被指数(NDVI)作为植被覆盖特征的指标,研究1982—2000年西藏地区地表植被覆盖的空间分布及时间序列变化,讨论了各季节植被活动的状况。并通过计算不同时段降水量与年最大NDVI、典型区域NDVI与降水、平均温度之间的相关系数,分析了不同植被类型下气象要素与植被覆盖的关系。     

降水气候强迫下非均匀地表区域平均径流的一种参数化方案

将任一中尺度区域的平均瞬间径流率考虑为区域平均降水量和地表土壤层水分渗透垦的余项.根据降水量在地理空间上分布的实测资料拟合其空间概率密度函数(PDF),并结合土壤入渗物理过程的数学描述及其经验公式,精确估计出地表土壤渗透率及其空间分布,由此建立区域地表径流率的统计-动力学估计方案.换言之,区域内地表产流率可视为区域平均降水量与区域平均的土壤下渗量之差值,而区域内土壤的平均下渗量又町分为非饱和区和饱和区两部分的下渗量来分别计算.就陆面水分循环的物理过程而言,地表入渗现象是在一定的下垫面特性基础上,由一定的水分供应源而形成的.根据大气降水向地表层输送水分的物理过程,在满足植被表层覆盖需水(截流水)和地表层土壤人渗水基础上,多余的降水量才会形成地表径流.凶此,推求地表产流率的主要关键在于地表土壤层需水量.为此奉文根据土壤水分通量方程推导出水分入渗公式.又从描述土壤水分和降水的空间PDF出发,推导出非均匀土壤含水量及降水气候强迫所形成的次网格尺度区域平均径流率计算公式.利用长江三角洲地区1996年降水量和土壤特性等实测资料建立区域平均地表径流率的估计公式,并对其影响凶素进行敏感性试验.结果表明,该方法与用Mosaic方法计算的区域径流率(或产流率)结果十分接近.由此可见,该文提出的降水气候强迫下非均匀地表区域平均径流的这种统计-动力参数化方案,具有相当的可靠性与可行性.     

用耦合模型模拟青弋江流域水文对气候变化的响应

用陆面模式SSiB与动态植被模型TRIFFID以及流域地形指数水文模型的耦合模型SSiB4T/TRIFFID模拟了长江下游的青弋江流域植被和水量平衡的动态过程,分析了气温和降水变化对流域径流和蒸发的影响。结果表明：（1）流域气温上升10C,径流减小6.7-9.7%;气温上升20C,径流减小11.7-17.4%;（2）降水增加5%,径流增加9.2-11.6%,降水减小5%,径流减小8.6-11.6%;（3）温度不变仅降水变化对流域蒸发影响很小,温度增加20C,流域总蒸发3-10月份增加8.0-10.7mm,其余月份增加5.4-7.1mm,1月和12月蒸发对温度增加最敏感;（4）气温上升20C,叶面积指数1月和12月增加,5-10月略有减小。降水和气温变化对青弋江流域径流影响明显且与植被类型有关,流域蒸发的变化主要受温度变化的控制。     

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.