冬季增城丘陵灌木林基本气象要素与湍流特征

介绍了广州增城丘陵灌木林通量站的地理环境、仪器布设,并利用该站2014年12月—2015年2月的观测资料,分析冬季华南丘陵灌木林下垫面近地层的风、温度、气压、湿度和大气稳定度等基本气象要素特征,以及湍流强度、风速归一化标准差的统计特征及通量特征。结果表明:（1）冬季增城站以北向小风为主,风速多小于4 m/s,日平均温度变化剧烈,平均相对湿度为67.94%;（2）大气稳定度集中于“±” 0.5之间,以近中性为主;冬季平均湍流强度I_u、I_v和I_w分别为0.37、0.31和0.19;在风速小于1 m/s时湍流发展旺盛,随风速增大,湍流强度减小并趋于稳定;（3）风速分量归一化标准差δ_u/u_*、δ_v/u_*、δ_w/u_*与稳定度Z/L之间满足莫宁-奥布霍夫相似理论,并呈1/3次方关系,且随|Z/L|的减小而减小,在中性层结（Z/L→0）均匀下垫面,其值为常数,分别为2.286、1.757和1.174;（4）冬季增城站的能量耗散以感热为主,晴天尤为显著;阴天感热与潜热大小相当;平均日间吸收CO₂极值为0.236 mg/(m2·s)。      

气象因素对常州市区PM2.5浓度影响

利用2016—2018年常州市区环境空气细颗粒物数据,结合同期地面气象资料,分析了常州市区PM_2.5以及气象因素的变化特征,并统计分析气象因素对PM_2.5浓度的影响。结果表明：常州市区PM_2.5、降水量、相对湿度和气温等具有明显季节性,呈夏季较高冬季较低,而气压夏季较低冬季较高的特征。相对湿度与PM_2.5呈正相关,即随着相对湿度的增加PM_2.5超标率和平均浓度均增加;降水对PM_2.5具有一定的清除作用,清除率与降水前PM_2.5浓度、降水量、降水强度有关,降水量、降水强度越大,则降水清除效果越好,而降水前PM_2.5浓度较小,则清除率不明显;常州市区偏西风时PM_2.5的超标率和平均浓度较其他风向较高;风速对常州市区PM_2.5的影响呈负相关,即风速越大PM_2.5超标率和平均浓度均减小;常州市区地面天气形势可以分为两种类型,第一种类型表现为气压较低气温较高,PM_2.5超标率以及平均浓度相对较低,而第二种类型表现为气压较高气温较低,PM_2.5超标率以及平均浓度相对较高。     

P-σ坐标系数值模式中气压梯度力的递推算法

在静力关系成立的条件下,推导出气压梯度力的递推算法,用以计算有地形P-σ坐标系模式中的气压梯度力.将递推算法与对应的6种直接计算方案相比后发现,在有气压梯度力的参考大气中,递推算法的计算精度相对于直接计算方案的精度有较大改进,特别是采用递推后的两种经典修正格式和Corby格式,能较大提高水平气压梯度力的计算精度,从而改善数值模式的预报效果和稳定性.     

一种依据初始场构造的参考大气廓线及对高原地区预报的影响

地形追随坐标系下的动力学方程组通常都要引入参考大气的概念来减少气压梯度力的计算误差.参考大气的一般取法是满足静止和静力平衡关系,实际大气被看成在参考大气上的偏差.实际大气与参考大气越接近,它们之间的偏差就越小,计算就越精确.参考大气的取法一般有等温大气、定常温度廓线、初始资料进行水平平均的垂直温度廓线等.本文在一个非静力的中尺度模式中选取经过水平平均的初始温度场的垂直廓线构造模式的参考大气廓线.但是初始资料在三维空间是离散的,这里根据三次样条函数的分段连续光滑的性质,构造了随起报时间不同而不同的,只随高度变化的解析形式的参考大气廓线.构造的参考大气廓线完全满足静力平衡关系,方程中的静力平衡部分可以得到最大限度地扣除.试验表明,与等温大气相比,这种根据初始状态选取的参考大气和实际资料更接近,在大地形附近气压梯度力的计算误差更小,500 hPa高原附近的24 h气压预报更精确.     

沙瓦特指数的一种迭代算法

为了提高沙瓦特指数的计算精度, 在沙瓦特指数计算方案中引入二分法:先用于计算抬升凝结高度上的要素, 进而计算假相当位温, 再依据假相当位温守恒性质用于计算500 hPa气块温度。与其他方案对比表明, 该迭代算法计算的沙瓦特指数精度较高; 与查表法所得的气块温度对比表明, 该迭代算法得到的气块温度的最大绝对误差为1.36 ℃, 平均误差 (-0.68 ℃) 与平均绝对误差 (0.69 ℃) 数值接近; 迭代算法与查表法之间以系统性偏差为主; 此外, 该迭代算法计算量小, 收敛速度快, 具有推广应用价值。     

包含冻结过程的广义位温及位涡特征分析

为了对比分析降水过程中不同表达形式热力学变量和位涡时空分布特点,本文针对2017年7月13～14日吉林省强降水过程,利用模式输出资料对常规位温（θ）、相当位温（θ_e）、包含凝结概率函数的广义位温（θ_Gao）、包含冻结概率函数的广义位温（θ_Wang）和同时涵盖凝结过程与冻结过程（θ_Gu）五种不同形式位温进行计算,并分析五种对应位涡[PV(θ)、PV(θ_e)、PV(θ_Gao)、PV(θ_Wang)、PV(θ_Gu)]与降水的关系。结果表明,引入冻结概率函数的广义位温（θ_Wang）和对应的广义湿位涡PV(θ_Wang)与强降水的对应性更好。θ_Wang与θ_Gao差异集中在降水区对流层中高层5～11 km,前者始终高于后者,最大差异达2.5 K,说明冻结概率函数的引入扩大了广义位温的适用范围,更适合描述降水区湿大气非均匀饱和热力状态。五种位涡的差异主要在降水区上空12 km以下,由θ_Gao和θ_Wang定义的位涡PV(θ_Gao)和PV(θ_Wang)的正负异常中心更为明显。相比于PV(θ_Gao)和PV(θ_Wang)异常值更大,差异可达±0.2 PVU,这主要是由于冻结概率函数的引入增大降水区上空广义位温,促使冻结区的湿位涡异常增强。     

2016—2021年呼和浩特市大气臭氧浓度变化特征分析

随着城市的发展,大气O₃的污染不容忽视。文章以呼和浩特市2016—2021年大气O₃浓度为研究对象,结合同期气象数据,采用相关分析等数理统计方法,分析了呼和浩特市大气O₃浓度年际、季节以及月份动态特征规律,探讨大气O₃浓度与关键影响气象要素之间的关系。结果表明：呼和浩特市大气O₃作为空气首要污染物的天数在逐年增加;6年间O₃浓度年际变化呈下降趋势,以-0.02μg·m^-3的年均速率变化;O₃浓度季节变化明显,夏季最高,其次是春季、秋季、冬季;O₃月平均浓度呈倒“V”形单峰变化,主要在6—8月浓度较高,12月浓度最低,O₃浓度超标日主要集中在6、7月;O₃浓度与气象要素呈密切关系,在非采暖期,O₃-8 h浓度与日平均气温和日照时数具有正相关趋势,与日平均气压、日平均相对湿度以及平均风速具有负相关趋势,降水对大气O₃     

桂林城区大气能见度与颗粒物浓度和气象因子关系研究

利用2015年1月至2017年6月桂林国家基本气象站能见度、相对湿度、气温、气压、降水等气象要素和PM₁₀、PM_2.5、PM_1.0颗粒物质量浓度资料,分析桂林城区大气能见度与颗粒物浓度和气象因子之间关系。结果表明：桂林城区大气能见度和PM₁₀、PM_2.5、PM_1.0呈对数关系,相关系数分别为-0.341、-0.461、-0.509,颗粒物对大气能见度影响在相对湿度为60%—70%时最为显著。在各气象因子中,大气能见度与风速的相关性最好,其次为相对湿度,与风速呈二次函数关系,与相对湿度呈幂指数关系,与温度相关性较小,与气压在秋冬季节呈正相关,相关系数冬季可达0.301,但在春、夏季节相关性不显著;利用颗粒物浓度和气象要素建立8种大气能见度非线性统计回归模型,比较后发现利用PM_1.0、风速、相对湿度、气温等因子建立的不同季节大气能见度拟合公式在实际检验中效果最优,能较好地模拟桂林地区大气能见度的变化。     

三种温湿参数下昆明准静止锋锋面位置及锋生函数诊断的对比分析

针对2008年1~2月昆明准静止锋天气过程,讨论了位温、相当位温和广义位温三个温湿参数及其对应的锋生函数,对比分析结果表明:（1）广义位温因其湿度因子权重过大,不适用于分析以温度梯度定义的昆明准静止锋锋面位置和强度;（2）理想状态下的干大气位温和饱和大气相当位温仅与气压和温度有关,均能较好地反映大气温度场特征,适用于昆明准静止锋的分析,比较而言相当位温对温度锋区的描述更为显著;（3）用位温或相当位温计算锋生函数所得结果直接表达了因温度场变化导致的锋面生消状况,而用广义位温计算锋生函数所得结果同时包含温度锋生（消）和湿度锋生（消）两种情况,不能简洁地判断昆明准静止锋锋面生消状况。     

青藏高原对流层顶高度与臭氧总量及上升运动的耦合关系

根据1979-2008年青藏高原地区14个探空站对流层项气压资料以及同期各标准等压面上的温度资料,分析了不同季节高原上空两类对流层顶高度与高空各层温度之间的关系;在此基础上,结合同期的NCEP／NCAR月平均再分析资料以及NASA提供的TOMS／SBUV月平均臭氧总量资料,分别讨论了高原上升运动以及高原臭氧总量与对流层顸高度的耦合关系。结果表明：高原第一（二）对流层顶高度全年处在300～200hPa（100hPa附近）高度,在季节变化、年际变化以及长期变化趋势上,两类对流层顸高度与各自对应高度层上的温度存在着密切的反相变化关系,当对流层顶高度偏高（低）时,相应高度上的温度偏低（高）。上升运动有助于两类对流层顶高度的抬升,尤其是当高空200（100）hPa附近有上升运动时,有利于第一（二）对流层项高度抬升。各季节高原臭氧总量与第二对流层顶高度均呈显著的负相关关系,当臭氧含量减少（增加）时,该对流层顶高度将偏高（偏低）,近年来伴随着高原臭氧总量的减少,高原第二对流层顸高度出现了明显的抬升。     

The development and initial tests of an atmospheric model based on a vertical coordinate with a smooth transition from terrain following to isentropic coordinates

An atmospheric model (η model) is developed by modifying the UW θ－σ hybrid model. In the η model, the vertical coordinate transforms smoothly from terrain following to isentropic coordinates. The model is developed to capitalize on the inherent advantage of numerical modeling in isentropic coordinates and to eliminate the interface between the sigma planetary boundary layer and isentropic free atmosphere present in the UW θ－σ model. This formulation provides the potential for the data assimilation and the application of higher order schemes. This paper describes the structure of the η model and presents results from initial numerical experiments. The first experiment tests the capability of the η model for simulating the baroclinic development process. In the 48-hr numerical weather forecast experiment, the η model produces reasonable precipitation and synoptic fields at all levels which are similar to those from the UW θ-σ model. The second and third experiments test the capability of the η model for conserving 1) the joint distribution of isentropic potential vorticity (IPV) and proxy ozone and 2) equivalent potential temperature under frictionless and isentropic conditions. These experiments show that distributions of IPV and proxy ozone in the pure isentropic domain and the distributions of prognostic and diagnostic equivalent potential temperature in the model domain remain highly correlated to day 10.     

On the Factors Controlling the Snow Surface and 2-m Air Temperatures Over the Arctic Sea Ice in Winter

Factors controlling the magnitudes of, and short-term variations in, the potential temperatures of the snow surface and the air at the height of 2 m θ_S and θ_{2 m} over Arctic sea ice in winter are analysed. The study addresses the winters of 1986–1987 and 1987–1988, and is based on the temperature, wind, and cloud observations made by Russian drifting ice stations. It also relies on the ERA40 re-analyses of the European Centre for Medium-Range Weather Forecasts, which were utilised to calculate the lateral heat advection at the sites of the ice stations. The cloud cover and wind speed were more important than the heat advection in controlling the magnitudes of θ_{2 m} and θ_S, while on a time scale of 24 h, during steady forcing conditions, the heat advection was the most important factor affecting the changes in θ_S and θ_{2 m}. During changing conditions, and considering individual factors separately, the monthly mean 24-h temperature changes were less than ± 5 °C: the effect of the cloud cover was the largest, and that of the heat advection was the smallest. When simultaneous changes in the three factors were analysed, the seasonal mean temperature changes were even of the order of ±15 °C, with the strongest warming events exceeding 35 K in a single day. The difference θ_S − θ_{2 m} reached its lowest seasonal mean values during conditions of clear skies (−1.3 °C), light winds (−1.3 °C) and warm-air advection (−0.8 °C). θ_S and θ_{2 m} followed each other closely, even during major synoptic-scale temperature variations.     

Spectra of temperature derivatives in the atmospheric surface layer

Measured spectra ofθ _x,θ _y,θ _z, the derivatives of temperature in streamwise, lateral and vertical directions, respectively, indicate that the spectral densities ofθ _z andθ _y are nearly similar but significantly different from the spectral density ofθ _x. The high-frequency parts of the three spectra satisfy, in a qualitative sense, local isotropy requirements. In the high-frequency end of the inertial subrange, the relative behaviour of spectra ofθ _x,θ _y andθ _z is also consistent with local isotropy.     

A Study in Search of Interconnection between Surface Parameters and Surrounding Synoptic And Subsynoptic Features

The paper reveals that the variations in parameters like u*, the scaling velocity and θ*. The scaling tempera-ture during the various phases of monsoon might be linked with subsynoptic features. The rise in u* is mainly connected with the presence of lower tropospheric cyclonic vorticity over a subsynoptic scale of the site. However the variations in θ* is mainly linked with the various phases of monsoon and θ* shows a sharp rise in presence of low level convective cloud.Besides the correlation studies of u and u*, θv and θv* , θv-θv0 and θv* are undertaken. The correlation be?tween θv and θv* is poor. In other two cases correlations are good. Besides u/u* , has shown good coefficient of variation values within the ζ range.     

The increased radius of the Aleutian eddies and their long-term evolution

The Alaskan Stream south of the Blizhnii Strait disintegrates into mesoscale Aleutian eddies, which provide a westward transport of warm water to the Kamchatka Current and upper Oyashio Current areas. Oceanographic observations from 1949 to 2009 are indicative of a substantial rise of temperature and salinity in the intermediate waters of the Aleutian and Kamchatka currents. A long-term temperature trend in the Aleutian Current in the intermediate layer (at isopicnal of 26.75σ_θ) amounted to ΔT = 0.013°C/year. A positive salinity trend was about 0.0014 psu/year. The ocean upperlayer salinity decreased with a rate of −0.0021 psu/year simultaneously with the increase in salinity of the intermediate layer. The decreased salinity of the upper layer and its increase in the intermediate layer result in the intensified halocline. The radius of the Aleutian eddies significantly increased (by 50–100%), which is equivalent, at least, to twofold increase in the volume of the water transported.     

Characteristics of microscale evapotranspiration: a comparative analysis

Summary ?Evapotranspiration characteristics on the point-scale (several hundred square meters) and the local scale (several square kilometers) are analysed by comparing a deterministic and a statistical – deterministic surface energy balance model. The vegetation surface variability is represented by both the surface heterogeneity and inhomogeneity. Heterogeneity means the mosaic of wet (wif ≠ 0) and dry (1-wif) fractions of vegetation surface, while inhomogeneity addresses small scale variations of soil moisture content. The microscale characteristics of evapotranspiration are considered in terms of analysing evapotranspiration E _v versus soil moisture content θ, relative frequency distribution characteristics of E _v (θ) and the aggregation algorithms for its estimation. The analyses are performed for loam soil type under different atmospheric forcing conditions. The main result is as follows: For dry vegetation surface (wif = 0), the relationship between the aggregated (θ_agg) and the area-averaged (θ _m ) soil moisture content is nonlinear and depends on both the states of the surface and the atmospheric forcing conditions. In the study, we assumed that there are no advective effects and mesoscale circulation patterns induced by surface discontinuities. Based on this fact it seems unlikely to be able to construct an aggregation algorithm for calculating θ_agg without inclusion of the atmospheric forcing conditions. This means that it will be difficult to construct a simple formula for calculating area-averaged transpiration, if it is possible at all. Received May 3, 2001; revised May 31, 2002; accepted June 3, 2002     

Flux-Variance Method for Latent Heat and Carbon Dioxide Fluxes in Unstable Conditions

Applied previously to momentum and heat fluxes, the present study extends the flux-variance method to latent heat and CO₂ fluxes in unstable conditions. Scalar similarity is also examined among temperature (θ), water vapour (q), and CO₂ (c). Temperature is adopted as the reference scalar, leading to two feasible strategies to estimate latent heat and CO₂ fluxes: the first one relies on flux-variance similarity relations for scalars, while the second is based on the parameterization of relative transport efficiency in terms of scalar correlation coefficient and a non-dimensional quantity. The relationship between the θ-to-q transport efficiency (λ _{θ q} ) and θ-q correlation coefficient (R _{θ q} ) is used to describe the intermediate hydrological conditions. We also parameterize the θ-to-c transport efficiency (λ _{θ c} ) as a function of the θ-c correlation coefficient (R _{θ c} ) by introducing a new non-dimensional ratio (α). The flux-variance method is a viable technique for flux gap-filling, when turbulence measurements of wind velocity are not available. It is worth noting that the extended method is not exempt from a correction for density effects when used for estimating water or carbon exchange.     

A Comparative Study of the Atmospheric Layers below First Lifting Condensation Level for Instantaneous Pre-Monsoon Thunderstorm Occurrence at Agartala (23o30’N, 91o15’E) and Ranchi (23o14’N, 85o14’E) of India

An attempt has been made to investigate the role of vertical wind shear, corrective instability and the thermodynamic parameter (θes - θe) below the first lifting condensation level (FLCL) in the occurrence of instanta-neous premonsoon thunderstorm over Agartala (AGT) and Ranchi (RNC) at 12 GMT Radiosonde data of 1988 have been utilized here. The study has however been confined to 1000 hPa-500 hPa range at most Here the convectively unstable layers with positive vertical wind shear upto 500 hPa have been termed as ‘Fa?vourable Layers’ (FL) and the level at which an initially stable layer turns out to be convectively unstable for the first time has been termed as ‘Transition Level’ (TL). It is observed that the changes in vertical wind shear are positive at TL at the time of occurrence of thunderstorm (TS) and the corresponding change is negative on fair-weather situa?tion Moreover, the 90% confidence interval for (θes - θe) reveals that for AGT the upper layer thermodynamic characteristic is important at the time of occurrence of TS whereas for RNC, the value of (θes - θe) at the surface is much more effective     

Quantitative reconstruction of the last interglacial vegetation and climate based on the pollen record from Lake Baikal, Russia

Changes in mean temperature of the coldest (T _c) and warmest month (T _w), annual precipitation (P _ann) and moisture index (α) were reconstructed from a continuous pollen record from Lake Baikal, Russia. The pollen sequence CON01-603-2 (53°57′N, 108°54′E) was recovered from a 386 m water depth in the Continent Ridge and dated to ca. 130–114.8 ky BP. This time interval covers the complete last interglacial (LI), corresponding to MIS 5e. Results of pollen analysis and pollen-based quantitative biome reconstruction show pronounced changes in the regional vegetation throughout the record. Shrubby tundra covered the area at the beginning of MIS 5e (ca. 130–128 ky), consistent with the end of the Middle Pleistocene glaciation. The late glacial climate was characterised by low winter and summer temperatures (T _c ~ −38 to −35°C and T _w~11–13°C) and low annual precipitation (P _ann~300 mm). However, the wide spread of tundra vegetation suggests rather moist environments associated with low temperatures and evaporation (reconstructed α~1). Tundra was replaced by boreal conifer forest (taiga) by ca. 128 ky BP, suggesting a transition to the interglacial. Taiga-dominant phase lasted until ca. 117.4 ky BP, e.g. about 10 ky. The most favourable climate conditions occurred during the first half of the LI. P _ann reached 500 mm soon after 128 ky BP. However, temperature changed more gradually. Maximum values of T _c ~ −20°C and T _w~16–17°C are reconstructed from about 126 ky BP. Conditions became gradually colder after ca. 121 ky BP. T _c dropped to ~ −27°C and T _w to ~15°C by 119.5 ky BP. The reconstructed increase in continentality was accompanied by a decrease in P _ann to ~400–420 mm. However, the climate was still humid enough (α~0.9) to support growth of boreal evergreen conifers. A sharp turn towards a dry climate is reconstructed after ca. 118 ky BP, causing retreat of forest and spread of cool grass-shrub communities. Cool steppe dominated the vegetation in the area between ca. 117.5 ky and 114.8 ky BP, suggesting the end of the interglacial and transition to the last glacial (MIS 5d). Shift to the new glaciation was characterised by cool and very dry conditions with T _c ~ −28 to −30°C, T _w~14–15°C, P _ann~250 mm and α~0.5.     

热带太平洋表面水中CO2对El Niño事件响应的数值模拟

文中用一个带生物泵的三维全球海洋碳循环模式模拟了热带太平洋表面水中CO2总量(TCO2)在ElNi