WRF-Hybrid背景误差协方差调整在台风同化及预报中的应用研究

利用雷达径向风单点试验和2006年8月超强台风"桑美"个例,首先研究了静态背景误差协方差的尺度化因子和方差在台风系统雷达资料同化中对台风路径和强度预报的影响。结果表明:在高时空分辨率的雷达资料同化中,较小的尺度化因子能显著改进对台风路径的预报;尺度化因子的影响比方差的影响更为显著。基于上述实验结果,进一步研究了WRFDA-Hybrid系统中"流依赖"控制变量的水平局地化和垂直局地化对台风预报的影响。试验结果表明:当"流依赖"的水平局地化距离与静态背景误差协方差的尺度化因子具有等效影响范围的时候,WRFDA-Hybrid能够得到比较合理的分析结果。同时针对雷达观测资料的空间分布特征,本文提出了一种新的基于雷达探测高度的垂直局地化方案,对台风的强度和路径预报均有显著的改进。     

西北干旱区荒漠生态系统通量贡献区模型研究

涡度相关仪通量值所代表的通量贡献区范围,对于通量观测塔的选址、仪器安装高度的确定以及通量观测数据的质量控制等具有重要的指导意义。利用通量贡献区模型对位于古尔班通古特沙漠试验场通量观测资料的空间代表性进行初步分析。结果表明：该荒漠区在大气稳定条件下90%的通量贡献区最远可以达到686.40 m,通量贡献函数最大点的位置在162.50 m;大气稳定时各风向的通量贡献区范围在生长末期均达到最大,生长初期和中期的源区变化因受到各风向风速和植被下垫面的影响而有差异;大气不稳定时不同生长时期各风向通量贡献区没有固定变化规律;通量源区大约有58.71%的信息来自于荒漠区通量观测塔西南至西北方,整个生长季生长末期通量贡献最多,所占比例为40.16%。由FSAM模型测得的通量贡献区范围可以较准确地反映荒漠生态系统下垫面的通量信息。     

大孔径闪烁仪测量戈壁地区感热通量

利用2008年6月11～30日在金塔开展的"绿洲系统非均匀下垫面能量水分交换和边界层过程观测与理论研究"期间第一阶段戈壁下垫面大孔径闪烁仪(LAS)的观测资料,用混合对流方法和自由对流方法分别计算了戈壁感热通量。结果表明,对于利用LAS资料计算地表感热通量的方法中,混合对流方法相对于自由对流方法更加适用,且混合对流方法中Andreas给出的参数相对于DeBruin的参数更加适用于戈壁下垫面。此外,LAS测得的感热通量相对涡动相关方法的值较大,提高了地表能量闭合度。     

Study on the processing method of nighttime CO_2 eddy covariance flux data in ChinaFLUX

At present, using Eddy Covariance (EC) method to estimate the "true value" of carbon sequestration in terrestrial ecosystem arrests more attention. However, one issue is how to solve the uncertainty of observations (especially the nighttime CO2 flux data) appearing in post-processing CO2 flux data. The ratio of effective and reliable nighttime EC CO2 flux data to all nighttime data is relatively low (commonly, less than 50%) for all the long-term and continuous observation stations in the world. Thus, the processing method of nighttime CO2 flux data and its effect analysis on estimating CO2 flux annual sums are very important. In this paper, the authors analyze and discuss the reasons for underestimating nighttime CO2 flux using EC method, and introduce the general theory and method for processing nighttime CO2 flux data. By analyzing the relationship between nighttime CO2 flux and air fraction velocity u., we present an alternate method, Average Values Test (AVT), to determine the thresholds of fraction velocity (u.c) for screening the effective nighttime CO2 flux data. Meanwhile, taking the data observed in Yucheng and Changbai Mountains stations for an example, we analyze and discuss the effects of different methods or parameters on nighttime CO2 flux estimations. Finally, based on the data of part ChinaFLUX stations and related literatures, empirical models of nighttime respiration at different sites in ChinaFLUX are summarized.     

On ensemble representation of the observation-error covariance in the Ensemble Kalman Filter

Evensen (2003) presents a modification of the Ensemble Kalman Filter (EnKF), in which the observation-error and background-error covariance matrices are both represented by ensembles, in contrast to the usual practice, where only the background error is so represented. It is shown that this modification can cause the ensemble to collapse to a single member, in the common situation where the number of observations is more than twice the number of ensemble members, and to be rank-deficient when the number of observations is greater than or equal to the ensemble size. It is also shown that some further modifications to the scheme, presented by Evensen as offering numerical efficiencies, can prevent this collapse. However, these latter modifications are shown in some simple numerical examples to require tuning to produce acceptable results, which are nevertheless inferior to those of the standard EnKF.Acknowledgements The author acknowledges useful discussions with Peter Steinle, and other participants at the EnKF workshop held in BMRC in November, 2003.     

Growing season evapotranspiration in boreal fens in the Athabasca Oil Sands Region: Variability and environmental controls

Current efforts to assess changes to the wetland hydrology caused by growing anthropogenic pressures in the Athabasca Oil Sands Region (AOSR) require well-founded spatial and temporal estimates of actual evapotranspiration (ET), which is the dominant component of the water budget in this region. This study assessed growing season (May–September) and peak growing season (July) ET variability at a treed moderate-rich fen and treed poor fen (in 2013–2018), open poor fen (in 2011–2014), and saline fen (in 2015–2018) using eddy covariance technique and a set of complementary environmental data. Seasonal fluctuations in ET were positively related to net radiation, air temperature and vapour pressure deficit and followed trends typical for the Boreal Plains (BP) and AOSR with highest rates in June–July. However, no strong effect of water table position on ET was found. Strong surface control on ET is evident from lower ET values than potential evapotranspiration (PET); the lowest ET/PET was observed at saline fen, followed by open fen, moderately treed fen, and heavily treed fen, suggesting a strong influence of vegetation on water loss. In most years PET exceeded precipitation (P), and positive relations between P/PET and ET were observed with the highest July ET rates occurring under P/PET ~1. However, during months with P/PET > 1, increased P/PET was associated with decreased July ET. With respect to 30-year mean values of air temperature and P in the area, both dry and wet, cool and warm growing seasons (GS) were observed. No clear trends between ET values and GS wetness/coldness were found, but all wet GS were characterized by peak growing seasons with high daily ET variability.     

鼎湖山针阔叶混交林生态系统能量平衡分析

能量平衡分析作为评价涡度相关法通量观测数据可靠性方法，备受学界重视。应用OLS（Ordinary least squares）和EBR（Energy Balance Ratio）2种方法，系统分析了广东省鼎湖山针阔叶混交林生态系统能量平衡特点，并分析各种涡度通量修正方法对能量平衡的影响，结果表明鼎湖山通量站平均能力平衡不闭合度为33％-47％，略高于普遍报道的不闭合度范围（10％-30％）。WPL修正、μ＊订正和坐标转换，使得能量平衡闭合度有所提高，但夜间特别是冬季能量平衡较差问题依然没有得到根本解决，表明夜间弱湍流并不是导致夜间能量平衡闭合度差的主要原因。为客观评价本通量站以及ChinaFLUX能量平衡状况和通量数据质量，确定涡度相关法CO2通量数据分析方法和改进策略提供依据。     

VARIATIONAL DATA ASSIMILATION USING WAVELET BACKGROUND ERROR COVARIANCE: INITIALIZATION OF TYPHOON KAEMI (2006)

Background error covariance plays an important role in any variational data assimilation system, because it determines how information from observations is spread in model space and between different model variables. In this paper, the use of orthogonal wavelets in representation of background error covariance over a limited area is studied. Based on the WRF model and its 3D-VAR system, an algorithm using orthogonal wavelets to model background error covariance is developed. Because each wavelet function contains information on both position and scale, using a diagonal correlation matrix in wavelet space gives the possibility to represent some anisotropic and inhomogeneous characteristics of background error covariance. The experiments show that local correlation functions are better modeled than spectral methods. The formulation of wavelet background error covariance is tested with the typhoon Kaemi (2006). The results of experiments indicate that the subsequent forecasts of typhoon Kaemi’s track and intensity are significantly improved by the new method.     

Characteristics of secondary circulations over an inhomogeneous surface simulated with large-eddy simulation

Large-eddy simulation is used to study secondary circulations in the convective boundary layer modulated as a result of horizontally varying surface properties and surface heat fluxes over flat terrain. The presence of heat flux heterogeneity and its alignment with respect to geostrophic wind influences the formation, strength and orientation of organized thermals. Results show boundary-attached roll formation along heat flux maxima in the streamwise direction. The streamwise organization of the updrafts and downdrafts formed downwind of heterogeneities leads to counter-rotating secondary circulations in the crosswind plane. The distribution of resolved-scale pressure deviations shows large pressure gradients in the crosswind plane. Spanwise and vertical velocity variances and heat flux profiles depict considerable spatial variability compared to a homogeneous forest simulation. Secondary circulations are observed for various ambient wind scenarios parallel and perpendicular to heterogeneities. In the presence of increased wind speed, thermals emerging from the heat flux heterogeneity are elongated, and organize along and downwind of large-scale heterogeneity in the streamwise direction. Simulation with a reduced heat flux shows a shallower circulation with a lower aspect ratio. Point measurements of heat flux inside the roll circulation could be overestimated by up to 15–25% compared to a homogeneous case.     

协方差阵二次型容许估计的一个必要条件

研究协方差阵Σ的二次型容许估计问题。设 y1,y2 ,… ,yniid,n≥ 2 ,y1与 p维正态分布N (β,Σ )有相同的前四阶矩。其中β =(β1,β2 ,… ,βp)′∈ Rp与Σ =(σij) p× p >0均未知。记 y =△ (y1,y2 ,… ,yn)′。在二次损失 L (d ,Σ ) =tr(d -Σ) 2下给出Σ的二次型估计 a S2 + nby-y-′是容许估计的必要条件为 :(n - 1) a + b + 2 max(a,b)≤ 1。此必要条件比张立振等协方差阵的二次型容许估计中的必要条件有了明显的加强