Consistency of Tropospheric Water Vapor between Reanalyses and Himawari-8/AHI Measurements over East Asia

High spatiotemporal resolution radiances from the advanced imagers onboard the new generation of geostationary weather satellites provide a unique opportunity to evaluate the abilities of various reanalysis datasets to depict multilayer tropospheric water vapor(WV),thereby enhancing our understanding of the deficiencies of WV in reanalysis datasets.Based on daily measurements from the Advanced Himawari Imager(AHI)onboard the Himawari-8 satellite in 2016,the bias features of multilayer WV from six reanalysis datasets over East Asia are thoroughly evaluated.The assessments show that wet biases exist in the upper troposphere in all six reanalysis datasets;in particular,these biases are much larger in summer.Overall,we find better depictions of WV in the middle troposphere than in the upper troposphere.The accuracy of WV in the ERA5 dataset is the highest,in terms of the bias magnitude,dispersion,and pattern similarity.The characteristics of the WV bias over the Tibetan Plateau are significantly different from those over other parts of East Asia.In addition,the reanalysis datasets all capture the shift of the subtropical high very well,with ERA5 performing better overall.

欧亚大陆积雪分布及其类型划分

利用1966-2012年欧亚大陆1 152个地面气象台站积雪深度资料,对欧亚大陆积雪深度、累计积雪天数和连续积雪天数的空间分布进行了分析,以连续积雪天数为标准对欧亚大陆季节性积雪类型进行了划分,并与应用累计积雪天数对积雪区类型的划分进行了比较研究.结果表明:欧亚大陆积雪分布具有显著纬度地带性特征,积雪深度、累计积雪天数和连续积雪天数的大值分布区均位于俄罗斯平原的东北部、科拉半岛、西西伯利亚平原、中西伯利亚高原以及俄罗斯远东北部大部分区域.与累计积雪天数划分方法相比,利用连续积雪天数对欧亚大陆季节性积雪分区,在前苏联地区积雪类型分区差异并不显著,但蒙古和中国的稳定积雪区明显缩减,青藏高原无稳定积雪区,中国大部分地区为非周期性不稳定积雪区.两种积雪分区划分方法比较结果显示,连续积雪天数划分方法更能体现积雪累积的连续性和持久性,更符合对稳定积雪和不稳定积雪的划分标准.

应用小波收缩方法剔除MODIS热红外波段数据条带噪声

采用多元并扫方式 (1km分辨率 1 0元并扫 ,5 0 0m分辨率 2 0元并扫 ,2 5 0m分辨率 4 0元并扫 )的MODIS传感器由于各探测单元在轨响应差异而引起的条带噪声对MODIS定量产品的反演计算精度造成一定影响。这种影响在MODIS的热红外波段尤其明显。为了尽可能减少这种影响 ,提高MODIS定量产品反演精度 ,提出应用小波收缩方法剔除MODIS数据条带噪声。研究首先使用连续小波变换方法分析MODIS条带数据 ,确定MODIS数据条带噪声在小波系数域中的尺度。其次 ,在分析条带噪声模式的基础上 ,使用小波收缩方法对MODIS 1B数据进行噪声剔除计算。最后比较了分别使用噪声剔除前后的两组MODIS 1B数据反演得到的部分云和大气定量遥感产品 ,结果表明 ,使用剔除噪声后的MODIS

中国寒区水文学研究的新阶段——记我国杰出寒区水文学家叶柏生研究员的创新与贡献

我国杰出的寒区水文学家叶柏生研究员不幸因公殉职.选取了叶柏生研究员若干代表性研究成果,包括与他人的合作研究成果,重点从冰川水文、冻土水文和区域水文变化三方面总结了其对寒区水文学发展所做的创新与贡献.文章列出了每项研究成果的核心内容,并给予了简要评述.所选成果中,冰川水文方面研究涉及冰川对河川径流的调节作用、冰川径流对气候变化的响应机理等;冻土水文研究着重介绍了多年冻土变化对流域径流过程及其变化影响方面的系统性成果;区域水文变化研究方面,选取了降水观测误差修正、气候变化对区域径流的影响等方面的创新成果.这些研究成果极大地提高了我国在世界寒区水文学研究的地位,对认识寒区水文过程及气候变化对水资源的影响具有重要科学意义.

背景场误差样本模拟对同化及数值预报效果的影响

合理估计背景场误差协方差矩阵(B)是做好变分同化的关键环节。利用控制变量随机扰动法(RandomCV)、增长模繁殖法(BGM)及NMC法等3种背景场样本模拟方法,基于WRFDA系统计算B矩阵,对B矩阵的特征及其对同化预报效果的影响进行了研究。B矩阵的特征分析和单点观测试验表明,NMC法与RandomCV法得到的B矩阵误差方差较大,在同化中观测的权重更大;RandomCV法得到的B矩阵,背景场误差中变量的长度尺度更大,说明同化中观测的水平影响范围更大。连续循环同化和预报试验表明:应用RandomCV法计算得到的B矩阵分析与预报的效果明显优于系统自带的以及BGM法得到的B矩阵,且效果与NMC法相当与NMC方法相比,采用RandomCV方法产生背景场样本具有时间和人力成本相对低的优点。

我国冰冻圈及其变化的脆弱性与评估方法

基于脆弱性研究文献,从脆弱性概念、评估内容、评估方法等方面系统总结了脆弱性评估研究的特点和趋势,开展脆弱性评估时应注意的系统、关注特征、灾害、尺度与时间等五个方面,剖析了区域尺度脆弱性评估涉及的一系列问题,介绍了脆弱性评估方法及其优缺点和适用性.探讨了我国冰冻圈及其变化的脆弱性概念,分析了二者之间的关系,论述了我国冰冻圈及其变化的脆弱性评估内容、评估尺度和评估方法.

利用1996年丽江地震序列反演震区应力状态

本文以1996年丽江M_S7.0地震为例,以地震目录和强震震源模型为输入参数,将力学计算与统计分析相结合,反演计算了1986～1996年间丽江震区2.5°×3.0°空间范围内的地壳应力状态.结果表明,丽江地区应力场最大主应力为近南北向,方位角为355°,最小主应力方位角为241°,二者倾角为近水平,应力形因子为O.398,最大主应力大小约为53 MPa,远小于完整岩石的强度.该方法可获得地壳深部绝对应力大小、有效摩擦系数,反演获得的震区地壳深部应力状态的方向性特征与震源机制解结果基本一致;该研究表明地震活动的空间分布可以为我们提供有关强震孕震的地壳深部应力状态信息.这些结果对了解震区应力场具有理论上和实践上的意义.

ETKF协方差膨胀方案对WRFDA混合同化及预报的影响

为讨论集合变换卡尔曼滤波不同协方差膨胀方案对集合-变分混合同化及预报的影响,开展了中国中东部区域一次连续大范围降水过程的连续10天的循环同化和预报试验。结果表明:4种不同的协方差膨胀方案相对于无协方差膨胀的方案,均有效地提高了混合同化和预报的效果。将同化时次之前所有膨胀系数平均值作为新膨胀系数的方案,同化和预报的效果均是最差的;其他3种协方差膨胀方案效果较为接近略有区别:对于风场,将预报误差协方差投影到集合子空间的方案和采用平均新息协方差信息的方案表现较好;对于温度场、湿度场和降水预报,采用平均新息协方差信息的方案和采用了同化时次前两次集合预报比率的方案较好。

雷暴云内电场力对起电和电荷结构的反馈作用

利用美国国家强风暴实验室(NSSL)发展的耦合了详细起电机制和放电过程的中尺度电耦合数值模式WRF(weather research forecasting)-Elec,在NSSL云微物理双参数化方案中增加了电场力对霰、雹粒子降落末速度的影响,完善了WRF-Elec模式的物理过程,建立了双向耦合WRF-Elec模式.利用改进后的WRF-Elec模式,通过敏感性数值实验,定量分析了雷暴云内电场力对起电和电荷结构的反馈作用.结果发现:雷暴云发展旺盛阶段,由于电场力作用,霰、雹粒子质量加权平均降落末速度的瞬时变化极值可以超过4 m/s,但这种情况仅出现在雷暴云内局部区域,并且维持时间较短;电场力对直径小且数浓度较低的霰和雹粒子影响较大,但这种影响不是由单一物理量决定,而是由电场强度和霰、雹粒子的电荷密度、极性以及粒子的直径与数浓度共同决定;电场力通过对霰、雹粒子降落末速度的调整,增强了雷暴云内感应、非感应起电率,且前者远大于后者,云内局部产生-0.6—1.2 n C/m^3总电荷密度的变化,从而使电荷结构重新分布,局部垂直电场强度增强5 k V/m,总闪电数增加,与此同时,雷暴云内降水粒子的微观增长过程也发生改变.总体上,电场力对雷暴云起电过程的作用为正反馈,电场力对雷暴云电荷结构的反馈作用不可忽略.

Temperature and Water Vapor Weighting Functions from Radiative Transfer Equation with Surface Emissivity and Solar Reflectivity

Linearization of Radiative Transfer Equation (RTE) is the key step in physical retrieval of atmospheric temperature and moisture profiles from InfRared (IR) sounder observations. In this paper, the successive forms of temperature and water vapor mixing ratio component weighting functions are derived by applying one term variation method to RTE with surface emissivity and solar reflectivity contained. Retrivals of temperature and water vapor mixing ratio profiles from simulated Atmospheric Infrared Sounder (AIRS) observations with surface emissivity and solar reflectivity are presented.

Atmospheric Profile Retrieval with AIRS Data and Validation at the ARM CART Site

The physical retrieval algorithm of atmospheric temperature and moisture distribution from the Atmospheric InfraRed Sounder （AIRS） radiances is presented. The retrieval algorithm is applied to AIRS clear-sky radiance measurements. The algorithm employs a statistical retrieval followed by a subsequent nonlinear physical retrieval. The regression coefficients for the statistical retrieval are derived from a dataset of global radiosonde observations （RAOBs） comprising atmospheric temperature, moisture, and ozone profiles. Evaluation of the retrieved profiles is performed by a comparison with RAOBs from the Atmospheric Radiation Measurement （ARM） Program Cloud And Radiation Testbed （CART） in Oklahoma, U. S. A.. Comparisons show that the physically-based AIRS retrievals agree with the RAOBs from the ARM CART site with a Root Mean Square Error （RMSE） of 1K on average for temperature profiles above 850 hPa, and approximately 10% on average for relative humidity profiles. With its improved spectral resolution, AIRS depicts more detailed structure than the current Geostationary Operational Environmental Satellite （GOES） sounder when comparing AIRS sounding retrievals with the operational GOES sounding products.

Analysis of Types and Changes of Village-level Economy in Rural Gongyi City, Henan Province Since 1990

Based on the research on the rural living standard in China in terms of annual net income per capita, we de-fine six types of village-level economy, i.e. "to be extremely poor", "to make a basic living", "to dress warmly and eat one's fill", "to try to enrich (to disengage poverty)", "to be well-off" and "to be affluent". The data of average annual net income of all the 292 villages between 1990 and 2004 in rural Gongyi City, Henan Province were collected, veri-fied and classified. By using standard deviation, coefficient of variation and regression analysis, it is found that the Gongyi's rural economy has boosted up remarkably from the relative-poverty and absolute-poverty stages in 1990 to the well-off in 2004. However, the absolute differences between villages present a trend of enlargement, while the rela-tive differences fluctuating. On the other hand, spatial analysis of village-level economy shows that most villages with relatively high economic development level were located along national expressway and most villages with abso-lute-poverty lay in remote mountainous areas in 1990. Since the 1990s, the rapid urbanization and industrialization have had strongly positive effects on rural economic growth. Initial economic foundation, natural resources and tradi-tional techniques also contribute to village economy. From the perspective of geography, villages with location advan-tages, such as near urban center or industrial parks, have more chances for their economic development and the "core-periphery" economic structure has been presented in the process of rural development.

Computing Streamfunction and Velocity Potential in a Limited Domain of Arbitrary Shape.Part I:Theory and Integral Formulae

The non-uniqueness of solution and compatibility between the coupled boundary conditions in computing velocity potential and streamfunction from horizontal velocity in a limited domain of arbitrary shape are revisited theoretically with rigorous mathematic treatments.Classic integral formulas and their variants are used to formulate solutions for the coupled problems.In the absence of data holes,the total solution is the sum of two integral solutions.One is the internally induced solution produced purely and uniquely by the domain internal divergence and vorticity,and its two components（velocity potential and streamfunction） can be constructed by applying Green＇s function for Poisson equation in unbounded domain to the divergence and vorticity inside the domain.The other is the externally induced solution produced purely but non-uniquely by the domain external divergence and vorticity,and the non-uniqueness is caused by the harmonic nature of the solution and the unknown divergence and vorticity distributions outside the domain.By setting either the velocity potential（or streamfunction） component to zero,the other component of the externally induced solution can be expressed by the imaginary（or real） part of the Cauchy integral constructed using the coupled boundary conditions and solvability conditions that exclude the internally induced solution.The streamfunction（or velocity potential） for the externally induced solution can also be expressed by the boundary integral of a double-layer（or singlelayer） density function.In the presence of data holes,the total solution includes a data-hole-induced solution in addition to the above internally and externally induced solutions.

Continuous-wave mid-infrared intracavity singly resonant optical parametric oscillator based on periodically poled lithium niobate

This paper reports a continuous-wave （CW） mid-infrared intracavity singly resonant optical parametric oscillator based on periodically poled lithium niobate （PPLN） pumped by a diode-end-pumped CW Nd：YVO4 laser. Considering the thermal lens effects, it adopted an optical ballast lens and the near-concentric cavity for better operation. At the PPLN＇s grating period of 28.5 μm and the temperature of 140℃, the maximum idler output power of 155 mW at 3.86 μm has been achieved when the 808 nm pump power is 8.5 W, leading to an optical-to-optical conversion efficiency of 1.82%.

Use of Total Precipitable Water Classification of A Priori Error and Quality Control in Atmospheric Temperature and Water Vapor Sounding Retrieval

This study investigates the use of dynamic a priori error information according to atmospheric moistness and the use of quality controls in temperature and water vapor profile retrievals from hyperspectral infrared （IR） sounders. Temperature and water vapor profiles are retrieved from Atmospheric InfraRed Sounder （AIRS） radiance measurements by applying a physical iterative method using regression retrieval as the first guess. Based on the dependency of first-guess errors on the degree of atmospheric moistness, the a priori first-guess errors classified by total precipitable water （TPW） are applied in the AIRS physical retrieval procedure. Compared to the retrieval results from a fixed a priori error, boundary layer moisture retrievals appear to be improved via TPW classification of a priori first-guess errors. Six quality control （QC） tests, which check non-converged or bad retrievals, large residuals, high terrain and desert areas, and large temperature and moisture deviations from the first guess regression retrieval, are also applied in the AIRS physical retrievals. Significantly large errors are found for the retrievals rejected by these six QCs, and the retrieval errors are substantially reduced via QC over land, which suggest the usefulness and high impact of the QCs, especially over land. In conclusion, the use of dynamic a priori error information according to atmospheric moistness, and the use of appropriate QCs dealing with the geographical information and the deviation from the first-guess as well as the conventional inverse performance are suggested to improve temperature and moisture retrievals and their applications.

Water Vapor and Cloud Radiative Forcings over the Pacific Ocean Simulated by the LASG/IAP AGCM:Sensitivity to Convection Schemes

Characteristics of the total clear-sky greenhouse effect （GA） and cloud radiative forcings （CRFs）, along with the radiative-related water vapor and cloud properties simulated by the Spectral Atmospheric Model developed by LASGIAP （SAMIL） are evaluated. Impacts of the convection scheme on the simulation of CRFs are discussed by using two AMIP （Atmospheric Model Inter-comparison Project） type simulations employing different convection schemes： the new Zhang-McFarlane （NZH） and Tiedtke （TDK） convection schemes. It shows that both the climatological GA and its response to El Nio warming are simulated well, both in terms of spatial pattern and magnitude. The impact of the convection scheme on GA is not significant. The climatological longwave CRF （LWCRF） and its response to El Nio warming are simulated well, but with a prominently weaker magnitude. The simulation of the climatology （response） of LWCRF in the NZH （TDK） run is slightly more realistic than in the TDK （NZH） simulation, indicating significant impacts of the convection scheme. The shortwave CRF （SWCRF） shows large biases in both spatial pattern and magnitude, and the results from the TDK run are better than those from the NZH run. A spuriously excessive negative climatological SWCRF over the southeastern Pacific and an insufficient response of SWCRF to El Nio warming over the tropical Pacific are seen in the NZH run. These two biases are alleviated in the TDK run, since it produces vigorous convection, which is related to the low threshold for convection to take place. Also, impacts of the convection scheme on the cloud profile are discussed.

Weak ENSO Asymmetry Due to Weak Nonlinear Air–Sea Interaction in CMIP5 Climate Models

State-of-the-art climate models have long-standing intrinsic biases that limit their simulation and projection capabilities. Significantly weak ENSO asymmetry and weakly nonlinear air-sea interaction over the tropical Pacific was found in CMIP5 （Coupled Model Intercomparison Project, Phase 5） climate models compared with observation. The results suggest that a weak nonlinear air-sea interaction may play a role in the weak ENSO asymmetry. Moreover, a weak nonlinearity in air-sea interaction in the models may be associated with the biases in the mean climate--the cold biases in the equatorial central Pacific. The excessive cold tongue bias pushes the deep convection far west to the western Pacific warm pool region and suppresses its development in the central equatorial Pacific. The deep convection has difficulties in further moving to the eastern equatorial Pacific, especially during extreme El Nifio events, which confines the westerly wind anomaly to the western Pacific. This weakens the eastern Pacific El Nifio events, especially the extreme E1 Nifio events, and thus leads to the weakened ENSO asymmetry in climate models. An accurate mean state structure （especially a realistic cold tongue and deep convection） is critical to reproducing ENSO events in climate models. Our evaluation also revealed that ENSO statistics in CMIP5 climate models are slightly improved compared with those of CMIP3. The weak ENSO asymmetry in CMIP5 is closer to the observation. It is more evident in CMIP5 that strong ENSO activities are usually accompanied by strong ENSO asymmetry, and the diversity of ENSO amplitude is reduced.

Inventory of Atmospheric Pollutant Emissions from Burning of Crop Residues in China Based on Satellite-retrieved Farmland Data

The burning of crop residues emits large quantities of atmospheric aerosols.Published studies have developed inventories of emissions from crop residue burning based on statistical data.In contrast,this study used satellite-retrieved land-cover data(1 km×1 km)as activity data to compile an inventory of atmospheric pollutants emitted from the burning of crop residues in China in 2015.The emissions of PM10,PM2.5,VOCs,NOx,SO2,CO,and NH3 from burning crop straw on nonirrigated farmland in China in 2015 were 610.5,598.4,584.4,230.6,35.4,3329.3,and 36.1 Gg(1 Gg=109 g),respectively;the corresponding emissions from burning paddy rice residues were 234.1,229.7,342.3,57.5,57.5,1122.1,and 21.5 Gg,respectively.The emissions from crop residue burning showed large spatial and temporal variations.The emissions of particulate matter and gaseous pollutants from crop residue burning in nonirrigated farmland were highest in east China,particularly in Shandong,Henan,Anhui,and Sichuan provinces.Emissions from burning paddy rice residue were highest in east and central China,with particularly high levels in Shandong,Jiangsu,Zhejiang,and Hunan provinces.The monthly variations in atmospheric pollutant emissions were similar among different regions,with the highest levels observed in October in north,northeast,northwest,east,and southwest China and in June and July in central and south China.The developed inventory of emissions from crop residue burning is expected to help improve air quality models by providing high-resolution spatial and temporal data.

Application and comparison of coaxial correlation diagram and hydrological model for reconstructing flood series under human disturbance

Intense human activities have greatly changed the flood generation conditions in most areas of the world, and have destroyed the consistency in the annual flood peak and volume series. For design flood estimation, coaxial correlation diagram and conceptual hydrological model are two frequently used tools to adjust and reconstruct the flood series under human disturbance. This study took a typical mountain catchment of the Haihe River Basin as an example to investigate the effects of human activities on flood regime and to compare and assess the two adjustment methods. The main purpose is to construct a conceptual hydrological model which can incorporate the effects of human activities. The results show that the coaxial correlation diagram is simple and widely-used, but can only adjust the time series of total flood volumes. Therefore, it is only applicable under certain conditions(e.g. There is a strong link between the flood peaks and volumes and the link is not significantly affected by human activities). The conceptual model is a powerful tool to adjust the time series of both flood peak flows and flood volumes over different durations provided that it is closely related to the catchment hydrological characteristics, specifically accounting for the effects of human activities, and incorporating expert knowledge when estimating or calibrating parameters. It is suggested that the two methods should be used together to cross check each other.

High-Power Continuous-Wave Diode-End-Pumped Intracavity Frequency Doubled Nd：YVO4 Laser at 671 nm with a Compact Three-Element Cavity

We report a high-power high-efficient continuous-wave （cw） diode-end-pumped Nd：YVO4 1342-nm laser with a short plane-parallel cavity and an efficient cw intracavity frequency-doubled red laser at 671 nm with a compact three-element cavity. At incident pump power of 20.6 W, a maximum output power of 7 W at 1342 nm is obtained with a slope effciency of 37.3%. By inserting a type-Ⅰ critical phase-matched LBO crystal as intracavity frequencydoubler, a cw red output as much as 2.85-W is achieved with an incident pump power of 16.9 W, inducing an optical-to-optical conversion efficiency of 16.9%. To the best of our knowledge, this is the highest output of diodepumped solid-state Nd：YVO4 red laser. During half an hour, the red output is very stable, and the instability of output power is less than 1%.