Development of the Surface Urban Energy and Water balance Scheme (SUEWS) for cold climate cities

The Surface Urban Energy and Water Balance Scheme (SUEWS) is developed to include snow. The processes addressed include accumulation of snow on the different urban surface types: snow albedo and density aging, snow melting and re-freezing of meltwater. Individual model parameters are assessed and independently evaluated using long-term observations in the two cold climate cities of Helsinki and Montreal. Eddy covariance sensible and latent heat fluxes and snow depth observations are available for two sites in Montreal and one in Helsinki. Surface runoff from two catchments (24 and 45 ha) in Helsinki and snow properties (albedo and density) from two sites in Montreal are also analysed. As multiple observation sites with different land-cover characteristics are available in both cities, model development is conducted independent of evaluation. The developed model simulates snowmelt related runoff well (within 19% and 3% for the two catchments in Helsinki when there is snow on the ground), with the springtime peak estimated correctly. However, the observed runoff peaks tend to be smoother than the simulated ones, likely due to the water holding capacity of the catchments and the missing time lag between the catchment and the observation point in the model. For all three sites the model simulates the timing of the snow accumulation and melt events well, but underestimates the total snow depth by 18–20% in Helsinki and 29–33% in Montreal. The model is able to reproduce the diurnal pattern of net radiation and turbulent fluxes of sensible and latent heat during cold snow, melting snow and snow-free periods. The largest model uncertainties are related to the timing of the melting period and the parameterization of the snowmelt. The results show that the enhanced model can simulate correctly the exchange of energy and water in cold climate cities at sites with varying surface cover.

An overview on the urban boundary-layer atmosphere network in Helsinki

The Helsinki Urban Boundary-Layer Atmosphere Network (UrBAN: http://urban.fmi.fi) is a dedicated research-grade observational network where the physical processes in the atmosphere above the city are studied. Helsinki UrBAN is the most poleward intensive urban research observation network in the world and thus will allow studying some unique features such as strong seasonality. The network's key purpose is for the understanding of the physical processes in the urban boundary layer and associated fluxes of heat, momentum, moisture, and other gases. A further purpose is to secure a research-grade database, which can be used internationally to validate and develop numerical models of air quality and weather prediction. Scintillometers, a scanning Doppler lidar, ceilometers, a sodar, eddy-covariance stations, and radiometers are used. This equipment is supplemented by auxiliary measurements, which were primarily set up for general weather and/or air-quality mandatory purposes, such as vertical soundings and the operational Doppler radar network. Examples are presented as a testimony to the potential of the network for urban studies, such as (i) evidence of a stable boundary layer possibly coupled to an urban surface, (ii) the comparison of scintillometer data with sonic anemometry above an urban surface, (iii) the application of scanning lidar over a city, and (iv) combination of sodar and lidar to give a fuller range of sampling heights for boundary layer profiling.

Infrared and millimetre-wave scintillometry in the suburban environment – Part 1: Structure parameters

Scintillometry, a form of ground-based remote sensing, provides the capability to estimate surface heat fluxes over scales of a few hundred metres to kilometres. Measurements are spatial averages, making this technique particularly valuable over areas with moderate heterogeneity such as mixed agricultural or urban environments. In this study, we present the structure parameters of temperature and humidity, which can be related to the sensible and latent heat fluxes through similarity theory, for a suburban area in the UK. The fluxes are provided in the second paper of this two-part series. A millimetre-wave scintillometer was combined with an infrared scintillometer along a 5.5 km path over northern Swindon. The pairing of these two wavelengths offers sensitivity to both temperature and humidity fluctuations, and the correlation between wavelengths is also used to retrieve the path-averaged temperature–humidity correlation. Comparison is made with structure parameters calculated from an eddy covariance station located close to the centre of the scintillometer path. The performance of the measurement techniques under different conditions is discussed. Similar behaviour is seen between the two data sets at sub-daily timescales. For the two summer-to-winter periods presented here, similar evolution is displayed across the seasons. A higher vegetation fraction within the scintillometer source area is consistent with the lower Bowen ratio observed (midday Bowen ratio < 1) compared with more built-up areas around the eddy covariance station. The energy partitioning is further explored in the companion paper.

Infrared and millimetre-wave scintillometry in the suburban environment – Part 2: Large-area sensible and latent heat fluxes

A millimetre-wave scintillometer was paired with an infrared scintillometer, enabling estimation of large-area evapotranspiration across northern Swindon, a suburban area in the UK. Both sensible and latent heat fluxes can be obtained using this "two-wavelength" technique, as it is able to provide both temperature and humidity structure parameters, offering a major advantage over conventional single-wavelength scintillometry. The first paper of this two-part series presented the measurement theory and structure parameters. In this second paper, heat fluxes are obtained and analysed. These fluxes, estimated using two-wavelength scintillometry over an urban area, are the first of their kind. Source area modelling suggests the scintillometric fluxes are representative of 5–10 km2. For comparison, local-scale (0.05–0.5 km2) fluxes were measured by an eddy covariance station. Similar responses to seasonal changes are evident at the different scales but the energy partitioning varies between source areas. The response to moisture availability is explored using data from 2 consecutive years with contrasting rainfall patterns (2011–2012). This extensive data set offers insight into urban surface-atmosphere interactions and demonstrates the potential for two-wavelength scintillometry to deliver fluxes over mixed land cover, typically representative of an area 1–2 orders of magnitude greater than for eddy covariance measurements. Fluxes at this scale are extremely valuable for hydro-meteorological model evaluation and assessment of satellite data products

Quantifying particle size and turbulent scale dependence of dust uplift in the Sahara using aircraft measurements

The first size-resolved airborne measurements of dust fluxes and the first dust flux measurements from the central Sahara are presented and compared with a parameterization by Kok (2011a). High-frequency measurements of dust size distribution were obtained from 0.16 to 300 µm diameter, and eddy covariance fluxes were derived. This is more than an order of magnitude larger size range than previous flux estimates. Links to surface emission are provided by analysis of particle drift velocities. Number flux is described by a −2 power law between 1 and 144 µm diameter, significantly larger than the 12 µm upper limit suggested by Kok (2011a). For small particles, the deviation from a power law varies with terrain type and the large size cutoff is correlated with atmospheric vertical turbulent kinetic energy, suggesting control by vertical transport rather than emission processes. The measured mass flux mode is in the range 30–100 µm. The turbulent scales important for dust flux are from 0.1 km to 1–10 km. The upper scale increases during the morning as boundary layer depth and eddy size increase. All locations where large dust fluxes were measured had large topographical variations. These features are often linked with highly erodible surface features, such as wadis or dunes. We also hypothesize that upslope flow and flow separation over such features enhance the dust flux by transporting large particles out of the saltation layer. The tendency to locate surface flux measurements in open, flat terrain means these favored dust sources have been neglected in previous studies.

Joint data assimilation of satellite reflectance and net ecosystem exchange data constrains ecosystem carbon fluxes at a high-elevation subalpine forest

We utilized an ecosystem process model (SIPNET, simplified photosynthesis and evapotranspiration model) to estimate carbon fluxes of gross primary productivity and total ecosystem respiration of a high-elevation coniferous forest. The data assimilation routine incorporated aggregated twice-daily measurements of the net ecosystem exchange of CO2 (NEE) and satellite-based reflectance measurements of the fraction of absorbed photosynthetically active radiation (fAPAR) on an eight-day timescale. From these data we conducted a data assimilation experiment with fifteen different combinations of available data using twice-daily NEE, aggregated annual NEE, eight-day f AP AR, and average annual fAPAR. Model parameters were conditioned on three years of NEE and fAPAR data and results were evaluated to determine the information content from the different combinations of data streams. Across the data assimilation experiments conducted, model selection metrics such as the Bayesian Information Criterion and Deviance Information Criterion obtained minimum values when assimilating average annual fAPAR and twice-daily NEE data. Application of wavelet coherence analyses showed higher correlations between measured and modeled fAPAR on longer timescales ranging from 9 to 12 months. There were strong correlations between measured and modeled NEE (R2, coefficient of determination, 0.86), but correlations between measured and modeled eight-day fAPAR were quite poor (R2 = −0.94). We conclude that this inability to determine fAPAR on eight-day timescale would improve with the considerations of the radiative transfer through the plant canopy. Modeled fluxes when assimilating average annual fAPAR and annual NEE were comparable to corresponding results when assimilating twice-daily NEE, albeit at a greater uncertainty. Our results support the conclusion that for this coniferous forest twice-daily NEE data are a critical measurement stream for the data assimilation. The results from this modeling exercise indicate that for this coniferous forest, average annuals for satellite-based fAPAR measurements paired with annual NEE estimates may provide spatial detail to components of ecosystem carbon fluxes in proximity of eddy covariance towers. Inclusion of other independent data streams in the assimilation will also reduce uncertainty on modeled values.

CO2 emission estimation in the urban environment: measurement of the CO2 storage term

Eddy covariance has been used in urban areas to evaluate the net exchange of CO2 between the surface and the atmosphere. Typically, only the vertical flux is measured at a height 2–3 times that of the local roughness elements; however, under conditions of relatively low instability, CO2 may accumulate in the airspace below the measurement height. This can result in inaccurate emissions estimates if the accumulated CO2 drains away or is flushed upwards during thermal expansion of the boundary layer. Some studies apply a single height storage correction; however, this requires the assumption that the response of the CO2 concentration profile to forcing is constant with height. Here a full seasonal cycle (7th June 2012 to 3rd June 2013) of single height CO2 storage data calculated from concentrations measured at 10 Hz by open path gas analyser are compared to a data set calculated from a concurrent switched vertical profile measured (2 Hz, closed path gas analyser) at 10 heights within and above a street canyon in central London. The assumption required for the former storage determination is shown to be invalid. For approximately regular street canyons at least one other measurement is required. Continuous measurements at fewer locations are shown to be preferable to a spatially dense, switched profile, as temporal interpolation is ineffective. The majority of the spectral energy of the CO2 storage time series was found to be between 0.001 and 0.2 Hz (500 and 5 s respectively); however, sampling frequencies of 2 Hz and below still result in significantly lower CO2 storage values. An empirical method of correcting CO2 storage values from under-sampled time series is proposed.

Seasonal surface urban energy balance and wintertime stability simulated using three land-surface models in the high-latitude city Helsinki

The performance of three urban land surface models, run in offline mode, with their default external parameters, is evaluated for two distinctly different sites in Helsinki: Torni and Kumpula. The former is a dense city centre site with 22% vegetation, while the latter is a suburban site with over 50% vegetation. At both locations the models are compared against sensible and latent heat fluxes measured using the eddy covariance technique, along with snow depth observations. The cold climate experienced by the city causes strong seasonal variations that include snow cover and stable atmospheric conditions. Most of the time the three models are able to account for the differences between the study areas as well as the seasonal and diurnal variability of the energy balance components. However, the performances are not systematic across the modelled components, season and surface type. The net all-wave radiation is well simulated, with the greatest uncertainties related to snowmelt timing, when the fraction of snow cover has a key role, particularly in determining the surface albedo. For the turbulent fluxes, more variation between the models is seen which can partly be explained by the different methods in their calculation and partly by surface parameter values. For the sensible heat flux, simulation of wintertime values was the main problem, which also leads to issues in predicting near-surface stabilities particularly at the dense city centre site. All models have the most difficulties in simulating latent heat flux. This study particularly emphasizes that improvements are needed in the parameterization of anthropogenic heat flux and thermal parameters in winter, snow cover in spring and evapotranspiration in order to improve the surface energy balance modelling in cold climate cities.

Surface Urban Energy and Water Balance Scheme (SUEWS): development and evaluation at two UK sites

The Surface Urban Energy and Water Balance Scheme (SUEWS) is evaluated at two locations in the UK: a dense urban site in the centre of London and a residential suburban site in Swindon. Eddy covariance observations of the turbulent fluxes are used to assess model performance over a twoyear period (2011-2013). The distinct characteristics of the sites mean their surface energy exchanges differ considerably. The model suggests the largest differences can be attributed to surface cover (notably the proportion of vegetated versus impervious area) and the additional energy supplied by human activities. SUEWS performs better in summer than winter, and better at the suburban site than the dense urban site. One reason for this is the bias towards suburban summer field campaigns in observational data used to parameterise this (and other) model(s). The suitability of model parameters (such as albedo, energy use and water use) for the UK sites is considered and, where appropriate, alternative values are suggested. An alternative parameterisation for the surface conductance is implemented, which permits greater soil moisture deficits before evaporation is restricted at non-irrigated sites. Accounting for seasonal variation in the estimation of storage heat flux is necessary to obtain realistic wintertime fluxes.

Fluxes of CO(2) above a plantation of Eucalyptus in southeast Brazil

Eddy-covariance measurements of net ecosystem exchange of CO(2) (NEE) and estimates of gross ecosystem productivity (GEP) and ecosystem respiration (R(E)) were obtained in a 2-4 year old Eucalyptus plantation during two years with very different winter rainfall In the first (drier) year the annual NEE GEP and RE were lower than the sums in the second (normal) year and conversely the total respiratory costs of assimilated carbon were higher in the dry year than in the normal year Although the net primary production (NPP) in the first year was 23% lower than that of the second year the decrease in the carbon use efficiency (CUE = NPP/GEP) was 11% and autotrophic respiration utilized more resources in the first dry year than in the second normal year The time variations in NEE were followed by NPP because in these young Eucalyptus plantations NEE is very largely dominated by NPP and heterotrophic respiration plays only a relatively minor role During the dry season a pronounced hysteresis was observed in the relationship between NEE and photosynthetically active radiation and NEE fluxes were inversely proportional to humidity saturation deficit values greater than 0 8 kPa Nighttime fluxes of CO(2) during calm conditions when the friction velocity (u) was below the threshold (0 25 ms(-1)) were estimated based on a Q(10) temperature-dependence relationship adjusted separately for different classes of soil moisture content which regulated the temperature sensitivity of ecosystem respiration (C) 2010 Elsevier B V All rights reserved

The energy and water balance of a Eucalyptus plantation in southeast Brazil

The eddy covariance method was used to measure energy and water balance of a plantation of Eucalyptus (grandis x urophylla) hybrids over a 2 year period. The average daily evaporation rates were 5.4 (+/- 2.0) mm day(-1) in summer, but fell to 1.2 (+/- 0.3) mm day(-1) in winter. In contrast, the sensible heat flux was relatively low in summer but dominated the energy balance in winter. Evaporation accounted for 80% and 26% of the available energy, in summer and winter respectively. The annual evaporation was 82% (1124 mm) and 96% (1235 mm) of the annual rainfall recorded during the first and second year, respectively. Daily average canopy and aerodynamic conductance to water vapour were in the summer 51.9 (+/- 38.4) mm s(-1) 84.1 (+/- 25.6) mm s(-1), respectively; and in the winter 6.0 (+/- 10.5) mm s(-1) and 111.6 (+/- 24.6) mm s(-1), respectively. (C) 2010 Elsevier B.V. All rights reserved.

The land-atmosphere water flux in the tropics

Tropical vegetation is a major source of global land surface evapotranspiration, and can thus play a major role in global hydrological cycles and global atmospheric circulation. Accurate prediction of tropical evapotranspiration is critical to our understanding of these processes under changing climate. We examined the controls on evapotranspiration in tropical vegetation at 21 pan-tropical eddy covariance sites, conducted a comprehensive and systematic evaluation of 13 evapotranspiration models at these sites, and assessed the ability to scale up model estimates of evapotranspiration for the test region of Amazonia. Net radiation was the strongest determinant of evapotranspiration (mean evaporative fraction was 0.72) and explained 87% of the variance in monthly evapotranspiration across the sites. Vapor pressure deficit was the strongest residual predictor (14%), followed by normalized difference vegetation index (9%), precipitation (6%) and wind speed (4%). The radiation-based evapotranspiration models performed best overall for three reasons: (1) the vegetation was largely decoupled from atmospheric turbulent transfer (calculated from X decoupling factor), especially at the wetter sites; (2) the resistance-based models were hindered by difficulty in consistently characterizing canopy (and stomatal) resistance in the highly diverse vegetation; (3) the temperature-based models inadequately captured the variability in tropical evapotranspiration. We evaluated the potential to predict regional evapotranspiration for one test region: Amazonia. We estimated an Amazonia-wide evapotranspiration of 1370 mm yr(-1), but this value is dependent on assumptions about energy balance closure for the tropical eddy covariance sites; a lower value (1096 mm yr(-1)) is considered in discussion on the use of flux data to validate and interpolate models.

Aerosol properties, in-canopy gradients, turbulent fluxes and VOC concentrations at a pristine forest site in Amazonia

Aerosol physical and chemical properties were measured in a forest site in central Amazonia (Cuieiras reservation, 2.61S; 60.21W) during the dry season of 2004 (Aug-Oct). Aerosol light scattering and absorption, mass concentration, elemental composition and size distributions were measured at three tower levels (Ground: 2 m; Canopy: 28 m, and Top: 40 m). For the first time, simultaneous eddy covariance fluxes of fine mode particles and volatile organic compounds (VOC) were measured above the Amazonian forest canopy. Aerosol fluxes were measured by eddy covariance using a Condensation Particle Counter (CPC) and a sonic anemometer. VOC fluxes were measured by disjunct eddy covariance using a Proton Transfer Reaction Mass Spectrometer (PTR-MS). At nighttime, a strong vertical gradient of phosphorus and potassium in the aerosol coarse mode was observed, with higher concentrations at Ground level. This suggests a source of primary biogenic particles below the canopy. Equivalent black carbon measurements indicate the presence of light-absorbing aerosols from biogenic origin. Aerosol number size distributions typically consisted of superimposed Aitken (76 nm) and accumulation modes (144 nm), without clear events of new particle formation. Isoprene and monoterpene fluxes reached respectively 7.4 and 0.82 mg m(-2) s(-1) around noon. An average fine particle flux of 0.05 +/- 0.10 10(6) m(-2) s(-1) was calculated, denoting an equilibrium between emission and deposition fluxes of fine mode particles at daytime. No significant correlations were found between VOC and fine mode aerosol concentrations or fluxes. (C) 2009 Elsevier Ltd. All rights reserved.

Atmosphere and hydrological controls of the evapotranspiration over a floodplain forest in the Bananal Island region, Amazonia

This article discusses seasonal and interannual variations of the evapotranspiration (ET) rates in Bananal Island floodplain, Brazil. Measurements included ET and sensible heat flux using the eddy covariance method, atmospheric forcings (net radiation, Rn, vapor pressure deficit, VPD, wind speed and air temperature), soil moisture profiles, groundwater level and flood height, taken from November 2003 to December 2006. For the hydrological years (October-September) of 2003/2004, 2004/2005 and 2005/2006, the accumulated precipitation was 1692, 1471, 1914 mm and the accumulated ET was 1361, 1318 and 1317 mm, respectively. Seasonal analyses indicated that ET decreased in the dry season (average 3.7 mm day(-1)), despite the simultaneous increase in Rn, air temperature and VPD. The increase of ET in the wet season and particularly in the flood period (average 4.1 mm day(-1)) showed that the free water surface evaporation strongly influenced the energy exchange. Soil moisture, which was substantially depleted during the dry season, and adaptative vegetation mechanisms such as leaf senescence contributed to limit the dry season ET. Strong drainage within permeable sandy soils helped to explain the soil moisture depletion. These results suggest that the Bananal flooding area shows a different pattern in relation to the upland Amazon forests, being more similar to the savanna strictu senso areas in central Brazil. For example, seasonal ET variation was not in phase with Rn; the wet season ET was higher than the dry season ET; and the system stored only a tiny memory of the flooding period, being sensitive to extended drought periods.

Análise e estimativa dos componentes do balanço de energia em ecossistema de manguezal amazônico

As medições e estimativas dos componentes do balanço de energia foram feitos acima da copa das árvores no ecossistema de manguezal natural, localizada a 30 km da cidade de Bragança-PA, entre novembro de 2002 e agosto de 2003. Os dados foram utilizados para a análise das variações sazonais e horárias do fluxo de calor sensível e calor latente, bem como a avaliação da partição de energia. Os dados meteorológicos foram coletados pela estação meteorológica automática (EMA) e os fluxos foram calculados utilizando-se a técnica de covariância de vórtices turbulentos. Os modelos de Penman-Monteith e Shuttleworth foram usados para estimar o fluxo de calor sensível e calor latente. O objetivo deste estudo foi analisar o equilíbrio e a partição de energia no manguezal, assim como fazer uma avaliação do comportamento de modelos empíricos para estimar os fluxos de energia. O saldo de radiação apresentou valores mais elevados no período menos chuvoso. A razão de Bowen mostrou valor geralmente baixo, o que indica que uma proporção maior de energia foi utilizada sob a forma de calor latente. O modelo Shuttleworth é mais eficiente na estimativa de fluxos de calor sensível. Para estimar o fluxo de calor latente do modelo de Penman-Monteith foi mais eficiente durante a estação seca e o modelo Shuttleworth durante a estação chuvosa.