稳态与湍流特征测试对通量数据质量的评价

应用CO2通量、感热和潜热通量稳态测试及垂直速度积分测试相结合的湍流数据质量评价方法对锡林浩特国家气候观象台2009年全年的湍流通量观测数据进行了质量评价。结果表明,锡林浩特通量观测数据质量总体较好;80％以上数据通过质量评价为较好数据;质量评价较好的数据分布呈白天大于夜间、夏季大于冬季的特点。经过对比分析,前期数据处理能剔除大部分质量不好的数据,但应用稳态测试与湍流积分测试筛选后,数据质量仍有明显提高,能量闭合率较前期数据分析前提高2％以上。     

基于涡度相关的春玉米逐日作物系数及蒸散模拟

作物系数是计算作物蒸散量的关键参数。利用2006—2008年和2011年辽宁锦州玉米农田生态系统的涡度相关、气象、作物发育期及叶面积指数观测数据,分析不受水分胁迫条件下玉米逐日作物系数特征及其与叶面积指数的关系。研究表明:作物系数与玉米农田实际蒸散均呈单峰型变化,约在7月末至8月初达到最大值 (玉米开花吐丝期)。在此基础上,建立了不受水分胁迫条件下玉米逐日作物系数与叶面积指数关系 (达到0.01显著性水平), 同时,采用积温表示的标准化生育期方法模拟相对叶面积指数,并建立了逐日作物系数与相对叶面积指数关系 (达到0.01显著性水平),解决了无叶面积观测地区玉米逐日实际蒸散量的计算。研究结果可为玉米农田用水管理以及灌溉措施的制定提供参考。     

LES Study of the Energy Imbalance Problem with Eddy Covariance Fluxes

The spatial representativeness of heat fluxes on the basis of single-tower measurements, and the mechanism of the so-called energy imbalance problem, are investigated through numerical experiments using large-eddy simulation (LES). LES experiments are done for the daytime atmospheric boundary layer heated over a flat surface, as a best-case scenario completely free of sensor errors and the uncertainties of field conditions. Imbalance is defined as the deviation of the `turbulent' heat flux at a grid point from the horizontally averaged `total' heat flux. Both the theoretical and numerical results of the present study suggest the limitation of single-tower measurements and the necessity of horizontally-distributed observation networks.The temporally averaged `turbulent' flux based on a point measurement systematically underestimates the `total' flux (negative imbalance). This is attributed to local advection effects caused by the existence of turbulent organized structures (TOS), whose time scale is much longer than that of thermal plumes. The temporal and spatial change of TOS patterns causes low-frequency trends in the velocity and temperature data resulting in large scatter of the flux estimates. The influences of geostrophic wind speed, averaging time, observation height, computational domain size and resolution on tower-measured fluxes are also discussed. Finally, it is suggested that a weak inhomogenity in surface heating may reduce the negative bias of flux estimates.     

A Re-examination of Density Effects in Eddy Covariance Measurements of CO2 Fluxes

Corrections of density effects resulting from air-parcel expansion/compression are important in interpreting eddy covariance fluxes of water vapor and CO2 when open-path systems are used. To account for these effects, mean vertical velocity and perturbation of the density of dry air are two critical parameters in treating those physical processes responsible for density variations. Based on various underlying assumptions, different studies have obtained different formulas for the mean vertical velocity and perturbation of the density of dry air, leading to a number of approaches to correct density effects. In this study, we re-examine physical processes related to different assumptions that are made to formulate the density effects. Specifically, we re-examine the assumptions of a zero dry air flux and a zero moist air flux in the surface layer, used for treating density variations, and their implications for correcting density effects. It is found that physical processes in relation to the assumption of a zero dry air flux account for the influence of dry air expansion/compression on density variations. Meanwhile, physical processes in relation to the assumption of a zero moist air flux account for the influence of moist air expansion/compression on density variations. In this study, we also re-examine mixing ratio issues. Our results indicate that the assumption of a zero dry air flux favors the use of the mixing ratio relative to dry air, while the assumption of a zero moist air flux favors the use of the mixing ratio relative to the total moist air. Additionally, we compare different formula for the mean vertical velocity, generated by air-parcel expansion/compression, and for density effect corrections using eddy covariance data measured over three boreal ecosystems.     

Multi-Scale Sensible Heat Fluxes in the Suburban Environment from Large-Aperture Scintillometry and Eddy Covariance

Sensible heat fluxes ( \(Q_{H}\) ) are determined using scintillometry and eddy covariance over a suburban area. Two large-aperture scintillometers provide spatially integrated fluxes across path lengths of 2.8 and 5.5 km over Swindon, UK. The shorter scintillometer path spans newly built residential areas and has an approximate source area of 2–4 \(\text {km}^{2}\) , whilst the long path extends from the rural outskirts to the town centre and has a source area of around 5–10 \(\text {km}^{2}\) . These large-scale heat fluxes are compared with local-scale eddy-covariance measurements. Clear seasonal trends are revealed by the long duration of this dataset and variability in monthly \(Q_{H}\) is related to the meteorological conditions. At shorter time scales the response of \(Q_{H}\) to solar radiation often gives rise to close agreement between the measurements, but during times of rapidly changing cloud cover spatial differences in the net radiation ( \(Q^{*}\) ) coincide with greater differences between heat fluxes. For clear days \(Q_{H}\) lags \(Q^{*}\) , thus the ratio of \(Q_{H}\) to \(Q^{*}\) increases throughout the day. In summer the observed energy partitioning is related to the vegetation fraction through use of a footprint model. The results demonstrate the value of scintillometry for integrating surface heterogeneity and offer improved understanding of the influence of anthropogenic materials on surface-atmosphere interactions.     

A Re-examination of Density Effects in Eddy Covariance Measurements of CO_2 Fluxes

Corrections of density effects resulting from air-parcel expansion/compression are important in interpreting eddy covariance fluxes of water vapor and CO2 when open-path systems are used. To account for these effects, mean vertical velocity and perturbation of the density of dry air are two critical parameters in treating those physical processes responsible for density variations. Based on various underlying assumptions, different studies have obtained different formulas for the mean vertical velocity and ...     

北京325米气象塔上CO2梯度观测数据质量控制与评价

近些年涡度相关系统在城市通量研究中得到了广泛应用,因城市下垫面的特殊性和复杂性,以及系统观测原理和仪器精度存在着局限性,使得结果存在10%以上偏差,故必须对原始观测数据进行质量控制和评价。本文对北京325 m气象塔七层高度上的CO₂(二氧化碳)、能量等通量进行了长期观测,研究了涡度相关技术在城市环境通量观测中的适用性,设计出了一套适合于复杂城市下垫面上的涡度相关系统资料前处理和质量控制方案,并对比了不同修正方案效果。分析结果表明,城市下垫面环境中计算CO₂通量的最优周期为30 min;二次坐标旋转法优于平面拟合法;频率响应修正后的感热通量提高5.21%,潜热通量和CO₂通量均提高9.42%。空气密度脉动修正Liu法(Liu,2005)优于WPL法(Webb et al.,1980);湍流谱在惯性副区满足-2/3次方定律,协方差谱满足-4/3次方定律。经过这套前处理和质量控制方案,原始数据中有79%能够用于基础研究,该质量控制与评价体系可为复杂城市下垫面通量研究提供参考。     

Comparison of Ozone Fluxes over a Maize Field Measured with Gradient Methods and the Eddy Covariance Technique

Ozone(O_3) fluxes were measured over a maize field using the eddy covariance(EC) technique and gradient methods.The main objective was to evaluate the performance of the gradient methods for measuring the O_3 flux by comparing them with the EC O_3 flux.In this study,turbulent exchange coefficients(K) calculated with three methods were compared.These methods were the aerodynamic gradient(AG) method(in which K is calculated by using wind speed and temperature gradients),the aerodynamic gradient combined with EC(AGEC) method,in which the friction velocity and other variables are based on EC measurements,and the modified Bowen ratio using the EC sensible heat flux and temperature gradient(MBR) method.Meanwhile,the effects of the measurement and calculation methods of the O_3 concentration gradient were analyzed.The results showed that:(1) on average,the transfer coefficient computed by the MBR method was 40% lower,and the coefficient determined with the AG method was 25% higher,than that determined with the AGEC method.(2) The gradient method's O_3 fluxes with the MBR,AGEC,and AG methods were 30.4% lower,11.7% higher,and 45.6% higher than the EC O_3 flux,respectively.(3) The effect of asynchronous O_3 concentration measurements on the O_3 gradient must be eliminated when using one analyzer to cyclically measure two-level O_3 concentrations.The accuracy of gradient methods for O_3 flux is related to the exchange coefficient calculation method,and its precision mainly depends on the quality of the O_3 gradient.     

Evaluation of Density Corrections to Methane Fluxes Measured by Open-Path Eddy Covariance over Contrasting Landscapes

Corrections accounting for air density fluctuations due to heat and water vapour fluxes must be applied to the measurement of eddy-covariance fluxes when using open-path sensors. Experimental tests and ecosystem observations have demonstrated the important role density corrections play in accurately quantifying carbon dioxide \((\hbox {CO}_{2})\) fluxes, but less attention has been paid to evaluating these corrections for methane \((\hbox {CH}_{4})\) fluxes. We measured \(\hbox {CH}_{4}\) fluxes with open-path sensors over a suite of sites with contrasting \(\hbox {CH}_{4}\) emissions and energy partitioning, including a pavement airfield, two negligible-flux ecosystems (drained alfalfa and pasture), and two high-flux ecosystems (flooded wetland and rice). We found that density corrections successfully re-zeroed fluxes in negligible-flux sites; however, slight overcorrection was observed above pavement. The primary impact of density corrections varied over negligible- and high-flux ecosystems. For negligible-flux sites, corrections led to greater than 100% adjustment in daily budgets, while these adjustments were only 3–10% in high-flux ecosystems. The primary impact to high-flux ecosystems was a change in flux diel patterns, which may affect the evaluation of relationships between biophysical drivers and fluxes if correction bias exists. Additionally, accounting for density effects to high-frequency \(\hbox {CH}_{4}\) fluctuations led to large differences in observed \(\hbox {CH}_{4}\) flux cospectra above negligible-flux sites, demonstrating that similar adjustments should be made before interpreting \(\hbox {CH}_{4}\) cospectra for comparable ecosystems. These results give us confidence in \(\hbox {CH}_{4}\) fluxes measured by open-path sensors, and demonstrate that density corrections play an important role in adjusting flux budgets and diel patterns across a range of ecosystems.     

仪器表面加热效应对临泽站开路涡动相关系统CO2通量的影响

LI-7500分析仪仪器表面加热效应对开路式涡动相关系统CO2通量观测结果影响显著,利用Burba校正方法对提高观测站CO2通量观测精度、年净生态系统碳交换量(NEE)估算、全球CO2交换量估算和气候变化模拟等具有重要意义.基于临泽站绿洲玉米农田的开路式涡动相关系统和小气候观测系统所获得的一年数据,利用Burba校正方法分别对LI-7500分析仪光路中的感热通量、大气潜热通量、CO2通量以及NEE的季节变化过程进行了修正.结果表明,LI-7500分析仪底部窗口和支杆部分热量交换是光路中感热通量校正量的主要贡献者,平均分别为6.81 W·m-2和2.68 W·m-2,其加热效应主要来源于太阳辐射和电子元件运行产生的热量;加热效应对潜热通量影响最小,平均校正量仅为0.24 W·m-2;Burba校正对CO2通量和NEE的季节变化影响显著,其平均校正量分别为19.14μg· CO2·m-2·s-1和313.21 mg·C·m-2·d-1,而且气温越低加热效应对通量的影响越显著;除生长季空气中水汽浓度显著高于非生长季而导致潜热校正量较大外,生长季其他各通量的校正量明显小于非生长季,生长季光路中感热、周围大气潜热、CO2通量和NEE日校正量分别为6.94 W·m-2、0.33 W·m-2、12.86 μg·CO2·m-2·s-1和161.58mg·C· r-2·d-1,分别是非生长季的60.4％,220.0％,50.6％和37.4％.若未进行仪器表面加热效应的Burba校正,在生长季和非生长季累计高估的绿洲玉米农田生态系统碳吸收量分别为25.85 g·C·m-2和88.47 g·C·m-2.     

涡动相关仪和大孔径闪烁仪观测通量的空间代表性

在对涡动相关仪和大孔径闪烁仪足迹模型进行敏感性分析的基础上,利用北京密云站2006年8月至2007年12月期间的观测资料,应用足迹模型对观测通量的空间代表性做了初步的分析.结果表明:涡动相关仪和大孔径闪烁仪的源区对风向、Monin-Obukhov长度,空气动力学粗糙度和观测高度/有效高度等因子比较敏感.密云站涡动相关仪的源区白天主要分布在仪器的西南与南面,而夜间则在东北与北面.大孔径闪烁仪的源区为西南一东北向分布.涡动相关仪各月源区形状不同,但大致分布在南北长1000 m、东两宽850 m的范嗣内,而LAS各月源区为西南一东北向分布,长约2060 m,最宽处约为620 m.对涡动相关仪通量有贡献的下垫面主要为园地(67%)和耕地(19%).其中园地的通量贡献比例在夏、秋季比较大,冬、春季稍小,而耕地则相反.大孔径闪烁仪的主要通量贡献源区为园地、耕地和居民地,通量贡献比例分别为49%,28%和11%.其中园地和耕地通量贡献率的变化趋势与涡动相关仪的观测结果一致,但没有涡动相关仪的变化明显.     

青藏高原东部草甸下垫面涡旋相关观测数据的质量控制及评价研究

目前利用涡旋相关仪观测地表通量的方法已被广泛采用, 但由于涡旋相关法是建立在方程假设简化的基础上, 这对观测数据的质量提出了一定的要求, 所以对观测数据的处理及通量结果的质量状况分析就显得尤为重要。本文以理塘站2006年7月中到8月中的数据为例, 对其进行质量控制, 并将质量控制后的通量计算结果与原始值进行比较。质量控制后感热通量和潜热通量总和有所增长, 动量通量值则受质量控制影响较大。质量评价的结果为:动量通量、感热通量、潜热通量中高质量数据所占比例为68.2%、60.6%、63.3%。     

Energy Balance Closure Using Eddy Covariance Above Two Different Land Surfaces and Implications for CO<Subscript>2</Subscript> Flux Measurements

Components of the surface energy balance of a mature boreal jack pine forest and a jack pine clearcut were analysed to determine the causes of the imbalance that is commonly observed in micrometeorological measurements. At the clearcut site (HJP02), a significant portion of the imbalance was caused by: (i) the overestimation of net radiation (R _n ) due to the inclusion of the tower in the field of view of the downward facing radiometers, and (ii) the underestimation of the latent heat flux (λE) due to the damping of high frequency fluctuations in the water vapour mixing ratio by the sample tube of the closed-path infrared gas analyzer. Loss of low-frequency covariance induced by insufficient averaging time as well as systematic advection of fluxes away from the eddy-covariance (EC) tower were discounted as significant issues. Spatial and temporal distributions of the total surface-layer heat flux (T), i.e. the sum of sensible heat flux (H) and λE, were well behaved and differences between the relative magnitudes of the turbulent fluxes for several investigated energy balance closure (C) classes were observed. Therefore, it can be assumed that micrometeorological processes that affected all turbulent fluxes similarly did not cause the variation in C. Turbulent fluxes measured at the clearcut site should not be forced to close the energy balance. However, at the mature forest site (OJP), loss of low-frequency covariance contributed significantly to the systematic imbalance when a 30-min averaging time was used, but the application of averaging times that were long enough to capture all of the low-frequency covariance was inadequate to resolve all of the high-frequency covariance. Although we found qualitative similarity between T and the net ecosystem exchange (NEE) of carbon dioxide (CO₂), forcing T to closure while retaining the Bowen ratio and applying the same factor to CO₂ fluxes (F _C ) cannot be generally recommended since it remains uncertain to what extent long wavelength contributions affect the relationship between T, F _C and C.     

基于两种潜在蒸散发算法的SPEI对中国干湿变化的分析

利用美国普林斯顿大学高分辨率的全球陆面同化数据集和美国国家环境预测中心的辐射再分析数据,根据Thornthwaite和Penman-Monteith公式分别计算了1948～2008年中国区域潜在蒸散发量;而后,使用降水和两套潜在蒸散发数据分别计算得到标准化降水蒸散发指数SPEI(Standardized Precipitation Evapotranspiration Index),并以此研究了1949～2008年中国区域干湿变化时空特征以及两种SPEI结果之间的差异;最后,给出了两种SPEI在中国的适用区域。结果表明:两种SPEI均显示中国地区整体上存在变干趋势,季节上以春季的变干趋势最为显著;空间上表现为以长江为界的南涝北旱,显著变干的区域有内蒙古中部、华北、东北以及四川东部地区,显著变湿的地区主要位于新疆北部和西部。同时,各种不同等级干旱也呈增加趋势,其中以中等干旱增加最为显著。1990年代中后期以来是中等和极端干旱发生最多的时期,空间上与SPEI显著减小的区域相对应。两种SPEI在 冬、春季差异最大,这主要是由于期间两种潜在蒸散发的计算结果之间存在很大差异。在Penman-Monteith公式中,由于空气动力项对冬、春季北方潜在蒸散发的贡献显著增加,基于该公式的SPEI相对而言能更合理地描述干湿变化特征。     

Oxygen Isotope Composition Of Human Teeth And The Record Of Climate Changes In France (Lorraine) During The Last 1700 Years

In order to study climate variations during the last 1700 years in eastern France, fifty-eight oxygen isotope compositions of phosphate were measured in human tooth enamel. The individuals, who lived in Lorraine, are assumed to have drunk local water derived directly from rainfall. According to previous work, drinking water is the main source of oxygen that sets the isotopic composition of phosphatic tissues in humans. The empirical fractionation equation determined from our data combined with those of Longinelli’s one [Geochim. Cosmochim. Acta 48 (1984) 385] was used to calculate the oxygen isotope composition of meteoric waters. The mean air temperature was inferred from these isotope ratios and the Von Grafenstein et al.’s [Geochim. Cosmochim. Acta 60 (1996) 4025] relationship between δ¹⁸O and air temperature. Oxygen isotope composition of present-day individuals yields a mean air temperature of 9.9± 1.7 ^∘C which is consistent with meteorological data. Application of this method to historical individuals results in mean air temperatures estimates 0 to 3 ^∘C higher than present-day air temperature. These warm air temperatures are not realistic during the so-called Little Ice Age for which an air-cooling of about 0.5 to 2 ^∘C has been documented. We propose that these relatively high δ¹⁸O values of human tooth enamel reflect higher mean δ¹⁸O values of meteoric water which can be attributed to an increased proportion of summer rainfall during the “Little Ice Age” time frame in Lorraine.     

An Analysis and Implications of Alternative Methods of Deriving the Density (WPL) Terms for Eddy Covariance Flux Measurements

We explore some of the underlying assumptions used to derive the density or WPL terms (Webb et al. (1980) Quart J Roy Meteorol Soc 106:85–100) required for estimating the surface exchange fluxes by eddy covariance. As part of this effort we recast the origin of the density terms as an assumption regarding the density fluctuations rather than as a (dry air) flux assumption. This new approach, which is similar to the expansion/compression approach of Liu (Boundary-Layer Meteorol 115:151–168, 2005), eliminates the dry-air mean advective vertical velocity from the development of the WPL terms and allows us to directly compare Liu’s assumptions for deriving the WPL terms with the analogous assumptions appropriate to the original expression of the WPL terms. We suggest, (i) that the main difference between these two approaches lies in the interpretation of the turbulent exchange flux, and (ii) that the original WPL formulation is the more appropriate approach. Given the importance of the WPL terms to accurate and reliable measurements of surface exchange fluxes, a careful analysis of their origins and their proper mathematical expression and interpretation is warranted.     

The Results of Precipitation Isotope Composition Monitoring in the Northern and Middle Urals

Russian Meteorology and Hydrology - The study presents the results of the monitoring of the isotope composition of precipitation which was performed using the GNIP (Global Network of Isotopes in...     

Evaluation of Cut-Off Frequency and Correction of Filter-Induced Phase Lag and Attenuation in Eddy Covariance Analysis of Turbulence Data

A recursive filter adopted for online eddy covariance analysis of turbulence data inthe atmospheric surface layer is revisited, and its properties and performanceare evaluated by means of concepts and methods developed in digital signal analysis.A rigorous estimate of effective cut-off frequency is derived along with an estimateof filter induced phase lag. Accordingly, filter design criteria are revised and variouseffects of parameter choice are assessed. Furthermore suitable corrections for compensating filter induced distortion, i.e., phase lag and attenuation, are proposed. The modified filter is tested on: (a) An artificial time series (random noise superimposed to `slow' sinusoidal signal); (b) turbulence data from field measurements. In case (a) the retrieval of input signal parameters is appreciably improved. In case (b) a better agreement with a Gaussian moving average is obtained, at lower computational cost.     

Estimating Global Distribution of Evapotranspiration and Water Balance Using Complementary Methods

ABSTRACT

The authors propose a modified complementary method to estimate regional evapotranspiration (ET) under different climatic and physical conditions using only meteorological data. The purpose of this study is to investigate the applicability of the modified complementary method for estimating global ET distribution and corresponding water balance. Gridded data from the Climate Research Unit, University of East Anglia, with 30 min spatial resolution and monthly time steps are used. Using the Thornthwaite water budget, monthly maps of global water surplus (precipitation minus ET) are produced. The results show good agreement with many previous studies. The average annual precipitation, ET, and water surplus are 690, 434, and 256?mm, respectively. The results show that the modified model can predict regional ET using meteorological data and can be used to assess global water resources. Consequently, the proposed method has strong potential for projecting water resource balance under future climate change.