# Harmonized gap-filled datasets from 20 urban flux tower sites

Mathew Lipson, Sue Grimmond, Martin Best, Winston T. L. Chow, Andreas Christen, Nektarios Chrysoulakis, Andrew Coutts, Ben Crawford, Stevan Earl, Jonathan Evans, et al.

► To cite this version:

Mathew Lipson, Sue Grimmond, Martin Best, Winston T. L. Chow, Andreas Christen, et al. Harmonized gap-filled datasets from 20 urban flux tower sites. Earth System Science Data, 2022, 14, pp.5157-5178. 10.5194/essd-14-5157-2022. insu-04471018

HAL Id: insu-04471018

https://insu.hal.science/insu-04471018

Submitted on 21 Feb 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_1.png" /></figure>

Harmonized gap-filled datasets from
20 urban flux tower sites

Mathew Lipson<fnref n="1" />,<fnref n="2" />, Sue Grimmond<fnref n="2" />, Martin Best<fnref n="3" />, Winston T. L. Chow<fnref n="4" />, Andreas Christen<fnref n="5" />,
Nektarios Chrysoulakis<fnref n="6" />, Andrew Coutts<fnref n="7" />, Ben Crawford<fnref n="8" />, Stevan Earl<fnref n="9" />, Jonathan Evans<fnref n="10" />,
Krzysztof Fortuniak<fnref n="11" />, Bert G. Heusinkveld<fnref n="12" />, Je-Woo Hong¹³, Jinkyu Hong<fnref n="14" />, Leena Järvi<fnref n="15" />,
Sungsoo Jo<fnref n="16" />, Yeon-Hee Kim<fnref n="17" />, Simone Kotthaus<fnref n="18" />, Keunmin Lee<fnref n="19" />, Valéry Masson<fnref n="20" />,
Joseph P. McFadden<fnref n="21" />, Oliver Michels<fnref n="22" />, Wlodzimierz Pawlak<fnref n="23" />, Matthias Roth<fnref n="24" />,
Hirofumi Sugawara<fnref n="25" />, Nigel Tapper<fnref n="26" />, Erik Velasco²⁷, and Helen Claire Ward²⁸

<footnote n="1">¹Australian Research Council (ARC) Centre of Excellence for Climate System Science, Climate Change Research Centre, Level 4, Mathews Building, UNSW Sydney, New South Wales 2052, Australia</footnote>

<footnote n="2">²Department of Meteorology, University of Reading, Reading, RG6 6ET, United Kingdom</footnote>

<footnote n="3">³Met Office, Fitzroy Road, Exeter, Devon, EX1 3PB, United Kingdom</footnote>

<footnote n="4">⁴College of Integrative Studies, Singapore Management University, Singapore, Singapore</footnote>

<footnote n="5">⁵Environmental Meteorology, Institute of Earth and Environmental Sciences, Faculty of Environment and Natural Resources, University of Freiburg, Freiburg, Germany</footnote>

<footnote n="6">⁶Foundation for Research and Technology Hellas, Institute of Applied and Computational Mathematics, Remote Sensing Lab, Heraklion, Greece</footnote>

<footnote n="7">⁷School of Earth, Atmosphere and Environment, Monash University, Melbourne, Australia</footnote>

<footnote n="8">⁸Geography and Environmental Sciences, University of Colorado, Denver, Colorado, USA</footnote>

<footnote n="9">⁹Global Institute of Sustainability and Innovation, Arizona State University, Tempe, Arizona, USA</footnote>

<footnote n="10">¹⁰UK Centre for Ecology & Hydrology, Crowmarsh Gifford, Wallingford, United Kingdom</footnote>

<footnote n="11">¹¹Department of Meteorology and Climatology, University of Lodz, Lodz, Poland</footnote>

<footnote n="12">¹²Department of Meteorology and Air Quality, Wageningen University, Wageningen, the Netherlands
¹³Korea Environment Institute, Sejong, Republic of Korea</footnote>

<footnote n="14">¹⁴Ecosystem-Atmosphere Process Lab, Department of Atmospheric Sciences,
Yonsei University, Seoul, Republic of Korea</footnote>

<footnote n="15">¹⁵Institute for Atmospheric and Earth System Research/Physics, Faculty of Science,
University of Helsinki, Helsinki, Finland</footnote>

<footnote n="16">¹⁶Ecosystem-Atmosphere Process Lab, Department of Atmospheric Sciences,
Yonsei University, Seoul, Republic of Korea</footnote>

<footnote n="17">¹⁷National Institute of Meteorological Sciences, Korea Meteorological Administration,
Jeju, Republic of Korea</footnote>

<footnote n="18">¹⁸Institut Pierre Simon Laplace (IPSL), CNRS, École Polytechnique, Institut Polytechnique de Paris,
91128 Palaiseau CEDEX, France</footnote>

<footnote n="19">¹⁹Ecosystem-Atmosphere Process Lab, Department of Atmospheric Sciences,
Yonsei University, Seoul, Republic of Korea</footnote>

<footnote n="20">²⁰Centre National de Recherches Météorologiques, University of Toulouse,
Météo-France and CNRS, Toulouse, France</footnote>

<footnote n="21">²¹Department of Geography and Earth Research Institute, University of California, Santa Barbara, USA</footnote>

<footnote n="22">²²Environmental Meteorology, Albert-Ludwigs-University, Freiburg, Germany</footnote>

<footnote n="23">²³Department of Meteorology and Climatology, University of Lodz, Lodz, Poland</footnote>

<footnote n="24">²⁴Department of Geography, National University of Singapore, Singapore, Singapore</footnote>

<footnote n="25">²⁵Department of Earth and Ocean Sciences, National Defense Academy of Japan, Yokosuka, Japan</footnote>

<footnote n="26">²⁶School of Earth, Atmosphere and Environment, Monash University, Melbourne, Australia</footnote>
<figure><img src="image_2.png" /></figure>

Harmonized gap-filled datasets from
20 urban flux tower sites

Mathew Lipson<fnref n="1" />,<fnref n="2" />, Sue Grimmond<fnref n="2" />, Martin Best<fnref n="3" />, Winston T. L. Chow<fnref n="4" />, Andreas Christen<fnref n="5" />,
Nektarios Chrysoulakis<fnref n="6" />, Andrew Coutts<fnref n="7" />, Ben Crawford<fnref n="8" />, Stevan Earl<fnref n="9" />, Jonathan Evans<fnref n="10" />,
Krzysztof Fortuniak<fnref n="11" />, Bert G. Heusinkveld<fnref n="12" />, Je-Woo Hong¹³, Jinkyu Hong<fnref n="14" />, Leena Järvi<fnref n="15" />,
Sungsoo Jo<fnref n="16" />, Yeon-Hee Kim<fnref n="17" />, Simone Kotthaus<fnref n="18" />, Keunmin Lee<fnref n="19" />, Valéry Masson<fnref n="20" />,
Joseph P. McFadden<fnref n="21" />, Oliver Michels<fnref n="22" />, Wlodzimierz Pawlak<fnref n="23" />, Matthias Roth<fnref n="24" />,
Hirofumi Sugawara<fnref n="25" />, Nigel Tapper<fnref n="26" />, Erik Velasco²⁷, and Helen Claire Ward²⁸

<footnote n="1">¹Australian Research Council (ARC) Centre of Excellence for Climate System Science, Climate Change Research Centre, Level 4, Mathews Building, UNSW Sydney, New South Wales 2052, Australia</footnote>

<footnote n="2">²Department of Meteorology, University of Reading, Reading, RG6 6ET, United Kingdom</footnote>

<footnote n="3">³Met Office, Fitzroy Road, Exeter, Devon, EX1 3PB, United Kingdom</footnote>

<footnote n="4">⁴College of Integrative Studies, Singapore Management University, Singapore, Singapore</footnote>

<footnote n="5">⁵Environmental Meteorology, Institute of Earth and Environmental Sciences, Faculty of Environment and Natural Resources, University of Freiburg, Freiburg, Germany</footnote>

<footnote n="6">⁶Foundation for Research and Technology Hellas, Institute of Applied and Computational Mathematics, Remote Sensing Lab, Heraklion, Greece</footnote>

<footnote n="7">⁷School of Earth, Atmosphere and Environment, Monash University, Melbourne, Australia</footnote>

<footnote n="8">⁸Geography and Environmental Sciences, University of Colorado, Denver, Colorado, USA</footnote>

<footnote n="9">⁹Global Institute of Sustainability and Innovation, Arizona State University, Tempe, Arizona, USA</footnote>

<footnote n="10">¹⁰UK Centre for Ecology & Hydrology, Crowmarsh Gifford, Wallingford, United Kingdom</footnote>

<footnote n="11">¹¹Department of Meteorology and Climatology, University of Lodz, Lodz, Poland</footnote>

<footnote n="12">¹²Department of Meteorology and Air Quality, Wageningen University, Wageningen, the Netherlands
¹³Korea Environment Institute, Sejong, Republic of Korea</footnote>

<footnote n="14">¹⁴Ecosystem-Atmosphere Process Lab, Department of Atmospheric Sciences,
Yonsei University, Seoul, Republic of Korea</footnote>

<footnote n="15">¹⁵Institute for Atmospheric and Earth System Research/Physics, Faculty of Science,
University of Helsinki, Helsinki, Finland</footnote>

<footnote n="16">¹⁶Ecosystem-Atmosphere Process Lab, Department of Atmospheric Sciences,
Yonsei University, Seoul, Republic of Korea</footnote>

<footnote n="17">¹⁷National Institute of Meteorological Sciences, Korea Meteorological Administration,
Jeju, Republic of Korea</footnote>

<footnote n="18">¹⁸Institut Pierre Simon Laplace (IPSL), CNRS, École Polytechnique, Institut Polytechnique de Paris,
91128 Palaiseau CEDEX, France</footnote>

<footnote n="19">¹⁹Ecosystem-Atmosphere Process Lab, Department of Atmospheric Sciences,
Yonsei University, Seoul, Republic of Korea</footnote>

<footnote n="20">²⁰Centre National de Recherches Météorologiques, University of Toulouse,
Météo-France and CNRS, Toulouse, France</footnote>

<footnote n="21">²¹Department of Geography and Earth Research Institute, University of California, Santa Barbara, USA</footnote>

<footnote n="22">²²Environmental Meteorology, Albert-Ludwigs-University, Freiburg, Germany</footnote>

<footnote n="23">²³Department of Meteorology and Climatology, University of Lodz, Lodz, Poland</footnote>

<footnote n="24">²⁴Department of Geography, National University of Singapore, Singapore, Singapore</footnote>

<footnote n="25">²⁵Department of Earth and Ocean Sciences, National Defense Academy of Japan, Yokosuka, Japan</footnote>

<footnote n="26">²⁶School of Earth, Atmosphere and Environment, Monash University, Melbourne, Australia</footnote>
<sup>27</sup>Molina Center for Energy and the Environment (MCE2), Boston, USA

<sup>28</sup>Department of Atmospheric and Cryospheric Sciences, University of Innsbruck, Innsbruck, Austria

**Correspondence:** Mathew Lipson (m.lipson@unsw.edu.au)

Received: 15 February 2022 – Discussion started: 3 June 2022

Revised: 23 September 2022 – Accepted: 22 October 2022 – Published: 22 November 2022

**Abstract.** A total of 20 urban neighbourhood-scale eddy covariance flux tower datasets are made openly available after being harmonized to create a 50 site-year collection with broad diversity in climate and urban surface characteristics. Variables needed as inputs for land surface models (incoming radiation, temperature, humidity, air pressure, wind and precipitation) are quality controlled, gap-filled and prepended with 10 years of reanalysis-derived local data, enabling an extended spin up to equilibrate models with local climate conditions. For both gap filling and spin up, ERA5 reanalysis meteorological data are bias corrected using tower-based observations, accounting for diurnal, seasonal and local urban effects not modelled in ERA5. The bias correction methods developed perform well compared to methods used in other datasets (e.g. WFDE5 or FLUXNET2015). Other variables (turbulent and upwelling radiation fluxes) are harmonized and quality controlled without gap filling. Site description metadata include local land cover fractions (buildings, roads, trees, grass etc.), building height and morphology, aerodynamic roughness estimates, population density and satellite imagery. This open collection can help extend our understanding of urban environmental processes through observational synthesis studies or in the evaluation of land surface environmental models in a wide range of urban settings. These data can be accessed from https://doi.org/10.5281/zenodo.7104984 (Lipson et al., 2022).

# 1 Background

Tower mounted instruments allow the measurement of land-atmosphere fluxes (e.g. energy, momentum, water, carbon) and local meteorological conditions. These observations are one of the fundamental ways of improving both our understanding and ability to predict biogeophysical and weather-related processes at local scales. Regional and global networks of flux tower sites have helped extend our knowledge of ecosystem and climate science (Novick et al., 2018; Beringer et al., 2016; Yamamoto et al., 2005; Valentini, 2003). Over the last 25 years, networks such as FLUXNET have progressively increased access to flux data through open-source collections (Pastorello et al., 2020), extending the reach and impact of individual site observations through synthesis studies (Baldocchi, 2020) and multi-site environmental modelling and model evaluation projects (Best et al., 2015; Ukkola et al., 2022). However, with few urban sites included, urban areas have not benefited from the improved understanding or more extensive model evaluations that these collections can facilitate.

Urban areas are unique ecosystems, distinct from natural or rural landscapes. First, most people live in cities (UN, 2019) and infrastructure is concentrated within them. Therefore, climate-related health and economic impacts fall disproportionately within urban areas. Second, urban infrastructure (e.g. buildings and roads) along with transient human activities (e.g. energy consumption and irrigation) fundamentally alter surface energy, water and mass exchanges with the atmosphere, modifying local and larger-scale environmental

conditions (Oke et al., 2017). Third, as built environments, urban areas are uniquely capable of actively mitigating and adapting to climate change.

Establishing and maintaining long-term flux sites in cities is particularly challenging because of the rarity of appropriate sites with homogenous fetch, the difficulty in gaining approval to access existing towers (e.g. for telecommunications), the cost of constructing tall towers over an aerodynamically rough surface and extremely limited long-term funding opportunities (Arnfield, 2003; Grimmond, 2006; Velasco and Roth, 2010; Feigenwinter et al., 2012; Grimmond and Ward, 2021). Thus, despite the diversity and importance of urban areas across the globe, urban flux tower data are relatively scarce, generally of short duration and rarely open source. Databases identifying urban observational programmes exist (e.g. the Urban Flux Network (Grimmond and Christen, 2012)), however urban flux tower datasets have not previously been brought together into a harmonized, gap-filled, open access collection.

We bring together quality-controlled data from 20 urban sites in an open collection that includes 50 observation years (Lipson et al., 2022). The sites are chosen to be diverse in both regional climates and urban characteristics. As evaluating land surface models is one key application for these data, we create continuous forcing datasets (i.e. with incoming radiation fluxes and other meteorological data) that are gap filled using site specific, bias corrected reanalysis data. Observations are also prepended with 10 years of site-specific reanalysis-derived meteorological data to allow modelled soil moisture and other conditions to equilibrate
<sup>27</sup>Molina Center for Energy and the Environment (MCE2), Boston, USA

<sup>28</sup>Department of Atmospheric and Cryospheric Sciences, University of Innsbruck, Innsbruck, Austria

**Correspondence:** Mathew Lipson (m.lipson@unsw.edu.au)

Received: 15 February 2022 – Discussion started: 3 June 2022

Revised: 23 September 2022 – Accepted: 22 October 2022 – Published: 22 November 2022

**Abstract.** A total of 20 urban neighbourhood-scale eddy covariance flux tower datasets are made openly available after being harmonized to create a 50 site-year collection with broad diversity in climate and urban surface characteristics. Variables needed as inputs for land surface models (incoming radiation, temperature, humidity, air pressure, wind and precipitation) are quality controlled, gap-filled and prepended with 10 years of reanalysis-derived local data, enabling an extended spin up to equilibrate models with local climate conditions. For both gap filling and spin up, ERA5 reanalysis meteorological data are bias corrected using tower-based observations, accounting for diurnal, seasonal and local urban effects not modelled in ERA5. The bias correction methods developed perform well compared to methods used in other datasets (e.g. WFDE5 or FLUXNET2015). Other variables (turbulent and upwelling radiation fluxes) are harmonized and quality controlled without gap filling. Site description metadata include local land cover fractions (buildings, roads, trees, grass etc.), building height and morphology, aerodynamic roughness estimates, population density and satellite imagery. This open collection can help extend our understanding of urban environmental processes through observational synthesis studies or in the evaluation of land surface environmental models in a wide range of urban settings. These data can be accessed from https://doi.org/10.5281/zenodo.7104984 (Lipson et al., 2022).

# 1 Background

Tower mounted instruments allow the measurement of land-atmosphere fluxes (e.g. energy, momentum, water, carbon) and local meteorological conditions. These observations are one of the fundamental ways of improving both our understanding and ability to predict biogeophysical and weather-related processes at local scales. Regional and global networks of flux tower sites have helped extend our knowledge of ecosystem and climate science (Novick et al., 2018; Beringer et al., 2016; Yamamoto et al., 2005; Valentini, 2003). Over the last 25 years, networks such as FLUXNET have progressively increased access to flux data through open-source collections (Pastorello et al., 2020), extending the reach and impact of individual site observations through synthesis studies (Baldocchi, 2020) and multi-site environmental modelling and model evaluation projects (Best et al., 2015; Ukkola et al., 2022). However, with few urban sites included, urban areas have not benefited from the improved understanding or more extensive model evaluations that these collections can facilitate.

Urban areas are unique ecosystems, distinct from natural or rural landscapes. First, most people live in cities (UN, 2019) and infrastructure is concentrated within them. Therefore, climate-related health and economic impacts fall disproportionately within urban areas. Second, urban infrastructure (e.g. buildings and roads) along with transient human activities (e.g. energy consumption and irrigation) fundamentally alter surface energy, water and mass exchanges with the atmosphere, modifying local and larger-scale environmental

conditions (Oke et al., 2017). Third, as built environments, urban areas are uniquely capable of actively mitigating and adapting to climate change.

Establishing and maintaining long-term flux sites in cities is particularly challenging because of the rarity of appropriate sites with homogenous fetch, the difficulty in gaining approval to access existing towers (e.g. for telecommunications), the cost of constructing tall towers over an aerodynamically rough surface and extremely limited long-term funding opportunities (Arnfield, 2003; Grimmond, 2006; Velasco and Roth, 2010; Feigenwinter et al., 2012; Grimmond and Ward, 2021). Thus, despite the diversity and importance of urban areas across the globe, urban flux tower data are relatively scarce, generally of short duration and rarely open source. Databases identifying urban observational programmes exist (e.g. the Urban Flux Network (Grimmond and Christen, 2012)), however urban flux tower datasets have not previously been brought together into a harmonized, gap-filled, open access collection.

We bring together quality-controlled data from 20 urban sites in an open collection that includes 50 observation years (Lipson et al., 2022). The sites are chosen to be diverse in both regional climates and urban characteristics. As evaluating land surface models is one key application for these data, we create continuous forcing datasets (i.e. with incoming radiation fluxes and other meteorological data) that are gap filled using site specific, bias corrected reanalysis data. Observations are also prepended with 10 years of site-specific reanalysis-derived meteorological data to allow modelled soil moisture and other conditions to equilibrate
with local climate conditions during model spin up. These
data can be used to drive land surface models offline at a
single grid point. Other variables (turbulent and upwelling
radiation fluxes) can be used to evaluative models in simulat-
ing land-atmosphere energy exchanges, or in observational
synthesis studies.

Along with the meteorological data, site characteristics and metadata are provided in a common format. The metadata include tower location, land cover fractions, building heights and morphology, aerodynamic roughness parameter estimates, population density, estimated anthropogenic heat fluxes, site photos and satellite imagery. This collection can help extend our ability to model and understanding of environmental processes in different urban settings.

## 2 Methods

### 2.1 Site selection

The initial motivation for collating these flux tower and site data is for use in the Urban-PLUMBER multi-site model evaluation project, currently underway (Lipson, 2021). Urban-PLUMBER draws on methods from the first international urban land surface model comparison (Grimmond et al., 2010, 2011) and the Protocol for the Analysis of Land Surface Models Benchmarking Evaluation Project (PLUMBER (Best et al., 2015)). The latter evaluated land surface models in non-urban (vegetated) areas, while Urban-PLUMBER evaluates land surface models at 20 urban sites (Fig. 1, Table 1).

There is a two-fold use of these observational data in the
model evaluation of Urban-PLUMBER:

1. To provide local-scale meteorological input forcing to drive land surface models;

2. To evaluate the performance of models, primarily assessing the local-scale exchange of radiant and turbulent heat fluxes between the surface and lower atmosphere.

With these objectives in mind, the following criteria are used
to select flux tower sites:

- appropriately sited for neighbourhood-scale conditions – i.e. within the inertial sub-layer, typically 2–5 times above the average building height and with relatively homogenous fetch (Grimmond, 2006; Barlow, 2014; Grimmond and Ward, 2021);

- requested observations available at 30 or 60 min resolution (Table 2);

- local site characteristics available for description and configuring models;

- a preference for longer datasets (as this allows seasonal and inter-annual variability to be included);

- collectively represent a diverse range of site characteristics and climates.

Potential sites are identified from published site lists (Grimmond and Christen, 2012; Oke et al., 2017) and open calls for data (e.g. community newsletters (Lipson et al., 2020a), international conferences ((Lipson et al., 2020b, c) and social media professional networks). We deemed 20 sites sufficient for the evaluation project (Table 1), together covering 50 site-years. Included sites have built fractions (i.e. plan area fraction of all impervious surfaces including roofs, roads, other paving etc.) from 0.05 to 0.965, and are located in four major Köppen-Geiger (Beck et al., 2018) climate classes (Fig. 1). Eleven sites are in temperate climates, eight in cold (or continental) climates, and one in each of tropical and arid climates.

Sites are reasonably distributed across mean temperature and precipitation for global urban locations, but gaps remain, particularly in warm, wet and very cold climates (Fig. 2). Some urban flux observations in understudied regions were not included (e.g. Ouagadougou (Offerle et al., 2005), São Paulo (Ferreira et al., 2013), Guangzhou (Shi et al., 2019), Beijing (Dou et al., 2019)) because they do not meet the model evaluation project needs because of the relatively short observed periods for the available data. These regions and climates have large urban populations with significant environmental challenges and have few urban flux tower sites compared with Northern Hemisphere temperate or continental locations (Grimmond, 2006; Roth et al., 2017). Understudied regions and climates should be included in future collections when appropriate time series become available.

### 2.2 Flux tower data

The observed data are provided in 30 or 60 min periods (Table 1), processed from high-frequency samples by individual observing groups. In the harmonized collection, timestamps are in coordinated universal time (UTC) indicating the end of the measurement period. Variables names and units use ALMA (Assistance for Land-surface Modelling Activities) conventions, a format used in previous land surface model comparisons.

Data are cleaned (Sect. 2.3: *Quality control*), forcing variables are gap filled (Sect. 2.4) and prepended with data derived from ERA5 (Sect. 2.5), after site-specific corrections (Sect. 2.6). An example of final prepended and gap filled data is shown in Fig. 3 for one site (UK-KingsCollege). Plots for all other variables and sites are also available in the collection (Lipson et al., 2022). Turbulent fluxes and upwelling radiation fluxes are not gap filled (Table 2).

Data are split into forcing and analysis variable sets (Table 2) to allow the forcing variables to be provided to modelling groups as input to run their models. The withheld analysis data are used by the coordinating group to assess the model outputs.
with local climate conditions during model spin up. These data can be used to drive land surface models offline at a single grid point. Other variables (turbulent and upwelling radiation fluxes) can be used to evaluative models in simulating land-atmosphere energy exchanges, or in observational synthesis studies.

Along with the meteorological data, site characteristics and metadata are provided in a common format. The metadata include tower location, land cover fractions, building heights and morphology, aerodynamic roughness parameter estimates, population density, estimated anthropogenic heat fluxes, site photos and satellite imagery. This collection can help extend our ability to model and understanding of environmental processes in different urban settings.

# 2 Methods

## 2.1 Site selection

The initial motivation for collating these flux tower and site data is for use in the Urban-PLUMBER multi-site model evaluation project, currently underway (Lipson, 2021). Urban-PLUMBER draws on methods from the first international urban land surface model comparison (Grimmond et al., 2010, 2011) and the Protocol for the Analysis of Land Surface Models Benchmarking Evaluation Project (PLUMBER (Best et al., 2015)). The latter evaluated land surface models in non-urban (vegetated) areas, while Urban-PLUMBER evaluates land surface models at 20 urban sites (Fig. 1, Table 1).

There is a two-fold use of these observational data in the model evaluation of Urban-PLUMBER:

1. To provide local-scale meteorological input forcing to drive land surface models;

2. To evaluate the performance of models, primarily assessing the local-scale exchange of radiant and turbulent heat fluxes between the surface and lower atmosphere.

With these objectives in mind, the following criteria are used to select flux tower sites:

- appropriately sited for neighbourhood-scale conditions – i.e. within the inertial sub-layer, typically 2–5 times above the average building height and with relatively homogenous fetch (Grimmond, 2006; Barlow, 2014; Grimmond and Ward, 2021);

- requested observations available at 30 or 60 min resolution (Table 2);

- local site characteristics available for description and configuring models;

- a preference for longer datasets (as this allows seasonal and inter-annual variability to be included);

- collectively represent a diverse range of site characteristics and climates.

Potential sites are identified from published site lists (Grimmond and Christen, 2012; Oke et al., 2017) and open calls for data (e.g. community newsletters (Lipson et al., 2020a), international conferences ((Lipson et al., 2020b, c) and social media professional networks). We deemed 20 sites sufficient for the evaluation project (Table 1), together covering 50 site-years. Included sites have built fractions (i.e. plan area fraction of all impervious surfaces including roofs, roads, other paving etc.) from 0.05 to 0.965, and are located in four major Köppen-Geiger (Beck et al., 2018) climate classes (Fig. 1). Eleven sites are in temperate climates, eight in cold (or continental) climates, and one in each of tropical and arid climates.

Sites are reasonably distributed across mean temperature and precipitation for global urban locations, but gaps remain, particularly in warm, wet and very cold climates (Fig. 2). Some urban flux observations in understudied regions were not included (e.g. Ouagadougou (Offerle et al., 2005), São Paulo (Ferreira et al., 2013), Guangzhou (Shi et al., 2019), Beijing (Dou et al., 2019)) because they do not meet the model evaluation project needs because of the relatively short observed periods for the available data. These regions and climates have large urban populations with significant environmental challenges and have few urban flux tower sites compared with Northern Hemisphere temperate or continental locations (Grimmond, 2006; Roth et al., 2017). Understudied regions and climates should be included in future collections when appropriate time series become available.

## 2.2 Flux tower data

The observed data are provided in 30 or 60 min periods (Table 1), processed from high-frequency samples by individual observing groups. In the harmonized collection, timestamps are in coordinated universal time (UTC) indicating the end of the measurement period. Variables names and units use ALMA (Assistance for Land-surface Modelling Activities) conventions, a format used in previous land surface model comparisons.

Data are cleaned (Sect. 2.3: *Quality control*), forcing variables are gap filled (Sect. 2.4) and prepended with data derived from ERA5 (Sect. 2.5), after site-specific corrections (Sect. 2.6). An example of final prepended and gap filled data is shown in Fig. 3 for one site (UK-KingsCollege). Plots for all other variables and sites are also available in the collection (Lipson et al., 2022). Turbulent fluxes and upwelling radiation fluxes are not gap filled (Table 2).

Data are split into forcing and analysis variable sets (Table 2) to allow the forcing variables to be provided to modelling groups as input to run their models. The withheld analysis data are used by the coordinating group to assess the model outputs.
<table><caption>**Table 1.** Site location and included observation (focus) period. Data providers may have longer observation periods available than are in this collection. Resolution is 30 min (or 60 min if denoted by *). All periods in universal time coordinated (UTC). US-Minneapolis data are split based on wind direction and fetch (Sect. 4.5).</caption><thead><tr><th>Site name</th><th>City</th><th>Country</th><th>Observed period</th><th>Latitude</th><th>Longitude</th><th>References</th></tr></thead><tbody><tr><td>AU-Preston</td><td>Melbourne</td><td>Australia</td><td>Aug 2003–Nov 2004</td><td>-37.7306</td><td>145.0145</td><td>Coutts et al. (2007a, b)</td></tr><tr><td>AU-SurreyHills</td><td>Melbourne</td><td>Australia</td><td>Feb 2004–Jul 2004</td><td>-37.8265</td><td>145.099</td><td>Coutts et al. (2007a, b)</td></tr><tr><td>CA-Sunset</td><td>Vancouver</td><td>Canada</td><td>Jan 2012–Dec 2016</td><td>49.2261</td><td>-123.078</td><td>Christen et al. (2011); Crawford and Christen (2015)</td></tr><tr><td>FI-Kumpula</td><td>Helsinki</td><td>Finland</td><td>Dec 2010–Dec 2013</td><td>60.2028</td><td>24.9611</td><td>Karsisto et al. (2016)</td></tr><tr><td>FI-Torni</td><td>Helsinki</td><td>Finland</td><td>Dec 2010–Dec 2013</td><td>60.1678</td><td>24.9387</td><td>Järvi et al. (2018); Nordbo et al. (2013)</td></tr><tr><td>FR-Capitole</td><td>Toulouse</td><td>France</td><td>Feb 2004–Mar 2005</td><td>43.6035</td><td>1.4454</td><td>Masson et al. (2008); Goret et al. (2019)</td></tr><tr><td>GR-HECKOR</td><td>Heraklion</td><td>Greece</td><td>Jun 2019–Jun 2020</td><td>35.3361</td><td>25.1328</td><td>Stagakis et al. (2019)</td></tr><tr><td>JP-Yoyogi</td><td>Tokyo</td><td>Japan</td><td>Mar 2016–Mar 2020*</td><td>35.6645</td><td>139.6845</td><td>Hirano et al. (2015); Ishidoya et al. (2020)</td></tr><tr><td>KR-Jungnang</td><td>Seoul</td><td>South Korea</td><td>Jan 2017–Apr 2019</td><td>37.5907</td><td>127.0794</td><td>Hong et al. (2020); Hong et al. (2022)</td></tr><tr><td>KR-Ochang</td><td>Ochang</td><td>South Korea</td><td>Jun 2015–Jul 2017</td><td>36.7197</td><td>127.4344</td><td>Hong et al. (2019, 2020)</td></tr><tr><td>MX-Escandon</td><td>Mexico City</td><td>Mexico</td><td>Jun 2011–Sep 2012</td><td>19.4042</td><td>-99.1761</td><td>Velasco et al. (2011, 2014)</td></tr><tr><td>NL-Amsterdam</td><td>Amsterdam</td><td>Netherlands</td><td>Jan 2019–Oct 2020</td><td>52.3665</td><td>4.8929</td><td>Steeneveld et al. (2020)</td></tr><tr><td>PL-Lipowa</td><td>Łódź</td><td>Poland</td><td>Jan 2008–Dec 2012*</td><td>51.7625</td><td>19.4453</td><td>Fortuniak et al. (2013); Pawlak et al. (2011)</td></tr><tr><td>PL-Narutowicza</td><td>Łódź</td><td>Poland</td><td>Jan 2008–Dec 2012*</td><td>51.7733</td><td>19.4811</td><td>Fortuniak et al. (2013, 2006)</td></tr><tr><td>SG-TelokKurau06</td><td>Singapore</td><td>Singapore</td><td>Apr 2006–Mar 2007</td><td>1.3143</td><td>103.9112</td><td>Roth et al. (2017)</td></tr><tr><td>UK-KingsCollege</td><td>London</td><td>UK</td><td>Apr 2012–Jan 2014</td><td>51.5118</td><td>-0.1167</td><td>Bjorkegren et al. (2015); Kotthaus and Grimmond (2014a, b)</td></tr><tr><td>UK-Swindon</td><td>Swindon</td><td>UK</td><td>May 2011–Apr 2013</td><td>51.5846</td><td>-1.7981</td><td>Ward et al. (2013)</td></tr><tr><td>US-Baltimore</td><td>Baltimore</td><td>USA</td><td>Jan 2002–Jan 2007*</td><td>39.4128</td><td>-76.5215</td><td>Crawford et al. (2011)</td></tr><tr><td>US-Minneapolis</td><td>Minneapolis</td><td>USA</td><td>Jun 2006–May 2009</td><td>44.9984</td><td>-93.1884</td><td>Peters et al. (2011); Menzer and McFadden (2017)</td></tr><tr><td>US-WestPhoenix</td><td>Phoenix</td><td>USA</td><td>Dec 2011–Jan 2013</td><td>44.9984</td><td>-93.1884</td><td>Chow (2017); Chow et al. (2014)</td></tr></tbody></table>
<table><caption>**Table 1.** Site location and included observation (focus) period. Data providers may have longer observation periods available than are in this collection. Resolution is 30 min (or 60 min if denoted by *). All periods in universal time coordinated (UTC). US-Minneapolis data are split based on wind direction and fetch (Sect. 4.5).</caption><thead><tr><th>Site name</th><th>City</th><th>Country</th><th>Observed period</th><th>Latitude</th><th>Longitude</th><th>References</th></tr></thead><tbody><tr><td>AU-Preston</td><td>Melbourne</td><td>Australia</td><td>Aug 2003–Nov 2004</td><td>-37.7306</td><td>145.0145</td><td>Coutts et al. (2007a, b)</td></tr><tr><td>AU-SurreyHills</td><td>Melbourne</td><td>Australia</td><td>Feb 2004–Jul 2004</td><td>-37.8265</td><td>145.099</td><td>Coutts et al. (2007a, b)</td></tr><tr><td>CA-Sunset</td><td>Vancouver</td><td>Canada</td><td>Jan 2012–Dec 2016</td><td>49.2261</td><td>-123.078</td><td>Christen et al. (2011); Crawford and Christen (2015)</td></tr><tr><td>FI-Kumpula</td><td>Helsinki</td><td>Finland</td><td>Dec 2010–Dec 2013</td><td>60.2028</td><td>24.9611</td><td>Karsisto et al. (2016)</td></tr><tr><td>FI-Torni</td><td>Helsinki</td><td>Finland</td><td>Dec 2010–Dec 2013</td><td>60.1678</td><td>24.9387</td><td>Järvi et al. (2018); Nordbo et al. (2013)</td></tr><tr><td>FR-Capitole</td><td>Toulouse</td><td>France</td><td>Feb 2004–Mar 2005</td><td>43.6035</td><td>1.4454</td><td>Masson et al. (2008); Goret et al. (2019)</td></tr><tr><td>GR-HECKOR</td><td>Heraklion</td><td>Greece</td><td>Jun 2019–Jun 2020</td><td>35.3361</td><td>25.1328</td><td>Stagakis et al. (2019)</td></tr><tr><td>JP-Yoyogi</td><td>Tokyo</td><td>Japan</td><td>Mar 2016–Mar 2020*</td><td>35.6645</td><td>139.6845</td><td>Hirano et al. (2015); Ishidoya et al. (2020)</td></tr><tr><td>KR-Jungnang</td><td>Seoul</td><td>South Korea</td><td>Jan 2017–Apr 2019</td><td>37.5907</td><td>127.0794</td><td>Hong et al. (2020); Hong et al. (2022)</td></tr><tr><td>KR-Ochang</td><td>Ochang</td><td>South Korea</td><td>Jun 2015–Jul 2017</td><td>36.7197</td><td>127.4344</td><td>Hong et al. (2019, 2020)</td></tr><tr><td>MX-Escandon</td><td>Mexico City</td><td>Mexico</td><td>Jun 2011–Sep 2012</td><td>19.4042</td><td>-99.1761</td><td>Velasco et al. (2011, 2014)</td></tr><tr><td>NL-Amsterdam</td><td>Amsterdam</td><td>Netherlands</td><td>Jan 2019–Oct 2020</td><td>52.3665</td><td>4.8929</td><td>Steeneveld et al. (2020)</td></tr><tr><td>PL-Lipowa</td><td>Łódź</td><td>Poland</td><td>Jan 2008–Dec 2012*</td><td>51.7625</td><td>19.4453</td><td>Fortuniak et al. (2013); Pawlak et al. (2011)</td></tr><tr><td>PL-Narutowicza</td><td>Łódź</td><td>Poland</td><td>Jan 2008–Dec 2012*</td><td>51.7733</td><td>19.4811</td><td>Fortuniak et al. (2013, 2006)</td></tr><tr><td>SG-TelokKurau06</td><td>Singapore</td><td>Singapore</td><td>Apr 2006–Mar 2007</td><td>1.3143</td><td>103.9112</td><td>Roth et al. (2017)</td></tr><tr><td>UK-KingsCollege</td><td>London</td><td>UK</td><td>Apr 2012–Jan 2014</td><td>51.5118</td><td>-0.1167</td><td>Bjorkegren et al. (2015); Kotthaus and Grimmond (2014a, b)</td></tr><tr><td>UK-Swindon</td><td>Swindon</td><td>UK</td><td>May 2011–Apr 2013</td><td>51.5846</td><td>-1.7981</td><td>Ward et al. (2013)</td></tr><tr><td>US-Baltimore</td><td>Baltimore</td><td>USA</td><td>Jan 2002–Jan 2007*</td><td>39.4128</td><td>-76.5215</td><td>Crawford et al. (2011)</td></tr><tr><td>US-Minneapolis</td><td>Minneapolis</td><td>USA</td><td>Jun 2006–May 2009</td><td>44.9984</td><td>-93.1884</td><td>Peters et al. (2011); Menzer and McFadden (2017)</td></tr><tr><td>US-WestPhoenix</td><td>Phoenix</td><td>USA</td><td>Dec 2011–Jan 2013</td><td>44.9984</td><td>-93.1884</td><td>Chow (2017); Chow et al. (2014)</td></tr></tbody></table>
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Location of flux tower sites in this collection. Each site Köppen–Geiger climate classification (Beck et al., 2018) and the built land fraction around the tower are indicated at the bottom of the figure.</figcaption></figure>

<figure><img src="image_3.png" /><figcaption>**Figure 2.** Climatology of included sites compared with more than 70 000 global urban areas. Mean temperature and annual precipitation at the 20 tower sites (red, truncated site name, Table 1) from tower observations; global urban locations (grey) from ERA5 surface data (Hersbach et al., 2020, 2018) (2000–2010) from grid nearest to locations identified in the Global Rural–Urban Mapping Project (GRUMP) (Center for International Earth Science Information Network – CIESIN – Columbia University et al., 2017). Locations with rainfall above 3000 mm yr⁻¹ (13 % of locations) and mean temperature below –3 °C (0.2 % of locations) are not shown.</figcaption></figure>

Some additional observed variables (Table 3) have, where practical, been included in the datasets after passing through the quality control steps. Missing forcing variables are obtained using bias corrected reanalysis data (Sect. 2.4). No gap filling is applied to analysis data or additional variables.

## 2.3 Quality control and assurance

For each site the 30 or 60 min variables are calculated by data providers from high-frequency samples after applying their own quality control measures (e.g. Aubinet et al., 2012; Feigenwinter et al., 2012; Kotthaus and Grimmond, 2012; Vitale et al., 2020). The harmonized collection consists of the data retained after undergoing five additional quality control steps, in the following order:

1. *Out-of-range.* Removal of unphysical values (e.g. negative shortwave radiation) using the ALMA expected range protocol (Bowling and Polcher, 2001).

2. *Night.* Nocturnal shortwave radiation set to zero, based on civil twilight (when the sun is 6° below the horizon (Forsythe et al., 1995)).

3. *Constant.* Four or more time steps with identical values (excluding zero values for shortwave radiation, rainfall and snowfall) are removed as suspicious.

4. *Outlier.* Values outside ±4 standard deviations for each hour in a rolling 30 d window (to account for diurnal and seasonal variations) removed. Repeat with a larger tolerance (±5 SD, standard deviations) until no outliers remain (Schmid et al., 2000; Vickers and Mahrt, 1997). The outlier test is not applied to precipitation.
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Location of flux tower sites in this collection. Each site Köppen–Geiger climate classification (Beck et al., 2018) and the built land fraction around the tower are indicated at the bottom of the figure.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Climatology of included sites compared with more than 70 000 global urban areas. Mean temperature and annual precipitation at the 20 tower sites (red, truncated site name, Table 1) from tower observations; global urban locations (grey) from ERA5 surface data (Hersbach et al., 2020, 2018) (2000–2010) from grid nearest to locations identified in the Global Rural–Urban Mapping Project (GRUMP) (Center for International Earth Science Information Network – CIESIN – Columbia University et al., 2017). Locations with rainfall above 3000 mm yr⁻¹ (13 % of locations) and mean temperature below –3 °C (0.2 % of locations) are not shown.</figcaption></figure>

Some additional observed variables (Table 3) have, where practical, been included in the datasets after passing through the quality control steps. Missing forcing variables are obtained using bias corrected reanalysis data (Sect. 2.4). No gap filling is applied to analysis data or additional variables.

## 2.3 Quality control and assurance

For each site the 30 or 60 min variables are calculated by data providers from high-frequency samples after applying their own quality control measures (e.g. Aubinet et al., 2012; Feigenwinter et al., 2012; Kotthaus and Grimmond, 2012; Vitale et al., 2020). The harmonized collection consists of the data retained after undergoing five additional quality control steps, in the following order:

1. *Out-of-range.* Removal of unphysical values (e.g. negative shortwave radiation) using the ALMA expected range protocol (Bowling and Polcher, 2001).

2. *Night.* Nocturnal shortwave radiation set to zero, based on civil twilight (when the sun is 6° below the horizon (Forsythe et al., 1995)).

3. *Constant.* Four or more time steps with identical values (excluding zero values for shortwave radiation, rainfall and snowfall) are removed as suspicious.

4. *Outlier.* Values outside ±4 standard deviations for each hour in a rolling 30 d window (to account for diurnal and seasonal variations) removed. Repeat with a larger tolerance (±5 SD, standard deviations) until no outliers remain (Schmid et al., 2000; Vickers and Mahrt, 1997). The outlier test is not applied to precipitation.
<figure><img src="image_4.png" /></figure>

**Figure 3.** Forcing time series with gap filling. Example shown for air temperature (Tair) at Kings College, London (UK-KingsCollege); (a) full forcing period, including 10 years of ERA5-derived data (red) prior to observations (black) used for model spin up, and (b) focus period used for analysis. Gaps are first filled from nearby tower measurements where available and short gaps (<= 2 h) are linearly interpolated (blue). Remaining gaps are filled using the ERA5-derived time series which is seasonally and diurnally bias corrected using site observations. White lines show 7 d mean values. Similar plots are available for other sites within the site data collection (Lipson et al., 2022).

5. *Visual.* Remaining suspect readings are removed manually via visual inspection.

These steps are undertaken in the processing script `qc_observations.py` (see Sect. 5: Code availability), including periods identified through visual inspection (21 instances across all data). Data removed through quality control are indicated in plots of each variable at each site included in the data collection (Lipson et al., 2022). Note that quality control steps which eliminate observations at particular times (e.g. at night or after rainfall) can introduce biases (Grimmond, 2006). In addition, the outlier check values (step 4) are somewhat arbitrary (as noted in Vickers and Mahrt, 1997). Therefore, we also provide the “raw” observations (prior to quality control discussed here) in the collection as they may be more appropriate for some types of analyses.

Communication, or human errors, also have the potential to degrade or invalidate data (Menard et al., 2021). As part of quality assurance, project coordinators prepared an observational data protocol (Lipson et al., 2021) to explicitly set out requirements for data providers prior to submission of their data. The protocol documented instrument siting requirements, variables and data formats, dataset length and resolution, necessary site characteristic information and metadata, as well as the expectations for data handling, use and au-

thorship. On receiving data, coordinators undertook further checks and identified errors that were not picked up by automated quality control. Identified errors included mislabelled variables and metadata, inconsistent timestamps and unit discrepancies. Many of the errors were identified by comparing provided data with secondary sources such as ERA5, nearby meteorological stations or previous publications. Errors were corrected collaboratively with data providers, some leading to corrections in primary data sources.

## 2.4 Gap filling

Three gap filling methods are used to create a continuous dataset for forcing variables, in the following order:

- contemporaneous and nearby flux tower or weather observing sites (where available from data providers);

- small gaps (<= 2 h) are linearly interpolated from the adjoining observations;

- larger gaps and a 10-year spin up period are filled with bias corrected ERA5 data (Sect. 2.6).

As only one site provided observed snowfall rate (JP-Yoyogi), ERA5 snowfall rates are used for all periods at other sites. At those sites the additional water equivalent
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Forcing time series with gap filling. Example shown for air temperature (Tair) at Kings College, London (UK-KingsCollege); (a) full forcing period, including 10 years of ERA5-derived data (red) prior to observations (black) used for model spin up, and (b) focus period used for analysis. Gaps are first filled from nearby tower measurements where available and short gaps (<= 2 h) are linearly interpolated (blue). Remaining gaps are filled using the ERA5-derived time series which is seasonally and diurnally bias corrected using site observations. White lines show 7 d mean values. Similar plots are available for other sites within the site data collection (Lipson et al., 2022).</figcaption></figure>

5. *Visual*. Remaining suspect readings are removed manually via visual inspection.

These steps are undertaken in the processing script `qc_observations.py` (see Sect. 5: Code availability), including periods identified through visual inspection (21 instances across all data). Data removed through quality control are indicated in plots of each variable at each site included in the data collection (Lipson et al., 2022). Note that quality control steps which eliminate observations at particular times (e.g. at night or after rainfall) can introduce biases (Grimmond, 2006). In addition, the outlier check values (step 4) are somewhat arbitrary (as noted in Vickers and Mahrt, 1997). Therefore, we also provide the “raw” observations (prior to quality control discussed here) in the collection as they may be more appropriate for some types of analyses.

Communication, or human errors, also have the potential to degrade or invalidate data (Menard et al., 2021). As part of quality assurance, project coordinators prepared an observational data protocol (Lipson et al., 2021) to explicitly set out requirements for data providers prior to submission of their data. The protocol documented instrument siting requirements, variables and data formats, dataset length and resolution, necessary site characteristic information and metadata, as well as the expectations for data handling, use and au-

thorship. On receiving data, coordinators undertook further checks and identified errors that were not picked up by automated quality control. Identified errors included mislabelled variables and metadata, inconsistent timestamps and unit discrepancies. Many of the errors were identified by comparing provided data with secondary sources such as ERA5, nearby meteorological stations or previous publications. Errors were corrected collaboratively with data providers, some leading to corrections in primary data sources.

## 2.4 Gap filling

Three gap filling methods are used to create a continuous dataset for forcing variables, in the following order:

- contemporaneous and nearby flux tower or weather observing sites (where available from data providers);

- small gaps (<= 2 h) are linearly interpolated from the adjoining observations;

- larger gaps and a 10-year spin up period are filled with bias corrected ERA5 data (Sect. 2.6).

As only one site provided observed snowfall rate (JP-Yoyogi), ERA5 snowfall rates are used for all periods at other sites. At those sites the additional water equivalent
<table><caption>**Table 2.** Forcing and analysis flux tower data variables. Short name description, units and positive direction use ALMA data conventions. Mean annual estimates of anthropogenic heat flux are included as site metadata. Analysis and additional data are not gap filled. Ground heat flux (Qg) is the heat flux into soil rather than total storage heat flux which is difficult to measure in urban areas (Grimmond and Oke, 1999).</caption><thead><tr><th>Variable</th><th>Description</th><th>Units</th><th>Positive direction</th><th>Gap filled</th><th>Bias correction</th></tr></thead><tbody><tr><td colspan="6"><strong>Forcing data</strong></td></tr><tr><td>SWdown</td><td>Downward shortwave radiation</td><td>W m<sup>-2</sup></td><td>Downward</td><td>yes</td><td>none</td></tr><tr><td>LWdown</td><td>Downward longwave radiation</td><td>W m<sup>-2</sup></td><td>Downward</td><td>yes</td><td>hourly and daily</td></tr><tr><td>Tair</td><td>Air temperature</td><td>K</td><td>–</td><td>yes</td><td>hourly and daily</td></tr><tr><td>Qair</td><td>Specific humidity</td><td>kg kg<sup>-1</sup></td><td>–</td><td>yes</td><td>hourly and daily</td></tr><tr><td>PSurf</td><td>Station air pressure</td><td>Pa</td><td>–</td><td>yes</td><td>hourly and daily</td></tr><tr><td>Wind_N</td><td>Northward wind component</td><td>m s<sup>-1</sup></td><td>Northward</td><td>yes</td><td>logarithmic law</td></tr><tr><td>Wind_E</td><td>Eastward wind component</td><td>m s<sup>-1</sup></td><td>Eastward</td><td>yes</td><td>logarithmic law</td></tr><tr><td>Rainf</td><td>Rainfall rate</td><td>kg m<sup>-2</sup> s<sup>-1</sup></td><td>Downward</td><td>yes</td><td>long-term precipitation</td></tr><tr><td>Snowf</td><td>Snowfall rate</td><td>kg m<sup>-2</sup> s<sup>-1</sup></td><td>Downward</td><td>yes</td><td>long term precipitation</td></tr><tr><td colspan="6"><strong>Analysis data</strong></td></tr><tr><td>SWup</td><td>Upward shortwave radiation</td><td>W m<sup>-2</sup></td><td>Upward</td><td>no</td><td>none</td></tr><tr><td>LWup</td><td>Upward longwave radiation</td><td>W m<sup>-2</sup></td><td>Upward</td><td>no</td><td>none</td></tr><tr><td>Qle</td><td>Latent heat flux</td><td>W m<sup>-2</sup></td><td>Upward</td><td>no</td><td>none</td></tr><tr><td>Qh</td><td>Sensible heat flux</td><td>W m<sup>-2</sup></td><td>Upward</td><td>no</td><td>none</td></tr><tr><td colspan="6"><strong>Additional data (optional)</strong></td></tr><tr><td>Qg</td><td>Ground heat flux into soil</td><td>W m<sup>-2</sup></td><td>Downward</td><td>no</td><td>none</td></tr><tr><td>Qtau</td><td>Momentum flux</td><td>N m<sup>-2</sup></td><td>Downward</td><td>no</td><td>none</td></tr><tr><td>Tair2m</td><td>Near surface air temperature (2 m)</td><td>K</td><td>–</td><td>no</td><td>none</td></tr><tr><td>SoilTemp</td><td>Soil temperature (depth in metadata)</td><td>K</td><td>–</td><td>no</td><td>none</td></tr></tbody></table>

<table><caption>**Table 3.** Site climate classification missing and additional variables. Climate classification from Köppen–Geiger global dataset (Beck et al., 2018). Table 2 gives variable definitions.</caption><thead><tr><th>Site name</th><th>Class</th><th>Climate description</th><th>Missing variables</th><th>Additional variables</th></tr></thead><tbody><tr><td>AU-Preston</td><td>Cfb</td><td>Temperate, no dry season, warm summer</td><td>Snowf</td><td>Qtau</td></tr><tr><td>AU-SurreyHills</td><td>Cfb</td><td>Temperate, no dry season, warm summer</td><td>Snowf</td><td>Qtau</td></tr><tr><td>CA-Sunset</td><td>Csb</td><td>Temperate, dry summer, warm summer</td><td>Snowf</td><td>Qtau, SoilTemp</td></tr><tr><td>FI-Kumpula</td><td>Dfb</td><td>Cold, no dry season, warm summer</td><td>Snowf</td><td></td></tr><tr><td>FI-Torni</td><td>Dfb</td><td>Cold, no dry season, warm summer</td><td>Snowf</td><td></td></tr><tr><td>FR-Capitole</td><td>Cfa</td><td>Temperate, no dry season, hot summer</td><td>Snowf</td><td>Qtau</td></tr><tr><td>GR-HECKOR</td><td>Csa</td><td>Temperate, dry summer, hot summer</td><td>Snowf</td><td>Qtau</td></tr><tr><td>JP-Yoyogi</td><td>Cfa</td><td>Temperate, no dry season, hot summer</td><td></td><td></td></tr><tr><td>KR-Jungnang</td><td>Dwa</td><td>Cold, dry winter, hot summer</td><td>Snowf</td><td></td></tr><tr><td>KR-Ochang</td><td>Dwa</td><td>Cold, dry winter, hot summer</td><td>Snowf</td><td></td></tr><tr><td>MX-Escandon</td><td>Cwb</td><td>Temperate, dry winter, warm summer</td><td>Snowf, LWdown*</td><td>Qtau</td></tr><tr><td>NL-Amsterdam</td><td>Cfb</td><td>Temperate, no dry season, warm summer</td><td>Snowf</td><td>Qtau</td></tr><tr><td>PL-Lipowa</td><td>Dfb</td><td>Cold, no dry season, warm summer</td><td>Snowf</td><td></td></tr><tr><td>PL-Narutowicza</td><td>Dfb</td><td>Cold, no dry season, warm summer</td><td>Snowf</td><td></td></tr><tr><td>SG-TelokKurau06</td><td>Af</td><td>Tropical, rainforest</td><td>Snowf</td><td></td></tr><tr><td>UK-KingsCollege</td><td>Cfb</td><td>Temperate, no dry season, warm summer</td><td>Snowf</td><td>Qtau</td></tr><tr><td>UK-Swindon</td><td>Cfb</td><td>Temperate, no dry season, warm summer</td><td>Snowf</td><td>Qtau</td></tr><tr><td>US-Baltimore</td><td>Cfa</td><td>Temperate, no dry season, hot summer</td><td>Snowf</td><td>Qtau, SoilTemp</td></tr><tr><td>US-Minneapolis</td><td>Dfa</td><td>Cold, no dry season, hot summer</td><td>Snowf</td><td>Qtau, SoilTemp, Qg</td></tr><tr><td>US-WestPhoenix</td><td>BWh</td><td>Arid, desert, hot</td><td>Snowf</td><td>Qtau, SoilTemp</td></tr></tbody></table>

<footnote>* Note: as MX-Escandon LWdown data are unavailable during the 2011–2012 focus period, but were available in 2006, this earlier period is used to determine bias correction for ERA5 LWdown data.</footnote>
<table><caption>**Table 2.** Forcing and analysis flux tower data variables. Short name description, units and positive direction use ALMA data conventions. Mean annual estimates of anthropogenic heat flux are included as site metadata. Analysis and additional data are not gap filled. Ground heat flux (Qg) is the heat flux into soil rather than total storage heat flux which is difficult to measure in urban areas (Grimmond and Oke, 1999).</caption><thead><tr><th>Variable</th><th>Description</th><th>Units</th><th>Positive direction</th><th>Gap filled</th><th>Bias correction</th></tr></thead><tbody><tr><td colspan="6"><strong>Forcing data</strong></td></tr><tr><td>SWdown</td><td>Downward shortwave radiation</td><td>W m<sup>-2</sup></td><td>Downward</td><td>yes</td><td>none</td></tr><tr><td>LWdown</td><td>Downward longwave radiation</td><td>W m<sup>-2</sup></td><td>Downward</td><td>yes</td><td>hourly and daily</td></tr><tr><td>Tair</td><td>Air temperature</td><td>K</td><td>–</td><td>yes</td><td>hourly and daily</td></tr><tr><td>Qair</td><td>Specific humidity</td><td>kg kg<sup>-1</sup></td><td>–</td><td>yes</td><td>hourly and daily</td></tr><tr><td>PSurf</td><td>Station air pressure</td><td>Pa</td><td>–</td><td>yes</td><td>hourly and daily</td></tr><tr><td>Wind_N</td><td>Northward wind component</td><td>m s<sup>-1</sup></td><td>Northward</td><td>yes</td><td>logarithmic law</td></tr><tr><td>Wind_E</td><td>Eastward wind component</td><td>m s<sup>-1</sup></td><td>Eastward</td><td>yes</td><td>logarithmic law</td></tr><tr><td>Rainf</td><td>Rainfall rate</td><td>kg m<sup>-2</sup> s<sup>-1</sup></td><td>Downward</td><td>yes</td><td>long-term precipitation</td></tr><tr><td>Snowf</td><td>Snowfall rate</td><td>kg m<sup>-2</sup> s<sup>-1</sup></td><td>Downward</td><td>yes</td><td>long term precipitation</td></tr><tr><td colspan="6"><strong>Analysis data</strong></td></tr><tr><td>SWup</td><td>Upward shortwave radiation</td><td>W m<sup>-2</sup></td><td>Upward</td><td>no</td><td>none</td></tr><tr><td>LWup</td><td>Upward longwave radiation</td><td>W m<sup>-2</sup></td><td>Upward</td><td>no</td><td>none</td></tr><tr><td>Qle</td><td>Latent heat flux</td><td>W m<sup>-2</sup></td><td>Upward</td><td>no</td><td>none</td></tr><tr><td>Qh</td><td>Sensible heat flux</td><td>W m<sup>-2</sup></td><td>Upward</td><td>no</td><td>none</td></tr><tr><td colspan="6"><strong>Additional data (optional)</strong></td></tr><tr><td>Qg</td><td>Ground heat flux into soil</td><td>W m<sup>-2</sup></td><td>Downward</td><td>no</td><td>none</td></tr><tr><td>Qtau</td><td>Momentum flux</td><td>N m<sup>-2</sup></td><td>Downward</td><td>no</td><td>none</td></tr><tr><td>Tair2m</td><td>Near surface air temperature (2 m)</td><td>K</td><td>–</td><td>no</td><td>none</td></tr><tr><td>SoilTemp</td><td>Soil temperature (depth in metadata)</td><td>K</td><td>–</td><td>no</td><td>none</td></tr></tbody></table>

<table><caption>**Table 3.** Site climate classification missing and additional variables. Climate classification from Köppen–Geiger global dataset (Beck et al., 2018). Table 2 gives variable definitions.</caption><thead><tr><th>Site name</th><th>Class</th><th>Climate description</th><th>Missing variables</th><th>Additional variables</th></tr></thead><tbody><tr><td>AU-Preston</td><td>Cfb</td><td>Temperate, no dry season, warm summer</td><td>Snowf</td><td>Qtau</td></tr><tr><td>AU-SurreyHills</td><td>Cfb</td><td>Temperate, no dry season, warm summer</td><td>Snowf</td><td>Qtau</td></tr><tr><td>CA-Sunset</td><td>Csb</td><td>Temperate, dry summer, warm summer</td><td>Snowf</td><td>Qtau, SoilTemp</td></tr><tr><td>FI-Kumpula</td><td>Dfb</td><td>Cold, no dry season, warm summer</td><td>Snowf</td><td></td></tr><tr><td>FI-Torni</td><td>Dfb</td><td>Cold, no dry season, warm summer</td><td>Snowf</td><td></td></tr><tr><td>FR-Capitole</td><td>Cfa</td><td>Temperate, no dry season, hot summer</td><td>Snowf</td><td>Qtau</td></tr><tr><td>GR-HECKOR</td><td>Csa</td><td>Temperate, dry summer, hot summer</td><td>Snowf</td><td>Qtau</td></tr><tr><td>JP-Yoyogi</td><td>Cfa</td><td>Temperate, no dry season, hot summer</td><td></td><td></td></tr><tr><td>KR-Jungnang</td><td>Dwa</td><td>Cold, dry winter, hot summer</td><td>Snowf</td><td></td></tr><tr><td>KR-Ochang</td><td>Dwa</td><td>Cold, dry winter, hot summer</td><td>Snowf</td><td></td></tr><tr><td>MX-Escandon</td><td>Cwb</td><td>Temperate, dry winter, warm summer</td><td>Snowf, LWdown*</td><td>Qtau</td></tr><tr><td>NL-Amsterdam</td><td>Cfb</td><td>Temperate, no dry season, warm summer</td><td>Snowf</td><td>Qtau</td></tr><tr><td>PL-Lipowa</td><td>Dfb</td><td>Cold, no dry season, warm summer</td><td>Snowf</td><td></td></tr><tr><td>PL-Narutowicza</td><td>Dfb</td><td>Cold, no dry season, warm summer</td><td>Snowf</td><td></td></tr><tr><td>SG-TelokKurau06</td><td>Af</td><td>Tropical, rainforest</td><td>Snowf</td><td></td></tr><tr><td>UK-KingsCollege</td><td>Cfb</td><td>Temperate, no dry season, warm summer</td><td>Snowf</td><td>Qtau</td></tr><tr><td>UK-Swindon</td><td>Cfb</td><td>Temperate, no dry season, warm summer</td><td>Snowf</td><td>Qtau</td></tr><tr><td>US-Baltimore</td><td>Cfa</td><td>Temperate, no dry season, hot summer</td><td>Snowf</td><td>Qtau, SoilTemp</td></tr><tr><td>US-Minneapolis</td><td>Dfa</td><td>Cold, no dry season, hot summer</td><td>Snowf</td><td>Qtau, SoilTemp, Qg</td></tr><tr><td>US-WestPhoenix</td><td>BWh</td><td>Arid, desert, hot</td><td>Snowf</td><td>Qtau, SoilTemp</td></tr></tbody></table>

<footnote>* Note: as MX-Escandon LWdown data are unavailable during the 2011–2012 focus period, but were available in 2006, this earlier period is used to determine bias correction for ERA5 LWdown data.</footnote>
from ERA5 snowfall is removed from subsequent observed rainfall until mass balance of observed total precipitation is achieved. This corrects melting snow being recorded as rainfall.

## 2.5 ERA5 reanalysis data

The ERA5 reanalysis product (Hersbach et al., 2020) assimilates global satellite, atmospheric and ground-based observations to constrain numerical weather prediction simulations, producing global output at 0.25° spatial and hourly temporal resolutions from 1959 to the present. It is therefore useful as a globally consistent and accessible source of meteorological data across space and time. ERA5, and its lower resolution predecessor ERA-Interim (Dee et al., 2011), have been used extensively to provide meteorological forcing data to drive land surface models and gap fill flux tower observations (Vuichard and Papale, 2015; Kokkonen et al., 2018; Pastorello et al., 2020; Ukkola et al., 2017, 2022).

The ERA5 hourly single level (Hersbach et al., 2018) dataset (retrieved from NCI Australia (Druken, 2020)) is used for gap filling missing observations within the focus periods (Table 1) and for the 10-year model spin up period. However, combining ERA5 data directly with urban flux tower observations has several deficiencies.

Grid-scale ERA5 data are not directly compatible with point-scale urban flux tower observations. This incompatibility is three-fold:

1. **Horizontally.** The ERA5 grid cell area (of the order of 500 km²) does not match the flux footprint from tower observations (of order 1 km²). The ERA5 surface characteristics, including elevation, are based on an average description for the grid which may differ from surface characteristics around the observing tower, particularly in coastal or mountainous regions (Martens et al., 2020) in which many cities are located.

2. **Vertically.** ERA5 provide near-surface variables (2 or 10 m above ground level), aligning with World Meteorological Organization (WMO) guidelines for standard regional observations taken over short grass (World Meteorological Organization, 2008). As the urban roughness elements (e.g. buildings) are much taller than grass, instruments are mounted on towers at heights greater than 2–5 times average building height in order to be located within the inertial sub layer or constant flux layer (Velasco and Roth, 2010; Barlow, 2014; Grimmond and Ward, 2021).

3. **Land surface.** As the current operational ERA5 modelling systems do not include an urban land surface scheme (Boussetta et al., 2013; McNorton et al., 2021), other land types (grass, crops, shrubs, trees etc.) are used to characterize the grid cell (Table 4). Urban land surfaces are well known to alter local meteorological

conditions (Oke et al., 2017), therefore ERA5 output will likely differ from locally observed conditions.

Outside of cities there are known diurnal and seasonal biases between the ERA5 near-surface variables and observations (Haiden et al., 2018; Betts et al., 2019; Nogueira, 2020; Martens et al., 2020; Jiang et al., 2021). These biases are an outcome of simplifying assumptions made in model parameterizations and inadequacies of modelling frameworks in general (Cucchi et al., 2020). Various approaches to reduce ERA5 biases in non-urban areas have been proposed. For example, the Water and Global Change (WATCH) Forcing Data (WFD) project use gridded observations to bias correct ERA-Interim data (Weedon et al., 2011), and more recently ERA5 data, creating the global WFDE5 dataset for impact studies (Cucchi et al., 2020). WFDE5 relies on the Climate Research Unit (CRU) monthly time series of gridded observations with resolution courser than ERA5 (New et al., 1999), requiring ERA5 to be regridded to a lower resolution. This may reduce the representativeness of the ERA5 data, particularly in heterogeneous or complex terrain.

Alternatively, local observations can be used to bias correct ERA data, e.g. the linear regression corrections using tower observations applied to FLUXNET datasets (Vuichard and Papale, 2015; Pastorello et al., 2020). However, linear methods neither conserve the variability of observations (Vuichard and Papale, 2015) (Sect. 3), nor can they correct diurnal timing differences within ERA5 data (e.g. out of phase from urban temporal profiles, which are typically delayed compared with non-urban surfaces used in ERA5, Fig. 4).

To account for the mischaracterization of sites (Table 4) and other listed deficiencies in ERA5, we develop a novel set of methods to bias correct ERA5 data to better represent observed urban conditions (Sect. 2.6).

## 2.6 Bias correction methods

Bias correction approaches used in the collection depend on the forcing variable (Table 2) and are described below.

## 2.7 Hourly and daily corrections

For incoming longwave radiation, air temperature, specific humidity and air pressure, the mean bias between ERA5 and local flux tower observations are calculated for each hour (h) and each day of a year (D) in a 60 d rolling window of a representative year (Fig. 4a). The calculated bias $\eta_{bias}(D, h)$ is subtracted from the complete ERA5 time series $\eta_{ERA5}(t)$ to create a new corrected time series:

$$ \eta(t) = \eta_{ERA5}(t) - \eta_{bias}(D, h). \qquad (1) $$

The ERA5 data are from the grid nearest the observation site with at least 50% land. The resulting corrected time series (e.g. Fig. 4b) is used for gap
from ERA5 snowfall is removed from subsequent observed rainfall until mass balance of observed total precipitation is achieved. This corrects melting snow being recorded as rainfall.

## 2.5 ERA5 reanalysis data

The ERA5 reanalysis product (Hersbach et al., 2020) assimilates global satellite, atmospheric and ground-based observations to constrain numerical weather prediction simulations, producing global output at 0.25° spatial and hourly temporal resolutions from 1959 to the present. It is therefore useful as a globally consistent and accessible source of meteorological data across space and time. ERA5, and its lower resolution predecessor ERA-Interim (Dee et al., 2011), have been used extensively to provide meteorological forcing data to drive land surface models and gap fill flux tower observations (Vuichard and Papale, 2015; Kokkonen et al., 2018; Pastorello et al., 2020; Ukkola et al., 2017, 2022).

The ERA5 hourly single level (Hersbach et al., 2018) dataset (retrieved from NCI Australia (Druken, 2020)) is used for gap filling missing observations within the focus periods (Table 1) and for the 10-year model spin up period. However, combining ERA5 data directly with urban flux tower observations has several deficiencies.

Grid-scale ERA5 data are not directly compatible with point-scale urban flux tower observations. This incompatibility is three-fold:

1. **Horizontally.** The ERA5 grid cell area (of the order of 500 km²) does not match the flux footprint from tower observations (of order 1 km²). The ERA5 surface characteristics, including elevation, are based on an average description for the grid which may differ from surface characteristics around the observing tower, particularly in coastal or mountainous regions (Martens et al., 2020) in which many cities are located.

2. **Vertically.** ERA5 provide near-surface variables (2 or 10 m above ground level), aligning with World Meteorological Organization (WMO) guidelines for standard regional observations taken over short grass (World Meteorological Organization, 2008). As the urban roughness elements (e.g. buildings) are much taller than grass, instruments are mounted on towers at heights greater than 2–5 times average building height in order to be located within the inertial sub layer or constant flux layer (Velasco and Roth, 2010; Barlow, 2014; Grimmond and Ward, 2021).

3. **Land surface.** As the current operational ERA5 modelling systems do not include an urban land surface scheme (Boussetta et al., 2013; McNorton et al., 2021), other land types (grass, crops, shrubs, trees etc.) are used to characterize the grid cell (Table 4). Urban land surfaces are well known to alter local meteorological

conditions (Oke et al., 2017), therefore ERA5 output will likely differ from locally observed conditions.

Outside of cities there are known diurnal and seasonal biases between the ERA5 near-surface variables and observations (Haiden et al., 2018; Betts et al., 2019; Nogueira, 2020; Martens et al., 2020; Jiang et al., 2021). These biases are an outcome of simplifying assumptions made in model parameterizations and inadequacies of modelling frameworks in general (Cucchi et al., 2020). Various approaches to reduce ERA5 biases in non-urban areas have been proposed. For example, the Water and Global Change (WATCH) Forcing Data (WFD) project use gridded observations to bias correct ERA-Interim data (Weedon et al., 2011), and more recently ERA5 data, creating the global WFDE5 dataset for impact studies (Cucchi et al., 2020). WFDE5 relies on the Climate Research Unit (CRU) monthly time series of gridded observations with resolution courser than ERA5 (New et al., 1999), requiring ERA5 to be regridded to a lower resolution. This may reduce the representativeness of the ERA5 data, particularly in heterogeneous or complex terrain.

Alternatively, local observations can be used to bias correct ERA data, e.g. the linear regression corrections using tower observations applied to FLUXNET datasets (Vuichard and Papale, 2015; Pastorello et al., 2020). However, linear methods neither conserve the variability of observations (Vuichard and Papale, 2015) (Sect. 3), nor can they correct diurnal timing differences within ERA5 data (e.g. out of phase from urban temporal profiles, which are typically delayed compared with non-urban surfaces used in ERA5, Fig. 4).

To account for the mischaracterization of sites (Table 4) and other listed deficiencies in ERA5, we develop a novel set of methods to bias correct ERA5 data to better represent observed urban conditions (Sect. 2.6).

## 2.6 Bias correction methods

Bias correction approaches used in the collection depend on the forcing variable (Table 2) and are described below.

## 2.7 Hourly and daily corrections

For incoming longwave radiation, air temperature, specific humidity and air pressure, the mean bias between ERA5 and local flux tower observations are calculated for each hour (h) and each day of a year (D) in a 60d rolling window of a representative year (Fig. 4a). The calculated bias $\eta_{bias}(D, h)$ is subtracted from the complete ERA5 time series $\eta_{ERA5}(t)$ to create a new corrected time series:

$$ \eta(t) = \eta_{ERA5}(t) - \eta_{bias}(D, h). \qquad (1) $$

The ERA5 data are from the grid nearest the observation site with at least 50% land. The resulting corrected time series (e.g. Fig. 4b) is used for gap
<table><caption>**Table 4.** Surface cover information as specified in ERA5 differs from actual tower site characteristics (see Table 6), and so ERA5 data are corrected (Sect. 2.6). Given the ERA5 surface roughness values vary slightly through time, the values listed are indicative (from 1 January 2000). Effective roughness is our correction accounting for observed urban mean wind speeds.</caption><thead><tr><th>Site</th><th>ERA5 low vegetation</th><th>ERA5 high vegetation</th><th>Low vegetation fraction</th><th>High vegetation fraction</th><th>Lake (or sea) fraction</th><th>Bare soil fraction</th><th>ERA5 surface roughness [m]</th><th>Effective roughness [m]</th></tr></thead><tbody><tr><td>AU-Preston</td><td>tall grass</td><td>interrupted forest</td><td>0.484</td><td>0.407</td><td>0.088</td><td>0.021</td><td>0.514</td><td>0.289</td></tr><tr><td>AU-SurreyHills</td><td>tall grass</td><td>interrupted forest</td><td>0.484</td><td>0.407</td><td>0.088</td><td>0.021</td><td>0.514</td><td>0.368</td></tr><tr><td>CA-Sunset</td><td>crops, mixed farming</td><td>evergreen needleleaf trees</td><td>0.205</td><td>0.723</td><td>0.071</td><td>0.000</td><td>1.077</td><td>1.508</td></tr><tr><td>FI-Kumpula</td><td>crops, mixed farming</td><td>evergreen needleleaf trees</td><td>0.296</td><td>0.352</td><td>0.137</td><td>0.215</td><td>0.708</td><td>0.703</td></tr><tr><td>FI-Torni</td><td>crops, mixed farming</td><td>evergreen needleleaf trees</td><td>0.296</td><td>0.352</td><td>0.137</td><td>0.215</td><td>0.708</td><td>0.424</td></tr><tr><td>FR-Capitole</td><td>crops, mixed farming</td><td>interrupted forest</td><td>0.920</td><td>0.050</td><td>0.004</td><td>0.025</td><td>0.291</td><td>0.519</td></tr><tr><td>GR-HECKOR</td><td>crops, mixed farming</td><td>interrupted forest</td><td>0.172</td><td>0.463</td><td>0.158</td><td>0.207</td><td>0.505</td><td>1.187</td></tr><tr><td>JP-Yoyogi</td><td>semidesert</td><td>no vegetation recorded</td><td>0.943</td><td>0.000</td><td>0.010</td><td>0.047</td><td>0.015</td><td>0.649</td></tr><tr><td>KR-Jungnang</td><td>crops, mixed farming</td><td>evergreen needleleaf trees</td><td>0.781</td><td>0.168</td><td>0.051</td><td>0.000</td><td>0.516</td><td>0.074</td></tr><tr><td>KR-Ochang</td><td>irrigated crops</td><td>interrupted forest</td><td>0.281</td><td>0.716</td><td>0.003</td><td>0.000</td><td>0.844</td><td>0.181</td></tr><tr><td>MX-Escandon</td><td>evergreen shrubs</td><td>mixed forest/woodland</td><td>0.743</td><td>0.216</td><td>0.006</td><td>0.035</td><td>0.404</td><td>0.229</td></tr><tr><td>NL-Amsterdam</td><td>crops, mixed farming</td><td>interrupted forest</td><td>0.867</td><td>0.061</td><td>0.056</td><td>0.015</td><td>0.248</td><td>0.254</td></tr><tr><td>PL-Lipowa</td><td>crops, mixed farming</td><td>interrupted forest</td><td>0.855</td><td>0.144</td><td>0.001</td><td>0.000</td><td>0.250</td><td>0.306</td></tr><tr><td>PL-Narutowicza</td><td>crops, mixed farming</td><td>interrupted forest</td><td>0.855</td><td>0.144</td><td>0.001</td><td>0.000</td><td>0.250</td><td>0.558</td></tr><tr><td>SG-TelokKurau06</td><td>irrigated crops</td><td>interrupted forest</td><td>0.905</td><td>0.021</td><td>0.074</td><td>0.000</td><td>0.335</td><td>0.309</td></tr><tr><td>UK-KingsCollege</td><td>crops, mixed farming</td><td>interrupted forest</td><td>0.609</td><td>0.372</td><td>0.020</td><td>0.000</td><td>0.504</td><td>0.315</td></tr><tr><td>UK-Swindon</td><td>crops, mixed farming</td><td>interrupted forest</td><td>0.727</td><td>0.251</td><td>0.001</td><td>0.021</td><td>0.397</td><td>0.146</td></tr><tr><td>US-Baltimore</td><td>crops, mixed farming</td><td>deciduous broadleaf trees</td><td>0.044</td><td>0.908</td><td>0.048</td><td>0.000</td><td>1.675</td><td>1.076</td></tr><tr><td>US-Minneapolis1</td><td>crops, mixed farming</td><td>interrupted forest</td><td>0.228</td><td>0.706</td><td>0.059</td><td>0.006</td><td>0.814</td><td>0.242</td></tr><tr><td>US-Minneapolis2</td><td>crops, mixed farming</td><td>interrupted forest</td><td>0.228</td><td>0.706</td><td>0.059</td><td>0.006</td><td>0.814</td><td>0.406</td></tr><tr><td>US-WestPhoenix</td><td>semidesert</td><td>evergreen needleleaf trees</td><td>0.949</td><td>0.051</td><td>0.000</td><td>0.000</td><td>0.084</td><td>0.404</td></tr></tbody></table>

filling site observations and for the spin up period. The subroutine *rolling_hourly_bias_correction* in the file *pipeline_functions.py* (Sect. 5) undertakes corrections with the following steps:

1. If observations have a 30 min resolution, average to 60 min to match ERA5 periods;

2. Remove ERA5 data for time periods where site observations are missing;

3. Calculate mean for each hour with both data sets to create a “representative” year of 366 d;

4. Extend to a 3-year period by duplication to provide smoother transitions at year end;

5. Calculate hourly means in a 60 d rolling window across the repeating time series, excluding data in windows with less than 30 observations. Repeat mean calculation for greater smoothing;

6. Calculate a time series of the bias between observed and ERA5 rolling means;

7. Fill gaps in the bias time series by linear interpolation through each hour separately;

8. Remove first and last year in the bias time series, using only the central year to bias correct each hour in the original ERA5 time series.

A 60 d rolling window is selected to smooth out individual weather events while still capturing seasonal variation. Repeating the representative year three times prior to smoothing ensures bias corrections match at the start and end of the

year. The resulting set of bias correction curves (Fig. 4a) has greater robustness when multiple years are available.

### 2.7.1 Logarithmic wind profile correction

Wind speed differences between ERA5 and site observations can result from errors in modelled synoptic-scale speeds, differences in representative heights, and differences in surface aerodynamic properties like roughness and displacement height. To correct bias and maintain standard deviations of the wind components (U, V) of observations at sensor height (z<sub>site</sub>), the following correction to ERA5 data is undertaken assuming both a logarithmic wind profile and neutral conditions (Goret et al., 2019):

$$u_{\text{corr}} = u_{\text{ERA5}} \frac{\ln\left(\frac{z_{\text{site}} - d_{\text{site}}}{z_{0,\text{site}}}\right)}{\ln\left(\frac{z_{\text{grid}} - d_{\text{grid}}}{z_{0,\text{grid}}}\right)}, \quad (2)$$

where $u_{corr}$ is the corrected wind speed at z<sub>site</sub>. The ERA5 wind ($u_{ERA5}$) at 10 m (z<sub>grid</sub>) is used with the site surface roughness (z<sub>0,site</sub>) and displacement height (d<sub>site</sub>), and grid roughness (z<sub>0,grid</sub>) (Table 4) while assuming grid displacement height (d<sub>grid</sub>) is zero for simplicity. If the resulting mean value of $u_{corr}$ differs from observed mean value by more than 0.01 m s<sup>-1</sup>, then z<sub>0,grid</sub> is iteratively adapted until this threshold accuracy in mean wind speed is achieved. Derived z<sub>0,grid</sub> values are given in Table 4 (last column). Note this approach ignores seasonal effects from vegetation phenology and directional effects but ensures mean wind speeds are appropriate at the urban z<sub>site</sub> while conserving variability within the ERA5 derived wind data.
<table><caption>**Table 4.** Surface cover information as specified in ERA5 differs from actual tower site characteristics (see Table 6), and so ERA5 data are corrected (Sect. 2.6). Given the ERA5 surface roughness values vary slightly through time, the values listed are indicative (from 1 January 2000). Effective roughness is our correction accounting for observed urban mean wind speeds.</caption><thead><tr><th>Site</th><th>ERA5 low vegetation</th><th>ERA5 high vegetation</th><th>Low vegetation fraction</th><th>High vegetation fraction</th><th>Lake (or sea) fraction</th><th>Bare soil fraction</th><th>ERA5 surface roughness [m]</th><th>Effective roughness [m]</th></tr></thead><tbody><tr><td>AU-Preston</td><td>tall grass</td><td>interrupted forest</td><td>0.484</td><td>0.407</td><td>0.088</td><td>0.021</td><td>0.514</td><td>0.289</td></tr><tr><td>AU-SurreyHills</td><td>tall grass</td><td>interrupted forest</td><td>0.484</td><td>0.407</td><td>0.088</td><td>0.021</td><td>0.514</td><td>0.368</td></tr><tr><td>CA-Sunset</td><td>crops, mixed farming</td><td>evergreen needleleaf trees</td><td>0.205</td><td>0.723</td><td>0.071</td><td>0.000</td><td>1.077</td><td>1.508</td></tr><tr><td>FI-Kumpula</td><td>crops, mixed farming</td><td>evergreen needleleaf trees</td><td>0.296</td><td>0.352</td><td>0.137</td><td>0.215</td><td>0.708</td><td>0.703</td></tr><tr><td>FI-Torni</td><td>crops, mixed farming</td><td>evergreen needleleaf trees</td><td>0.296</td><td>0.352</td><td>0.137</td><td>0.215</td><td>0.708</td><td>0.424</td></tr><tr><td>FR-Capitole</td><td>crops, mixed farming</td><td>interrupted forest</td><td>0.920</td><td>0.050</td><td>0.004</td><td>0.025</td><td>0.291</td><td>0.519</td></tr><tr><td>GR-HECKOR</td><td>crops, mixed farming</td><td>interrupted forest</td><td>0.172</td><td>0.463</td><td>0.158</td><td>0.207</td><td>0.505</td><td>1.187</td></tr><tr><td>JP-Yoyogi</td><td>semidesert</td><td>no vegetation recorded</td><td>0.943</td><td>0.000</td><td>0.010</td><td>0.047</td><td>0.015</td><td>0.649</td></tr><tr><td>KR-Jungnang</td><td>crops, mixed farming</td><td>evergreen needleleaf trees</td><td>0.781</td><td>0.168</td><td>0.051</td><td>0.000</td><td>0.516</td><td>0.074</td></tr><tr><td>KR-Ochang</td><td>irrigated crops</td><td>interrupted forest</td><td>0.281</td><td>0.716</td><td>0.003</td><td>0.000</td><td>0.844</td><td>0.181</td></tr><tr><td>MX-Escandon</td><td>evergreen shrubs</td><td>mixed forest/woodland</td><td>0.743</td><td>0.216</td><td>0.006</td><td>0.035</td><td>0.404</td><td>0.229</td></tr><tr><td>NL-Amsterdam</td><td>crops, mixed farming</td><td>interrupted forest</td><td>0.867</td><td>0.061</td><td>0.056</td><td>0.015</td><td>0.248</td><td>0.254</td></tr><tr><td>PL-Lipowa</td><td>crops, mixed farming</td><td>interrupted forest</td><td>0.855</td><td>0.144</td><td>0.001</td><td>0.000</td><td>0.250</td><td>0.306</td></tr><tr><td>PL-Narutowicza</td><td>crops, mixed farming</td><td>interrupted forest</td><td>0.855</td><td>0.144</td><td>0.001</td><td>0.000</td><td>0.250</td><td>0.558</td></tr><tr><td>SG-TelokKurau06</td><td>irrigated crops</td><td>interrupted forest</td><td>0.905</td><td>0.021</td><td>0.074</td><td>0.000</td><td>0.335</td><td>0.309</td></tr><tr><td>UK-KingsCollege</td><td>crops, mixed farming</td><td>interrupted forest</td><td>0.609</td><td>0.372</td><td>0.020</td><td>0.000</td><td>0.504</td><td>0.315</td></tr><tr><td>UK-Swindon</td><td>crops, mixed farming</td><td>interrupted forest</td><td>0.727</td><td>0.251</td><td>0.001</td><td>0.021</td><td>0.397</td><td>0.146</td></tr><tr><td>US-Baltimore</td><td>crops, mixed farming</td><td>deciduous broadleaf trees</td><td>0.044</td><td>0.908</td><td>0.048</td><td>0.000</td><td>1.675</td><td>1.076</td></tr><tr><td>US-Minneapolis1</td><td>crops, mixed farming</td><td>interrupted forest</td><td>0.228</td><td>0.706</td><td>0.059</td><td>0.006</td><td>0.814</td><td>0.242</td></tr><tr><td>US-Minneapolis2</td><td>crops, mixed farming</td><td>interrupted forest</td><td>0.228</td><td>0.706</td><td>0.059</td><td>0.006</td><td>0.814</td><td>0.406</td></tr><tr><td>US-WestPhoenix</td><td>semidesert</td><td>evergreen needleleaf trees</td><td>0.949</td><td>0.051</td><td>0.000</td><td>0.000</td><td>0.084</td><td>0.404</td></tr></tbody></table>

filling site observations and for the spin up period. The subroutine *rolling_hourly_bias_correction* in the file *pipeline_functions.py* (Sect. 5) undertakes corrections with the following steps:

1. If observations have a 30 min resolution, average to 60 min to match ERA5 periods;

2. Remove ERA5 data for time periods where site observations are missing;

3. Calculate mean for each hour with both data sets to create a “representative” year of 366 d;

4. Extend to a 3-year period by duplication to provide smoother transitions at year end;

5. Calculate hourly means in a 60 d rolling window across the repeating time series, excluding data in windows with less than 30 observations. Repeat mean calculation for greater smoothing;

6. Calculate a time series of the bias between observed and ERA5 rolling means;

7. Fill gaps in the bias time series by linear interpolation through each hour separately;

8. Remove first and last year in the bias time series, using only the central year to bias correct each hour in the original ERA5 time series.

A 60 d rolling window is selected to smooth out individual weather events while still capturing seasonal variation. Repeating the representative year three times prior to smoothing ensures bias corrections match at the start and end of the

year. The resulting set of bias correction curves (Fig. 4a) has greater robustness when multiple years are available.

## 2.7.1 Logarithmic wind profile correction

Wind speed differences between ERA5 and site observations can result from errors in modelled synoptic-scale speeds, differences in representative heights, and differences in surface aerodynamic properties like roughness and displacement height. To correct bias and maintain standard deviations of the wind components (U, V) of observations at sensor height (z<sub>site</sub>), the following correction to ERA5 data is undertaken assuming both a logarithmic wind profile and neutral conditions (Goret et al., 2019):

$$u_{\text{corr}} = u_{\text{ERA5}} \frac{\ln\left(\frac{z_{\text{site}} - d_{\text{site}}}{z_{0, \text{site}}}\right)}{\ln\left(\frac{z_{\text{grid}} - d_{\text{grid}}}{z_{0, \text{grid}}}\right)}, \quad (2)$$

where $u_{corr}$ is the corrected wind speed at z<sub>site</sub>. The ERA5 wind (u<sub>ERA5</sub>) at 10 m (z<sub>grid</sub>) is used with the site surface roughness (z<sub>0,site</sub>) and displacement height (d<sub>site</sub>), and grid roughness (z<sub>0,grid</sub>) (Table 4) while assuming grid displacement height (d<sub>grid</sub>) is zero for simplicity. If the resulting mean value of u<sub>corr</sub> differs from observed mean value by more than 0.01 m s<sup>-1</sup>, then z<sub>0,grid</sub> is iteratively adapted until this threshold accuracy in mean wind speed is achieved. Derived z<sub>0,grid</sub> values are given in Table 4 (last column). Note this approach ignores seasonal effects from vegetation phenology and directional effects but ensures mean wind speeds are appropriate at the urban z<sub>site</sub> while conserving variability within the ERA5 derived wind data.
<figure><img src="image_5.png" /><figcaption>Figure 4. Urban-PLUMBER reanalysis bias correction methods. Demonstrated using air temperature (Tair) for the grid containing the King's College London site (UK-KingsCollege). (a) Hourly (colour) bias calculated for each day of a “representative” year and applied to entire ERA5 time series; (b) diurnal hourly mean Urban-PLUMBER correction (UP, red), observations (black), original ERA5 data (blue), WFDE5 bias corrected data (green) and linear bias correction method used in FLUXNET (LN, yellow). Our new UP method has smaller mean absolute errors (MAE) overall, and can correct both pattern and phase errors of ERA5 (Sect. 3).</figcaption></figure>

## 2.7.2 Long-term precipitation correction

Total precipitation is an important variable in urban land surface models because of the effect on soil moisture which evolves over multi-year periods (Best and Grimmond, 2014). Most of the observational datasets included are not long enough to capture interannual variations. We therefore use longer term total precipitation (P) from nearby stations from the Global Historical Climatology Network – Daily (GHCND) (Menne et al., 2012) over a 10 year period to correct ERA5 rain and snow fluxes ($\phi$) at each time step:

$$ \phi_{\text{corr}} = \frac{\sum_{i=1}^{10 \text{ yrs}} P_{\text{GHCND}}}{\sum_{i=1}^{10 \text{ yrs}} P_{\text{ERA5}}} \phi_{\text{ERA5}(t)}. \quad (3) $$

Gaps in the nearest GHCND station data are progressively filled by the next nearest station until no gaps are present. If gaps could not be filled with GHNCD stations within 2° of latitude and longitude from the flux tower, and no alternative records are found (e.g. from national meteorological and hydrological services), then ERA5 rates are used unadjusted (viz., KR-Ochang and KR-Jungnang). This assumes precipitation occurs on the same dates and at the same times in ERA5 and observed datasets, which may become less valid under increasingly convective conditions.

## 2.7.3 Linear bias correction

The FLUXNET2015 collection of 212 flux tower sites (Pastorello et al., 2020) bias correction method uses linear regression between site observations and reanalysis data to derive one slope (s) and intercept (b) per site, hence “unbiasing” all ERA time steps (i):

$$ \text{LN}_i = s \cdot \text{ERA5}_i + b. \quad (4) $$

Following FLUXNET2015, global radiation and wind fields are assigned an intercept of zero, and precipitation is not linearly modified (Vuichard and Papale, 2015). FLUXNET2015 uses the coarser resolution ERA-Interim (spatial: 0.5° cf. 0.25°; temporal 3 h cf. 1 h) than ERA5. In this evaluation we use ERA5, which is found to be a consistent improvement over ERA-Interim (Albergel et al., 2018). However, after assessment we chose not to use a linear method (LN) for correcting variables in this collection. Its description is retained here for comparison purposes (Sect. 3).

# 3 Gap filling evaluation

Site observations are quality controlled by individual data providers and collectively for this project (Sect. 2.3). The observed data required for forcing land surface models are then gap filled using a novel method of bias correcting reanalysis data. In this section, four methods which draw on ERA5 data are evaluated:

1. ERA5. Nearest land-based 0.25° resolution ERA5 (Hersbach et al., 2018) grid without bias correction.

2. W5. Nearest WFDE5 (Cucchi et al., 2020) grid (which uses bias correction from 0.5° CRU monthly gridded observations).

3. UP. The Urban-PLUMBER methods described here (using site observations for bias correction).

4. LN. Linear methods based on FLUXNET2015 (Vuichard and Papale, 2015; Pastorello et al., 2020) (using site observations for bias correction).

To evaluate the methods available for gap filling, quality-controlled tower site observations ($O_i$) are used to assess the calculated value ($\eta$) at timestep i using three metrics:
<figure><img src="image_6.png" /><figcaption>Figure 4. Urban-PLUMBER reanalysis bias correction methods. Demonstrated using air temperature (Tair) for the grid containing the King's College London site (UK-KingsCollege). (a) Hourly (colour) bias calculated for each day of a “representative” year and applied to entire ERA5 time series; (b) diurnal hourly mean Urban-PLUMBER correction (UP, red), observations (black), original ERA5 data (blue), WFDE5 bias corrected data (green) and linear bias correction method used in FLUXNET (LN, yellow). Our new UP method has smaller mean absolute errors (MAE) overall, and can correct both pattern and phase errors of ERA5 (Sect. 3).</figcaption></figure>

## 2.7.2 Long-term precipitation correction

Total precipitation is an important variable in urban land surface models because of the effect on soil moisture which evolves over multi-year periods (Best and Grimmond, 2014). Most of the observational datasets included are not long enough to capture interannual variations. We therefore use longer term total precipitation (P) from nearby stations from the Global Historical Climatology Network – Daily (GHCND) (Menne et al., 2012) over a 10 year period to correct ERA5 rain and snow fluxes ($\phi$) at each time step:

$$ \phi_{corr} = \frac{\sum_{i=1}^{10 \text{ yrs}} P_{\text{GHCND}}}{\sum_{i=1}^{10 \text{ yrs}} P_{\text{ERA5}}} \phi_{\text{ERA5}(t)}. \quad (3) $$

Gaps in the nearest GHCND station data are progressively filled by the next nearest station until no gaps are present. If gaps could not be filled with GHNCD stations within 2° of latitude and longitude from the flux tower, and no alternative records are found (e.g. from national meteorological and hydrological services), then ERA5 rates are used unadjusted (viz., KR-Ochang and KR-Jungnang). This assumes precipitation occurs on the same dates and at the same times in ERA5 and observed datasets, which may become less valid under increasingly convective conditions.

## 2.7.3 Linear bias correction

The FLUXNET2015 collection of 212 flux tower sites (Pastorello et al., 2020) bias correction method uses linear regression between site observations and reanalysis data to derive one slope (s) and intercept (b) per site, hence “unbiasing” all ERA time steps (i):

$$ LN_i = s \cdot \text{ERA5}_i + b. \quad (4) $$

Following FLUXNET2015, global radiation and wind fields are assigned an intercept of zero, and precipitation is not linearly modified (Vuichard and Papale, 2015). FLUXNET2015 uses the coarser resolution ERA-Interim (spatial: 0.5° cf. 0.25°; temporal 3 h cf. 1 h) than ERA5. In this evaluation we use ERA5, which is found to be a consistent improvement over ERA-Interim (Albergel et al., 2018). However, after assessment we chose not to use a linear method (LN) for correcting variables in this collection. Its description is retained here for comparison purposes (Sect. 3).

# 3 Gap filling evaluation

Site observations are quality controlled by individual data providers and collectively for this project (Sect. 2.3). The observed data required for forcing land surface models are then gap filled using a novel method of bias correcting reanalysis data. In this section, four methods which draw on ERA5 data are evaluated:

1. ERA5. Nearest land-based 0.25° resolution ERA5 (Hersbach et al., 2018) grid without bias correction.

2. W5. Nearest WFDE5 (Cucchi et al., 2020) grid (which uses bias correction from 0.5° CRU monthly gridded observations).

3. UP. The Urban-PLUMBER methods described here (using site observations for bias correction).

4. LN. Linear methods based on FLUXNET2015 (Vuichard and Papale, 2015; Pastorello et al., 2020) (using site observations for bias correction).

To evaluate the methods available for gap filling, quality-controlled tower site observations ($O_i$) are used to assess the calculated value ($\eta$) at timestep $i$ using three metrics:
a. Mean bias error (MBE): $\frac{\sum_{i=1}^{n} \eta_i - O_i}{n}$

b. Mean absolute error (MAE): $\frac{\sum_{i=1}^{n} |\eta_i - O_i|}{n}$

c. Normalized standard deviation (nSD):

$$ \sqrt{\frac{\sum_{i=1}^{n} (\eta_i - \bar{\eta})^2}{n-1}} / \sqrt{\frac{\sum_{i=1}^{n} (O_i - \bar{O})^2}{n-1}} $$

Where $\bar{O}$ and $\bar{\eta}$ are time-averaged over n data points. The timestamps for each variable are made consistent between the data sources: ERA5 are hourly time ending or instantaneous data (Hersbach et al., 2018); 60 min site observations are hour ending; 30 min site observations are converted to 60 min time ending by averaging; whereas as the WFDE5 SWdown, LWdown and Rainf are natively 60 min time beginning (Cucchi et al., 2020), they are shifted forward to match the time ending timestamps; and WFDE5 Tair, Qair, Psurf and Wind are instantaneous samples on the hour so their timestamp remains unchanged.

To summarize each metric for the 20 sites, we use box-plots (Fig. 5) for the seven forcing variables. The evaluation inherently considers the net differences associated with both the spatial (vertical and surface cover) differences and errors (model and observation) from the two datasets. Uncorrected ERA5 (blue, Fig. 5) biases are generally negative for Tair, LWdown and Wind, and generally positive for Qair and SWdown. These biases can be partly explained by the ERA5 framework not including an urban surface model. For example, the well-documented warmer air temperature in cities (urban heat island) are not modelled in ERA5 because natural land surfaces are assumed in simulations (Table 4), although ERA5 can include an urban signal if the data assimilated are from within urban areas (Tang et al., 2021). Qair shows a general positive bias as evapotranspiration will be overestimated in ERA5 without an urban land surface representation. Likewise, ERA5 SWdown are overestimated and LWdown underestimated possibly because urban air pollution effects are not included (Oke, 1988).

Other discrepancies between ERA5 data and site data arise from elevation differences in height above sea level (a.s.l.) of the ERA5 grid and site. For example, the MX-Escandon tower in Mexico City measurement height is 2277 m a.s.l., whereas the ERA5 grid cell is assigned a surface elevation of 2540 m a.s.l. because the cell includes nearby mountains. This 263 m difference causes a negative bias to Psurf of 2594 Pa and contributes to a Tair difference of -2.15 K. Additionally, orographic uplift increases the grid cell rainfall, leading to positive ERA5 bias of $2.25 \times 10^{-5} \text{ kg m}^{-2} \text{ s}^{-1}$ (+710 mm yr$^{-1}$) compared with the MX-Escandon site observations. The ERA5 rainfall bias is even more pronounced at CA-Sunset in Vancouver, Canada, with a +1178 mm yr$^{-1}$ bias. These results are consistent with other studies highlighting discrepancies between reanalysis and local data in mountainous regions (Kokkonen et al., 2018).

The other three methods (WFDE5 (W5), linear debiasing (LN) and the Urban-PLUMBER corrections (UP)) apply bias

corrections to ERA5 data, and so may be expected to reduce ERA5 errors. However, W5 does not reduce errors at these sites, most likely because the observations used for W5 bias correction are at very different spatial scales to the flux tower footprints (2500 km$^2$ versus 1 km$^2$) and include observations from non-urban locations. The LN and UP methods eliminate spatial mismatches by drawing on local site or nearby rain gauge observations. As such, they do reduce the mean bias error to near zero for most variables. Notably, the UP methods outperform LN methods in normalized standard deviation. As Vuichard and Papale (2015) noted, linear methods do not conserve the variability of observations, nor can they correct for diurnal phase shifts of some variables observed in urban areas (Fig. 4). Therefore, we consider the UP methods to be the most appropriate at these urban sites. However, we apply no bias corrections to SWdown because the hourly and daily corrections (i.e. UP methods applied to other variables) adversely impact the standard deviation errors, as does the LN method for SWdown.

## 4 Data availability

Data described in this manuscript can be accessed from https://doi.org/10.5281/zenodo.7104984 (Lipson et al., 2022) under a Creative Commons Attribution licence (CC-BY-4.0).

We recommend data users consult with site contributing authors and/or the coordination team early (i.e. planning stage) in projects that plan to use these data. Relevant contacts are included in site metadata.

## 5 Code availability

Code used to process datasets are available at https://doi.org/10.5281/zenodo.7108466 (Lipson, 2022a).

Code used to create paper figures are available at https://doi.org/10.5281/zenodo.6590941 (Lipson, 2022b).

## 6 Data records

### 6.1 Data format

Time series data and site descriptive metadata are recorded in both plain text and netCDF4 (Rew et al., 1989) formats. Each site folder contains the following time series:

- [sitename]_raw_observations_[version]: site observed before project-wide quality control and gap filling (Table 1 gives period)

- [sitename]_clean_observations_[version]: after project-wide quality control and gap filling (Table 1 gives period)

- [sitename]_metforcing_[version]: continuous observations with reanalysis-derived data after quality control
a. Mean bias error (MBE): $\frac{\sum_{i=1}^{n} \eta_i - O_i}{n}$

b. Mean absolute error (MAE): $\frac{\sum_{i=1}^{n} |\eta_i - O_i|}{n}$

c. Normalized standard deviation (nSD):

$$ \sqrt{\frac{\sum_{i=1}^{n} (\eta_i - \bar{\eta})^2}{n-1}} / \sqrt{\frac{\sum_{i=1}^{n} (O_i - \bar{O})^2}{n-1}} $$

Where $\bar{O}$ and $\bar{\eta}$ are time-averaged over n data points. The timestamps for each variable are made consistent between the data sources: ERA5 are hourly time ending or instantaneous data (Hersbach et al., 2018); 60 min site observations are hour ending; 30 min site observations are converted to 60 min time ending by averaging; whereas as the WFDE5 SWdown, LWdown and Rainf are natively 60 min time beginning (Cucchi et al., 2020), they are shifted forward to match the time ending timestamps; and WFDE5 Tair, Qair, Psurf and Wind are instantaneous samples on the hour so their timestamp remains unchanged.

To summarize each metric for the 20 sites, we use box-plots (Fig. 5) for the seven forcing variables. The evaluation inherently considers the net differences associated with both the spatial (vertical and surface cover) differences and errors (model and observation) from the two datasets. Uncorrected ERA5 (blue, Fig. 5) biases are generally negative for Tair, LWdown and Wind, and generally positive for Qair and SWdown. These biases can be partly explained by the ERA5 framework not including an urban surface model. For example, the well-documented warmer air temperature in cities (urban heat island) are not modelled in ERA5 because natural land surfaces are assumed in simulations (Table 4), although ERA5 can include an urban signal if the data assimilated are from within urban areas (Tang et al., 2021). Qair shows a general positive bias as evapotranspiration will be overestimated in ERA5 without an urban land surface representation. Likewise, ERA5 SWdown are overestimated and LWdown underestimated possibly because urban air pollution effects are not included (Oke, 1988).

Other discrepancies between ERA5 data and site data arise from elevation differences in height above sea level (a.s.l.) of the ERA5 grid and site. For example, the MX-Escandon tower in Mexico City measurement height is 2277 m a.s.l., whereas the ERA5 grid cell is assigned a surface elevation of 2540 m a.s.l. because the cell includes nearby mountains. This 263 m difference causes a negative bias to Psurf of 2594 Pa and contributes to a Tair difference of -2.15 K. Additionally, orographic uplift increases the grid cell rainfall, leading to positive ERA5 bias of $2.25 \times 10^{-5} \text{ kg m}^{-2} \text{ s}^{-1}$ (+710 mm yr$^{-1}$) compared with the MX-Escandon site observations. The ERA5 rainfall bias is even more pronounced at CA-Sunset in Vancouver, Canada, with a +1178 mm yr$^{-1}$ bias. These results are consistent with other studies highlighting discrepancies between reanalysis and local data in mountainous regions (Kokkonen et al., 2018).

The other three methods (WFDE5 (W5), linear debiasing (LN) and the Urban-PLUMBER corrections (UP)) apply bias

corrections to ERA5 data, and so may be expected to reduce ERA5 errors. However, W5 does not reduce errors at these sites, most likely because the observations used for W5 bias correction are at very different spatial scales to the flux tower footprints (2500 km$^2$ versus 1 km$^2$) and include observations from non-urban locations. The LN and UP methods eliminate spatial mismatches by drawing on local site or nearby rain gauge observations. As such, they do reduce the mean bias error to near zero for most variables. Notably, the UP methods outperform LN methods in normalized standard deviation. As Vuichard and Papale (2015) noted, linear methods do not conserve the variability of observations, nor can they correct for diurnal phase shifts of some variables observed in urban areas (Fig. 4). Therefore, we consider the UP methods to be the most appropriate at these urban sites. However, we apply no bias corrections to SWdown because the hourly and daily corrections (i.e. UP methods applied to other variables) adversely impact the standard deviation errors, as does the LN method for SWdown.

## 4 Data availability

Data described in this manuscript can be accessed from https://doi.org/10.5281/zenodo.7104984 (Lipson et al., 2022) under a Creative Commons Attribution licence (CC-BY-4.0).

We recommend data users consult with site contributing authors and/or the coordination team early (i.e. planning stage) in projects that plan to use these data. Relevant contacts are included in site metadata.

## 5 Code availability

Code used to process datasets are available at https://doi.org/10.5281/zenodo.7108466 (Lipson, 2022a).

Code used to create paper figures are available at https://doi.org/10.5281/zenodo.6590941 (Lipson, 2022b).

## 6 Data records

### 6.1 Data format

Time series data and site descriptive metadata are recorded in both plain text and netCDF4 (Rew et al., 1989) formats. Each site folder contains the following time series:

- [sitename]_raw_observations_[version]: site observed before project-wide quality control and gap filling (Table 1 gives period)

- [sitename]_clean_observations_[version]: after project-wide quality control and gap filling (Table 1 gives period)

- [sitename]_metforcing_[version]: continuous observations with reanalysis-derived data after quality control
<figure><img src="image_6.png" /><figcaption>**Figure 5.** Evaluation of bias correction methods. Four methods (colour) to create gap filled observed time series data: ERA5 (blue), WFDE5 (W5, green), linear debiasing (LN, orange), UP (red, this study) using (row 1) mean bias error, (row 2) mean absolute error, (row 3) normalized standard deviation, with the 20 individual sites (dots), and ideal agreement with observations (red line) and boxplot showing distribution. The UP corrections (selected for use in this study) have lower overall errors (cf. other methods) except SWdown, where no corrections to ERA5 are applied.</figcaption></figure>

and gap filling (Table 2; forcing data for 10-year spin up, then Table 1 periods)

- `[sitename]_era5_corrected_[version]`: continuous time series (1990–2020) of bias corrected ERA5 reanalysis meteorological data (as used for gap filling and prepending metforcing observations).

Each site folder also contains the following site metadata:

- `[sitename]_sitedata_[version].csv`: comma separated text file for site characteristics metadata e.g. latitude, longitude, surface cover fraction and morphology (Tables 5, 6). This site characteristic data are also included within the metforcing netcdf (for convenience)

- index.html: a summary page of site information in html format, including site characteristics, site images, time series, gap filling, quality control and diurnal plots.

## 6.2 Time series metadata

The time series files include the following metadata:

- *title*: short description of the file

- *summary*: longer description of the file

- *sitename*: site code (e.g. AU-Preston)

- *long_sitename*: site long name, including city and country information

- *version*: version of current file

- *timesequencing*: start of time series in UTC (includes spin up) with period ending timestamps

- *timecoverage*: end of time series in UTC with period ending timestamps
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Evaluation of bias correction methods. Four methods (colour) to create gap filled observed time series data: ERA5 (blue), WFDE5 (W5, green), linear debiasing (LN, orange), UP (red, this study) using (row 1) mean bias error, (row 2) mean absolute error, (row 3) normalized standard deviation, with the 20 individual sites (dots), and ideal agreement with observations (red line) and boxplot showing distribution. The UP corrections (selected for use in this study) have lower overall errors (cf. other methods) except SWdown, where no corrections to ERA5 are applied.</figcaption></figure>

and gap filling (Table 2; forcing data for 10-year spin up, then Table 1 periods)

- `[sitename]_era5_corrected_[version]`: continuous time series (1990–2020) of bias corrected ERA5 reanalysis meteorological data (as used for gap filling and prepending metforcing observations).

Each site folder also contains the following site metadata:

- `[sitename]_sitedata_[version].csv`: comma separated text file for site characteristics metadata e.g. latitude, longitude, surface cover fraction and morphology (Tables 5, 6). This site characteristic data are also included within the metforcing netcdf (for convenience)

- index.html: a summary page of site information in html format, including site characteristics, site images, time series, gap filling, quality control and diurnal plots.

## 6.2 Time series metadata

The time series files include the following metadata:

- *title*: short description of the file

- *summary*: longer description of the file

- *sitename*: site code (e.g. AU-Preston)

- *long_sitename*: site long name, including city and country information

- *version*: version of current file

- *timesequencing*: start of time series in UTC (includes spin up) with period ending timestamps

- *timecoverage*: end of time series in UTC with period ending timestamps
- `time_analysis_start`: start of observed (focus) period in UTC with period ending timestamps

- `time_shown_in`: time standard (always UTC)

- `local_utf_offsetHours`: offset in hours of local time from UTC

- `timestep_interval_seconds`: period of block averaging in seconds (time step)

- `timestep_number_spinup`: number of time steps prior to observed focus period

- `timestep_number_analysis`: number of time steps in observed focus period

- `projectContact`: contact details for the Urban-PLUMBER project coordinators

- `observationsContact`: contact details of the observational site data providers

- `observations_reference`: published references associated with the observations

- `date_created`: date and time of creation of this file

- `source`: repository for processing code

- `comment`: additional comments associated with this dataset (e.g. excluded wind sectors).

Text file time series include metadata headers indicated with a hash (#) at line beginnings. Columns are headed by variable names in ALMA format (Table 2). NetCDF4 files include identical data, with additional attributes for each variable:

- `long_name`: plain language description of variable

- `standard_name`: equivalent variable name under the CF (climate and forecast) conventions

- `units`: SI (international system) units

- `ancillary_variables`: name of the associated quality control flag variable.

NetCDF files also include site characteristics parameter values and descriptions (Table 5). Times in all datasets are UTC. Python programmes are provided (Lipson, 2022a) to convert UTC times to local standard time (`convert_UTF_to_local_time.py`), and netCDF to text (`convert_nc_to_text.py`).

## 6.3 Site characteristics metadata

Site characteristics (Tables 5, 6) are essential for any use of these data, and fundamental to application of land surface models. These metadata are provided in two machine readable forms (plain text in csv files and netCDF4). The metadata are primarily drawn from published sources or as advised by the data providers. If local parameters are not known, values are estimated from high-resolution global datasets or derived from empirical relations. The sources for each parameter are included within the site characteristic metadata.

There are numerous methods to estimate the probable extent and weighting of turbulent fluxes footprints relative to the eddy covariance sensors located on a flux tower (Velasco and Roth, 2010). The eddy covariance flux footprint provides a basis to identify which area (and weighting) should be used to estimate the land surface fractions impacting the measurements. Some studies in this collection determine the footprint and resulting land cover fractions dynamically (e.g. for each 30 min period based on that period's observed atmospheric variables such as stability and wind direction), whereas others used a constant radius (e.g. based on the footprint climatology or rule of thumb) (Table 6). Standardizing the method to determine land cover fractions across sites is beyond the scope of this work, so users of metadata should be mindful of these differences.

Different methods for estimating surface roughness length and zero-plane displacement height can give significantly different values (Kent et al., 2017). Given that different data providers have derived values using different methods, we also provide values using two consistent morphometric methods (Macdonald (Macdonald et al., 1998) and Kanda (Kanda et al., 2013); Table 5, parameters 26–29) derived from surface fraction and building height parameters within the measurement footprint (Table 6). The Kanda modification to the Macdonald method accounts for the variability in roughness element height, resulting in larger displacement heights which are closer to estimates made with anemometric methods (Kent et al., 2017). The Macdonald method assumes that all the buildings have the same average roughness element height. However, care must be taken when using Kanda values as some urban land surface models expect displacement height to be always lower than average building height (Hertwig et al., 2020).

Where not known, building height standard deviation ($\sigma_H$) is estimated from an empirical relation to building mean height ($H_{ave}$) (Kanda et al., 2013):

$$ \sigma_H = 1.05 H_{ave} - 3.7. \qquad (5) $$

Similarly, unknown local wall to plan area ratios ($\lambda_w$) are derived from roof area fraction ($\lambda_p$) and canyon height to width ratio ($H/W$) assuming an infinite canyon geometry (Masson et al., 2020):

$$ \lambda_w = 2(1 - \lambda_p) H/W. \qquad (6) $$
- `time_analysis_start`: start of observed (focus) period in UTC with period ending timestamps

- `time_shown_in`: time standard (always UTC)

- `local_utf_offsetHours`: offset in hours of local time from UTC

- `timestep_interval_seconds`: period of block averaging in seconds (time step)

- `timestep_number_spinup`: number of time steps prior to observed focus period

- `timestep_number_analysis`: number of time steps in observed focus period

- `projectContact`: contact details for the Urban-PLUMBER project coordinators

- `observationsContact`: contact details of the observational site data providers

- `observations_reference`: published references associated with the observations

- `date_created`: date and time of creation of this file

- `source`: repository for processing code

- `comment`: additional comments associated with this dataset (e.g. excluded wind sectors).

Text file time series include metadata headers indicated with a hash (#) at line beginnings. Columns are headed by variable names in ALMA format (Table 2). NetCDF4 files include identical data, with additional attributes for each variable:

- `long_name`: plain language description of variable

- `standard_name`: equivalent variable name under the CF (climate and forecast) conventions

- `units`: SI (international system) units

- `ancillary_variables`: name of the associated quality control flag variable.

NetCDF files also include site characteristics parameter values and descriptions (Table 5). Times in all datasets are UTC. Python programmes are provided (Lipson, 2022a) to convert UTC times to local standard time (`convert_UTF_to_local_time.py`), and netCDF to text (`convert_nc_to_text.py`).

## 6.3 Site characteristics metadata

Site characteristics (Tables 5, 6) are essential for any use of these data, and fundamental to application of land surface models. These metadata are provided in two machine readable forms (plain text in csv files and netCDF4). The metadata are primarily drawn from published sources or as advised by the data providers. If local parameters are not known, values are estimated from high-resolution global datasets or derived from empirical relations. The sources for each parameter are included within the site characteristic metadata.

There are numerous methods to estimate the probable extent and weighting of turbulent fluxes footprints relative to the eddy covariance sensors located on a flux tower (Velasco and Roth, 2010). The eddy covariance flux footprint provides a basis to identify which area (and weighting) should be used to estimate the land surface fractions impacting the measurements. Some studies in this collection determine the footprint and resulting land cover fractions dynamically (e.g. for each 30 min period based on that period's observed atmospheric variables such as stability and wind direction), whereas others used a constant radius (e.g. based on the footprint climatology or rule of thumb) (Table 6). Standardizing the method to determine land cover fractions across sites is beyond the scope of this work, so users of metadata should be mindful of these differences.

Different methods for estimating surface roughness length and zero-plane displacement height can give significantly different values (Kent et al., 2017). Given that different data providers have derived values using different methods, we also provide values using two consistent morphometric methods (Macdonald (Macdonald et al., 1998) and Kanda (Kanda et al., 2013); Table 5, parameters 26–29) derived from surface fraction and building height parameters within the measurement footprint (Table 6). The Kanda modification to the Macdonald method accounts for the variability in roughness element height, resulting in larger displacement heights which are closer to estimates made with anemometric methods (Kent et al., 2017). The Macdonald method assumes that all the buildings have the same average roughness element height. However, care must be taken when using Kanda values as some urban land surface models expect displacement height to be always lower than average building height (Hertwig et al., 2020).

Where not known, building height standard deviation ($\sigma_H$) is estimated from an empirical relation to building mean height ($H_{ave}$) (Kanda et al., 2013):

$$ \sigma_H = 1.05 H_{ave} - 3.7. \qquad (5) $$

Similarly, unknown local wall to plan area ratios ($\lambda_w$) are derived from roof area fraction ($\lambda_p$) and canyon height to width ratio ($H/W$) assuming an infinite canyon geometry (Masson et al., 2020):

$$ \lambda_w = 2(1 - \lambda_p) H/W. \qquad (6) $$
<table><caption>**Table 5.** Site characteristic metadata description and units. Parameters are determined for the turbulent flux footprint extent (Table 6), except 1–4 which are applicable to the tower itself, and 19 which is a function of the radiometer field of view (Offerle et al., 2003) and differs from the turbulent flux footprint (Schmid et al., 1991).</caption><thead><tr><th>ID</th><th>Parameter</th><th>Units</th><th>Description</th></tr></thead><tbody><tr><td>1</td><td>latitude</td><td>degrees_north</td><td>Latitude of tower</td></tr><tr><td>2</td><td>longitude</td><td>degrees_east</td><td>Longitude of tower</td></tr><tr><td>3</td><td>ground_height</td><td>m</td><td>Height above sea level of base of tower</td></tr><tr><td>4</td><td>measurement_heightAbove</td><td>m</td><td>Height above ground level (agl) of eddy covariance equipment on tower</td></tr><tr><td>5</td><td>impervious_area_fractions</td><td>1</td><td>Plan area fraction of all impervious (hard) surfaces, including roofs, roads, paths and paved areas</td></tr><tr><td>6</td><td>tree_area_fractions</td><td>1</td><td>Plan area fraction of tree canopy (&gt;2 m)</td></tr><tr><td>7</td><td>grass_area_fractions</td><td>1</td><td>Plan area fraction of grass or other vegetation (&lt;2 m)</td></tr><tr><td>8</td><td>bare_soil_area_fractions</td><td>1</td><td>Plan area fraction of bare soil</td></tr><tr><td>9</td><td>water_area_fractions</td><td>1</td><td>Plan area fraction of water</td></tr><tr><td>10</td><td>roof_area_fractions</td><td>1</td><td>Plan area fraction of roofs (λ<sub>p</sub>)</td></tr><tr><td>11</td><td>road_area_fractions</td><td>1</td><td>Plan area fraction of roads</td></tr><tr><td>12</td><td>other_paved_area_fractions</td><td>1</td><td>Plan area fraction of hard surfaces on ground excluding roads (e.g. paths, plazas, car parks)</td></tr><tr><td>13</td><td>building_mean_height</td><td>m</td><td>Mean height above ground of buildings (H<sub>ave</sub>)</td></tr><tr><td>14</td><td>tree_mean_height</td><td>m</td><td>Mean height above ground of trees</td></tr><tr><td>15</td><td>roughness_length_momentum</td><td>m</td><td>Aerodynamic roughness length for momentum as reported in the literature or advised by data providers</td></tr><tr><td>16</td><td>displacement_height</td><td>m</td><td>Zero-plane displacement height as reported in literature or advised by data providers</td></tr><tr><td>17</td><td>canyon_height_width_ratio</td><td>1</td><td>Mean building height to mean street canyon width (distance between buildings) ratio (H/W)</td></tr><tr><td>18</td><td>wall_to_plan_area_ratio</td><td>1</td><td>Sum of wall surface area to plan area ratio (λ<sub>w</sub>)</td></tr><tr><td>19</td><td>average_albedo_at_midday</td><td>1</td><td>Median site albedo at midday (local standard time) for available observations</td></tr><tr><td>20</td><td>residentPopulation_density</td><td>person km<sup>-2</sup></td><td>Resident (night) population density</td></tr><tr><td>21</td><td>anthropogenic_heat flux_mean</td><td>W m<sup>-2</sup></td><td>Anthropogenic heat flux annual mean</td></tr><tr><td>22</td><td>topsoil_clay_fractions</td><td>1</td><td>Clay fraction of topsoil</td></tr><tr><td>23</td><td>topsoil_sand_fractions</td><td>1</td><td>Sand fraction of topsoil</td></tr><tr><td>24</td><td>topsoil_bulk_density</td><td>kg m<sup>-3</sup></td><td>Bulk (dry) density of topsoil</td></tr><tr><td>25</td><td>building_height_standard_deviation</td><td>m</td><td>Aerodynamic roughness length for momentum calculated by the Macdonald morphometric method</td></tr><tr><td>26</td><td>roughness_length_momentumMac</td><td>m</td><td>Zero-plane displacement height calculated by the Macdonald morphometric method</td></tr><tr><td>27</td><td>displacement_heightMac</td><td>m</td><td>Zero-plane displacement height calculated by the Macdonald morphometric method</td></tr><tr><td>28</td><td>roughness_length_momentumKanda</td><td>m</td><td>Aerodynamic roughness length for momentum calculated by the Kanda morphometric method</td></tr><tr><td>29</td><td>displacement_heightKanda</td><td>m</td><td>Zero-plane displacement height calculated by the Kanda morphometric method</td></tr></tbody></table>

Frontal area index ($\lambda_f$) is sometimes reported in site literature without $H/W$ or $\lambda_w$, in which case these are estimated (again assuming an infinite canyon geometry) with (Porson et al., 2010):

$$ \lambda_f = \frac{2}{\pi} (1 - \lambda_p) H/W. \qquad (7) $$

Where not known or provided, mean annual anthropogenic heat flux (Varquez et al., 2021) or soil characteristics (Hengl, 2018a, b, c) are estimated from global datasets at 1 km or lower resolutions.

## 6.4 Data flags

Each variable for each time step has a quality control (qc) flag. For example, LWdown_qc lists qc flags for LWdown at each time step. Flag numbers are consistent across all variables:

0. observed by measurement at site and passes project quality control tests
<table><caption>**Table 5.** Site characteristic metadata description and units. Parameters are determined for the turbulent flux footprint extent (Table 6), except 1–4 which are applicable to the tower itself, and 19 which is a function of the radiometer field of view (Offerle et al., 2003) and differs from the turbulent flux footprint (Schmid et al., 1991).</caption><thead><tr><th>ID</th><th>Parameter</th><th>Units</th><th>Description</th></tr></thead><tbody><tr><td>1</td><td>latitude</td><td>degrees_north</td><td>Latitude of tower</td></tr><tr><td>2</td><td>longitude</td><td>degrees_east</td><td>Longitude of tower</td></tr><tr><td>3</td><td>ground_height</td><td>m</td><td>Height above sea level of base of tower</td></tr><tr><td>4</td><td>measurement_heightAbove-ground</td><td>m</td><td>Height above ground level (agl) of eddy covariance equipment on tower</td></tr><tr><td>5</td><td>impervious_area_fraction</td><td>1</td><td>Plan area fraction of all impervious (hard) surfaces, including roofs, roads, paths and paved areas</td></tr><tr><td>6</td><td>tree_area_fraction</td><td>1</td><td>Plan area fraction of tree canopy (&gt;2 m)</td></tr><tr><td>7</td><td>grass_area_fraction</td><td>1</td><td>Plan area fraction of grass or other vegetation (&lt;2 m)</td></tr><tr><td>8</td><td>bare_soil_area_fraction</td><td>1</td><td>Plan area fraction of bare soil</td></tr><tr><td>9</td><td>water_area_fraction</td><td>1</td><td>Plan area fraction of water</td></tr><tr><td>10</td><td>roof_area_fraction</td><td>1</td><td>Plan area fraction of roofs (λ<sub>p</sub>)</td></tr><tr><td>11</td><td>road_area_fraction</td><td>1</td><td>Plan area fraction of roads</td></tr><tr><td>12</td><td>other_paved_area_fraction</td><td>1</td><td>Plan area fraction of hard surfaces on ground excluding roads (e.g. paths, plazas, car parks)</td></tr><tr><td>13</td><td>building_mean_height</td><td>m</td><td>Mean height above ground of buildings (H<sub>ave</sub>)</td></tr><tr><td>14</td><td>tree_mean_height</td><td>m</td><td>Mean height above ground of trees</td></tr><tr><td>15</td><td>roughness_length_momentum</td><td>m</td><td>Aerodynamic roughness length for momentum as reported in the literature or advised by data providers</td></tr><tr><td>16</td><td>displacement_height</td><td>m</td><td>Zero-plane displacement height as reported in literature or advised by data providers</td></tr><tr><td>17</td><td>canyon_height_width_ratio</td><td>1</td><td>Mean building height to mean street canyon width (distance between buildings) ratio (H/W)</td></tr><tr><td>18</td><td>wall_to_plan_area_ratio</td><td>1</td><td>Sum of wall surface area to plan area ratio (λ<sub>w</sub>)</td></tr><tr><td>19</td><td>average_albedo_at_midday</td><td>1</td><td>Median site albedo at midday (local standard time) for available observations</td></tr><tr><td>20</td><td>residentPopulation_density</td><td>person km<sup>-2</sup></td><td>Resident (night) population density</td></tr><tr><td>21</td><td>anthropogenic_heat flux_mean</td><td>W m<sup>-2</sup></td><td>Anthropogenic heat flux annual mean</td></tr><tr><td>22</td><td>topsoil_clay_fraction</td><td>1</td><td>Clay fraction of topsoil</td></tr><tr><td>23</td><td>topsoil_sand_fraction</td><td>1</td><td>Sand fraction of topsoil</td></tr><tr><td>24</td><td>topsoilBulk_density</td><td>kg m<sup>-3</sup></td><td>Bulk (dry) density of topsoil</td></tr><tr><td>25</td><td>building_height_standard_deviation</td><td>m</td><td>Aerodynamic roughness length for momentum calculated by the Macdonald morphometric method</td></tr><tr><td>26</td><td>roughness_length_momentum Mac</td><td>m</td><td>Zero-plane displacement height calculated by the Macdonald morphometric method</td></tr><tr><td>27</td><td>displacement_height Mac</td><td>m</td><td>Zero-plane displacement height calculated by the Macdonald morphometric method</td></tr><tr><td>28</td><td>roughness_length_momentum_kanda</td><td>m</td><td>Aerodynamic roughness length for momentum calculated by the Kanda morphometric method</td></tr><tr><td>29</td><td>displacement_height_kanda</td><td>m</td><td>Zero-plane displacement height calculated by the Kanda morphometric method</td></tr></tbody></table>

Frontal area index ($\lambda_f$) is sometimes reported in site literature without $H/W$ or $\lambda_w$, in which case these are estimated (again assuming an infinite canyon geometry) with (Porson et al., 2010):

$$ \lambda_f = \frac{2}{\pi} (1 - \lambda_p) H/W. \qquad (7) $$

Where not known or provided, mean annual anthropogenic heat flux (Varquez et al., 2021) or soil characteristics (Hengl, 2018a, b, c) are estimated from global datasets at 1 km or lower resolutions.

## 6.4 Data flags

Each variable for each time step has a quality control (qc) flag. For example, LWdown_qc lists qc flags for LWdown at each time step. Flag numbers are consistent across all variables:

0. observed by measurement at site and passes project quality control tests
<table><caption>**Table 6.** Select site characteristic values (see Table 5 for definitions). Other site characteristic values and sources are provided within the collection (Lipson et al., 2022). Areas analysed for land cover fractions and roughness parameters are based on either a static radius around the flux tower (value given) or a dynamic footprint model (fpm). For the latter, the spatial extents are the order of a few hundred metres but are dynamic varying for example with atmospheric stability and wind direction (Grimmond and Ward, 2021).</caption><thead><tr><th>Parameter</th><th>flux footprint extent</th><th>ground height</th><th>measurement height above ground</th><th>impervious area fraction</th><th>tree_area fraction</th><th>grass_area fraction</th><th>bare_soil_area fraction</th><th>water Area fraction</th><th>roof Area fraction</th><th>road Area fraction</th><th>other paved Area fraction</th><th>building mean height</th></tr></thead><tbody><tr><td>AU-Preston</td><td>500 m</td><td>93</td><td>40</td><td>0.620</td><td>0.225</td><td>0.150</td><td>0.005</td><td>0</td><td>0.445</td><td>0.130</td><td>0.045</td><td>6.4</td></tr><tr><td>AU-SurreyHills</td><td>500 m</td><td>97</td><td>38</td><td>0.54</td><td>0.29</td><td>0.15</td><td>0.01</td><td>0.01</td><td>0.39</td><td>0.09</td><td>0.06</td><td>7.2</td></tr><tr><td>CA-Sunset</td><td>fpm</td><td>78</td><td>24.8</td><td>0.68</td><td>0.12</td><td>0.20</td><td>0</td><td>0</td><td>0.23</td><td>0.20</td><td>0.25</td><td>4.9</td></tr><tr><td>FI-Kumpula</td><td>1000 m</td><td>29</td><td>31</td><td>0.46</td><td>0.30</td><td>0.24</td><td>0</td><td>0</td><td>0.14</td><td>0.32</td><td>0</td><td>12.6</td></tr><tr><td>FI-Torni</td><td>1000 m</td><td>15.2</td><td>60</td><td>0.77</td><td>0.15</td><td>0.07</td><td>0</td><td>0.01</td><td>0.37</td><td>0.25</td><td>0.15</td><td>17.9</td></tr><tr><td>FR-Capitole</td><td>500 m</td><td>143</td><td>48.05</td><td>0.90</td><td>0.08</td><td>0.02</td><td>0</td><td>0</td><td>0.62</td><td>0.28</td><td>0</td><td>15</td></tr><tr><td>GR-HECKOR</td><td>fpm</td><td>30</td><td>27</td><td>0.916</td><td>0.040</td><td>0.016</td><td>0.010</td><td>0.019</td><td>0.516</td><td>0.201</td><td>0.199</td><td>11.3</td></tr><tr><td>JP-Yoyogi</td><td>500 m</td><td>39</td><td>52</td><td>0.92</td><td>0.06</td><td>0.01</td><td>0.01</td><td>0</td><td>0.41</td><td>0.32</td><td>0.19</td><td>9.0</td></tr><tr><td>KR-Jungnang</td><td>500 m</td><td>22</td><td>41.5</td><td>0.965</td><td>0</td><td>0.019</td><td>0.016</td><td>0</td><td>0.588</td><td>0.377</td><td>0</td><td>8.648</td></tr><tr><td>KR-Ochang</td><td>500 m</td><td>60</td><td>19</td><td>0.470</td><td>0.184</td><td>0.333</td><td>0.013</td><td>0</td><td>0.133</td><td>0.337</td><td>0</td><td>7.384</td></tr><tr><td>MX-Escandon</td><td>fpm</td><td>22.40</td><td>37</td><td>0.94</td><td>0.06</td><td>0</td><td>0</td><td>0</td><td>0.57</td><td>0.37</td><td>0</td><td>9.69</td></tr><tr><td>NL-Amsterdam</td><td>500 m</td><td>0</td><td>40</td><td>0.68</td><td>0.15</td><td>0</td><td>0</td><td>0.17</td><td>0.44</td><td>0.07</td><td>0.17</td><td>14.2</td></tr><tr><td>PL-Lipowa</td><td>fpm</td><td>204</td><td>37</td><td>0.76</td><td>0.16</td><td>0.08</td><td>0</td><td>0</td><td>0.35</td><td>0.21</td><td>0.20</td><td>10.2</td></tr><tr><td>PL-Narutowicza</td><td>500 m</td><td>221</td><td>42</td><td>0.65</td><td>0.22</td><td>0.09</td><td>0.04</td><td>0</td><td>0.29</td><td>0.19</td><td>0.17</td><td>16</td></tr><tr><td>SG-TelokKurau06</td><td>1000 m</td><td>5</td><td>20.7</td><td>0.85</td><td>0.11</td><td>0.04</td><td>0</td><td>0</td><td>0.39</td><td>0.12</td><td>0.34</td><td>9.9</td></tr><tr><td>UK-KingsCollege</td><td>fpm</td><td>14.5</td><td>50.3</td><td>0.79</td><td>0.03</td><td>0.04</td><td>0</td><td>0.14</td><td>0.40</td><td>0.39</td><td>0</td><td>21.3</td></tr><tr><td>UK-Swindon</td><td>500 m</td><td>108</td><td>12.5</td><td>0.49</td><td>0.09</td><td>0.36</td><td>0.06</td><td>0</td><td>0.16</td><td>0.15</td><td>0.18</td><td>4.5</td></tr><tr><td>US-Baltimore</td><td>1000 m</td><td>157</td><td>37.2</td><td>0.313</td><td>0.536</td><td>0.138</td><td>0.007</td><td>0.006</td><td>0.160</td><td>0.153</td><td>0</td><td>5.6</td></tr><tr><td>US-Minneapolis1</td><td>fpm</td><td>301</td><td>40</td><td>0.21</td><td>0.38</td><td>0.36</td><td>0</td><td>0.05</td><td>0.12</td><td>0.05</td><td>0.04</td><td>5.05</td></tr><tr><td>US-Minneapolis2</td><td>fpm</td><td>301</td><td>40</td><td>0.05</td><td>0.2</td><td>0.73</td><td>0</td><td>0.02</td><td>0.01</td><td>0</td><td>0.04</td><td>5.05</td></tr><tr><td>US-WestPhoenix</td><td>fpm</td><td>340</td><td>22.1</td><td>0.48</td><td>0.05</td><td>0.10</td><td>0.37</td><td>0</td><td>0.26</td><td>0.22</td><td>0</td><td>4.5</td></tr></tbody></table>
<table><caption>**Table 6.** Select site characteristic values (see Table 5 for definitions). Other site characteristic values and sources are provided within the collection (Lipson et al., 2022). Areas analysed for land cover fractions and roughness parameters are based on either a static radius around the flux tower (value given) or a dynamic footprint model (fpm). For the latter, the spatial extents are the order of a few hundred metres but are dynamic varying for example with atmospheric stability and wind direction (Grimmond and Ward, 2021).</caption><thead><tr><th>Parameter</th><th>flux footprint extent</th><th>ground height</th><th>measurement height above ground</th><th>impervious area fraction</th><th>tree_area fraction</th><th>grass_area fraction</th><th>bare_soil_area fraction</th><th>water Area fraction</th><th>roof Area fraction</th><th>road Area fraction</th><th>other paved Area fraction</th><th>building mean height</th></tr></thead><tbody><tr><td>AU-Preston</td><td>500 m</td><td>93</td><td>40</td><td>0.620</td><td>0.225</td><td>0.150</td><td>0.005</td><td>0</td><td>0.445</td><td>0.130</td><td>0.045</td><td>6.4</td></tr><tr><td>AU-Surrey Hills</td><td>500 m</td><td>97</td><td>38</td><td>0.54</td><td>0.29</td><td>0.15</td><td>0.01</td><td>0.01</td><td>0.39</td><td>0.09</td><td>0.06</td><td>7.2</td></tr><tr><td>CA-Sunset</td><td>fpm</td><td>78</td><td>24.8</td><td>0.68</td><td>0.12</td><td>0.20</td><td>0</td><td>0</td><td>0.23</td><td>0.20</td><td>0.25</td><td>4.9</td></tr><tr><td>FI-Kumpula</td><td>1000 m</td><td>29</td><td>31</td><td>0.46</td><td>0.30</td><td>0.24</td><td>0</td><td>0</td><td>0.14</td><td>0.32</td><td>0</td><td>12.6</td></tr><tr><td>FI-Torni</td><td>1000 m</td><td>15.2</td><td>60</td><td>0.77</td><td>0.15</td><td>0.07</td><td>0</td><td>0.01</td><td>0.37</td><td>0.25</td><td>0.15</td><td>17.9</td></tr><tr><td>FR-Capitole</td><td>500 m</td><td>143</td><td>48.05</td><td>0.90</td><td>0.08</td><td>0.02</td><td>0</td><td>0</td><td>0.62</td><td>0.28</td><td>0</td><td>15</td></tr><tr><td>GR-HECKOR</td><td>fpm</td><td>30</td><td>27</td><td>0.916</td><td>0.040</td><td>0.016</td><td>0.010</td><td>0.019</td><td>0.516</td><td>0.201</td><td>0.199</td><td>11.3</td></tr><tr><td>JP-Yoyogi</td><td>500 m</td><td>39</td><td>52</td><td>0.92</td><td>0.06</td><td>0.01</td><td>0.01</td><td>0</td><td>0.41</td><td>0.32</td><td>0.19</td><td>9.0</td></tr><tr><td>KR-Jungnang</td><td>500 m</td><td>22</td><td>41.5</td><td>0.965</td><td>0</td><td>0.019</td><td>0.016</td><td>0</td><td>0.588</td><td>0.377</td><td>0</td><td>8.648</td></tr><tr><td>KR-Ochang</td><td>500 m</td><td>60</td><td>19</td><td>0.470</td><td>0.184</td><td>0.333</td><td>0.013</td><td>0</td><td>0.133</td><td>0.337</td><td>0</td><td>7.384</td></tr><tr><td>MX-Escandon</td><td>fpm</td><td>22.40</td><td>37</td><td>0.94</td><td>0.06</td><td>0</td><td>0</td><td>0</td><td>0.57</td><td>0.37</td><td>0</td><td>9.69</td></tr><tr><td>NL-Amsterdam</td><td>500 m</td><td>0</td><td>40</td><td>0.68</td><td>0.15</td><td>0</td><td>0</td><td>0.17</td><td>0.44</td><td>0.07</td><td>0.17</td><td>14.2</td></tr><tr><td>PL-Lipowa</td><td>fpm</td><td>204</td><td>37</td><td>0.76</td><td>0.16</td><td>0.08</td><td>0</td><td>0</td><td>0.35</td><td>0.21</td><td>0.20</td><td>10.2</td></tr><tr><td>PL-Narutowicza</td><td>500 m</td><td>221</td><td>42</td><td>0.65</td><td>0.22</td><td>0.09</td><td>0.04</td><td>0</td><td>0.29</td><td>0.19</td><td>0.17</td><td>16</td></tr><tr><td>SG-TelokKurau06</td><td>1000 m</td><td>5</td><td>20.7</td><td>0.85</td><td>0.11</td><td>0.04</td><td>0</td><td>0</td><td>0.39</td><td>0.12</td><td>0.34</td><td>9.9</td></tr><tr><td>UK-KingsCollege</td><td>fpm</td><td>14.5</td><td>50.3</td><td>0.79</td><td>0.03</td><td>0.04</td><td>0</td><td>0.14</td><td>0.40</td><td>0.39</td><td>0</td><td>21.3</td></tr><tr><td>UK-Swindon</td><td>500 m</td><td>108</td><td>12.5</td><td>0.49</td><td>0.09</td><td>0.36</td><td>0.06</td><td>0</td><td>0.16</td><td>0.15</td><td>0.18</td><td>4.5</td></tr><tr><td>US-Baltimore</td><td>1000 m</td><td>157</td><td>37.2</td><td>0.313</td><td>0.536</td><td>0.138</td><td>0.007</td><td>0.006</td><td>0.160</td><td>0.153</td><td>0</td><td>5.6</td></tr><tr><td>US-Minneapolis1</td><td>fpm</td><td>301</td><td>40</td><td>0.21</td><td>0.38</td><td>0.36</td><td>0</td><td>0.05</td><td>0.12</td><td>0.05</td><td>0.04</td><td>5.05</td></tr><tr><td>US-Minneapolis2</td><td>fpm</td><td>301</td><td>40</td><td>0.05</td><td>0.2</td><td>0.73</td><td>0</td><td>0.02</td><td>0.01</td><td>0</td><td>0.04</td><td>5.05</td></tr><tr><td>US-WestPhoenix</td><td>fpm</td><td>340</td><td>22.1</td><td>0.48</td><td>0.05</td><td>0.10</td><td>0.37</td><td>0</td><td>0.26</td><td>0.22</td><td>0</td><td>4.5</td></tr></tbody></table>
<table><caption>**Table 7.** Site wind sector exclusions. Sites with sensible and latent heat fluxes excluded because of land cover or land use differences by wind sectors as described in the reference provided. Maps of these sectors are provided in the site data collection (Lipson et al., 2022).</caption><thead><tr><th>Site name</th><th>Sectors excluded</th><th>Reason</th><th>Reference</th></tr></thead><tbody><tr><td>FI-Kumpula</td><td>0–180°, 320–360°</td><td>surface inhomogeneity</td><td>Karsisto et al. (2016)</td></tr><tr><td>FI-Torni</td><td>40–150°</td><td>flow interference from tower</td><td>Järvi et al. (2018)</td></tr><tr><td>JP-Yoyogi</td><td>170–260°</td><td>surface inhomogeneity</td><td>Ishidoya et al. (2020)</td></tr><tr><td>US-Minneapolis1</td><td>75–285°</td><td>surface inhomogeneity</td><td>Menzer and McFadden (2017)</td></tr><tr><td>US-Minneapolis2</td><td>0–180°, 270–360°</td><td>120–180°: flow interference from tower<br>270–360°, 0–120°: surface inhomogeneity</td><td>Menzer and McFadden (2017)</td></tr></tbody></table>

<table><caption>**Table 8.** Funding acknowledgements for individual sites.</caption><thead><tr><th>Site</th><th>Contributing author</th><th>Site funding acknowledgements</th></tr></thead><tbody><tr><td>AU-Preston</td><td>Andrew Coutts, Nigel Tapper</td><td>–</td></tr><tr><td>AU-SurreyHills</td><td>Andrew Coutts, Nigel Tapper</td><td>–</td></tr><tr><td>CA-Sunset</td><td>Andreas Christen, Oliver Michels</td><td>Canadian Foundation for Climate and Atmospheric Sciences (CFCAS, Project “Environmental Prediction in Canadian Cities (EpiCC)”) and the Natural and Engineering Research Council of Canada (NSERC, RGPIN-03958, RGPAS-507854). Some instruments were supported by the Canada Foundation for Innovation (CFI, IF 2015, grant no. 33600) and BCKDF. We acknowledge the support of BC Hydro to operate the tower.</td></tr><tr><td>FI-Kumpula</td><td>Leena Järvi</td><td>ICOS Finland</td></tr><tr><td>FI-Torni</td><td>Leena Järvi</td><td>ICOS Finland</td></tr><tr><td>FR-Capitole</td><td>Valéry Masson</td><td>Météo-France and CNRS</td></tr><tr><td>GR-HECKOR</td><td>Nektarios Chrysoulakis</td><td>EU Horizon 2020 Research and Innovation Programme, under grant agreement no. 870337 project CURE (http://cure-copernicus.eu, last access: 7 November 2022)</td></tr><tr><td>JP-Yoyogi</td><td>Hirofumi Sugawara</td><td>Japan Society for the Promotion of Science KAKENHI grants (nos. 24241008, 15H02814, 18K01129 and 19H01975), and the Environment Research and Technology Development Fund (JPMEERF20191009) of the Environmental Restoration and Conservation Agency of Japan</td></tr><tr><td>KR-Jungnang</td><td>Jinkyu Hong, Sungsoo Jo, Yeon-Hee</td><td>Korea Meteorological Administration Research and Development Program “Development of Production Techniques on User-Customized Weather information” under grant (KMA2018-00622) and National Research Foundation of Korea (NRF) grants funded by the Korean government (NRF-2018R1A5A1024958)</td></tr><tr><td>KR-Ochang</td><td>Jinkyu Hong, Je-Woo Hong, Keunmin Lee</td><td>Korea Meteorological Administration Research and Development Program under grant KMI2021-01610</td></tr><tr><td>MX-Escandon</td><td>Erik Velasco</td><td>National Institute of Ecology and Climate Change (INECC) and the Mexico City’s Secretariat for the Environment (SEDEMA) through the Molina Center for Energy and the Environment (MCE2)</td></tr><tr><td>NL-Amsterdam</td><td>Bert Heusinkveld</td><td>Netherlands Organisation for Scientific Research (NWO) project 864.14.007 and the Amsterdam Institute for Advanced Metropolitan Solutions (AMS) project VIR16002.</td></tr><tr><td>PL-Lipowa</td><td>Wlodzimierz Pawlak, Krzysztof Fortuniak</td><td>University of Lodz, Poland, IDUB grant – number of decision 57/2021</td></tr><tr><td>PL-Narutowicza</td><td>Wlodzimierz Pawlak, Krzysztof Fortuniak</td><td>University of Lodz, Poland, IDUB grant – number of decision 57/2021</td></tr><tr><td>SG-TelokKurau06</td><td>Matthias Roth</td><td>Ministry of Education, Singapore</td></tr><tr><td>UK-KingsCollege</td><td>Simone Kotthaus, Sue Grimmond</td><td>EUfp7 grant agreement no. 211345 (BRIDGE), NERC ClearfLo (NE/H003231/1), NERC ARSF (GB08/19), EPSRC (EP/I00159X/1, EP/I00159X/2) and KCL</td></tr><tr><td>UK-Swindon</td><td>Helen Ward, Jonathan Evans, Sue Grimmond</td><td>NERC NE/H52479X/1</td></tr><tr><td>US-Baltimore</td><td>Sue Grimmond, Ben Crawford</td><td>National Science Foundation (BCS-0095284, DEB-9714835) and USDA Forest Service</td></tr><tr><td>US-Minneapolis</td><td>Joeseph McFadden</td><td>NASA Earth Science Division (NNG04GN80G)</td></tr><tr><td>US-WestPhoenix</td><td>Stevan Earl, Winston Chow</td><td>National Science Foundation (DEB-1832016) Central Arizona–Phoenix Long-Term Ecological Research Program (CAP LTER)</td></tr></tbody></table>
<table><caption>**Table 7.** Site wind sector exclusions. Sites with sensible and latent heat fluxes excluded because of land cover or land use differences by wind sectors as described in the reference provided. Maps of these sectors are provided in the site data collection (Lipson et al., 2022).</caption><thead><tr><th>Site name</th><th>Sectors excluded</th><th>Reason</th><th>Reference</th></tr></thead><tbody><tr><td>FI-Kumpula</td><td>0–180°, 320–360°</td><td>surface inhomogeneity</td><td>Karsisto et al. (2016)</td></tr><tr><td>FI-Torni</td><td>40–150°</td><td>flow interference from tower</td><td>Järvi et al. (2018)</td></tr><tr><td>JP-Yoyogi</td><td>170–260°</td><td>surface inhomogeneity</td><td>Ishidoya et al. (2020)</td></tr><tr><td>US-Minneapolis1</td><td>75–285°</td><td>surface inhomogeneity</td><td>Menzer and McFadden (2017)</td></tr><tr><td>US-Minneapolis2</td><td>0–180°, 270–360°</td><td>120–180°: flow interference from tower<br>270–360°, 0–120°: surface inhomogeneity</td><td>Menzer and McFadden (2017)</td></tr></tbody></table>

<table><caption>**Table 8.** Funding acknowledgements for individual sites.</caption><thead><tr><th>Site</th><th>Contributing author</th><th>Site funding acknowledgements</th></tr></thead><tbody><tr><td>AU-Preston</td><td>Andrew Coutts, Nigel Tapper</td><td>–</td></tr><tr><td>AU-SurreyHills</td><td>Andrew Coutts, Nigel Tapper</td><td>–</td></tr><tr><td>CA-Sunset</td><td>Andreas Christen, Oliver Michels</td><td>Canadian Foundation for Climate and Atmospheric Sciences (CFCAS, Project “Environmental Prediction in Canadian Cities (EpiCC)”) and the Natural and Engineering Research Council of Canada (NSERC, RGPIN-03958, RGPAS-507854). Some instruments were supported by the Canada Foundation for Innovation (CFI, IF 2015, grant no. 33600) and BCKDF. We acknowledge the support of BC Hydro to operate the tower.</td></tr><tr><td>FI-Kumpula</td><td>Leena Järvi</td><td>ICOS Finland</td></tr><tr><td>FI-Torni</td><td>Leena Järvi</td><td>ICOS Finland</td></tr><tr><td>FR-Capitole</td><td>Valéry Masson</td><td>Météo-France and CNRS</td></tr><tr><td>GR-HECKOR</td><td>Nektarios Chrysoulakis</td><td>EU Horizon 2020 Research and Innovation Programme, under grant agreement no. 870337 project CURE (http://cure-copernicus.eu, last access: 7 November 2022)</td></tr><tr><td>JP-Yoyogi</td><td>Hirofumi Sugawara</td><td>Japan Society for the Promotion of Science KAKENHI grants (nos. 24241008, 15H02814, 18K01129 and 19H01975), and the Environment Research and Technology Development Fund (JPMEERF20191009) of the Environmental Restoration and Conservation Agency of Japan</td></tr><tr><td>KR-Jungnang</td><td>Jinkyu Hong, Sungsoo Jo, Yeon-Hee</td><td>Korea Meteorological Administration Research and Development Program “Development of Production Techniques on User-Customized Weather information” under grant (KMA2018-00622) and National Research Foundation of Korea (NRF) grants funded by the Korean government (NRF-2018R1A5A1024958)</td></tr><tr><td>KR-Ochang</td><td>Jinkyu Hong, Je-Woo Hong, Keunmin Lee</td><td>Korea Meteorological Administration Research and Development Program under grant KMI2021-01610</td></tr><tr><td>MX-Escandon</td><td>Erik Velasco</td><td>National Institute of Ecology and Climate Change (INECC) and the Mexico City’s Secretariat for the Environment (SEDEMA) through the Molina Center for Energy and the Environment (MCE2)</td></tr><tr><td>NL-Amsterdam</td><td>Bert Heusinkveld</td><td>Netherlands Organisation for Scientific Research (NWO) project 864.14.007 and the Amsterdam Institute for Advanced Metropolitan Solutions (AMS) project VIR16002.</td></tr><tr><td>PL-Lipowa</td><td>Wlodzimierz Pawlak, Krzysztof Fortuniak</td><td>University of Lodz, Poland, IDUB grant – number of decision 57/2021</td></tr><tr><td>PL-Narutowicza</td><td>Wlodzimierz Pawlak, Krzysztof Fortuniak</td><td>University of Lodz, Poland, IDUB grant – number of decision 57/2021</td></tr><tr><td>SG-TelokKurau06</td><td>Matthias Roth</td><td>Ministry of Education, Singapore</td></tr><tr><td>UK-KingsCollege</td><td>Simone Kotthaus, Sue Grimmond</td><td>EUfp7 grant agreement no. 211345 (BRIDGE), NERC ClearfLo (NE/H003231/1), NERC ARSF (GB08/19), EPSRC (EP/I00159X/1, EP/I00159X/2) and KCL</td></tr><tr><td>UK-Swindon</td><td>Helen Ward, Jonathan Evans, Sue Grimmond</td><td>NERC NE/H52479X/1</td></tr><tr><td>US-Baltimore</td><td>Sue Grimmond, Ben Crawford</td><td>National Science Foundation (BCS-0095284, DEB-9714835) and USDA Forest Service</td></tr><tr><td>US-Minneapolis</td><td>Joeseph McFadden</td><td>NASA Earth Science Division (NNG04GN80G)</td></tr><tr><td>US-WestPhoenix</td><td>Stevan Earl, Winston Chow</td><td>National Science Foundation (DEB-1832016) Central Arizona–Phoenix Long-Term Ecological Research Program (CAP LTER)</td></tr></tbody></table>
1. filled by observation; interpolated from site observations over short (2 h) periods OR filled by observations from nearby (<10 km) stations over longer periods

2. filled by ERA5; derived from ERA5 with site specific bias correction

3. missing or removed through quality control (occurs only in time series without gap filling).

## 6.5 Wind sector exclusions

Turbulent flux data are excluded from certain wind directions (Table 7) because of

- interference on flow from tower structure (as identified by data providers);

- markedly different land cover characteristics from sectors of interest (with guidance from data providers).

The US-Minneapolis site has different surface cover by wind direction but is retained in the collection because of both a long observation period and its distinct land cover characteristics. Following previous studies (Menzer and McFadden, 2017) we subdivide these data into low density residential area (northern sectors, US-Minneapolis2) and irrigated grassland with few built structures (south, US-Minneapolis1). Each are given their own site time series and metadata, resulting in 21 datasets.

**Author contributions.** ML, SG and MB conceived and coordinated the project, and prepared the protocols for contributing authors. ML collated site datasets, wrote processing code, developed and undertook analysis of bias correction methods, prepared figures, prepared datasets and drafted the paper with guidance from SG and MB. All other authors (listed alphabetically) collected primary data, prepared site information, processed datasets for inclusion in the collection and contributed to the paper. Table 8 lists contributing author site affiliation.

**Competing interests.** The contact author has declared that none of the authors has any competing interests.

**Disclaimer.** Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Acknowledgements.** We would like to thank the vast number of people involved in day-to-day running of these sites that have been involved in instrument and tower installations (permitting, purchasing and site installation), routine (and unexpected event) maintenance, data collection, routine data processing and final data processing. We thank all those that have provided sites for the towers to be located and sometimes power and internet access. We ac-

knowledge the essential funding for the instrumentation and other infrastructure, for staff (administrative, technical and scientific) and students for these activities. We also thank those who offered data for use in this project which are not included at this time.

**Financial support.** The project coordinating team has been supported by UNSW Sydney and the Australian Research Council (ARC) Centre of Excellence for Climate System Science (grant no. CE110001028), University of Reading, the Met Office and ERC urbisphere 855005. Computation support has been provided by the ARC Centre of Excellence for Climate Extremes (grant no. CE170100023) and National Computational Infrastructure (NCI) Australia. It contains modified Copernicus Climate Change Service Information. Site-affiliated acknowledgements are listed in Table 8.

**Review statement.** This paper was edited by Hanqin Tian and reviewed by two anonymous referees.

## References

Albergel, C., Dutra, E., Munier, S., Calvet, J.-C., Munoz-Sabater, J., de Rosnay, P., and Balsamo, G.: ERA-5 and ERA-Interim driven ISBA land surface model simulations: which one performs better?, Hydrol. Earth Syst. Sci., 22, 3515–3532, https://doi.org/10.5194/hess-22-3515-2018, 2018.

Arnfield, A. J.: Two decades of urban climate research: a review of turbulence, exchanges of energy and water, and the urban heat island, Int. J. Climatol., 23, 1–26, https://doi.org/10.1002/joc.859, 2003.

Aubinet, M., Vesala, T., and Papale, D.: Eddy Covariance: A Practical Guide to Measurement and Data Analysis, 1st edn., edited by: Aubinet, M., Vesala, T., and Papale, D., Springer Science & Business Media, 451 pp., ISBN 978-94-007-9630-0, 2012.

Baldocchi, D. D.: How eddy covariance flux measurements have contributed to our understanding of *Global Change Biology*, Glob. Change Biol., 26, 242–260, https://doi.org/10.1111/gcb.14807, 2020.

Barlow, J. F.: Progress in observing and modelling the urban boundary layer, Urban Climate, 10, Part 2, 216–240, https://doi.org/10.1016/j.uclim.2014.03.011, 2014.

Beck, H. E., Zimmermann, N. E., McVicar, T. R., Vergopolan, N., Berg, A., and Wood, E. F.: Present and future Köppen-Geiger climate classification maps at 1-km resolution, Scientific Data, 5, 180214, https://doi.org/10.1038/sdata.2018.214, 2018.

Beringer, J., Hutley, L. B., McHugh, I., Arndt, S. K., Campbell, D., Cleugh, H. A., Cleverly, J., Resco de Dios, V., Eamus, D., Evans, B., Ewenz, C., Grace, P., Griebel, A., Haverd, V., Hinko-Najera, N., Huete, A., Isaac, P., Kanniah, K., Leuning, R., Liddell, M. J., Macfarlane, C., Meyer, W., Moore, C., Pendall, E., Phillips, A., Phillips, R. L., Prober, S. M., Restrepo-Coupe, N., Rutledge, S., Schroder, I., Silberstein, R., Southall, P., Yee, M. S., Tapper, N. J., van Gorsel, E., Vote, C., Walker, J., and Wardlaw, T.: An introduction to the Australian and New Zealand flux tower network – OzFlux, Biogeosciences, 13, 5895–5916, https://doi.org/10.5194/bg-13-5895-2016, 2016.
1. filled by observation; interpolated from site observations over short (2 h) periods OR filled by observations from nearby (<10 km) stations over longer periods

2. filled by ERA5; derived from ERA5 with site specific bias correction

3. missing or removed through quality control (occurs only in time series without gap filling).

## 6.5 Wind sector exclusions

Turbulent flux data are excluded from certain wind directions (Table 7) because of

- interference on flow from tower structure (as identified by data providers);

- markedly different land cover characteristics from sectors of interest (with guidance from data providers).

The US-Minneapolis site has different surface cover by wind direction but is retained in the collection because of both a long observation period and its distinct land cover characteristics. Following previous studies (Menzer and McFadden, 2017) we subdivide these data into low density residential area (northern sectors, US-Minneapolis2) and irrigated grassland with few built structures (south, US-Minneapolis1). Each are given their own site time series and metadata, resulting in 21 datasets.

**Author contributions.** ML, SG and MB conceived and coordinated the project, and prepared the protocols for contributing authors. ML collated site datasets, wrote processing code, developed and undertook analysis of bias correction methods, prepared figures, prepared datasets and drafted the paper with guidance from SG and MB. All other authors (listed alphabetically) collected primary data, prepared site information, processed datasets for inclusion in the collection and contributed to the paper. Table 8 lists contributing author site affiliation.

**Competing interests.** The contact author has declared that none of the authors has any competing interests.

**Disclaimer.** Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Acknowledgements.** We would like to thank the vast number of people involved in day-to-day running of these sites that have been involved in instrument and tower installations (permitting, purchasing and site installation), routine (and unexpected event) maintenance, data collection, routine data processing and final data processing. We thank all those that have provided sites for the towers to be located and sometimes power and internet access. We ac-

knowledge the essential funding for the instrumentation and other infrastructure, for staff (administrative, technical and scientific) and students for these activities. We also thank those who offered data for use in this project which are not included at this time.

**Financial support.** The project coordinating team has been supported by UNSW Sydney and the Australian Research Council (ARC) Centre of Excellence for Climate System Science (grant no. CE110001028), University of Reading, the Met Office and ERC urbisphere 855005. Computation support has been provided by the ARC Centre of Excellence for Climate Extremes (grant no. CE170100023) and National Computational Infrastructure (NCI) Australia. It contains modified Copernicus Climate Change Service Information. Site-affiliated acknowledgements are listed in Table 8.

**Review statement.** This paper was edited by Hanqin Tian and reviewed by two anonymous referees.

## References

Albergel, C., Dutra, E., Munier, S., Calvet, J.-C., Munoz-Sabater, J., de Rosnay, P., and Balsamo, G.: ERA-5 and ERA-Interim driven ISBA land surface model simulations: which one performs better?, Hydrol. Earth Syst. Sci., 22, 3515–3532, https://doi.org/10.5194/hess-22-3515-2018, 2018.

Arnfield, A. J.: Two decades of urban climate research: a review of turbulence, exchanges of energy and water, and the urban heat island, Int. J. Climatol., 23, 1–26, https://doi.org/10.1002/joc.859, 2003.

Aubinet, M., Vesala, T., and Papale, D.: Eddy Covariance: A Practical Guide to Measurement and Data Analysis, 1st edn., edited by: Aubinet, M., Vesala, T., and Papale, D., Springer Science & Business Media, 451 pp., ISBN 978-94-007-9630-0, 2012.

Baldocchi, D. D.: How eddy covariance flux measurements have contributed to our understanding of *Global Change Biology*, Glob. Change Biol., 26, 242–260, https://doi.org/10.1111/gcb.14807, 2020.

Barlow, J. F.: Progress in observing and modelling the urban boundary layer, Urban Climate, 10, Part 2, 216–240, https://doi.org/10.1016/j.uclim.2014.03.011, 2014.

Beck, H. E., Zimmermann, N. E., McVicar, T. R., Vergopolan, N., Berg, A., and Wood, E. F.: Present and future Köppen-Geiger climate classification maps at 1-km resolution, Scientific Data, 5, 180214, https://doi.org/10.1038/sdata.2018.214, 2018.

Beringer, J., Hutley, L. B., McHugh, I., Arndt, S. K., Campbell, D., Cleugh, H. A., Cleverly, J., Resco de Dios, V., Eamus, D., Evans, B., Ewenz, C., Grace, P., Griebel, A., Haverd, V., Hinko-Najera, N., Huete, A., Isaac, P., Kanniah, K., Leuning, R., Liddell, M. J., Macfarlane, C., Meyer, W., Moore, C., Pendall, E., Phillips, A., Phillips, R. L., Prober, S. M., Restrepo-Coupe, N., Rutledge, S., Schroder, I., Silberstein, R., Southall, P., Yee, M. S., Tapper, N. J., van Gorsel, E., Vote, C., Walker, J., and Wardlaw, T.: An introduction to the Australian and New Zealand flux tower network – OzFlux, Biogeosciences, 13, 5895–5916, https://doi.org/10.5194/bg-13-5895-2016, 2016.
Best, M. J. and Grimmond, C. S. B.: Importance of initial state and atmospheric conditions for urban land surface models' performance, Urban Climate, 10, 387–406, https://doi.org/10.1016/j.uclim.2013.10.006, 2014.

Best, M. J., Abramowitz, G., Johnson, H. R., Pitman, A. J., Balsamo, G., Boone, A., Cuntz, M., Decharme, B., Dirmeyer, P. A., Dong, J., Ek, M., Guo, Z., Haverd, V., Hurk, B. J. J. van den Nearing, G. S., Pak, B., Peters-Lidard, C., Santanello, J. A., Stevens, L., and Vuichard, N.: The Plumbing of Land Surface Models: Benchmarking Model Performance, J. Hydrometeorol., 16, 1425–1442, https://doi.org/10.1175/JHM-D-14-0158.1, 2015.

Betts, A. K., Chan, D. Z., and Desjardins, R. L.: Near-Surface Bases in ERA5 Over the Canadian Prairies, Front. Environ. Sci., 7, https://doi.org/10.3389/fenvs.2019.00129, 2019.

Bjorkegren, A. B., Grimmond, C. S. B., Kotthaus, S., and Malamud, B. D.: CO₂ emission estimation in the urban environment: Measurement of the CO₂ storage term, Atmos. Environ., 122, 775–790, https://doi.org/10.1016/j.atmosenv.2015.10.012, 2015.

Boussetta, S., Balsamo, G., Beljaars, A., Panareda, A.-A., Calvet, J.-C., Jacobs, C., Hurk, B. van den, Viterbo, P., Lafont, S., Dutra, E., Jarlan, L., Balzarolo, M., Papale, D., and Werf, G. van der: Natural land carbon dioxide exchanges in the ECMWF integrated forecasting system: Implementation and offline validation, J. Geophys. Res.-Atmos., 118, 5923–5946, https://doi.org/10.1002/jgrd.50488, 2013.

Bowling, L. and Polcher, J.: The ALMA data exchange convention, https://web.lmd.jussieu.fr/~polcher/ALMA/ (last access: 3 June 2021), 2001.

Center for International Earth Science Information Network – CIESIN – Columbia University, International Food Policy Research Institute – IFPRI, The World Bank, and Centro Internacional de Agricultura Tropical – CIAT: Global Rural-Urban Mapping Project, Version 1 (GRUMPv1): Settlement Points, Revision 01, NASA Socioeconomic Data and Applications Center (SEDAC) [data set], Palisades, NY, https://doi.org/10.7927/H4BC3WG1, 2017.

Chow, W.: Eddy covariance data measured at the CAP LTER flux tower located in the west Phoenix, AZ neighborhood of Maryvale from 2011-12-16 through 2012-12-31, EDI Data Portal [data set], https://doi.org/10.6073/PASTA/FED17D67583EDA16C439216CA40B0669, 2017.

Chow, W. T. L., Volo, T. J., Vivoni, E. R., Jenerette, G. D., and Ruddell, B. L.: Seasonal dynamics of a suburban energy balance in Phoenix, Arizona, Int. J. Climatol., 34, 3863–3880, https://doi.org/10.1002/joc.3947, 2014.

Christen, A., Coops, N. C., Crawford, B. R., Kellett, R., Liss, K. N., Olchovski, I., Tooke, T. R., van der Laan, M., and Voogt, J. A.: Validation of modeled carbon-dioxide emissions from an urban neighborhood with direct eddy-covariance measurements, Atmos. Environ., 45, 6057–6069, https://doi.org/10.1016/j.atmosenv.2011.07.040, 2011.

Coutts, A. M., Beringer, J., and Tapper, N. J.: Characteristics influencing the variability of urban CO₂ fluxes in Melbourne, Australia, Atmos. Environ., 41, 51–62, https://doi.org/10.1016/j.atmosenv.2006.08.030, 2007a.

Coutts, A. M., Beringer, J., and Tapper, N. J.: Impact of Increasing Urban Density on Local Climate: Spatial and Temporal Variations in the Surface Energy Balance in Mel-

bourne, Australia, J. Appl. Meteor. Climatol., 46, 477–493, https://doi.org/10.1175/JAM2462.1, 2007b.

Crawford, B. and Christen, A.: Spatial source attribution of measured urban eddy covariance CO₂ fluxes, Theor. Appl. Climatol., 119, 733–755, https://doi.org/10.1007/s00704-014-1124-0, 2015.

Crawford, B., Grimmond, C. S. B., and Christen, A.: Five years of carbon dioxide fluxes measurements in a highly vegetated suburban area, Atmos. Environ., 45, 896–905, https://doi.org/10.1016/j.atmosenv.2010.11.017, 2011.

Cucchi, M., Weedon, G. P., Amici, A., Bellouin, N., Lange, S., Müller Schmied, H., Hersbach, H., and Buontempo, C.: WFDE5: bias-adjusted ERA5 reanalysis data for impact studies, Earth Syst. Sci. Data, 12, 2097–2120, https://doi.org/10.5194/essd-12-2097-2020, 2020.

Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bidlot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Hólm, E. V., Isaksen, L., Källberg, P., Köhler, M., Matricardi, M., McNally, A. P., Monge-Sanz, B. M., Morcrette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J.-N., and Vitart, F.: The ERA-Interim reanalysis: configuration and performance of the data assimilation system, Q. J. Roy. Meteor. Soc., 137, 553–597, https://doi.org/10.1002/qj.828, 2011.

Dou, J., Grimmond, S., Cheng, Z., Miao, S., Feng, D., and Liao, M.: Summertime surface energy balance fluxes at two Beijing sites, Int. J. Climatol., 39, 2793–2810, https://doi.org/10.1002/joc.5989, 2019.

Druken, K.: ERA5 Replicated Datasets, NCI Data Catalogue [data set], https://doi.org/10.25914/5F48874388857, 2020.

Feigenwinter, C., Vogt, R., and Christen, A.: Eddy Covariance Measurements Over Urban Areas, in: Eddy Covariance: A Practical Guide to Measurement and Data Analysis, edited by: Aubinet, M., Vesala, T., and Papale, D., Springer Netherlands, Dordrecht, 377–397, https://doi.org/10.1007/978-94-007-2351-1_16, 2012.

Ferreira, M. J., de Oliveira, A. P., and Soares, J.: Diurnal variation in stored energy flux in São Paulo city, Brazil, Urban Climate, 5, 36–51, https://doi.org/10.1016/j.uclim.2013.06.001, 2013.

Forsythe, W. C., Rykiel, E. J., Stahl, R. S., Wu, H., and Schoolfield, R. M.: A model comparison for daylength as a function of latitude and day of year, Ecol. Model., 80, 87–95, https://doi.org/10.1016/0304-3800(94)00034-F, 1995.

Fortuniak, K., Kłysik, K., and Siedlecki, M.: New measurements of the energy balance components in Łódź, in: Preprints, sixth International Conference on Urban Climate, Göteborg, Sweden, Sixth International Conference On Urban Climate, Göteborg, Sweden, 12–16 June, 2006, 64–67, http://meteo.geo.uni.lodz.pl/kf/publikacje_kf_PDF/r2006_ICUC6_p64_Fortuniak_etal.pdf (last access: 14 November 2022), 2006.

Fortuniak, K., Pawlak, W., and Siedlecki, M.: Integral Turbulence Statistics Over a Central European City Centre, Bound.-Lay. Meteorol., 146, 257–276, https://doi.org/10.1007/s10546-012-9762-1, 2013.

Goret, M., Masson, V., Schoetter, R., and Moine, M.-P.: Inclusion of CO₂ flux modelling in an urban canopy layer model and an evaluation over an old European
Best, M. J. and Grimmond, C. S. B.: Importance of initial state and atmospheric conditions for urban land surface models' performance, Urban Climate, 10, 387–406, https://doi.org/10.1016/j.uclim.2013.10.006, 2014.

Best, M. J., Abramowitz, G., Johnson, H. R., Pitman, A. J., Balsamo, G., Boone, A., Cuntz, M., Decharme, B., Dirmeyer, P. A., Dong, J., Ek, M., Guo, Z., Haverd, V., Hurk, B. J. J. van den, Nearing, G. S., Pak, B., Peters-Lidard, C., Santanello, J. A., Stevens, L., and Vuichard, N.: The Plumbing of Land Surface Models: Benchmarking Model Performance, J. Hydrometeorol., 16, 1425–1442, https://doi.org/10.1175/JHM-D-14-0158.1, 2015.

Betts, A. K., Chan, D. Z., and Desjardins, R. L.: Near-Surface Bases in ERA5 Over the Canadian Prairies, Front. Environ. Sci., 7, https://doi.org/10.3389/fenvs.2019.00129, 2019.

Bjorkegren, A. B., Grimmond, C. S. B., Kotthaus, S., and Malamud, B. D.: CO₂ emission estimation in the urban environment: Measurement of the CO₂ storage term, Atmos. Environ., 122, 775–790, https://doi.org/10.1016/j.atmosenv.2015.10.012, 2015.

Boussetta, S., Balsamo, G., Beljaars, A., Panareda, A.-A., Calvet, J.-C., Jacobs, C., Hurk, B. van den, Viterbo, P., Lafont, S., Dutra, E., Jarlan, L., Balzarolo, M., Papale, D., and Werf, G. van der: Natural land carbon dioxide exchanges in the ECMWF integrated forecasting system: Implementation and offline validation, J. Geophys. Res.-Atmos., 118, 5923–5946, https://doi.org/10.1002/jgrd.50488, 2013.

Bowling, L. and Polcher, J.: The ALMA data exchange convention, https://web.lmd.jussieu.fr/~polcher/ALMA/ (last access: 3 June 2021), 2001.

Center for International Earth Science Information Network – CIESIN – Columbia University, International Food Policy Research Institute – IFPRI, The World Bank, and Centro Internacional de Agricultura Tropical – CIAT: Global Rural-Urban Mapping Project, Version 1 (GRUMPv1): Settlement Points, Revision 01, NASA Socioeconomic Data and Applications Center (SEDAC) [data set], Palisades, NY, https://doi.org/10.7927/H4BC3WG1, 2017.

Chow, W.: Eddy covariance data measured at the CAP LTER flux tower located in the west Phoenix, AZ neighborhood of Maryvale from 2011-12-16 through 2012-12-31, EDI Data Portal [data set], https://doi.org/10.6073/PASTA/FED17D67583EDA16C439216CA40B0669, 2017.

Chow, W. T. L., Volo, T. J., Vivoni, E. R., Jenerette, G. D., and Ruddell, B. L.: Seasonal dynamics of a suburban energy balance in Phoenix, Arizona, Int. J. Climatol., 34, 3863–3880, https://doi.org/10.1002/joc.3947, 2014.

Christen, A., Coops, N. C., Crawford, B. R., Kellett, R., Liss, K. N., Olchovski, I., Tooke, T. R., van der Laan, M., and Voogt, J. A.: Validation of modeled carbon-dioxide emissions from an urban neighborhood with direct eddy-covariance measurements, Atmos. Environ., 45, 6057–6069, https://doi.org/10.1016/j.atmosenv.2011.07.040, 2011.

Coutts, A. M., Beringer, J., and Tapper, N. J.: Characteristics influencing the variability of urban CO₂ fluxes in Melbourne, Australia, Atmos. Environ., 41, 51–62, https://doi.org/10.1016/j.atmosenv.2006.08.030, 2007a.

Coutts, A. M., Beringer, J., and Tapper, N. J.: Impact of Increasing Urban Density on Local Climate: Spatial and Temporal Variations in the Surface Energy Balance in Mel-

bourne, Australia, J. Appl. Meteor. Climatol., 46, 477–493, https://doi.org/10.1175/JAM2462.1, 2007b.

Crawford, B. and Christen, A.: Spatial source attribution of measured urban eddy covariance CO₂ fluxes, Theor. Appl. Climatol., 119, 733–755, https://doi.org/10.1007/s00704-014-1124-0, 2015.

Crawford, B., Grimmond, C. S. B., and Christen, A.: Five years of carbon dioxide fluxes measurements in a highly vegetated suburban area, Atmos. Environ., 45, 896–905, https://doi.org/10.1016/j.atmosenv.2010.11.017, 2011.

Cucchi, M., Weedon, G. P., Amici, A., Bellouin, N., Lange, S., Müller Schmied, H., Hersbach, H., and Buontempo, C.: WFDE5: bias-adjusted ERA5 reanalysis data for impact studies, Earth Syst. Sci. Data, 12, 2097–2120, https://doi.org/10.5194/essd-12-2097-2020, 2020.

Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bidlot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Hólm, E. V., Isaksen, L., Källberg, P., Köhler, M., Matricardi, M., McNally, A. P., Monge-Sanz, B. M., Morcrette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J.-N., and Vitart, F.: The ERA-Interim reanalysis: configuration and performance of the data assimilation system, Q. J. Roy. Meteor. Soc., 137, 553–597, https://doi.org/10.1002/qj.828, 2011.

Dou, J., Grimmond, S., Cheng, Z., Miao, S., Feng, D., and Liao, M.: Summertime surface energy balance fluxes at two Beijing sites, Int. J. Climatol., 39, 2793–2810, https://doi.org/10.1002/joc.5989, 2019.

Druken, K.: ERA5 Replicated Datasets, NCI Data Catalogue [data set], https://doi.org/10.25914/5F48874388857, 2020.

Feigenwinter, C., Vogt, R., and Christen, A.: Eddy Covariance Measurements Over Urban Areas, in: Eddy Covariance: A Practical Guide to Measurement and Data Analysis, edited by: Aubinet, M., Vesala, T., and Papale, D., Springer Netherlands, Dordrecht, 377–397, https://doi.org/10.1007/978-94-007-2351-1_16, 2012.

Ferreira, M. J., de Oliveira, A. P., and Soares, J.: Diurnal variation in stored energy flux in São Paulo city, Brazil, Urban Climate, 5, 36–51, https://doi.org/10.1016/j.uclim.2013.06.001, 2013.

Forsythe, W. C., Rykiel, E. J., Stahl, R. S., Wu, H., and Schoolfield, R. M.: A model comparison for daylength as a function of latitude and day of year, Ecol. Model., 80, 87–95, https://doi.org/10.1016/0304-3800(94)00034-F, 1995.

Fortuniak, K., Kłysik, K., and Siedlecki, M.: New measurements of the energy balance components in Łódź, in: Preprints, sixth International Conference on Urban Climate, Göteborg, Sweden, Sixth International Conference On Urban Climate, Göteborg, Sweden, 12–16 June, 2006, 64–67, http://meteo.geo.uni.lodz.pl/kf/publikacje_kf_PDF/r2006_ICUC6_p64_Fortuniak_etal.pdf (last access: 14 November 2022), 2006.

Fortuniak, K., Pawlak, W., and Siedlecki, M.: Integral Turbulence Statistics Over a Central European City Centre, Bound.-Lay. Meteorol., 146, 257–276, https://doi.org/10.1007/s10546-012-9762-1, 2013.

Goret, M., Masson, V., Schoetter, R., and Moine, M.-P.: Inclusion of CO₂ flux modelling in an urban canopy layer model and an evaluation over an old European
city centre, Atmospheric Environment: X, 3, 100042, https://doi.org/10.1016/j.aeaoa.2019.100042, 2019.

Grimmond, C. S. B.: Progress in measuring and observing the urban atmosphere, Theor. Appl. Climatol., 84, 3–22, https://doi.org/10.1007/s00704-005-0140-5, 2006.

Grimmond, C. S. B. and Oke, T. R.: Heat Storage in Urban Areas: Local-Scale Observations and Evaluation of a Simple Model, J. Appl. Meteor., 38, 922–940, https://doi.org/10.1175/1520-0450(1999)038<0922:HSIUAL>2.0.CO;2, 1999.

Grimmond, C. S. B., Blackett, M., Best, M. J., Barlow, J., Baik, J.-J., Belcher, S. E., Bohnenstengel, S. I., Calmet, I., Chen, F., Dandou, A., Fortuniak, K., Gouvea, M. L., Hamdi, R., Hendry, M., Kawai, T., Kawamoto, Y., Kondo, H., Krayenhoff, E. S., Lee, S.-H., and Loridan, T.: The International Urban Energy Balance Models Comparison Project: First Results from Phase 1, J. Appl. Meteorol. Clim., 49, 1268–1292, https://doi.org/10.1175/2010JAMC2354.1, 2010.

Grimmond, C. S. B., Blackett, M., Best, M. J., Baik, J.-J., Belcher, S. E., Beringer, J., Bohnenstengel, S. I., Calmet, I., Chen, F., Coutts, A., Dandou, A., Fortuniak, K., Gouvea, M. L., Hamdi, R., Hendry, M., Kanda, M., Kawai, T., Kawamoto, Y., Kondo, H., Krayenhoff, E. S., Lee, S.-H., Loridan, T., Martilli, A., Masson, V., Miao, S., Oleson, K., Ooka, R., Pigeon, G., Porson, A., Ryu, Y.-H., Salamanca, F., Steeneveld, G. J., Tombrou, M., Voogt, J. A., Young, D. T., and Zhang, N.: Initial results from Phase 2 of the international urban energy balance model comparison, Int. J. Climatol., 31, 244–272, https://doi.org/10.1002/joc.2227, 2011.

Grimmond, S. and Christen, A.: Flux measurements in urban ecosystems, FluxLetter – Newsletter of Fluxnet, 5, 1–8, https://fluxnet.org/wp-content/uploads/FluxLetter_Vol5_No1.pdf (last access: 14 November 2022), 2012.

Grimmond, S. and Ward, H. C.: Urban Measurements and Their Interpretation, in: Springer Handbook of Atmospheric Measurements, edited by: Foken, T., Springer International Publishing, Cham, 1407–1437, https://doi.org/10.1007/978-3-030-52171-4_52, 2021.

Haiden, T., Sandu, I., Balsamo, G., Arduini, G., and Beljaars, A.: Addressing biases in near-surface forecasts, ECMWF Newsletter, 157, 20–25, 2018.

Hengl, T.: Clay content in % (kg/kg) at 6 standard depths (0, 10, 30, 60, 100 and 200 cm) at 250 m resolution (v0.2), Zenodo [data set], https://doi.org/10.5281/zenodo.2525663, 2018a.

Hengl, T.: Sand content in % (kg/kg) at 6 standard depths (0, 10, 30, 60, 100 and 200 cm) at 250 m resolution (v0.2), Zenodo [data set], https://doi.org/10.5281/zenodo.2525662, 2018b.

Hengl, T.: Soil bulk density (fine earth) 10 × kg/mcubic at 6 standard depths (0, 10, 30, 60, 100 and 200 cm) at 250 m resolution (v0.2), Zenodo [data set], https://doi.org/10.5281/ZENODO.2525665, 2018c.

Hersbach, H., Bell, B., Berrisford, P., Biavati, G., Horányi, A., Muñoz Sabater, J., Nicolas, J., Peubey, C., Radu, R., Rozum, I., and others: ERA5 hourly data on single levels from 1979 to present, Copernicus Climate Change Service (C3S) Climate Data Store (CDS) [data set], https://doi.org/10.24381/cds.adbb2d47, 2018.

Hersbach, H., Bell, B., Berrisford, P., Hirahara, S., Horányi, A., Muñoz-Sabater, J., Nicolas, J., Peubey, C., Radu, R., Schepers, D., Simmons, A., Soci, C., Abdalla, S., Abellan, X., Balsamo, G., Bechtold, P., Biavati, G., Bidlot, J., Bonavita, M., Chiara, G.

D., Dahlgren, P., Dee, D., Diamantakis, M., Dragani, R., Fleming, J., Forbes, R., Fuentes, M., Geer, A., Haimberger, L., Healy, S., Hogan, R. J., Hólm, E., Janisková, M., Keeley, S., Laloyaux, P., Lopez, P., Lupu, C., Radnoti, G., Rosnay, P. de, Rozum, I., Vamborg, F., Villaume, S., and Thépaut, J.-N.: The ERA5 global reanalysis, Q. J. Roy. Meteor. Soc., 146, 1999–2049, https://doi.org/10.1002/qj.3803, 2020.

Hertwig, D., Grimmond, S., Hendry, M. A., Saunders, B., Wang, Z., Jeoffrion, M., Vidale, P. L., McGuire, P. C., Bohnenstengel, S. I., Ward, H. C., and Kotthaus, S.: Urban signals in high-resolution weather and climate simulations: role of urban land-surface characterisation, Theor. Appl. Climatol., 142, 701–728, https://doi.org/10.1007/s00704-020-03294-1, 2020.

Hirano, T., Sugawara, H., Murayama, S., and Kondo, H.: Diurnal Variation of CO₂ Flux in an Urban Area of Tokyo, Sola, 11, 100–103, https://doi.org/10.2151/sola.2015-024, 2015.

Hong, J., Lee, K., and Hong, J.-W.: Observational data of Ochang and Jungnang in Korea, EAPL at Yonsei University [data set], https://doi.org/10.22647/EAPL-OC_JN2021, 2020.

Hong, J.-W., Hong, J., Chun, J., Lee, Y. H., Chang, L.-S., Lee, J.-B., Yi, K., Park, Y.-S., Byun, Y.-H., and Joo, S.: Comparative assessment of net CO₂ exchange across an urbanization gradient in Korea based on eddy covariance measurements, Carbon Balance and Management, 14, 13, https://doi.org/10.1186/s13021-019-0128-6, 2019.

Hong, S., Kim, J., Byun, Y., Hong, J., Hong, J., Lee, K., Y. Park, Lee, S., and Kim, Y.: Intra-urban variations of the CO₂ fluxes at the surface-atmosphere interface in the Seoul metropolitan area, Carbon Balance and Management, in revision, 2022.

Ishidoya, S., Sugawara, H., Terao, Y., Kaneyasu, N., Aoki, N., Tsuboi, K., and Kondo, H.: O₂: CO₂ exchange ratio for net turbulent flux observed in an urban area of Tokyo, Japan, and its application to an evaluation of anthropogenic CO₂ emissions, Atmos. Chem. Phys., 20, 5293–5308, https://doi.org/10.5194/acp-20-5293-2020, 2020.

Järvi, L., Rannik, Ü., Kokkonen, T. V., Kurppa, M., Karppinen, A., Kouznetsov, R. D., Rantala, P., Vesala, T., and Wood, C. R.: Uncertainty of eddy covariance flux measurements over an urban area based on two towers, Atmos. Meas. Tech., 11, 5421–5438, https://doi.org/10.5194/amt-11-5421-2018, 2018.

Jiang, Q., Li, W., Fan, Z., He, X., Sun, W., Chen, S., Wen, J., Gao, J., and Wang, J.: Evaluation of the ERA5 reanalysis precipitation dataset over Chinese Mainland, J. Hydrol., 595, 125660, https://doi.org/10.1016/j.jhydrol.2020.125660, 2021.

Kanda, M., Inagaki, A., Miyamoto, T., Gryschka, M., and Raasch, S.: A New Aerodynamic Parametrization for Real Urban Surfaces, Bound.-Lay. Meteorol., 148, 357–377, https://doi.org/10.1007/s10546-013-9818-x, 2013.

Karsisto, P., Fortelius, C., Demuzere, M., Grimmond, C. S. B., W., O. K., Kouznetsov, R., Masson, V., and Järvi, L.: Seasonal surface urban energy balance and wintertime stability simulated using three land-surface models in the high-latitude city Helsinki, Q. J. Roy. Meteor. Soc., 142, 401–417, https://doi.org/10.1002/qj.2659, 2016.

Kent, C. W., Grimmond, S., Barlow, J., Gatey, D., Kotthaus, S., Lindberg, F., and Halios, C. H.: Evaluation of Urban Local-Scale Aerodynamic Parameters: Implications for the Vertical Profile of Wind Speed and for Source Areas, Bound.-Lay. Meteorol., 164, 183–213, https://doi.org/10.1007/s10546-017-0248-z, 2017.
city centre, Atmospheric Environment: X, 3, 100042, https://doi.org/10.1016/j.aeaoa.2019.100042, 2019.

Grimmond, C. S. B.: Progress in measuring and observing the urban atmosphere, Theor. Appl. Climatol., 84, 3–22, https://doi.org/10.1007/s00704-005-0140-5, 2006.

Grimmond, C. S. B. and Oke, T. R.: Heat Storage in Urban Areas: Local-Scale Observations and Evaluation of a Simple Model, J. Appl. Meteor., 38, 922–940, https://doi.org/10.1175/1520-0450(1999)038<0922:HSIUAL>2.0.CO;2, 1999.

Grimmond, C. S. B., Blackett, M., Best, M. J., Barlow, J., Baik, J.-J., Belcher, S. E., Bohnenstengel, S. I., Calmet, I., Chen, F., Dandou, A., Fortuniak, K., Gouvea, M. L., Hamdi, R., Hendry, M., Kawai, T., Kawamoto, Y., Kondo, H., Krayenhoff, E. S., Lee, S.-H., and Loridan, T.: The International Urban Energy Balance Models Comparison Project: First Results from Phase 1, J. Appl. Meteorol. Clim., 49, 1268–1292, https://doi.org/10.1175/2010JAMC2354.1, 2010.

Grimmond, C. S. B., Blackett, M., Best, M. J., Baik, J.-J., Belcher, S. E., Beringer, J., Bohnenstengel, S. I., Calmet, I., Chen, F., Coutts, A., Dandou, A., Fortuniak, K., Gouvea, M. L., Hamdi, R., Hendry, M., Kanda, M., Kawai, T., Kawamoto, Y., Kondo, H., Krayenhoff, E. S., Lee, S.-H., Loridan, T., Martilli, A., Masson, V., Miao, S., Oleson, K., Ooka, R., Pigeon, G., Porson, A., Ryu, Y.-H., Salamanca, F., Steeneveld, G. J., Tombrou, M., Voogt, J. A., Young, D. T., and Zhang, N.: Initial results from Phase 2 of the international urban energy balance model comparison, Int. J. Climatol., 31, 244–272, https://doi.org/10.1002/joc.2227, 2011.

Grimmond, S. and Christen, A.: Flux measurements in urban ecosystems, FluxLetter – Newsletter of Fluxnet, 5, 1–8, https://fluxnet.org/wp-content/uploads/FluxLetter_Vol5_No1.pdf (last access: 14 November 2022), 2012.

Grimmond, S. and Ward, H. C.: Urban Measurements and Their Interpretation, in: Springer Handbook of Atmospheric Measurements, edited by: Foken, T., Springer International Publishing, Cham, 1407–1437, https://doi.org/10.1007/978-3-030-52171-4_52, 2021.

Haiden, T., Sandu, I., Balsamo, G., Arduini, G., and Beljaars, A.: Addressing biases in near-surface forecasts, ECMWF Newsletter, 157, 20–25, 2018.

Hengl, T.: Clay content in % (kg/kg) at 6 standard depths (0, 10, 30, 60, 100 and 200 cm) at 250 m resolution (v0.2), Zenodo [data set], https://doi.org/10.5281/zenodo.2525663, 2018a.

Hengl, T.: Sand content in % (kg/kg) at 6 standard depths (0, 10, 30, 60, 100 and 200 cm) at 250 m resolution (v0.2), Zenodo [data set], https://doi.org/10.5281/zenodo.2525662, 2018b.

Hengl, T.: Soil bulk density (fine earth) 10 × kg/mcubic at 6 standard depths (0, 10, 30, 60, 100 and 200 cm) at 250 m resolution (v0.2), Zenodo [data set], https://doi.org/10.5281/ZENODO.2525665, 2018c.

Hersbach, H., Bell, B., Berrisford, P., Biavati, G., Horányi, A., Muñoz Sabater, J., Nicolas, J., Peubey, C., Radu, R., Rozum, I., and others: ERA5 hourly data on single levels from 1979 to present, Copernicus Climate Change Service (C3S) Climate Data Store (CDS) [data set], https://doi.org/10.24381/cds.adbb2d47, 2018.

Hersbach, H., Bell, B., Berrisford, P., Hirahara, S., Horányi, A., Muñoz-Sabater, J., Nicolas, J., Peubey, C., Radu, R., Schepers, D., Simmons, A., Soci, C., Abdalla, S., Abellan, X., Balsamo, G., Bechtold, P., Biavati, G., Bidlot, J., Bonavita, M., Chiara, G.

D., Dahlgren, P., Dee, D., Diamantakis, M., Dragani, R., Fleming, J., Forbes, R., Fuentes, M., Geer, A., Haimberger, L., Healy, S., Hogan, R. J., Hólm, E., Janisková, M., Keeley, S., Laloyaux, P., Lopez, P., Lupu, C., Radnoti, G., Rosnay, P. de, Rozum, I., Vamborg, F., Villaume, S., and Thépaut, J.-N.: The ERA5 global reanalysis, Q. J. Roy. Meteor. Soc., 146, 1999–2049, https://doi.org/10.1002/qj.3803, 2020.

Hertwig, D., Grimmond, S., Hendry, M. A., Saunders, B., Wang, Z., Jeoffrion, M., Vidale, P. L., McGuire, P. C., Bohnenstengel, S. I., Ward, H. C., and Kotthaus, S.: Urban signals in high-resolution weather and climate simulations: role of urban land-surface characterisation, Theor. Appl. Climatol., 142, 701–728, https://doi.org/10.1007/s00704-020-03294-1, 2020.

Hirano, T., Sugawara, H., Murayama, S., and Kondo, H.: Diurnal Variation of CO₂ Flux in an Urban Area of Tokyo, Sola, 11, 100–103, https://doi.org/10.2151/sola.2015-024, 2015.

Hong, J., Lee, K., and Hong, J.-W.: Observational data of Ochang and Jungnang in Korea, EAPL at Yonsei University [data set], https://doi.org/10.22647/EAPL-OC_JN2021, 2020.

Hong, J.-W., Hong, J., Chun, J., Lee, Y. H., Chang, L.-S., Lee, J.-B., Yi, K., Park, Y.-S., Byun, Y.-H., and Joo, S.: Comparative assessment of net CO₂ exchange across an urbanization gradient in Korea based on eddy covariance measurements, Carbon Balance and Management, 14, 13, https://doi.org/10.1186/s13021-019-0128-6, 2019.

Hong, S., Kim, J., Byun, Y., Hong, J., Hong., J., Lee, K., Y. Park, Lee, S., and Kim, Y.: Intra-urban variations of the CO₂ fluxes at the surface-atmosphere interface in the Seoul metropolitan area, Carbon Balance and Management, in revision, 2022.

Ishidoya, S., Sugawara, H., Terao, Y., Kaneyasu, N., Aoki, N., Tsuboi, K., and Kondo, H.: O₂: CO₂ exchange ratio for net turbulent flux observed in an urban area of Tokyo, Japan, and its application to an evaluation of anthropogenic CO₂ emissions, Atmos. Chem. Phys., 20, 5293–5308, https://doi.org/10.5194/acp-20-5293-2020, 2020.

Järvi, L., Rannik, Ü., Kokkonen, T. V., Kurppa, M., Karppinen, A., Kouznetsov, R. D., Rantala, P., Vesala, T., and Wood, C. R.: Uncertainty of eddy covariance flux measurements over an urban area based on two towers, Atmos. Meas. Tech., 11, 5421–5438, https://doi.org/10.5194/amt-11-5421-2018, 2018.

Jiang, Q., Li, W., Fan, Z., He, X., Sun, W., Chen, S., Wen, J., Gao, J., and Wang, J.: Evaluation of the ERA5 reanalysis precipitation dataset over Chinese Mainland, J. Hydrol., 595, 125660, https://doi.org/10.1016/j.jhydrol.2020.125660, 2021.

Kanda, M., Inagaki, A., Miyamoto, T., Gryschka, M., and Raasch, S.: A New Aerodynamic Parametrization for Real Urban Surfaces, Bound.-Lay. Meteorol., 148, 357–377, https://doi.org/10.1007/s10546-013-9818-x, 2013.

Karsisto, P., Fortelius, C., Demuzere, M., Grimmond, C. S. B., W., O. K., Kouznetsov, R., Masson, V., and Järvi, L.: Seasonal surface urban energy balance and wintertime stability simulated using three land-surface models in the high-latitude city Helsinki, Q. J. Roy. Meteor. Soc., 142, 401–417, https://doi.org/10.1002/qj.2659, 2016.

Kent, C. W., Grimmond, S., Barlow, J., Gatey, D., Kotthaus, S., Lindberg, F., and Halios, C. H.: Evaluation of Urban Local-Scale Aerodynamic Parameters: Implications for the Vertical Profile of Wind Speed and for Source Areas, Bound.-Lay. Meteorol., 164, 183–213, https://doi.org/10.1007/s10546-017-0248-z, 2017.
Kokkonen, T. V., Grimmond, C. S. B., Räty, O., Ward, H. C., Christen, A., Oke, T. R., Kotthaus, S., and Järvi, L.: Sensitivity of Surface Urban Energy and Water Balance Scheme (SUEWS) to downscaling of reanalysis forcing data, Urban Climate, 23, 36–52, https://doi.org/10.1016/j.uclim.2017.05.001, 2018.

Kotthaus, S. and Grimmond, C. S. B.: Identification of Micro-scale Anthropogenic CO₂, heat and moisture sources – Processing eddy covariance fluxes for a dense urban environment, Atmos. Environ., 57, 301–316, https://doi.org/10.1016/j.atmosenv.2012.04.024, 2012.

Kotthaus, S. and Grimmond, C. S. B.: Energy exchange in a dense urban environment – Part I: Temporal variability of long-term observations in central London, Urban Climate, 10, Part 2, 261–280, https://doi.org/10.1016/j.uclim.2013.10.002, 2014a.

Kotthaus, S. and Grimmond, C. S. B.: Energy exchange in a dense urban environment – Part II: Impact of spatial heterogeneity of the surface, Urban Climate, 10, Part 2, 281–307, https://doi.org/10.1016/j.uclim.2013.10.001, 2014b.

Lipson, M.: Urban-PLUMBER: A multi-site model evaluation project for urban areas – Project Home, https://urban-plumber.github.io/, last access: 1 July 2021.

Lipson, M.: Code to process “Harmonized gap-filled datasets from 20 urban flux tower sites” for the Urban-PLUMBER projects, Zenodo [code], https://doi.org/10.5281/zenodo.7108466, 2022a.

Lipson, M.: Code to produce figures in the manuscript: “Harmonized, gap-filled dataset from 20 urban flux tower sites”, Zenodo [code], https://doi.org/10.5281/zenodo.6590941, 2022b.

Lipson, M., Grimmond, S., and Best, M.: A new multi-site evaluation project for modelling in urban areas, Urban Climate News, 15–16, http://www.urban-climate.org/wp-content/uploads/IAUC075.pdf (last access: 14 November 2022), 2020a.

Lipson, M. J., Grimmond, S., Best, M. J., Abramowitz, G., Pitman, A. J., and Ward, H. C.: Urban-PLUMBER: A new evaluation and benchmarking project for land surface models in urban areas, EGU General Assembly 2020, Online, 4–8 May 2020, EGU2020-20987, https://doi.org/10.5194/egusphere-egu2020-20987, 2020b.

Lipson, M., Grimmond, S., Best, M. J., Abramowitz, G., and Pitman, A. J.: Urban-PLUMBER—Evaluation and Benchmarking of Land Surface Models in Urban Areas, 100th American Meteorological Society Annual Meeting, Boston, USA, 2020, https://ams.confex.com/ams/2020Annual/meetingapp.cgi/Paper/364914 (last access: 14 November 2022), 2020c.

Lipson, M., Grimmond, S., and Best, M.: Protocol for observational data used in Urban-PLUMBER, Zenodo, https://doi.org/10.5281/zenodo.5518278, 2021.

Lipson, M., Grimmond, S., Best, M., Chow, W., Christen, A., Chrysoulakis, N., Coutts, A., Crawford, B., Earl, S., Evans, J., Fortuniak, K., Heusinkveld, B. G., Hong, J.-W., Hong, J., Järvi, L., Jo, S., Kim, Y.-H., Kotthaus, S., Lee, K., Masson, V., McFadden, J. P., Michels, O., Pawlak, W., Roth, M., Sugawara, H., Tapper, N., Velasco, E., and Ward, H. C.: Data for “Harmonized gap-filled dataset from 20 urban flux tower sites” for the Urban-PLUMBER project, Zenodo [data set], https://doi.org/10.5281/zenodo.7104984, 2022.

Macdonald, R. W., Griffiths, R. F., and Hall, D. J.: An improved method for the estimation of surface roughness of obstacle arrays, Atmos. Environ., 32, 1857–1864, https://doi.org/10.1016/S1352-2310(97)00403-2, 1998.

Martens, B., Schumacher, D. L., Wouters, H., Muñoz-Sabater, J., Verhoest, N. E. C., and Miralles, D. G.: Evaluating the land-surface energy partitioning in ERA5, Geosci. Model Dev., 13, 4159–4181, https://doi.org/10.5194/gmd-13-4159-2020, 2020.

Masson, V., Gomes, L., Pigeon, G., Liousse, C., Pont, V., Lagouarde, J.-P., Voogt, J., Salmond, J., Oke, T. R., Hidalgo, J., Legain, D., Garrouste, O., Lac, C., Connan, O., Briottet, X., Lachérade, S., and Tulet, P.: The Canopy and Aerosol Particles Interactions in TOulouse Urban Layer (CAPI-TOUL) experiment, Meteorol. Atmos. Phys., 102, 135–157, https://doi.org/10.1007/s00703-008-0289-4, 2008.

Masson, V., Heldens, W., Bocher, E., Bonhomme, M., Bucher, B., Burmeister, C., de Munck, C., Esch, T., Hidalgo, J., Kanani-Sühring, F., Kwok, Y.-T., Lemonsu, A., Lévy, J.-P., Maronga, B., Pavlik, D., Petit, G., See, L., Schoetter, R., Tornay, N., Votsis, A., and Zeidler, J.: City-descriptive input data for urban climate models: Model requirements, data sources and challenges, Urban Climate, 31, 100536, https://doi.org/10.1016/j.uclim.2019.100536, 2020.

McNorton, J. R., Arduini, G., Bousserez, N., Agustí-Panareda, A., Balsamo, G., Boussetta, S., Choulga, M., Hadade, I., and Hogan, R. J.: An Urban Scheme for the ECMWF Integrated Forecasting System: Single-Column and Global Offline Application, J. Adv. Model. Earth Sy., 13, e2020MS002375, https://doi.org/10.1029/2020MS002375, 2021.

Menard, C. B., Essery, R., Krinner, G., Arduini, G., Bartlett, P., Boone, A., Brutel-Vuilmet, C., Burke, E., Cuntz, M., Dai, Y., Decharme, B., Dutra, E., Fang, X., Fierz, C., Gusev, Y., Hagemann, S., Haverd, V., Kim, H., Lafaysse, M., Marke, T., Na sonova, O., Nitta, T., Niwano, M., Pomeroy, J., Schädler, G., Se menov, V. A., Smirnova, T., Strasser, U., Swenson, S., Turkov, D., Wever, N., and Yuan, H.: Scientific and Human Errors in a Snow Model Intercomparison, B. Am. Meteorol. Soc., 102, E61– E79, https://doi.org/10.1175/BAMS-D-19-0329.1, 2021.

Menne, M. J., Durre, I., Vose, R. S., Gleason, B. E., and Houston, T. G.: An Overview of the Global Historical Climatology Network-Daily Database, J. Atmos. Ocean. Tech., 29, 897–910, https://doi.org/10.1175/JTECH-D-11-00103.1, 2012.

Menzer, O. and McFadden, J. P.: Statistical partitioning of a three-year time series of direct urban net CO₂ flux measurements into biogenic and anthropogenic components, Atmos. Environ., 170, 319–333, https://doi.org/10.1016/j.atmosenv.2017.09.049, 2017.

New, M., Hulme, M., and Jones, P.: Representing Twentieth-Century Space-Time Climate Variability. Part I: Development of a 1961–90 Mean Monthly Terrestrial Climatology, J. Climate, 12, 829–856, https://doi.org/10.1175/1520-0442(1999)012<0829:RTCSTC>2.0.CO;2, 1999.

Nogueira, M.: Inter-comparison of ERA-5, ERA-interim and GPCP rainfall over the last 40 years: Process-based analysis of systematic and random differences, J. Hydrol., 583, 124632, https://doi.org/10.1016/j.jhydrol.2020.124632, 2020.

Nordbo, A., Järvi, L., Haapanala, S., Moilanen, J., and Vesala, T.: Intra-City Variation in Urban Morphology and Turbulence Structure in Helsinki, Finland, Bound.-Lay. Meteorol., 146, 469–496, https://doi.org/10.1007/s10546-012-9773-y, 2013.

Novick, K. A., Biederman, J. A., Desai, A. R., Litvak, M. E., Moore, D. J. P., Scott, R. L., and Torn, M. S.: The AmeriFlux network: A coalition of the willing, Agr. Forest Meteorol., 249, 444–456, https://doi.org/10.1016/j.agrformet.2017.10.009, 2018.
Kokkonen, T. V., Grimmond, C. S. B., Räty, O., Ward, H. C., Christen, A., Oke, T. R., Kotthaus, S., and Järvi, L.: Sensitivity of Surface Urban Energy and Water Balance Scheme (SUEWS) to downscaling of reanalysis forcing data, Urban Climate, 23, 36–52, https://doi.org/10.1016/j.uclim.2017.05.001, 2018.

Kotthaus, S. and Grimmond, C. S. B.: Identification of Micro-scale Anthropogenic CO₂, heat and moisture sources – Processing eddy covariance fluxes for a dense urban environment, Atmos. Environ., 57, 301–316, https://doi.org/10.1016/j.atmosenv.2012.04.024, 2012.

Kotthaus, S. and Grimmond, C. S. B.: Energy exchange in a dense urban environment – Part I: Temporal variability of long-term observations in central London, Urban Climate, 10, Part 2, 261–280, https://doi.org/10.1016/j.uclim.2013.10.002, 2014a.

Kotthaus, S. and Grimmond, C. S. B.: Energy exchange in a dense urban environment – Part II: Impact of spatial heterogeneity of the surface, Urban Climate, 10, Part 2, 281–307, https://doi.org/10.1016/j.uclim.2013.10.001, 2014b.

Lipson, M.: Urban-PLUMBER: A multi-site model evaluation project for urban areas – Project Home, https://urban-plumber.github.io/, last access: 1 July 2021.

Lipson, M.: Code to process “Harmonized gap-filled datasets from 20 urban flux tower sites” for the Urban-PLUMBER projects, Zenodo [code], https://doi.org/10.5281/zenodo.7108466, 2022a.

Lipson, M.: Code to produce figures in the manuscript: “Harmonized, gap-filled dataset from 20 urban flux tower sites”, Zenodo [code], https://doi.org/10.5281/zenodo.6590941, 2022b.

Lipson, M., Grimmond, S., and Best, M.: A new multi-site evaluation project for modelling in urban areas, Urban Climate News, 15–16, http://www.urban-climate.org/wp-content/uploads/IAUC075.pdf (last access: 14 November 2022), 2020a.

Lipson, M. J., Grimmond, S., Best, M. J., Abramowitz, G., Pitman, A. J., and Ward, H. C.: Urban-PLUMBER: A new evaluation and benchmarking project for land surface models in urban areas, EGU General Assembly 2020, Online, 4–8 May 2020, EGU2020-20987, https://doi.org/10.5194/egusphere-egu2020-20987, 2020b.

Lipson, M., Grimmond, S., Best, M. J., Abramowitz, G., and Pitman, A. J.: Urban-PLUMBER—Evaluation and Benchmarking of Land Surface Models in Urban Areas, 100th American Meteorological Society Annual Meeting, Boston, USA, 2020, https://ams.confex.com/ams/2020Annual/meetingapp.cgi/Paper/364914 (last access: 14 November 2022), 2020c.

Lipson, M., Grimmond, S., and Best, M.: Protocol for observational data used in Urban-PLUMBER, Zenodo, https://doi.org/10.5281/zenodo.5518278, 2021.

Lipson, M., Grimmond, S., Best, M., Chow, W., Christen, A., Chrysoulakis, N., Coutts, A., Crawford, B., Earl, S., Evans, J., Fortuniak, K., Heusinkveld, B. G., Hong, J.-W., Hong, J., Järvi, L., Jo, S., Kim, Y.-H., Kotthaus, S., Lee, K., Masson, V., McFadden, J. P., Michels, O., Pawlak, W., Roth, M., Sugawara, H., Tapper, N., Velasco, E., and Ward, H. C.: Data for “Harmonized gap-filled dataset from 20 urban flux tower sites” for the Urban-PLUMBER project, Zenodo [data set], https://doi.org/10.5281/zenodo.7104984, 2022.

Macdonald, R. W., Griffiths, R. F., and Hall, D. J.: An improved method for the estimation of surface roughness of obstacle arrays, Atmos. Environ., 32, 1857–1864, https://doi.org/10.1016/S1352-2310(97)00403-2, 1998.

Martens, B., Schumacher, D. L., Wouters, H., Muñoz-Sabater, J., Verhoest, N. E. C., and Miralles, D. G.: Evaluating the land-surface energy partitioning in ERA5, Geosci. Model Dev., 13, 4159–4181, https://doi.org/10.5194/gmd-13-4159-2020, 2020.

Masson, V., Gomes, L., Pigeon, G., Liousse, C., Pont, V., Lagouarde, J.-P., Voogt, J., Salmond, J., Oke, T. R., Hidalgo, J., Legain, D., Garrouste, O., Lac, C., Connan, O., Briottet, X., Lachérade, S., and Tulet, P.: The Canopy and Aerosol Particles Interactions in TOulouse Urban Layer (CAPI-TOUL) experiment, Meteorol. Atmos. Phys., 102, 135–157, https://doi.org/10.1007/s00703-008-0289-4, 2008.

Masson, V., Heldens, W., Bocher, E., Bonhomme, M., Bucher, B., Burmeister, C., de Munck, C., Esch, T., Hidalgo, J., Kanani-Sühring, F., Kwok, Y.-T., Lemonsu, A., Lévy, J.-P., Maronga, B., Pavlik, D., Petit, G., See, L., Schoetter, R., Tornay, N., Votsis, A., and Zeidler, J.: City-descriptive input data for urban climate models: Model requirements, data sources and challenges, Urban Climate, 31, 100536, https://doi.org/10.1016/j.uclim.2019.100536, 2020.

McNorton, J. R., Arduini, G., Bousserez, N., Agustí-Panareda, A., Balsamo, G., Boussetta, S., Choulga, M., Hadade, I., and Hogan, R. J.: An Urban Scheme for the ECMWF Integrated Forecasting System: Single-Column and Global Offline Application, J. Adv. Model. Earth Sy., 13, e2020MS002375, https://doi.org/10.1029/2020MS002375, 2021.

Menard, C. B., Essery, R., Krinner, G., Arduini, G., Bartlett, P., Boone, A., Brutel-Vuilmet, C., Burke, E., Cuntz, M., Dai, Y., Decharme, B., Dutra, E., Fang, X., Fierz, C., Gusev, Y., Hagemann, S., Haverd, V., Kim, H., Lafaysse, M., Marke, T., Na sonova, O., Nitta, T., Niwano, M., Pomeroy, J., Schädler, G., Se menov, V. A., Smirnova, T., Strasser, U., Swenson, S., Turkov, D., Wever, N., and Yuan, H.: Scientific and Human Errors in a Snow Model Intercomparison, B. Am. Meteorol. Soc., 102, E61– E79, https://doi.org/10.1175/BAMS-D-19-0329.1, 2021.

Menne, M. J., Durre, I., Vose, R. S., Gleason, B. E., and Houston, T. G.: An Overview of the Global Historical Climatology Network-Daily Database, J. Atmos. Ocean. Tech., 29, 897–910, https://doi.org/10.1175/JTECH-D-11-00103.1, 2012.

Menzer, O. and McFadden, J. P.: Statistical partitioning of a three-year time series of direct urban net CO₂ flux measurements into biogenic and anthropogenic components, Atmos. Environ., 170, 319–333, https://doi.org/10.1016/j.atmosenv.2017.09.049, 2017.

New, M., Hulme, M., and Jones, P.: Representing Twentieth-Century Space-Time Climate Variability. Part I: Development of a 1961–90 Mean Monthly Terrestrial Climatology, J. Climate, 12, 829–856, https://doi.org/10.1175/1520-0442(1999)012<0829:RTCSTC>2.0.CO;2, 1999.

Nogueira, M.: Inter-comparison of ERA-5, ERA-interim and GPCP rainfall over the last 40 years: Process-based analysis of systematic and random differences, J. Hydrol., 583, 124632, https://doi.org/10.1016/j.jhydrol.2020.124632, 2020.

Nordbo, A., Järvi, L., Haapanala, S., Moilanen, J., and Vesala, T.: Intra-City Variation in Urban Morphology and Turbulence Structure in Helsinki, Finland, Bound.-Lay. Meteorol., 146, 469–496, https://doi.org/10.1007/s10546-012-9773-y, 2013.

Novick, K. A., Biederman, J. A., Desai, A. R., Litvak, M. E., Moore, D. J. P., Scott, R. L., and Torn, M. S.: The AmeriFlux network: A coalition of the willing, Agr. Forest Meteorol., 249, 444–456, https://doi.org/10.1016/j.agrformet.2017.10.009, 2018.