Estimation of Turbulent Sensible Heat and Momentum Fluxes over a Heterogeneous Urban Area Using a Large Aperture Scintillometer

Sang-Hyun LEE,Jun-Ho LEE,and Bo-Young KIM

1Department of Atmospheric Science,Kongju National University,Gongju 314-701,South Korea

2Department of Optical Engineering,Kongju National University,Gongju 314-701,South Korea

Estimation of Turbulent Sensible Heat and Momentum Fluxes over a Heterogeneous Urban Area Using a Large Aperture Scintillometer

Sang-Hyun LEE?1,Jun-Ho LEE2,and Bo-Young KIM2

1Department of Atmospheric Science,Kongju National University,Gongju 314-701,South Korea

2Department of Optical Engineering,Kongju National University,Gongju 314-701,South Korea

The accurate determination of surface-layer turbulent fluxes over urban areas is critical to understanding urban boundary layer(UBL)evolution.In this study,a remote-sensing technique using a large aperture scintillometer(LAS)was investigated to estimate surface-layerturbulentfluxes overa highly heterogeneous urban area.The LAS system,with an opticalpath length of 2.1 km,was deployed in an urban area characterized by a complicated land-use mix(residential houses,water body,bare ground,etc.).The turbulent sensible heat(QH)and momentum fluxes(τ)were estimated from the scintillation measurements obtained from the LAS system during the cold season.Three-dimensional LAS footprint modeling was introduced to identify the source areas(“footprint”)of the estimated turbulent fluxes.The analysis results showed that the LAS-derived turbulent fluxes for the highly heterogeneous urban area revealed reasonable temporal variation during daytime on clear days,in comparison to the land-surface process-resolving numerical modeling.A series of sensitivity tests indicated that the overall uncertainty in the LAS-derived daytime QHwas within 20%–30%in terms of the influence of input parameters and the nondimensional similarity function for the temperature structure function parameter,while the estimation errors inτwere less sensitive to the factors of influence,except aerodynamic roughness length.The 3D LAS footprint modeling characterized the source areas of the LAS-derived turbulent fluxes in the heterogeneous urban area,revealing that the representative spatial scales of the LAS system deployed with the 2.1 km optical path distance ranged from 0.2 to 2 km2(a“micro-αscale”), depending on local meteorological conditions.

heterogeneous surface,large aperture scintillometer,mesoscale modeling,surface energy balance,turbulent exchange

1.Introduction

Turbulent transfer of heat,moisture,and momentum continuously occurs at the interface between the atmosphere and underlying Earth surface.The atmospheric boundary layer (ABL)is formed as a consequence of the interactions on time scales of a few hours or less(Garratt,1992;Foken,2008). Therefore,accurate determination of the turbulent fluxes in the surface layer is of primary importance to understanding the spatial and temporal variability of atmospheric flow and thermal structure through the ABL,which sequentially influences the formation of the characteristic local climate and meteorology,as well as air quality.

Generally,eddy-covariance(EC)systems and scintillometers have been widely used to determine surface-layer turbulent fluxes.The EC system measures surface-layer turbulentfluxeson the orderofhundredsofmetersorless,which has been successfully used to provide the local-scale turbulent fluxes for various natural and urban surfaces(e.g.,Offerle et al.,2006).The EC system directly measures atmospheric turbulent fluctuations,and thus the EC technique has great merit in producing accurate turbulent fluxes for relatively homogeneous flat surfaces.However,it is difficult to determine the sensor’s location formeasuring turbulentfluxes over highly heterogeneous surfaces,due to the uncertainty of surface source areas of turbulent fluxes observed at a particular location.Unlike the EC technique,the scintillometry technique uses Monin–Obukhov similarity theory to estimate area-averaged surface-layer turbulent fluxes.Thus,the turbulent fluxes can be determined at various spatial scales ranging from 100 m to a few kilometers,depending on the scintillometer type and optical path distance.For example,a small aperture scintillometer(SAS)uses a short optical path length (a few hundred meters),while a large aperture scintillometer (LAS)has a longer path length of up to 5 km.

Many field experiments based on the scintillometry technique have been conducted to determine surface-layer turbulent fluxes over natural surfaces(e.g.,De Bruin et al.,1995, 2002;Chehbouni et al.,2000;Beyrich et al.,2002;Zeweldi et al.,2010;Evans et al.,2012)and urban areas(e.g.,Kanda et al.,2002;Lagouarde et al.,2006;Pauscher,2010;Ward et al.,2014),and showed the usefulness of LAS-derived turbulent fluxes through direct comparison with EC measurements(e.g.,McAneney et al.,1995;Hartogensis et al.,2003; Lagouarde et al.,2006;Zeweldi et al.,2010;Liu et al.,2011, 2013;Xu et al.,2013;Ward et al.,2014).For example, Lagouarde et al.(2006)showed the LAS-derived turbulent fluxes with a path length of 1.8 km compared well to the EC fluxes over the city of Marseille,France.Ward et al.(2014) compared two LAS systems with different path lengths of 2.8 and 5.5 km and an EC instrument for surface-layer turbulent flux estimation over a residential area in Swindon,UK.These studies demonstrated the applicability of the LAS system over homogeneous residential areas in determining surfacelayer turbulent fluxes in the urban environment.

Most urban areas at spatial scales of a few km or less are,morphologically,highly heterogeneous due to the diversity of buildings,natural/artificial land-use mix,and topography.Even though previous studies have successfully determined the surface-layer turbulent fluxes over homogeneous urban areas using the EC system and/or LAS system(e.g., Lagouarde et al.,2006;Ward et al.,2014),there are few studies concerning turbulent flux estimation over heterogeneous urban surfaces.As discussed,the EC technique is difficult to apply for highly heterogeneous surfaces,whereas the scintillation method using the LAS system can provide surfacelayer turbulent fluxes integrated over the LAS-path distance, suffering less due to the surface heterogeneity than the EC system.In the present study,the aim was to investigate the ability of the LAS technique to estimate surface-layer turbulent fluxes over an urban area with a significant mix of different land-use types.First,surface-layer turbulent sensible heat and momentum fluxes in a highly heterogeneous urban area were estimated using the LAS system.Then,for verification of the LAS-derived turbulent fluxes,the results were compared to those of a land-surface process-resolving numerical model and radiative flux measurements.Three-dimensional LAS footprint modeling was also performed to determine the source areas of the estimated turbulent fluxes over the heterogeneous urban area.

The remainder of the paper is structured as follows:Section 2 briefly describes the scintillation methodology for estimating surface-layer turbulent fluxes and the 3D LAS footprint modeling.The experimental site,LAS instrumentation, and land-surface process-resolving modeling are described in section 3.Section 4 presents the LAS-derived turbulent fluxes and footprints over the heterogeneous urban area,and discusses the uncertainty of the LAS-derived turbulent fluxes through a series of sensitivity tests.A summary and conclusions follow in section 5.

2.Methodology

2.1.Scintillation method

When electromagnetic radiation propagates through the atmosphere,its intensity is affected by the atmospheric turbulence,leading to fluctuations in the refractive index of the air (e.g.,Frehlich and Ochs,1990).The path-integrated structure parameter of the refractive index〈C2n〉can be given(Wang et al.,1978)by

where r(=x/LLAS)is the normalized distance along the LAS path(LLAS)between transmitter and receiver,and C2n(r)is the structure parameter at distance r.Here,x and LLASindicate distance from the transmitter and optical path length, respectively.W(r)is the LAS weighting function,which has a mathematical bell-shape curve with a maximum value in the middle of the optical path(e.g.,Chehbouni et al.,2000; Meijninger and De Bruin,2000;Timmermans et al.,2009; Geli et al.,2012):

where K(=2π/λ)is the optical wavenumber,k is the turbulent spatial wavenumber,?n(k)(=0.033k-11/3)is the 3D Kolmogorov spectrum of the refractive index,J1is a Bessel function of the first kind for order one,with x1=K Dr/2 and x2=[K D(1-r)]/2,where D is the aperture diameter of the LAS transmitter/receiver.

whereσ2lnIis the statistical variance of the logarithm of intensity fluctuation obtained from measured radiation intensity. Hereafter,〈C2n〉is denoted as C2nfor simplicity.

Because the intensity fluctuation changes with air temperature and humidity,C2ncan be expressed as a combination of the structure parameters of temperature(C2T),humidity(C2Q), and their covariance fluctuations(C2TQ).The influence by atmospheric pressure fluctuation can be negligible(Hill et al., 1980).The scintillation is predominantly influenced by temperature fluctuation forvisible and infrared radiation(Wesely, 1976;Moene,2003),leading to a simplified relation between C2n(m-2/3)and C2T(K2m-2/3)as follows:

where T is the air temperature(K),P is the atmosphericpressure(Pa),andβis the Bowen ratio,which represents a humidity correction in scintillation.

where zeffis the LAS effective height(m),d is the zero-plane displacement height(m),and fTis a dimensionless universal function for the temperature structure function parameter. The L is the Obukhov length,defined by

where k(=0.4)is the von Karman constant and g(=9.8 m s-2)is the gravitational acceleration.The non-dimensional universal function for the temperature structure function parameter can be formulated for unstable atmospheric conditions following Andreas(1988):

The friction velocity(u?)is also estimated by the following similarity relation(e.g.,Panofsky and Dutton,1984;Foken, 2008):

where U is the mean wind speed at zeff,andψmis the integrated universal function for momentum flux(e.g.,Paulson, 1970;Businger et al.,1971).The integrated universal function used in this study for the unstable atmospheric surface layer is(H¨ogstr¨om,1988)

Given the similarity relations,the temperature scale(T?), friction velocity(u?),and the Obukhov length(L)can be solved iteratively.The turbulent sensible heat flux(QH; W m-2)and momentum flux(τ;kg m-1s-2)can be calculated as follows:

whereρis the air density(kg m-3)and cpis the specific heat of the air at a constant pressure(J kg-1K-1).

Additional parameters required for estimation of the turbulent fluxes are the LAS effective height(zeff),surface aerodynamic roughness length(z0),and zero-plane displacement height(d).First,the LAS effective height can be estimated iteratively from the following relation(Hartogensis et al., 2003):

where z(r)is the optical beam height at distance r.This relation indicates that zeffis dependent on atmospheric stability, and thus it was estimated at every time increment(30 min)of the available data in this study.

2.2.Three-dimensional LAS footprint modeling

Many flux footprint models have been suggested based on analytical and numerical approaches and applied in the interpretation of EC flux measurements(e.g.,Schuepp et al., 1990;Horst and Weil,1992,1994;Schmid,1994;Horst, 1999;Hsieh et al.,2000;Kormann and Meixner,2001;Kljun et al.,2004).The 3D LAS footprint can be calculated as an extension of the flux footprint model;that is,a combination of the flux footprint calculated at each segment location along the optical path and the LAS spatial weighting function(e.g., Hoedjes et al.,2002;Meijninger et al.,2002;Hartogensis et al.,2003).The turbulent flux footprint relation used for the LAS footprint modeling in this study is briefly introduced below.

The turbulent flux F(x,y,zm)measured at height zmcan be related to the footprint function f(x,y,zm)and the spatial distribution of the surface fluxes F0(x′,y′,z′=0)as follows (Horst and Weil,1992):

Here,the flux footprint function(units:m-2)is determined by the product of the crosswind distribution function Dy(x,y) (units:m-1)and the crosswind-integrated footprint function fy(x,zm)(units:m-1)(van Ulden,1978;Horst and Weil, 1992)as

where x is the downwind distance from the measurement location,and y is the crosswind distance from the centerline. The dispersion to the crosswind direction is assumed to be symmetric using a Gaussian distribution function(Pasquill, 1974):

whereσyis the standard deviation of the plume in the y dimension,which relies on atmospheric turbulence intensity, upwind distance(x),and mean wind speed[u(x)]asσvx/u(x) for a short-range dispersion(Schmid,1994).Here,the standard deviation of lateral wind fluctuations,σv(=1.9u?),was calculated following Panofsky and Dutton(1984).On the other hand,the footprint model by Hsieh et al.(2000)was used to calculate the crosswind-integrated footprint:

where D and P are specific constants that depend on atmospheric stability,and zuis a length scale defined as

The constant parameters are D=0.28 and P=0.59 for unstable conditions(Hsieh et al.,2000).

3.Measurements and land-surface processresolving numerical modeling

3.1.Site description

The field experiment was conducted in a small urban area in Gongju,South Korea,where there is a population of 117 000 people.The city is located in a basin surrounded by extensive agricultural and hilly regions.It is also located inland, about 50 km west of the Yellow Sea.Figure 1a shows the experimental site within the city.The urbanized area is largely composed of single(1–3 story)residential houses,and a river (named“Geum River”)crosses over the city.Relatively new buildings are located alongside the northern stretches of the river,showing the directions of the city’s sprawl.The vegetation fraction within the built-up region ranges from 5% to 10%,and hilly mountains surround the city in the south and east directions.The measurement site was highly heterogeneous in terms of its land cover and topography,which is a typical morphological configuration that can be found in many cities of South Korea.

3.2.LAS instrumentation

The LAS instruments were deployed at two buildings of the Kongju National University,located within the city,to investigate the ability in estimating turbulent sensible heat and momentum fluxes over the heterogeneous urban area(Fig. 1a).The receiver(marked A in Fig.1a)and transmitter (marked B in Fig.1a)were installed on the rooftops of buildings that are 49 m and 48 m above mean sea level,and the optical path extends 2.1 km from southwest to northeast,crossing the Geum River(A–B in Fig.1a).This configuration of the LAS instruments measured area-averaged turbulent sensible heat and momentum fluxes that were representative of the heterogeneous urban area(residential area,bare ground, vegetated area,water body,etc.).

The LAS system used was the LAS MkII manufactured by Kipp&Zonen,with an aperture diameter of D=0.149 m.The light source of the LAS MkII transmitter operates at a near-infrared wavelength of 850 nm,at which the scintillations measured by C2nare primarily influenced by turbulent temperature fluctuations.Based on Kleissl et al.(2010), the threshold C2nvalue for signal saturation was estimated as 4.06×10-14m-2/3for our LAS experimental configuration. The LAS sampled C2nat a 1 Hz frequency,which were averaged to 15 s intervals for turbulent flux estimation.The air temperature,wind speed,and atmospheric pressure were measured at the same temporal frequency from separate me-teorological sensors,but which were connected to the LAS system at the same height of the receiver.In addition,a 7 m meteorological tower on the rooftop of an eight-story building was also operated at the location of 320 m north of the receiver along the optical path.The meteorological data measured from the meteorological tower were used for the calculation of the LAS-derived turbulent fluxes.

3.3.Case selection and surface input parameters

The field measurements were conducted for 29 days during 5–19 November and 5–18 December 2013,in which adverse meteorological conditions(e.g.,fog)affected the measurements for many days.In this study,two sequential clear days on 11 and 12 November 2013 were selected to investigate the capability of the LAS in estimating the turbulent fluxes over the heterogeneous surfaces.Furthermore,since the LAS cannot discriminate atmospheric stability from the measured C2nvalues,the daytime period(0800–1700 LST) data were analyzed for estimation of the LAS-derived turbulent fluxes in unstable atmospheric conditions.

Figure 1b shows the topographical variation along the LAS beam path and the optical beam height relative to the ground.From Eq.(11),the LAS effective heights(zeff)were estimated as 25.1–27.2 m during the daytime periods of 11 and 12 November,which satisfied a“saturation condition”required for reasonable estimation of turbulent fluxes using the LAS system(Kipp&Zonen,2012).The Bowen ratio (β)was initially assigned as 1.The effective aerodynamic roughness length and displacement height were assigned as z0=0.3 m and d=2.1 m,respectively.The aerodynamic parameters were determined following the general rule of z0≈0.1h and d≈0.7h(Garratt,1992;Grimmond and Oke, 1999),in which h(～3 m)indicates the mean height of obstacles along the beam path.Generally,the determination of surface aerodynamic parameters can be obtained by applying a morphological method(e.g.,Raupach,1994;Macdonald et al.,1998)and/or surface-layer similarity relation using EC measurements taken in the atmospheric surface layer for homogeneous surfaces(e.g.,Grimmond and Oke,1999;Roth, 2000).However,accurate estimation of the surface aerodynamic parameters is a very difficult task for heterogeneous surfaces due to intrinsic limitations of the methodologies.No standard method that determines the surface aerodynamic parameters for heterogeneous surfaces exists,which inevitably leads to uncertainty in the LAS-derived surface-layer turbulent fluxes.In this study,a series of sensitivity tests were performed to quantify the influence of errors in the surface input parameters on the LAS-derived turbulent fluxes.

3.4.Land-surface process-resolving numerical modeling

The ability of the LAS equipment used in this study in estimating surface-layer turbulent fluxes has been validated against the EC system for homogeneous land surfaces (e.g.,Wilson et al.,2013).Due to the difference of the footprints of the EC and LAS-derived turbulent fluxes,direct comparison of the two turbulent fluxes is difficult for highly heterogeneous surfaces.In order to verify the LAS-derived turbulent fluxes over the heterogeneous field experimental area,land-surface process-resolving modeling was used in the present study.Specifically,we conducted very high resolution urban numerical modeling using the Weather Research and Forecasting(WRF)model(Skamarock et al., 2008).Given a proper experimental configuration,the WRF model has the capability to explicitly represent significant surface land-use heterogeneity over the field experiment region and to resolve associated land-surface physical(radiative/dynamic/hydrological)processes.

The WRF model was configured with five nested domains,covering East Asia with a 16.2 km grid spacing in domain 1 and nesting down to the Gongju region with a 200 m grid spacing in domain 5(Fig.2).To realistically represent the significant heterogeneity of land-use and topography over the experimental region,very high resolution static datasets of the Moderate Resolution Imaging Spectroradiometer (MODIS)-based land-use(Kang et al.,2010;spatial resolution of 7.5′′)and Shuttle Radar Topography Mission(SRTM) topography(http://www2.jpl.nasa.gov/srtm/index.html;spatial resolution of 3′′)were implemented as input data of the WRF model.Physical processes applied in the simulation included:Goddard shortwave radiation(Chou and Suarez, 1994),Rapid Radiative Transfer Model(RRTM)longwave radiation(Mlawer et al.,1997),Lin cloud microphysics(Lin et al.,1983),and a turbulent kinetic energy(TKE)–based turbulence scheme(Janji′c,2002).Land-surface processes were parameterized following the Noah land surface model(Chen and Dudhia,2001)with a single-layer urban canopy model (Kusaka et al.,2001)for urban patches.

The WRF simulation was conducted starting from 0000 UTC 9 November 2013 to 0000 UTC 13 November 2013. In order to realistically reproduce surface-layer meteorology and relevant turbulent fluxes,observed wind fields from the meteorological site were nudged during the simulation period.The WRF model calculates surface-layer turbulent sensible heat and momentum fluxes at every grid mesh.The WRF-simulated turbulent fluxes with a 30 min time interval were extracted along the LAS beam path,which were multiplied by the LAS weighting function for comparison with the LAS-derived fluxes.In this study,a WRF model sensitivity simulation was also conducted,by replacing the Geum River with a natural surface type(“shrubland”)to further interpret the LAS-derived turbulent fluxes(discussed in section 4.2).

4.Results

4.1.Meteorological conditions and LAS footprints

The weather conditions on 11 and 12 November 2013 were fairly good,with little influence by fog and clouds.As shown in Fig.3,the prevailing winds were northeasterly or northwesterly during the daytime on both days,except for the morning hours on 12 November.The wind speed changed within the range of0.5–4 m s-1during the daytime,being relatively low in the morning hours(0.5–2 m s-1)and relatively high in the afternoon(2–4 m s-1),while the air temperatureshowed a typical temporal variation for both the days,ranging from 0?C to 9?C.

Based on the meteorological conditions(Table 1),3D LAS footprint modeling was conducted to interpret the LAS-derived turbulentfluxes.In this study,the flux footprints were calculated every 10 m in discrete locations along the LAS optical path,and then the estimates were multiplied by the LAS weighting function to model the 3D LAS footprints.As shown in Fig.4,the calculated LAS footprints clearly revealed the source regions of the LAS-derived turbulent fluxes over the heterogeneous urban area.The source regions of the LAS-derived fluxes depended highly on the angle between the wind direction and the LAS beam path,as well as the meteorological conditions(e.g.,atmospheric stability,wind speed,friction velocity).When the wind direction was almost parallel to the LAS optical path(Figs.4a and c),the calculated footprint indicated that 90%of the daytime turbulent fluxes were originated from areas within a width of 100–300 m along the optical path(0.2–0.6 km2).It is also clear that the width of in fluence became narrower as the atmosphericstability decreased(more unstable).On the other hand,the daytime turbulent fluxes of 90%were attributed to larger areas,expanding up to about 1 km(～2 km2)upwind from the LAS optical path,when the wind direction was almost perpendicular to the optical path(Figs.4b and d).The results also show thatthe LAS-derived turbulentfluxeswere the least influenced by sources near the LAS transmitter/receiver,as discussed by Chehbouni et al.(2000).

Table 1.Time-averaged meteorological variables used for the representative LAS(large aperture scintillometer)footprint modeling.The wind direction and speed were obtained from the meteorologicaltower,and the Obukhov length(L)and friction velocity(u?)were estimated using the scintillation method.

By virtue of the 3D footprint analysis,it was found that the LAS system,operating with an optical path length of 2.1 km and effective height of about 26 m,can represent the areaaveraged turbulent fluxes over this heterogeneous urban area on a micro-αscale,i.e.,0.2–2 km.The results imply that the LAS’s spatial-integration capability will be very useful to evaluate modeled surface-layer turbulent fluxes over significantly heterogeneous surfaces on the spatial scale.

4.2.LAS-derived sensible heat and momentum fluxes

The measured Cn2ranged from 0.1 to 3.5×10-14m-2/3and 0.1 to 2.8×10-14m-2/3during the daytime periods (0800–1700 LST)on 11 and 12 November,respectively,thus satisfying the saturation-free condition described in section 3.2.The LAS-derived turbulent sensible heat and momentum fluxes were calculated every 15 s from the Cn2values following the procedure described in section 2,and then the fluxes were averaged over 30 min intervals.Figure 5 shows that the LAS-derived momentum fluxes(τ)on 11 and 12 November 2013,which were well correlated with the mean wind speed profiles shown in Fig.3,ranged from 0.05 to 0.35 N m-2during the daytime.The land-surface process-resolving WRF simulation also compared well to the estimation,indicating the LAS estimation is able to represent the spatial-integrated momentum fluxes over heterogeneous surfaces.The ratio of friction velocity to mean wind speed(u?/U)varied from 0.14 in the early morning to 0.11 at later times on 11 November, and from 0.22 during morning to 0.10 in the afternoon on 12 November.The u?/U difference between the morning and the afternoon is attributed to the temporal variation of atmospheric stability,following Eq.(8).

Figure 6 shows that the LAS-derived daytime sensible heat flux(QH)varied between 25 and 160 W m-2on 11 November,and between 20 and 120 W m-2on 12 November.The LAS-derived QHquantities were different in magnitude and temporal variation between the two days,with a dramatic reduction at around 1400 LST 12 November.Temporal fluctuations in the LAS-derived turbulent fluxes re flected the temporal variability of the measured scintillation intensity for the daytime periods.The calculated temporal variability(1σ) in QHwas estimated to be as much as about 21%relative to the 30 min averaged values on 11 November,and 22%on 12 November.The temporal variability inτand its footprint dependency were not as large as in QH(Fig.5)because the momentum fluxes were not directly associated with the scintillation intensity.

The LAS-derived heat fluxes were in good agreement with the land-surface process-resolving simulation results, especially during morning hours(0800–1200 LST).Further comparison to the net radiative fluxes measured at the meteorological site shows that the estimated daily mean QH(107 and 72 W m-2)corresponded to about 49%and 37%of the net radiation of 219 and 194 W m-2on 11 and 12 November 2013,respectively.Despite the intrinsic difficulty in the quantitative verification of the LAS-derived turbulent fluxes over heterogeneous surfaces,the comparison results of the LAS-derived QHto the WRF simulation and the flux ratio analysis show that the scintillation method can reasonably estimate area-averaged surface-layer turbulent fluxes.

The distinctive features of the LAS-derived QHtemporal variation are that the peak QH(especially on 11 November)occurred earlier in the day and that morning QHwas enhanced compared to that in the afternoon(asymmetric distribution).Generally,a peak QHvalue in the surface-layer en-ergy balance is found at 1300–1500 LST on clear days(e.g., Lee and Park,2008;Lee,2011;Loridan et al.,2013).In addition,the maximum values tend to lag in time over urban areas compared to natural surfaces(Lee and Baik,2011).To further investigate the features in the LAS-derived QH,another WRF simulation was conducted by replacing the Geum River with a natural surface type(“shrubland”).The WRF simulation revealed that the peak time lag and the morning enhancement in the LAS-derived QHcan be partly attributed to the existence of the river,where signi ficantly larger QHvalues(～200 W m-2)were found in morning hours(0800–1000 LST)due to the large temperature difference between the water body(～14.0?C)and the overlying atmosphere in the cold season.

4.3.Sensitivity analysis

The LAS system produced surface-layer turbulent fluxes over this heterogeneous urban area within a reasonable uncertainty range(1σ).In this section,the results from several sensitivity tests are discussed to further quantify the influence of input parameters(surface aerodynamic parameters and Bowen ratio)and surface-layer similarity function for the temperature structure function on the LAS-derived turbulent fluxes.

The surface aerodynamic parameters(z0and d)could be estimated by an average of the surface parameters of each homogeneous patch consisting of heterogeneous source areas, as described in section 2.1.However,in practice it is dif ficult to determine the surface parameters for highly heterogeneous surfaces because of the complexity of the morphology of buildings/obstacles and the topography.The uncertainty in the surface aerodynamic parameters also influences the determination of the LAS effective height(zeff).As a realistic ballpark guess of the surface parameters for the study area (e.g.,Garratt,1992;Foken,2008),an aerodynamic roughness length of 0.1 m and 1.0 m and an LAS effective height of about±20%(21 and 31 m)from the base estimation were applied in the sensitivity tests.Besides the surface aerodynamic parameters,the influence of the Bowen ratio was also considered in the sensitivity tests,with values ofβ=0.5 and β=2.0.

Figure 7 shows that the estimatedτwas very sensitive to z0,leading to a reduction ofabout-34%for z0=0.3 m butan increase of about+80%for z0=1.0 m,compared to the base estimation on both days.The influence of the Bowen ratio on the LAS-derivedτwas negligible,while the change in zeffled to momentum-flux difference within±5%.On the other hand,as shown in Fig.8,the estimated QHwas more variable to the change in zeff,ranging between±17%on a daytime average.The change in z0led to the flux difference of between -6%and 14%on 11 November,which was lower than that of between-4%and 9%on 12 November.The difference inβled to the smallest(＜4%)changes in QH.Overall,the sensitivity tests indicate that the uncertainty in the surface parameters may affect the LAS-derived QHwithin±20%.

To quantify the uncertainty by the choice of the temperature structure function parameter,the LAS turbulent fluxes were estimated using four different non-dimensional similarity functions for the temperature structure function parameter.When calculating turbulent sensible heat flux using thescintillation method,previous studies have used various nondimensional similarity functions(e.g.,Meijninger and De Bruin,2000;Lagouarde et al.,2006;Evans et al.,2012;Geli et al.,2012).Along with the surface-layer structure function suggested by Andreas(1988),the function by De Bruin et al. (1993)has been popularly used,which is given for unstable conditions(z/L＜0)as follows:

Foken and Kretschmer(1990)and Thiermann and Grassl (1992)proposed different functional forms as follows:

where k=0.4,and

respectively.By comparing the non-dimensional similarity functions for unstable atmospheric conditions,it was found that the function by Andreas(1988)had higher fTthan that by De Bruin et al.(1993),except for near neutral conditions (-z/L＜0.01),while it had relatively lower fTthan the functions proposed by Foken and Kretschmer(1990)and Thiermann and Grassl(1992)for a wide range of unstable atmospheric conditions(Fig.9).

As shown in Fig.10,the QHfrom the Andreas function was lower by～14%(range:8%–17%)in daytime average than that by the De Bruin function,but larger by about 9%and 10%than those from the functions of Foken and Kretschmer(1990)and Thiermann and Grassl(1992),respectively.The estimated QHquantities were well correlated in time with each other,with a correlation coefficient of 0.99.

The influence of the temperature similarity functions onτ were less than 2%in daytime mean flux(not shown).It is worth noting that Lagouarde et al.(2006)reported the EC-derived QHroughly ranged between the LAS-derived QHusing the De Bruin and Andreas functions for a reasonably homogeneous urban area.Overall,the results indicate that the uncertainty caused by the non-dimensional similarity function forthe temperature structure function parameteriswithin ±15%in the LAS-derived QH.

5.Summary and conclusions

This study tested the ability of the LAS technique in determining surface-layer turbulent sensible heat and momentum fluxes over a heterogeneous urban area,in comparison with land-surface process-resolving numerical modeling. The LAS system was deployed and operated over a highly heterogeneous urban area(Gongju,South Korea)during the cold season with an optical path length of 2.1 km.Threedimensional LAS footprint modeling was also performed to interpret the LAS-derived turbulent fluxes over the heterogeneous surfaces.The WRF model,which was configured to explicitly represent the land-use heterogeneity over the experimental site and associated physical processes,was introduced for verification of the LAS-derived turbulent fluxes.

The LAS-derived turbulentfluxes(QHandτ)showed reasonable temporal variation during the daytime on the two clear days of 11 and 12 November 2013 when compared to the results of the land-surface process-resolving modeling and radiative flux measurements.A series of sensitivity tests showed that the overall uncertainty in the daytimeQHestimate had 20%–30%as an upper limit:within±20% in surface aerodynamic parameters and±15%in the nondimensional similarity function for the temperature structure function parameter.In addition,the estimated error inτwas very sensitive to the determination of the aerodynamic roughness length.It is noted that only two days of data were used in the analysis,which may also have contributed uncertainties to some extent.On the other hand,the footprint modeling analysis showed that the LAS-derived surface-layer turbulent fluxes estimated with the 2.1 km optical path distance were contributed by source areas within 0.2–2 km2(“microαscale”)depending on local meteorological conditions.Unlike the flux footprint modeling of the EC system,the LAS footprint(or source area)was significantly influenced by an incident angle of wind to the LAS beam path,as well as by atmospheric conditions such as atmospheric stability and wind speed.

Recently,mesoscale meteorological and environmental models have undergone rapid development by the inclusion of urban canopy models(e.g.,Masson,2000;Kusaka et al., 2001;Martilli et al.,2002;Lee and Park,2008;Oleson et al., 2008;Porson et al.,2010;Lee,2011;Ryu et al.,2011)for better representation of urban land-surface processes(e.g., Lee et al.,2011).Along with numerical model development,quantitative validation of surface-layer turbulent fluxes modeled by 3D mesoscale meteorological models is very important to improving forecast skill,as well as understanding boundary layer structures.Even though previous studies have successfully determined the surface-layer turbulent fluxes over various homogeneous surfaces using the EC system and/orLAS system(e.g.,Lagouarde etal.,2006;Zeweldi et al.,2010;Evans et al.,2012;Ward et al.,2014),few studies have provided estimations over heterogeneous urban surfaces.The present study demonstrates that the LAS technique can provide realistic surface-layer turbulent fluxes over highly heterogeneous urban areas within a reasonable uncertainty range.In addition,it shows that LAS footprint modeling can provide valuable insight for model–measurement comparisons by the interpretation of sources’contributions to the LAS estimates.Furthermore,this study indicates that LAS-derived turbulent fluxes can be useful for the verification of mesoscale meteorological and environmental modeling conducted over heterogeneous surfaces.

Acknowledgements.The authors gratefully acknowledge the two anonymous reviewers and the editor for their invaluable comments.We also thank Doo-Il LEE for carrying out the WRF simulation.This work was supported by the Korea Meteorological Administration Research and Development Program(Grant No.CATER 2012-3081).

REFERENCES

Andreas,E.L.,1988:Atmospheric stability from scintillation measurements.Appl.Opt.,27,2241–2246.

Beyrich,F.,H.A.R.De Bruin,W.M.L.Meijninger,J.W.Schipper,and H.Lohse,2002:Results from one-year continuous operation of a large aperture scintillometer over a heterogeneous land surface.Bound.-Layer Meteor.,105,85–97.

Businger,J.A.,J.C.Wyngaard,Y.Izumi,and E.F.Bradley,1971: Flux-profile relationships in the atmospheric surface layer.J. Atmos.Sci.,28,181–189.

Chehbouni,A.,and Coauthors,2000:Estimation of heat and momentum fluxes over complex terrain using a large aperture scintillometer.Agricultural and Forest Meteorology,105, 215–226.

Chen,F.,and J.Dudhia,2001:Coupling an advanced land surfacehydrology model with the Penn State-NCAR MM5 modeling system.Part I:Model implementation and sensitivity.Mon. Wea.Rev.,129,569–585.

Chou,M.D.,and M.J.Suarez,1994:An Efficient Thermal Infrared Radiation Parameterization for Use in General Circulation Models.NASA Tech Memo,104606,Vol.3,85 pp.

De Bruin,H.A.R.,W.Kohsiek,and B.J.J.M.van den Hurk, 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.,B.J.J.M.van den Hurk,and W.Kohsiek, 1995:The scintillation method tested over a dry vineyard area.Bound.-Layer Meteor.,76,25–40.

De Bruin,H.A.R.,W.M.L.Meijninger,A.S.Smedman,and M.Magnusson,2002:Displaced-beam small aperture scintillometer test.Part I:The WINTEX data-set.Bound.-Layer Meteor.,105,129–148.

Evans,J.G.,D.D.McNeil,J.W.Finch,T.Murray,R.J.Harding,H.C.Ward,and A.Verhoef,2012:Determination of turbulent heat fluxes using a large aperture scintillometer over undulating mixed agricultural terrain.Agricultural and Forest Meteorology,166–167,221–233.

Foken,T.,2008:Micrometeorology.Springer,306 pp.

Foken,T.,and D.Kretschmer,1990:Stability dependence of the temperature structure parameter.Bound.-Layer Meteor.,53, 185–189.

Frehlich,R.G.,and G.R.Ochs,1990:Effects of saturation on the optical scintillometer.Appl.Opt.,29,548–553.

Garratt,J.R.,1992:The Atmospheric Boundary Layer.Cambridge University Press,316 pp.

Geli,H.M.E.,C.M.U.Neale,D.Watts,J.Osterberg,H.A. R.De Bruin,W.Kohsiek,R.T.Pack,and L.E.Hipps, 2012:Scintillometer-based estimates of sensible heat flux using lidar-derived surface roughness.Journal of Hydrometeorology,13,1317–1331.

Grimmond,C.S.B.,and T.R.Oke,1999:Aerodynamic properties of urban areas derived from analysis of surface form.J.Appl. Meteor.,38,1262–1292.

Hartogensis,O.K.,C.J.Watts,J.C.Rodriguez,and H.A.R. De Bruin,2003:Derivation of an effective height for scintillometers:La Poza experiment in northwest Mexico.Journal of Hydrometeorology,4,915–928.

Hill,R.J.,S.F.Clifford,and R.S.Lawrence,1980:Refractiveindex and absorption fluctuations in the infrared caused by temperature,humidity,and pressure fluctuations.Journal of the Optical Society of America,70,1192–1205.

Hoedjes,J.C.B.,R.M.Zuurbier,and C.J.Watts,2002:Large aperture scintillometer used over a homogeneous irrigated area,partly affected by regional advection.Bound.-Layer Meteor.,105,99–117.

H¨ogstr¨om,U.,1988:Non-dimensional wind and temperatureprofiles in the atmospheric surface layer:A re-evaluation. Bound.-Layer Meteor.,42,55–78.

Horst,T.W.,1999:The footprint for estimation of atmospheresurface exchange fluxes by profile techniques.Bound.-Layer Meteor.,90,171–188.

Horst,T.W.,and J.C.Weil,1992:Footprint estimation for scalar flux measurements in the atmospheric surface layer.Bound.-Layer Meteor.,59,279–296.

Horst,T.W.,and J.C.Weil,1994:How far is far enough?The fetch requirements for micrometeorological measurement of surface fluxes.J.Atmos.Oceanic Technol.,11,1018–1025.

Hsieh,C.I.,G.Katul,and T.W.Chi,2000:An approximate analytical model for footprint estimation of scalar fluxes in thermally stratified atmospheric flows.Advances in Water Resources,23,765–772.

Janji′c,Z.I.,2002:Nonsingular implementation of the Mellor-Yamada Level 2.5 Scheme in the NCEP Meso Model.NCEP Office Note,No.437,61 pp.

Kanda,M.,R.Moriwaki,M.Roth,and T.Oke,2002:Areaaveraged 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.

Kang,J.-H.,M.-S.Suh,and J.-H.Kwak,2010:Land cover classification over East Asian region using recent MODIS NDVI data(2006–2008).Atmosphere,20,415–426.(Korean with an English abstract)

Kipp&Zonen,2012:Instruction Manual.Delft,,Netherlands,86 pp.

Kleissl,J.,O.K.Hartogensis,and J.D.Gomez,2010:Test of scintillometer saturation correction methods using field experimental data.Bound.-Layer Meteor.,137,493–507.

Kljun,N.,P.Calanca,M.W.Rotach,and H.P.Schmid,2004: A simple parameterisation for flux footprint predictions. Bound.-Layer Meteor.,112,503–523.

Kormann,R.,and F.X.Meixner,2001:An analytical footprint model for non-neutral stratification.Bound.-Layer Meteor., 99,207–224.

Kusaka,H.,H.Kondo,Y.Kikegawa,and F.Kimura,2001:A simple single-layer urban canopy model for atmospheric models: Comparison with multi-layer and slab models.Bound.-Layer Meteor.,101,329–358.

Lagouarde,J.P.,M.Irvine,J.M.Bonnefond,C.S.B.Grimmond,N.Long,T.R.Oke,J.A.Salmond,and B.Offerle, 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.

Lee,S.-H.,2011:Further development of the vegetated urban canopy model including a grass-covered surface parametrization and photosynthesis effects.Bound.-Layer Meteor.,140, 315–342.

Lee,S.-H.,and S.-U.Park,2008:A vegetated urban canopy model for meteorological and environmental modelling.Bound.-Layer Meteor.,126,73–102.

Lee,S.-H.,and J.-J.Baik,2011:Evaluation of the vegetated urban canopy model(VUCM)and its impacts on urban boundary layer simulation.Asia-Pacific Journal of Atmospheric Sciences,47(2),151–165.

Lee,S.-H.,and Coauthors,2011:Evaluation of urban surface parameterizations in the WRF model using measurements during the Texas Air Quality Study 2006 field campaign.Atmos. Chem.Phys.,11,2127–2143.

Lin,Y.L.,R.D.Farley,and H.D.Orville,1983:Bulk parameterization of the snow field in a cloud model.J.Climate Appl. Meteor.,22,1065–1092.

Liu,S.M.,Z.W.Xu,W.Z.Wang,Z.Z.Jia,M.J.Zhu,and J. M.Wang,2011:A comparison of eddy-covariance and large aperture scintillometer measurements with respect to the energy balance closure problem.Hydrology and Earth System Sciences,15,1291–1306.

Liu,S.M.,Z.W.Xu,Z.L.Zhu,Z.Z.Jia,and M.J.Zhu,2013: Measurements of evapotranspiration from eddy-covariance systems and large aperture scintillometers in the Hai River Basin,China.J.Hydrol.,487,24–38.

Loridan,T.,F.Lindberg,O.Jorba,S.Kotthaus,S.Grossman-Clarke,and C.S.B.Grimmond,2013:High resolution simulation of the variability of surface energy balance fluxes across central London with urban zones for energy partitioning.Bound.-Layer Meteor.,147,493–523.

Martilli,A.,A.Clappier,and M.W.Rotach,2002:An urban surface exchange parameterisation for mesoscale models. Bound.-Layer Meteor.,104,261–304.

Macdonald,R.W.,R.F.Griffiths,and D.J.Hall,1998:An improved method for the estimation of surface roughness of obstacle arrays.Atmos.Environ.,32,1857–1864.

Masson,V.,2000:A physically-based scheme for the urban energy budget in atmospheric models.Bound.-Layer Meteor., 94,357–397.

McAneney,K.J.,A.E.Green,and M.S.Astill,1995:Largeaperture scintillometry:The homogeneous case.Agricultural and Forest Meteorology,76,149–162.

Meijninger,W.M.L.,and H.A.R.De Bruin,2000:The sensible heat fluxes over irrigated areas in western Turkey determined with a large aperture scintillometer.J.Hydrol.,229,42–49.

Meijninger,W.M.L.,A.E.Green,O.K.Hartogensis,W.Kohsiek, J.C.B.Hoedjes,R.M.Zuurbier,and H.A.R.De Bruin, 2002:Determination of area-averaged water vapour fluxes with large aperture and radio wave scintillometers over a heterogeneous surface-Flevoland field experiment.Bound.-Layer Meteor.,105,63–83.

Mlawer,E.J.,S.J.Taubman,P.D.Brown,M.J.Iacono,and S.A. Clough,1997:Radiative transfer for inhomogeneous atmospheres:RRTM,a validated correlated-k model for the longwave.J.Geophys.Res.,102,16 663–16 682.

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.

Offerle,B.,C.S.B.Grimmond,K.Fortuniak,and W.Pawlak, 2006:Intraurban differencesof surface energy fluxesin a central European city.J.Appl.Meteor.Climatol.,45,125–136.

Oleson,K.W.,G.B.Bonan,J.Feddema,M.Vertenstein,and C.S.B.Grimmond,2008:An urban parameterization for a global climate model.Part I:Formulation and evaluation for two cities.J.Appl.Meteor.Climatol.,47,1038–1060.

Panofsky,H.A.,and J.A.Dutton,1984:Atmospheric Turbulence: Models and Methods for Engineering Applications.John Wiley and Sons,New York,397 pp.

Pasquill,F.,1974:Atmospheric Diffusion.2nd ed.,John Wiley& Sons,New York,425 pp.

Paulson,C.A.,1970:The mathematical representation of wind speed and temperature profiles in the unstable atmospheric surface layer.J.Appl.Meteor.,9,857–861.

Pauscher,L.,2010:Scintillometer measurements above the urban area of London.Diploma thesis,Dept.of Micrometeorology,University of Bayreuth,95 pp.

Porson,A.,P.A.Clark,I.N.Harman,M.J.Best,and S.E.Belcher, 2010:Implementation of a new urban energy budget scheme in the MetUM.Part I:Description and idealized simulations. Quart.J.Roy.Meteor.Soc.,136,1514–1529.

Raupach,M.R.,1994:Simplified expressions for vegetation roughness length and zero-plane displacement as functions of canopy height and area index.Bound.-Layer Meteor.,71, 211–216.

Roth,M.W.,2000:Review of atmospheric turbulence over cities. Quart.J.Roy.Meteor.Soc.,126,941–990.

Ryu,Y.-H.,J.-J.Baik,and S.-H.Lee,2011:A new single-layer urban canopy model for use in mesoscale atmospheric models. J.Appl.Meteor.Climatol.,50,1773–1794.

Schmid,H.P.,1994:Source areas for scalars and scalar fluxes. Bound.-Layer Meteor.,67,293–318.

Schuepp,P.H.,M.Y.Leclerc,J.I.Macpherson,and R.L.Desjardins,1990:Footprint prediction of scalar fluxes from analytical solutions of the diffusion equation.Bound.-Layer Meteor.,50,355–373.

Skamarock,W.C.,and Coauthors,2008:A description of the advanced research WRF version 3.NCAR Technical Note, NCAR/TN-475+STR,113 pp.

Tatarskii,V.I.,1961:Wave Propagation in a Turbulent Medium. McGraw-Hill,285 pp.

Thiermann,V.,and H.Grassl,1992:The measurement of turbulent surface-layer fluxes by use of bichromatic scintillation. Bound.-Layer Meteor.,58,367–389.

Timmermans,W.J.,Z.Su,and A.Olioso,2009:Footprintissues in scintillometry over heterogeneous landscapes.Hydrol.Earth Syst.Sci.,13(11),2179–2190.

van Ulden,A.P.,1978:Simple estimates for vertical diffusion from sources near the ground.Atmos.Environ.,12,2125–2129.

Wang,T.I.,G.R.Ochs,and S.F.Clifford,1978:A saturationresistant optical scintillometer to measure C2n.Journal of the Optical Society of America,69,334–338.

Ward,H.C.,J.G.Evans,and C.S.B.Grimmond,2014:Multiscale sensible heat fluxes in the suburban environment from large-aperture scintillometry and eddy covariance.Bound.-Layer Meteor.,152,65–89.

Wesely,M.L.,1976:The combined effect of temperature and humidity fluctuations on refractive index.J.Appl.Meteor.,15, 43–49.

Wilson,K.M.,A.van Tol,and J.Mes,2013:The upgraded Kipp&Zonen LAS MkII large aperture scintillometer instrument specifications.Tubingen Atmospheric Physics Symposium“Scintillometers and Applications”,7–9 Oct 2013,Tubingen,Germany.

Wyngaard,J.C.,Y.Izumi,and S.A.Collins Jr.,1971:Behavior of the refractive-index-structure parameter near the ground. Journal of the Optical Society of America,61,1646–1650.

Xu,Z.W.,and Coauthors,2013:Intercomparison of surface energy flux measurement systems used during the HiWATERMUSOEXE.J.Geophys.Res.,118,13 140–13 157.

Zeweldi,D.A.,M.Gebremichael,J.M.Wang,T.Sammis,J. Kleissl,and D.Miller,2010:Intercomparison of sensible heat flux from large aperture scintillometer and eddy covariance methods:Field experiment over a homogeneous semi-arid region.Bound.-Layer Meteor.,135,151–159.

:Lee,S.-H.,J.-H.Lee,and B.-Y.Kim,2015:Estimation of turbulent sensible heat and momentum fluxes over a heterogeneous urban area using a large aperture scintillometer.Adv.Atmos.Sci.,32(8),1092–1105,

10.1007/s00376-015-4236-2.

25 October 2014;revised 16 December 2014;accepted 9 January 2015)

?Corresponding author:Sang-Hyun LEE

Email:sanghyun@kongju.ac.kr