Urban Land Surface Models (ULSMs) simulate energy and water exchanges between the urban surface and atmosphere. When part of numerical weather prediction, ULSMs provide a lower boundary for the atmosphere and improve the applicability of model results in the urban environment compared with non-urban land surface models. However, earlier systematic ULSM comparison projects assessed the energy balance but ignored the water balance which is coupled to the energy balance. Here, we analyze the water balance representation in 19 ULSMs participating in the Urban-PLUMBER project using results for 20 sites spread across a range of climates and urban form characteristics. As observations for most water fluxes are unavailable, we examine the water balance closure, flux timing, and magnitude with a score derived from seven indicators. We find that the water budget is only closed in 57% of the model-site combinations assuming closure when annual total incoming fluxes (precipitation and irrigation) fluxes are within 3% of the outgoing (all other) fluxes. Results show the timing is better captured than magnitude. No ULSM has passed all good water balance indicators for any site. Our results indicate models could be improved by explicitly verifying water balance closure and revising runoff parameterizations. By expanding ULSM evaluation to the water balance and related to latent heat flux performance, we demonstrate the benefits of evaluating processes with direct feedback mechanisms to the processes of interest.

Evapotranspiration (ET) is a key process in the hydrological cycle that can help mitigate urban heat. ET depends on the surface cover, as the surface affects the partitioning of precipitation between runoff and evapotranspiration. In urban neighborhoods, this surface cover is highly heterogeneous. The resulting neighborhood-scale ET is observed with eddy-covariance systems. However, these observations represent the signal from wind- and stability-dependent footprints resulting in a continuously changing surface cover composition. This continuous change prevents quantitative analysis of the separate types. Here, we disentangle this neighborhood-scale ET at two urban sites in Berlin attributing the ET dynamics to the four major surface cover types in the footprint: impervious surfaces, low vegetation, high vegetation, and open water. Starting from the surface, we reconstruct ET based on patch-scale observations and conceptual models. Alternatively, we start with the eddy-covariance observations and attribute ET to the surface cover types solving a system of equations for four eddy-covariance systems with different footprints. Although starting at the surface yields more robust results, both approaches indicate that vegetation is responsible for more ET than proportional to its surface fraction, and evaporation from impervious surfaces although less cannot be neglected. We confirm the intuitive relation between ET and the surface cover fractions based on a wide range of surface compositions.

Water storage plays an important role in mitigating heat and flooding in urban areas. Assessment of the water storage capacity of cities remains challenging due to the inherent heterogeneity of the urban surface. Traditionally, effective storage has been estimated from runoff. Here, we present a novel approach to estimate effective water storage capacity from recession rates of observed evaporation during precipitation-free periods. We test this approach for cities at neighborhood scale with eddy-covariance based latent heat flux observations from fourteen contrasting sites with different local climate zones, vegetation cover and characteristics, and climates. Based on analysis of 583 drydowns, we find storage capacities to vary between 1.3-28.4 mm, corresponding to e-folding timescales of 1.8-20.1 days. This makes the storage capacity at least one order of magnitude smaller than the observed values for natural ecosystems, reflecting an evaporation regime characterised by extreme water limitation.