Skip to main content
Have a personal or library account? Click to login
Revised eddy covariance flux calculation methodologies – effect on urban energy balance Cover

Revised eddy covariance flux calculation methodologies – effect on urban energy balance

Open Access
|Jan 2012

Figures & Tables

Fig. 1. 

An overview of the surrounding area of the EC tower: white height contour lines (m) superimposed on a Google Earth picture. The locations of the SMEAR III tower and the auxiliary meteorological instruments are marked with a red star and circle, respectively. Black lines indicate the borders of the three sectors (urban 320–40°, road 40–180° and vegetation 180–320°).

Table 1. Calculation procedure comparison. The steps in eddy covariance raw data post-processing are shown in the columns, and the affected fluxes are in brackets: sensible heat (H) and open- and closed-path latent heat, (LE OP, LE CP). The names of the different calculation method combinations and their ID numbers are referred to in the text as rows. The calculation steps that are within the iteration are marked with a grey background. See text for further details

ID number Averaging period (min) Despiking and linear detrending Coordinate rotation (all) Relative to dry air (LECP) Lag time determination (all) Cross-wind correction (H) Iteration of fluxes (all) Spectral correction (all)a Sonic heat correction (H) Density correction (LEOP) Surface heating correction (LEOP) 1 Standard 30 x 2-D x x x x exp x x – 2 No spec. 30 x 2-D x x x x – x x – 3 Theor. spec. 30 x 2-D x x x x theor x x – 4 No iteration 30 x 2-D x x x – exp x x – 5 PFb 30 x PF x x x x exp x x – 6 SHC 30 x 2-D x x x x exp x x x 7 Theor. lagc 30 x 2-D x theor x x exp x x – 8 5 min 5 x 2-D x x x x exp x x – 9 10 min 10 x 2-D x x x x exp x x – 10 20 min 20 x 2-D x x x x exp x x – 11 60 min 60 x 2-D x x x x exp x x –

[i] aExperimental (exp) or theoretical (theor) spectral correction for LE CP .

[ii] bPF=planar fitting somewhat as in Wilczak et al. (2001).

[iii] cTheoretical lag window for LE CP as shown in Fig. 7.

Table 2. Flux omittance percentages for momentum (τ), sensible heat (H) and closed- and open-path latent heat (LE CP, LE OP). Values are given for omittance due to all criteria for all data and for day- and night-time data, and for omittance only due to flux stationarity and intermittency. Also, the percentage of unacceptable skewness and kurtosis of the zonal wind (u), meridional wind (v), vertical wind (w), temperature (T), H2O concentration from the open-path analyser (q OP) and closed-path analyser (q CP) are given. Values are calculated for raw 10 Hz data after despiking. See text for further details.

% τ HLEOPLECP All criteria 21.74 38.11 46.20 46.30 All criteria, day-time 19.64 40.05 48.69 57.05 All criteria, night-time 24.07 33.29 38.78 33.79 Flux stationarity 18.30 30.98 28.15 33.08 Flux intermittency 10.53 21.49 32.43 30.07 uv w TqOPqCP Skewness 0.17 0.30 0.03 0.88 9.91 2.49 Kurtosis 0.55 0.97 0.25 2.42 12.8 7.69
Fig. 2. 

(a) 5-d running means of air temperature (T, °C) and wind speed (U, ms−1) at 31 m height. (b) Daily precipitation (mm d−1) and 5-d running mean of relative humidity (RH,%).

Table 3. Flux error analysis for momentum (τ), sensible heat (H) and closed- and open-path latent heat (LE CP, LE OP). Mean random flux error (units of flux), the ratio between the random flux error and the observed flux (%), the detection limit (units of flux) and the ratio between the detection limit and the observed flux (%) for all quality checked data in 2009.

τ HLECPLEOP Detection limit 0.0414 kg m−1 s−2 5.02 W m−2 4.73 W m−2 5.20 W m−2 Detection limit/flux 9.95% 10.60% 10.69% 10.81% Random uncertainty 0.0464 kg m−1 s−2 6.39 W m−2 6.70 W m−2 7.63 W m−2 Random uncertainty/flux (%) 12.36% 9.05% 16.14% 15.75%
Fig. 3. 

(a) The ratio of the flux detection limit to the observed flux for momentum (τ), sensible heat (H) and closed- and open-path latent heat (LE CP, LE OP). Values are given for three wind direction sectors: urban (320–40°), road (40–180°) and vegetation (180–320°). Bars indicate means, and errorbars the 25th and 75th percentiles. (b) Same as (a) but for the ratio between the random flux error and the observed flux. Data are quality checked (as per text) and covers 2009.

Fig. 4. 

Time series of net radiation (R n ), sensible heat flux (H), latent heat flux measured with a closed- and open-path analyser (LE CP , LE OP), and residuals of the energy balance (Res CP , Res OP), which is defined Res=R n −H−LE. Grey dots indicate 30 min average fluxes, and the continuous and dashed lines are the 5-d running mean averages of day-time and night-time data, respectively. Turbulent flux data have been quality screened according to the text.

Fig. 5. 

Average diurnal cycles for 20 degree wind direction sectors. Hours are labelled with 2-hour resolution and run from 1:00 (local time) closest to the centre to 23:00 (local time) on the outside of the subplots. Colours indicate the median of the energy flux divided by the absolute value of net radiation (R n); see colourbar. The energy fluxes are (a) sensible heat flux (H), (b) closed-path latent heat flux (LE CP), (c) open-path latent heat flux (LE OP), (d) residual of the energy balance using LE CP and (e) residual of the energy balance using LE OP. All turbulent fluxes have been quality checked (see text), ratios exceeding 10 have been omitted from the plot, and data cover the whole period from October 2007 to September 2010 inclusive.

Table 4. Heat fluxes in cities with a cold temperate climate (latitude >45°N). The ranges of mean diurnal courses are given for the sensible (H) and latent (LE) heat fluxes, the energy balance residual (Res=R n −H−LE=ΔQ S −Q F) in addition to the anthropogenic heat flux (Q F).

Reference City Latitude (°N) Data period and season H (W m−2) LE (W m−2) Res (W m−2) QF (W m−2) Oke (1988) Moscow 55.750 1 yr – – – 127a Montreal 45.503 99 Vancouver 49.291 19 Fairbanks 64.849 6 Grimmond and Oke (1999) Vancouver, industr. 49.267 1 month fall, 3 months summer 69b 17b −200–400c – Vancouver, suburb. 49.250 84 31 −80–250 Walsh et al. (2004) Vancouver 49.226 1 yr 0–190 −5–50 – – Lemonsu et al. (2008) Montreal 45.503 4 weeks spring 0–200 0–50 −100–330 – Bergeron and Strachan (2010) Montreal, rural 45.547 2 winters −10–25 0–35 −30–80 – Montreal, suburb. 45.501 −5–130 0–35 −60–140 7–13d Montreal, urban 45.328 0–140 0–20 −50–145 25–45 Steinecke (1999) Reykjavik 64.134 1 yr – – – 35a Christen and Vogt, 2004 Basel, urban 1 47.556 1 month summer 10–250 0–90 –100–200 20e Basel, urban 2 47.552 0–250 0–100 −80–200 10 Basel, urban 3 47.555 −10–200 0–40 −80–130 20 Basel, suburban 47.552 −5–180 5–150 −100–180 5 Vogt et al. (2006) Basel 47.556 1 month summer −20–400c 0–100c – – Offerle et al. (2005) Łódź 51.767 2 yr −5–200 10–105 −100–180f −30–50f Mårtensson et al. (2006) Stockholm 59.303 49 d spring 27.7g – – – Martin et al. (2009) London 51.518 ~1 month campaigns, different seasons −20–125 – – – Gothenburg 57.700 −5–60 Edinburgh 55.950 −25–170 Manchester 53.483 −10–205 Helfter et al. (2011) London 51. 518 1 yr − 0–60h – – Iamarino et al. (2011) Greater London 51.500 4 yr – – – 11a Current study Helsinki 60.210 3 yr −17–125 11–95I 10–88 −62–89I –53–69 13e

[i] aAnnual means obtained by summing up contributions of human metabolism, vehicular traffic and buildings, using, e.g., energy consumption and population statistics.

[ii] bDaily means.

[iii] cIndividual values.

[iv] dSame as footnote ‘a’ but for the range of values in a diurnal course.

[v] eAnnual average estimated assuming that ΔQ S has to be zero over a year.

[vi] fΔQ S is obtained using the Element surface temperature method, where element surface temperatures are measured and the heat transfer through the elements are modelled. Q F=H+LE+ΔQ S −R n .

[vii] gMean of the whole period.

[viii] hRange is from autumn night to summer day.

[ix] iValues for open-path and closed-path gas analysers.

Fig. 6. 

Calculation procedure effects on turbulent flux magnitudes for momentum (τ), sensible heat (H) and closed- and open-path latent heat (LE CP, LE OP), in addition to effects on energy balance residuals (Res CP, Res OP). Bars give the means over all stabilities; and means for different stability classes are denoted by ∇ (stable, ζ>0.01), o (neutral, −0.01<ζ<0.01), Δ (unstable, ζ<−0.01). Errorbars denote 25th and 75th percentiles, and the median is denoted by a filled circle (•). Red crosses show fluxes that are not affected by a certain calculation procedure. Deviations are calculated as (flux-fluxStandard) except for Res where an absolute value of the denominator is used since corrections can easily change the sign of Res.

Fig. 7. 

Lag time (s) of closed-path H2O measurements as a function of relative humidity (%) and an exponential with 95% confidence intervals. The theoretical lag time (t lag,theor) and the theoretical lag window (rectangle) are also given. The colourbar gives the latent heat flux magnitude (LE, W m−2). Data are from May and June 2010.

Fig. 8. 

Normalised cross-correlation functions of closed-path H2O measurements as a function of time shift (s) for different relative humidities. The time lag has been taken into account and negative time shifts correspond to H2O signal lagging the vertical wind speed signal. Thick lines are means and shaded patches are standard deviations. Data are from September 2009.

Fig. 9. 

Flux convergence for momentum (τ), sensible heat (H) and closed- and open-path latent heat (LE CP, LE OP). Proportion of flux Ogives that converge within a 30 min averaging period (black) or have an extreme value before converging within 219 min (white). No quality-checked data fell into the class of non-convergent Ogives.

Fig. 10. 

(a) Monthly mean residual of the surface energy balance (W m−2, eq. 2) with closed-path (Res CP) and open-path (Res OP) analyser evaporation, (b) difference between a modified Res CP and the standard Res CP, (c) same for Res OP. Values in (b) and (c) are given for different calculation procedures: no spectral correction, theoretical spectral correction of LE CP, no iteration, planar fitting, surface heat correction of LE OP and theoretical lag time of LE CP. Only times when the data for all calculation types are available are included, and some months are thus omitted due to a lack of data. See text for further details.

Language: English
Page range: 18184 - 18184
Submitted on: Jul 19, 2011
Accepted on: Jan 27, 2012
Published on: Jan 1, 2012
Published by: Stockholm University Press
In partnership with: Paradigm Publishing Services

© 2012 Annika Nordbo, Leena Järvi, Timo Vesala, published by Stockholm University Press
This work is licensed under the Creative Commons Attribution 4.0 License.