## 1. Introduction

The scintillometry method is based on atmospheric refraction. For optical waves, the refraction is dominated by atmospheric temperature fluctuations, so one can confidently obtain the structure parameter of temperature

The derivation of

While scintillometry has been widely applied over grassland and cropland, and some of those studies have been over heterogeneous terrain (e.g., Schüttemeyer et al. 2006; Ezzahar et al. 2007; Ward et al. 2011), only few scintillometry studies have been performed over cities (Kanda et al. 2002; Lagouarde et al. 2006; Roth et al. 2006; Masson et al. 2008; Pauscher 2010; Salmond et al. 2012; Zieliński et al. 2012; Wood et al. 2013b). Cities cause difficulties in the interpretation of scintillometer data, given the possible heterogeneous surface below the path—in particular, the difficulty in determining effective height (including knowledge of the zero-plane displacement height). Furthermore, all studies over cities have been conducted in midlatitude cities. Therefore, our aim is to understand

In the future, if a model–measurement intercomparison is to be done, then one desires a purely measured quantity to compare with a purely modeled one. For the case of scintillometer-derived fluxes, one compares a flux derived with MOST using experimental values of a stability parameter with a flux derived with MOST using a modeled stability parameter. It would be desirable to compare measured

In this paper, we report (i) the first results from many months' data from two scintillometers and compare them with data from two sonic anemometers over Helsinki, Finland, and (ii) the evaluation of an algorithm for obtaining

## 2. Materials and methods

### a. Instrumentation and site description

Two large-aperture scintillometers (BLS900, Scintec AG, Germany) and two sonic anemometers (USA-1, Metek GmbH, Germany) have been part of a suite of equipment installed across the city of Helsinki, Finland, on the coast of the Gulf of Finland (Fig. 1; Table 1). They are a part of the Helsinki Urban Boundary-Layer Atmosphere Network (URBAN; http://urban.fmi.fi; Wood et al. 2013a), with increasing activity in observing the urban boundary layer in particular since 2004. The site is characterized by the vicinity of the sea and the strong seasonality in climate, caused by the high latitude (>60°N) and semicontinental climate. Downtown Helsinki is located in a peninsula protruding southward (Fig. 1a), including the sites of Torni, Sitra, and Elisa. Most of the land area on the map is 5–15 m above mean sea level (MSL), with some hills up to about 30 m.

Instrument positions (see also Fig. 1). Effective heights are as per Eq. (14) for scintillometers using data in Figs. 1b and 1d, and directionally averaged *z*–*z _{d}* for sonic anemometers.

It is desirable to have near-horizontal scintillometer paths, given that

The 10-Hz sonic-anemometer time series, of the three components of wind and sonic temperature (i.e., virtual temperature), were used to calculate turbulent fluxes: friction velocity

### b. Obtaining structure parameter from scintillometers

From 500-Hz light intensity measurements, 1-min averages of the structure parameter of the refractive index of air

^{−5}K hPa

^{−1}(Scintec 2011). Atmospheric pressure measurements from SMEAR III–Kumpula (HMP243, Vaisala Oyj, Vantaa, Finland) were used for both scintillometers. Air temperature measurements were from the Elisa mast (HMP45D, Vaisala Oyj) for the downtown scintillometer, and an average temperature of Elisa and SMEAR III–Kumpula (homemade platinum resistance thermometer Pt-100) masts for the city-scale scintillometer.

### c. Obtaining structure parameter from sonic anemometers

We estimate the structure parameter of temperature from time series of virtual temperature (i.e., uncorrected 10-Hz temperature data directly from the sonic anemometer) by fitting the parameters of model high-frequency spectra proposed by Kouznetsov and Kallistratova (2010). Because that reference is not permanent, we provide further details here.

*T*), one can write a power spectral density in the temporal domain as

_{υ}^{−1}) and

For practical calculations of the spectra, we used a one-size-fits-all approach: a fixed segment length of about 10 s (128 data points at 10 Hz), which was generally well within the inertial subrange of turbulence for the heights of several tens of meters. The calculations with the data used for this study have shown small dependence (<20%–30%) of the resulting structure parameters on the segment length within a range from 32 to 1024 points. Too-small segments lead to the decrease of the structure parameter accuracy at low values, when the spectrum is strongly affected by measurement uncertainties. This effect was most pronounced for low values of the temperature structure parameters. If segments are too long, then the signal from beyond the inertial subrange appears in the estimated spectrum. This results in a degradation of fitting performance of the model spectrum [Eq. (6)] and a corresponding increase in chi square. The choice of a segment length could be performed with some approximation of turbulent spectra (Kaimal et al. 1972); however, a simpler way is to choose a fixed segment length and then discard the resulting structure parameters that do not meet the quality criteria for the accuracy and goodness of fit.

Unlike the earlier methods (Greenhut and Mastrantonio 1989; Beyrich et al. 2005), the model spectrum [Eq. (6)] does not require any procedure to select the high-frequency limit of the inertial subrange in measured spectra, but it can use the spectra up to the Nyquist frequency. Attenuation of spectra, due to insufficient sensor response at high frequencies, was not observed; however, it can be incorporated into the model spectrum. Were this attenuation to become substantial, it would be detected with the chi-square test.

In this study, the sonic-derived structure parameter data were rejected primarily based on the chi-square test. The main reason for the test failure is too-low wind speed that invalidates the Taylor hypothesis for temperature fluctuations. The rejected fractions are 7.6% and 14.2% of all 30-min-average values for the SMEAR III–Kumpula and downtown (Torni) sites, respectively.

### d. Estimating sensible heat flux from structure parameter

^{−2}is acceleration due to gravity,

^{−3}) and

^{−1}kg

^{−1}). Different expressions for

*H*, for example, as is used in situations where EC data are not available (e.g., Hartogensis et al. 2003).

*u*is along-beam distance) to represent different conditions along the beam. However, in the calculation of

*H*, one single height is needed to be representative of the single value of

*H*for unstable stratification. The scintillometer effective height takes into account the spatial averaging and the stability dependence. The effective height as defined in the eddy-covariance community is simply the difference between the sonic-anemometer height and the zero-plane displacement height in the flux footprint.

Estimates of sensible heat flux are subject to some uncertainty. For greater certainty a source-area-model technique (Kormann and Meixner 2001) could be performed to estimate the zero-plane displacement height from morphological techniques (Grimmond and Oke 1999) for different upwind sectors and stabilities (and thus different

## 3. Results

### a. Structure parameter from the scintillometers

The variability of ^{−3} to 10^{−1} K^{2} m^{−2/3}) occurred during daytime unstable conditions. But also, there were a few high values during November, January, and February nights, consistent with the large occurrence of negative sensible heat fluxes over Helsinki. Indeed, the sonic-derived sensible heat fluxes were negative for 45% of the time in those particular studied months, partly driven by the long nights and snow cover (Wood et al. 2013a). The lowest

Comparing

### b. Comparison of structure parameter from sonic and scintillometer

The ratios

The ratio of

A comparison between sonic and scintillometer downtown (Figs. 5a,b) for many days of data shows a good agreement for *r* = 0.85). This corroborates the usability of the spectral method to calculate

Two case days demonstrate the time evolution of

### c. Relationship with sensible heat flux

It is not possible to estimate

The quality of *H* from the sonic than the

These sensible heat flux comparisons are only slightly worse than other scintillometer studies (Lagouarde et al. 2006; Roth et al. 2006; Zieliński et al. 2012)—that is, giving high correlation coefficients of 0.85–0.95, but also rms errors of 20–100 W m^{−2}. Note, however, that the absolute scatter in Figs. 5c–f is mostly caused by larger values of

## 4. Summary

Measurements were performed in the city of Helsinki using two large-aperture scintillometers (Scintec BLS900) and two sonic anemometers (Metek USA-1). One scintillometer has been installed in relatively homogeneous terrain downtown with a 1.8-km path; the other one has a longer city-scale path and more heterogeneity underneath. Sonics were installed at the end points of the city-scale scintillometer path (ideal points near the center of scintillometer beams were not possible because of practical constraints in finding available measurement locations in cities). The values of

A robust method to derive

The challenges in estimating sensible heat flux from ^{−2}.

The climate conditions in Helsinki represent a challenge for scintillometry, given the urban surface and the high-latitude location giving rise to negative sensible heat fluxes even above an urban surface. The values of sensible heat flux are typically small, which poses a problem to MOST-based methods. Moreover, shallow boundary layers of a few tens of meters that often occur over Helsinki can invalidate the surface layer scaling, which is a prerequisite of MOST. In such cases, a surface layer might not exist.

A challenge now presents itself: how do we make best use of

## Acknowledgments

This work has been supported by the EC FP7 ERC Grant 227915 “Atmospheric planetary boundary layers: Physics, modelling and role in earth system,” Academy of Finland (Projects 138328, 1118615, and ICOS-Finland 263149), and the Russian Foundation for Basic Research (Project 13-05-00846). Kari Riikonen, Erkki Siivola, Petri Keronen, and Sami Haapanala provided technical support. We are grateful to the reviewers for their valuable comments and to Timo Vesala, Sylvain Joffre, Ari Karppinen, Lukas Pauscher, Helen Ward, Oscar Hartogensis, Daniëlle van Dinther, and Sue Grimmond for the fruitful discussions. The results and conclusions in this study were made at specific locations and with specific equipment configurations. They should not be used to judge the general performance of instruments or particular manufacturers.

## APPENDIX

### Saturation of the Scintillometers

To quantify the effect of saturation (Kohsiek et al. 2006) for our scintillometers, the histograms of ^{−2} K^{2} m^{−2/3} for the longer path with a very pronounced transition and (ii) at about 3 × 10^{−2} K^{2} m^{−2/3} for the shorter path with a less-pronounced transition. The Scintec manual (Scintec 2011) states that the maximum measurable values of ^{2} m^{−2/3} for 2- and 5-km paths, respectively. The manual's limit values are well above the values for which we observed saturation in our instruments. Note that only the data reported valid by the Scintec SRun 1.09 software are used for the histograms. Such behavior is likely caused by the inability of the implemented correction (Clifford et al. 1974) to correct all the data affected by saturation. It is not clear if this is an implementation problem or a problem of the correction itself. Most values of ^{−2} K^{2} m^{−2/3}, and for them the downtown histograms coincide quite well. We thus consider those data reliable.

## REFERENCES

Andreas, E. L, 2012: Two experiments on using a scintillometer to infer the surface fluxes of momentum and sensible heat.

,*J. Appl. Meteor. Climatol.***51**, 1685–1701.Aubinet, M., Vesala T. , and Papale D. , Eds., 2012:

*Eddy Covariance: A Practical Guide to Measurement and Data Analysis*. Springer Atmospheric Sciences Series, Vol. 12, Springer, 460 pp.Beyrich, F., Kouznetsov R. D. , Leps J.-P. , Lüdi A. , Meijninger W. M. L. , and Weisensee U. , 2005: Structure parameters for temperature and humidity from simultaneous eddy-covariance and scintillometer measurements.

,*Meteor. Z.***14**, 641–649.Beyrich, F., and Coauthors, 2012: Towards a validation of scintillometer measurements: The LITFASS-2009 Experiment.

,*Bound.-Layer Meteor.***144**, 83–112.Chehbouni, A., 2000: Estimation of heat and momentum fluxes over complex terrain using a large aperture scintillometer.

,*Agric. For. Meteor.***105**, 215–226.Clifford, S. F., Ochs G. R. , and Lawrence R. S. , 1974: Saturation of optical scintillation by strong turbulence.

,*J. Opt. Soc. Amer.***64**, 148.Coulter, R. L., and Doran J. C. , 2002: Spatial and temporal occurrences of intermittent turbulence during CASES99.

,*Bound.-Layer Meteor.***105**, 329–349.de Bruin, H. A. R., Kohsiek W. , and van den Hurk B. J. J. M. , 1993: A verification of some methods to determine the fluxes of momentum, sensible heat, and water vapour using standard deviation and structure parameter of scalar meteorological quantities.

,*Bound.-Layer Meteor.***63**, 231–257.de Bruin, H. A. R., van den Hurk B. J. J. M. , and Kohsiek W. , 1995: The scintillation method tested over a dry vineyard area.

,*Bound.-Layer Meteor.***76**, 25–40.Evans, J. G., McNeil D. D. , Finch J. W. , Murray T. , Harding R. J. , Ward H. C. , and Verhoef A. , 2012: Determination of turbulent heat fluxes using a large aperture scintillometer over undulating mixed agricultural terrain.

,*Agric. For. Meteor.***166–167**, 221–233.Ezzahar, J., Chehbouni A. , Hoedjes J. C. B. , and Chehbouni A. , 2007: On the application of scintillometry over heterogeneous grids.

,*J. Hydrol.***334**, 493–501.Foken, T., and Wichura B. , 1996: Tools for quality assessment of surface-based flux measurements.

,*Agric. For. Meteor.***78**, 83–105.Greenhut, G. K., and Mastrantonio G. , 1989: Turbulence kinetic energy budget profiles derived from Doppler sodar measurements.

,*J. Appl. Meteor.***28**, 99–106.Grimmond, C. S. B., and Oke T. , 1999: Aerodynamic properties of urban areas derived from analysis of surface form.

,*J. Appl. Meteor.***38**, 1262–1292.Hartogensis, O. K., Watts C. J. , Rodriguez J.-C. , and de Bruin H. A. R. , 2003: Derivation of an effective height for scintillometers: La Poza experiment in northwest Mexico.

,*J. Hydrometeor.***4**, 915–928.Hill, R. J., 1992: Review of optical scintillation methods of measuring the refractive-index spectrum, inner scale and surface fluxes.

,*Waves Random Media***2**, 179–201.Hollinger, D. Y., and Richardson A. D. , 2005: Uncertainty in eddy covariance measurements and its application to physiological models.

,*Tree Physiol.***25**, 873–885.HSY, 2008: SeutuCD. Helsinki Region Environmental Services Authority, CD-ROM.

Järvi, L., and Coauthors, 2009: The urban measurement station SMEAR 3: Continuous monitoring of air pollution and surface–atmosphere interactions in Helsinki, Finland.

,*Boreal Environ. Res.***14**, 86–109.Kaimal, J. C., Izumi Y. , Wyngaard J. C. , and Cote R. , 1972: Spectral characteristics of surface-layer turbulence.

,*Quart. J. Roy. Meteor. Soc.***98**, 563–589.Kanda, M., Moriwaki R. , Roth M. , and Oke T. , 2002: Area-averaged sensible heat flux and a new method to determine zero-plane displacement length over an urban surface using scintillometry.

,*Bound.-Layer Meteor.***105**, 177–193.Kleissl, J., Gomez J. , Hong S.-H. , Hendrickx J. M. H. , Rahn T. , and Defoor W. , 2008: Large aperture scintillometer intercomparison study

,*Bound.-Layer Meteor.***128**, 133–150.Kohsiek, W., Meijninger W. M. L. , Debruin H. A. R. , and Beyrich F. , 2006: Saturation of the large aperture scintillometer.

,*Bound.-Layer Meteor.***121**, 111–126.Kormann, R., and Meixner F. X. , 2001: An analytical footprint model for non-neutral stratification.

,*Bound.-Layer Meteor.***99**, 207–224.Kouznetsov, R. D., and Kallistratova M. , 2010: Anisotropy of a small-scale turbulence in the atmospheric boundary layer and its effect on acoustic backscattering.

*Extended Abstracts, 15th Int. Symp. for the Advancement of Boundary Layer Remote Sensing,*Paris, France, IPSL, O-TUR/02. [Available online at http://www.isars2010.uvsq.fr/images/stories/Presentations/O_TUR_02_Kouznetsov.pdf.]Kukkonen, J., and Coauthors, 2005: Analysis and evaluation of selected local-scale PM10 air pollution episodes in four European cities: Helsinki, London, Milan and Oslo.

,*Atmos. Environ.***39**, 2759–2773.Lagouarde, J.-P., Irvine M. , Bonnefond J.-M. , Grimmond C. S. B. , Long N. , Oke T. R. , Salmond J. A. , and Offerle B. , 2006: Monitoring the sensible heat flux over urban areas using large aperture scintillometry: Case study of Marseille city during the ESCOMPTE Experiment.

,*Bound.-Layer Meteor.***118**, 449–476.Marple, S. L., Jr., 1987:

*Digital Spectral Analysis with Applications*. Prentice Hall, 492 pp.Masson, V., and Coauthors, 2008: The Canopy and Aerosol Particles Interactions in Toulouse Urban Layer (CAPITOUL) Experiment.

,*Meteor. Atmos. Phys.***102**, 135–157.Meijninger, W. M. L., Beyrich F. , Lüdi A. , Kohsiek W. , and Bruin H. A. R. D. , 2005: Scintillometer-based turbulent fluxes of sensible and latent heat over a heterogeneous land surface—A contribution to LITFASS-2003.

,*Bound.-Layer Meteor.***121**, 89–110.Moene, A. F., 2003: Effects of water vapour on the structure parameter of the refractive index for near-infrared radiation.

,*Bound.-Layer Meteor.***107**, 635–653.Moene, A. F., Meijninger W. M. L. , Hartogensis O. K. , Kohsiek W. , and de Bruin H. A. R. , 2004: A review of the relationships describing the signal of a large aperture scintillometer. University of Wageningen Internal Rep. 2004/2, 40 pp.

Moene, A. F., Hartogensis O. K. , and Beyrich F. , 2009: Developments in scintillometry.

,*Bull. Amer. Meteor. Soc.***90**, 694–698.Monin, A. S., and Yaglom A. M. , 2007:

*Statistical Fluid Mechanics: Mechanics of Turbulence.*Vol. 2. Dover Publications, 896 pp.National Land Survey of Finland, cited 2010: Lidar scanning database PaITuli. [Available online at www.csc.fi/paituli.]

Nordbo, A., Järvi L. , and Vesala T. , 2012: Revised eddy covariance flux calculation methodologies—Effect on urban energy balance.

,*Tellus***64B**, 18 184.Nordbo, A., Järvi L. , Haapanala S. , Moilanen J. , and Vesala T. , 2013: Intra-city variation in urban morphology and turbulence structure in Helsinki, Finland.

,*Bound.-Layer Meteor.***146,**469–496, doi:10.1007/s10546-012-9773-y.Pauscher, L., 2010: Scintillometer measurements above the urban area of London. M.S. thesis, Department of Micrometeorology, University of Bayreuth, 104 pp.

Press, W. H., Teukolsky S. A. , Vetterling W. T. , and Flannery B. P. , 2007:

*Numerical Recipes: The Art of Scientific Computing*. 3rd ed. Cambridge University Press, 1256 pp.Roth, M., Salmond J. A. , and Satyanarayana A. N. V. , 2006: Methodological considerations regarding the measurement of turbulent fluxes in the urban roughness sublayer: The role of scintillometery.

,*Bound.-Layer Meteor.***121**, 351–375.Salmond, J. A., Roth M. , Oke T. R. , Christen A. , and Voogt J. A. , 2012: Can surface-cover tiles be summed to give neighborhood fluxes in cities?

,*J. Appl. Meteor. Climatol.***51**, 133–149.Schotanus, P., Nieuwstadt F. T. M. , and Bruin H. A. R. , 1983: Temperature measurement with a sonic anemometer and its application to heat and moisture fluxes.

,*Bound.-Layer Meteor.***26**, 81–93.Schüttemeyer, D., Moene A. F. , Holtslag A. A. M. , Bruin H. A. R. , and van Giesen N. , 2006: Surface fluxes and characteristics of drying semi-arid terrain in West Africa.

,*Bound.-Layer Meteor.***118**, 583–612.Scintec, 2011: Scintec boundary layer scintillometer hardware manual. Version 2.07, Scintec, 67 pp.

Tatarskii, V., 1971:

*The Effects of the Turbulent Atmosphere on Wave Propagation.*Israel Program for Scientific Translations, 472 pp.Thiermann, V., and Grassl H. , 1992: The measurement of turbulent surface-layer fluxes by use of bichromatic scintillation.

,*Bound.-Layer Meteor.***58**, 367–389.Travouillon, T., Ashley M. C. B. , Burton M. G. , Storey J. W. V. , and Loewenstein R. F. , 2003: Atmospheric turbulence at the South Pole and its implications for astronomy.

,*Astron. Astrophys.***400**, 1163–1172.Tunick, A., 2005: Toward increasing the accuracy and realism of future optical turbulence calculations.

,*Meteor. Atmos. Phys.***90**, 159–164.Ward, H. C., Evans J. G. , and Grimmond C. S. B. , 2011: Effects of non-uniform crosswind fields on scintillometry measurements.

,*Bound.-Layer Meteor.***141**, 143–163.Wesely, M. L., 1976: The combined effect of temperature and humidity fluctuations on refractive index.

,*J. Appl. Meteor.***15**, 43–49.Wood, C. R., and Coauthors, 2013a: An overview on the Urban Boundary-Layer Atmosphere Network in Helsinki.

, in press.*Bull. Amer. Meteor. Soc.*Wood, C. R., Pauscher L. , Kotthaus S. , Barlow J. F. , Ward H. , Lane S. , and Grimmond S. , 2013b: Wind observations over the River Thames in central London.

,*Sci. Total Environ.***442**, 527–533.Zieliński, M., Fortuniak K. , Pawlak W. , and Siedlecki M. , 2012: Sensible heat flux in Łódź—Scintillometer and eddy covariance measurements.

*Extended Abstracts, Eighth Int. Conf. on Urban Climates,*Dublin, Ireland, IAUC, 110, 4 pp.