Numerical framework for the computation of urban flux footprints employing large-eddy simulation and Lagrangian stochastic modeling

Abstract. Conventional footprint models cannot account for the heterogeneity of the urban landscape imposing a pronounced uncertainty on the spatial interpretation of eddy-covariance (EC) flux measurements in urban studies. This work introduces a computational methodology that enables the generation of detailed footprints in arbitrarily complex urban flux measurements sites. The methodology is based on conducting high-resolution large-eddy simulation (LES) and Lagrangian stochastic (LS) particle analysis on a model that features a detailed topographic description of a real urban environment. The approach utilizes an arbitrarily sized target volume set around the sensor in the LES domain, to collect a dataset of LS particles which are seeded from the potential source area of the measurement and captured at the sensor site. The urban footprint is generated from this dataset through a piecewise postprocessing procedure, which divides the footprint evaluation into multiple independent processes that each yield an intermediate result. These results are ultimately selectively combined to produce the final footprint. The strategy reduces the computational cost of the LES–LS simulation and incorporates techniques to account for the complications that arise when the EC sensor is mounted on a building instead of a conventional flux tower. The presented computational framework also introduces a result assessment strategy which utilizes the obtained urban footprint together with a detailed land cover type dataset to estimate the potential error that may arise if analytically derived footprint models were employed instead. The methodology is demonstrated with a case study that concentrates on generating the footprint for a building-mounted EC measurement station in downtown Helsinki, Finland, under the neutrally stratified atmospheric boundary layer.