Effects of Coordinate Rotation on Turbulent Flux Measurements during Wintertime Haze Pollution in Beijing, China

GUO Xiao-Feng

State Key Laboratory of Atmospheric Boundary Layer Physics and Atmospheric Chemistry (LAPC), Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100083, China

Effects of Coordinate Rotation on Turbulent Flux Measurements during Wintertime Haze Pollution in Beijing, China

GUO Xiao-Feng

State Key Laboratory of Atmospheric Boundary Layer Physics and Atmospheric Chemistry (LAPC), Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100083, China

Eddy-covariance observations from the Beijing 325-m meteorological tower are used to evaluate the effects of coordinate rotation on the turbulent exchange of momentum and scalars during wintertime haze pollution (January-February 2013). Two techniques are used in the present evaluation; namely, the natural wind coordinate (NWC) and the planar fit coordinate (PFC), with the latter being applied by means of two methods for linear regression (i.e., overall and sector-wise). The different techniques show a general agreement in both turbulent fluxes and transport efficiencies, especially evident at the lower, 140-m level above the ground (compared to the higher, 280-m level), perhaps implying that the selection of a technique for coordinate rotation (NWC or PFC) is less of a concern for a sufficiently low level, despite the complexities of urban terrain. Additionally, sector-wise regression is a recommended approach for practical application of the PFC in a complex urban environment subjected to particulate pollution, because this method is found to produce a better correlation between the mean vertical velocity at the 140- and 280-m heights.

coordinate rotation, eddy-covariance method, particulate air pollution, turbulent exchange, urban environment

1　Introduction

Over recent years, Chinese cities have been afflicted with severe particulate air pollution, which frequently occurs in the form of consecutive episodes of haze within the urban boundary layer (UBL; Lin et al., 2013; Sun et al., 2013; Guo et al., 2014; Ji et al., 2014). To advance mechanistic understanding of persistent haze pollution in Beijing, a multitude of observational and modeling studies have been carried out to investigate various aspects of the UBL. For instance, poor dispersal capacity has been found to be characteristic of the UBL under stagnant meteorological conditions, particularly during prolonged occurrences of haze. Thermally stable stratification tends to prevail under hazy conditions, which not only suppresses the dispersion of air pollutants but also reduces turbulent exchange between the urban terrain and the at-mosphere.

The eddy-covariance (EC) method has gained popularity in UBL studies, thanks to its capability in producing direct measurements of both turbulent fluxes and related parameters. Nonetheless, in practice, its deployment in an urban setting is complicated, owing to a number of physical complexities associated with urbanized surfaces, such as tall roughness elements, temporal unsteadiness of airflows, and the heterogeneous distribution of scalar sources/sinks. EC turbulence data thus necessitate a suite of theoretical and empirical corrections, without which calculated fluxes are insufficiently accurate for studies of urban air pollution. Existing techniques of flux correction have, in fact, been established with little regard for all the complexities of the city environment, and so their application to urban areas likely involves inevitable uncertainties (e.g., Feigenwinter et al., 2012). Given the critical importance of correcting turbulent flux measurements over undulating, heterogeneous terrain (e.g., Lee, 1998; Baldocchi et al., 2000), it is necessary to scrutinize the applicability of relevant techniques by using observational data from urban environments.

The topic of the present paper is the calculation of EC fluxes during wintertime haze pollution in Beijing; specifically, the period 6 January to 28 February 2013 (see, for example, Sun et al. (2014) for in-depth investigations). The work presented involves evaluating the effects of coordinate rotation on both the turbulent fluxes and transport efficiencies with regard to momentum, heat, water vapor, and CO2. The conclusions reached herein apply primarily to conditions of haze pollution; a separate investigation is needed to clarify to what extent they apply to clean-air conditions.

2　Tower measurements and analysis

The present analysis relies on turbulence data from the Beijing 325-m meteorological tower, which is located in a compact building development (39°58＇N, 116°22＇E; 49 m above sea level). The buildings nearby vary considerably in height (mostly 20-50 m), forming a “compact midrise”settlement. During 6 January to 28 February 2013, two identical EC systems were in operation at the heights of 140 and 280 m, composed of three-dimensional sonic anemometers (CSAT3) (Campbell Scientific Inc., Logan, Utah, USA) and LI-7500 open-path gas analyzers (Li-Cor Inc., Lincoln, Nebraska, USA). These instruments measured three components of turbulent wind velocity, virtualtemperature, and water vapor and CO2concentrations; all turbulence data were collected at 10 Hz. A total of 1253 hourly measurements were available simultaneously at the two heights, and a subset of 1096 half-hourly segments passed Foken and Wichura's (1996) stationarity test, which were involved in the subsequent evaluations.

The following variables serve the analyses of turbulent exchange and transport efficiency at each height (140 and 280 m):

When applying the EC method, two techniques of coordinate rotation are adopted to evaluate their effects on turbulent fluxes and transport efficiencies, which are commonly known as the natural wind coordinate (NWC) and planar fit coordinate (PFC). Details of the algorithm can be found in, for example, McMillen (1988), Lee et al. (2004), and Wilczak et al. (2001). In brief, the NWC involves a two-angle (double) rotation, based on half-hourly mean velocity components measured in the sonic anemometer's instrument coordinate system; consequential averages of both the lateral and the vertical components are equal to zero (), so that the mean velocity vector is rotated into the longitudinal direction (＞0). The PFC also involves a procedure of double rotation, whose angles depend on a multiple linear regression (Wilczak et al., 2001):

where U, V, and W denote half-hourly mean velocity components measured in the instrument coordinate system, which are illustrated in Fig. 1 for each height (140 and 280 m); b0,1,2are regression coefficients used to derive two rotation angles. The PFC has a working assumption that the mean velocity vectors, namely (++UVWijk), are aligned within a plane that defines the local streamline coordinates; i, j, and k denote the unit vectors. After rotation, the ‘vertical' velocity (i.e., w) is rendered perpendicular to the plane described by Eq. (7). Given the complexity of urban airflows (as implied in Fig. 1), two methods are presently trialed for the regression of Eq. (7), viz. for each height:

1) Overall regression is made by using all of the mean velocity vectors (Ui, Vj, Wk) available in the dataset, viz. without regard to wind direction.

2) Sector-wise regression is made separately for each of the four wind sectors in the U-V plane, as indicated in Fig. 1 by Roman numerals (I, II, III, and IV).

Figure 2 clarifies that, as referenced to the unrotated measurements, adopting the PFC produces a better correlation between the mean vertical velocity () at the 140-and 280-m heights, especially evident by means of sector-wise regression of Eq. (7); similar effects were reported previously in Lee (1998), based on observations from a tall forest site. Moreover, Figs. 2b and 2c show that the two methods of regression lead to marked differences in, particularly at 280 m where values ofare mostly positive (52%) and negative (77%), due to overall and sector-wise regressions respectively. It remains to be seen whether their effects are still discernible in the flux measurements. For continued evaluations, the NWC serves as a benchmark, against which results derived from the PFC are compared.

3　Results

An appropriate technique of coordinate rotation ought not only to produce accurate fluxes but to give accurate estimates of transport efficiency. Figure 3 compares flux measurements (i.e.,*u, H, LE, and FCO2) derived from the PFC against those derived from the NWC; accordingly, Fig. 4 compares their respective transport efficiencies (i.e., Rwm, RwT, Rwq, and RwCO2). All comparisons are made separately for the 140- and 280-m heights. Table 1 gives several comparative statistics, including the mean absolute difference (MAD), mean bias estimate (MBE), root-mean-square error (RMSE), and mean absolute percentage difference (MAPD).

As Fig. 3 shows, there is an overall agreement in the turbulent fluxes (i.e., u*, H, LE, and FCO2) calculated using different techniques of coordinate rotation. Such agreement is particularly clear at 140 m, as confirmed by comparatively low values of MAD, RMSE, and MAPD (see Table 1); relative differences (MAPD) in the fluxes between the PFC and NWC are mostly within 10% at the 140-m height, while those at 280 m frequently exceed 15%. Interestingly, the two regression methods (i.e., overall and sector-wise) involved with the PFC lead to insignificant differences in the fluxes, despite yielding markedly different vertical velocities (i.e.,, as seen in Figs. 2b and 2c). Noticeable differences are occasionally visible only in theu*measurements at 280 m (see Fig. 3a). It is therefore suggested that, when applying the PFC under conditions of haze pollution, the selection of regression method is of a concern primarily for calculatingu*at a sufficiently high level in the UBL, where airflows could be increasingly subjected to directionally varied (sector-dependent) aerodynamic disturbance.

Figure 4 shows that vertical transport efficiencies of both momentum (i.e., Rwm) and scalars (i.e., RwT, Rwq, and RwCO2) are in good agreement between different techniques of coordinate rotation. At both the 140- and 280-m heights, results derived from the PFC lead, overall, to a negligible bias relative to those derived from the NWC (viz. the MBE values in Table 1 are almost zero, for instance). Associated values of MAPD, moreover, suggest that the selection of a technique for coordinate rotation (NWC or PFC) is less of a concern for calculating transport efficiencies at a lower level in a haze-polluted UBL (i.e., 140 m herein); nonetheless, their average effects approach or even exceed 20% at a higher level (i.e., 280 m herein). There are two potential explanations: (1) comparedto 140 m, flux ‘source areas' for 280 m have a larger spatial extent and probably a more heterogeneous distribution of the source/sink, rendering turbulence parameters more sensitive to the orientation of a coordinate system; and (2) turbulence measurements are demanding with respect to instrument set-up, whereas the observation data at 280 m might be more significantly influenced by mechanical shaking of this tall tower.

Table 1 Height-specific comparative statistics for flux measurements (u*, H, LE, and FCO2) and transport efficiencies (Rwm, RwT, Rwq, and RwCO2) between the PFC and the NWC. MAD, MBE, and RMSE have the units of m s-1for u*, W m-2for H and LE, and mg m-2s-1for FCO2; all are non-dimensional for Rwm, RwT, Rwq, and RwCO2.

4　Summary

Using UBL turbulence data of the Beijing 325-m meteorological tower, this paper provides a brief evaluation of the coordinate rotation effects on turbulent momentum and scalar exchange during haze pollution in the boreal winter of 2013. Candidate techniques include the NWC and PFC, with the latter being applied by means of two regression methods (i.e., overall and sector-wise). A general agreement is found between the different techniques, with regard to both turbulent fluxes and transport efficiencies; in particular, at a lower level above the ground (i.e., 140 m, as compared with 280 m). Such agreement appears interesting for an extremely complex urban environment. Regarding applications of the PFC, sector-wise regression of Eq. (7) is recommended so as to improve the correlation between the mean vertical velocity () at the 140- and 280-m heights. The two regression methods differ slightly inu*measured at 280 m; and, at either height, they do not yield noticeable differences in scalar fluxes of sensible heat, latent heat, and CO2.

Acknowledgements. This study was funded by the National Basic Research Program of China (Grant No. 2014CB447900), National High Technology Research and Development Program (Grant No. 2014AA06A512), and Ministry of Environmental Protection of China through its Special Funds for Scientific Research on Public Welfare (Grant No. 201409001).

References

Baldocchi, D., J. Finnigan, K. Wilson, et al., 2000: On measuring net ecosystem carbon exchange over tall vegetation on complex terrain, Bound.-Layer Meteor., 96, 257-291.

Detto, M., and G. G. Katul, 2007: Simplified expressions for adjusting higher-order turbulent statistics obtained from open path gas analyzers, Bound.-Layer Meteor., 122, 205-216.

Feigenwinter, C., R. Vogt, and A. Christen, 2012: Eddy covariance measurements over urban areas, in: Eddy Covariance—A Practical Guide to Measurement and Data Analysis, Springer, Heidelberg, 377-397.

Foken, T., and B. Wichura, 1996: Tools for quality assessment of surface-based flux measurements, Agric. For. Meteor., 78, 83-105.

Guo, X.-F., T. Yang, S.-G. Miao, et al., 2014: Urban boundary-layer stability and turbulent exchange during consecutive episodes of particle air pollution in Beijing, China, Atmos. Oceanic Sci. Lett., 7, 62-66.

Guo, X., H. Zhang, X. Cai, et al., 2009: Flux-variance method for latent heat and carbon dioxide fluxes in unstable conditions, Bound.-Layer Meteor., 131, 363-384.

Ji, D., L. Li, Y. Wang, et al., 2014: The heaviest particulate air pollution episodes occurred in northern China in January, 2013: Insights gained from observation, Atmos. Environ., 92, 546-556.

Lee, X., 1998: On micrometeorological observations of surface-air exchange over tall vegetation, Agric. For. Meteor., 91, 39-49.

Lee, X., J. Finnigan, and K. T. Paw U, 2004: Coordinate systems and flux bias error, in: Handbook of micrometeorology: A Guide for Surface Flux Measurements and Analysis, Kluwer Academic Publishers, Dordrecht, 33-66.

Lin, J.-T., D. Pan, and R.-X. Zhang, 2013: Trend and interannual variability of Chinese air pollution since 2000 in association with socioeconomic development: A brief overview, Atmos. Oceanic. Sci. Lett., 6, 84-89.

McMillen, R. T., 1988: An eddy correlation technique with extended applicability to non-simple terrain, Bound.-Layer Meteor., 43, 231-245.

Rajewski, D. A., E. S. Takle, J. K. Lundquist, et al., 2014: Changes in fluxes of heat, H2O, and CO2caused by a large wind farm, Agric. For. Meteor., 194, 175-187.

Sun, Y., Q. Jiang, Z. Wang, et al., 2014: Investigation of the sources and evolution processes of severe haze pollution in Beijing in January 2013, J. Geophys. Res. Atmos., 119, 4380-4398.

Sun, Y., T. Song, G. Tang, et al., 2013: The vertical distribution of PM2.5and boundary-layer structure during summer haze in Beijing, Atmos. Environ., 74, 413-421.

Wang, L., D. Li, Z. Gao, et al., 2014: Turbulent transport of momentum and scalars above an urban canopy, Bound.-Layer Meteor., 150, 485-511.

Webb, E. K., G. I. Pearman, and R. Leuning, 1980: Correction of flux measurements for density effects due to heat and water vapour transfer, Quart. J. Roy. Meteor. Soc., 106, 85-100.

Wilczak, J. M., S. P. Oncley, and S. A. Stage, 2001: Sonic anemometer tilt correction algorithms, Bound.-Layer Meteor., 99, 127-150.

Guo, X.-F., 2015: Effects of coordinate rotation on turbulent flux measurements during wintertime haze pollution in Beijing, China, Atmos. Oceanic. Sci. Lett., 8, 67-71,

10.3878/AOSL20140079.

27 September 2014; revised 17 November 2014; accepted 5 December 2014; published 16 March 2015

GUO Xiao-Feng, hsiaofengguo@gmail.com