A simple step heating approach for wall surface temperature estimation in the SOlar and LongWave Environmental Irradiance Geometry (SOLWEIG) model
Abstract. The urban climate is highly influenced by its building geometry, material characteristics, street orientation and high fraction of impermeable surfaces. All of these influence the microclimate and the resulting outdoor thermal comfort. Mean radiant temperature (Tmrt) is often used as an estimator for heat exposure as it is one of the most important variables governing outdoor human thermal comfort on clear, calm and warm days. The highest values of Tmrt are commonly found in front of sunlit facades where a human is exposed to high levels of direct and reflected shortwave radiation from the sun, as well as high levels of longwave radiation emitted from surrounding sunlit walls. As a consequence, outdoor thermal comfort modelling requires accurate simulation of wall surface temperatures (Ts).
The aim of this study is to present a step heating approach for calculating wall Ts in the SOlar and LongWave Environmental Irradiance Geometry model (SOLWEIG) and quantifying how it influences Tmrt. This method requires information on material characteristics, i.e. specific heat capacity, density, thermal conductivity, albedo and thickness of the outer layer of the wall, as well as radiation balance at the wall surface, and ambient air temperature. Simulated Ts is compared to observed Ts of two white walls (albedo = 0.5) in Gothenburg, Sweden; one wooden wall and one plaster brick wall. The simulations show high agreement with the 15,394 observations, with R2 = 0.93 and RMSE = 2.09 °C for the wooden wall and R2 = 0.94 and RMSE = 1.94 °C for the plaster brick wall. For the walls presented here, this new parameterization scheme results in differences in Tmrt of up to 2.5 °C compared to the previous version of SOLWEIG.
With this new approach SOLWEIG can be used to evaluate the effect of building materials on outdoor thermal comfort. The speed and accuracy of this approach suggests that it also could be applied in other areas where Ts of walls are important, for example building energy models and urban energy balance models.
General Comments
In this paper, a new parameterization scheme for the SOLWEIG model is introduced and tested for a dataset from two walls with similar aspect and albedo, but different thermal properties. The evaluation shows good performance for the wall temperatures and a demonstration is made of the impact on SOLWEIG-predicted MRT. Overall, I think the work provides the basis for a useful contribution and the introduction of the scheme will broaden the capabilities of the SOLWEIG model. The step heating approach chosen seems to be based in part on calculation efficiency, although tests are not made to demonstrate its performance in this respect and I thought there was scope to include some sensitivity testing of the choice of thermal parameters. More generally, it was not clear to me if there are any assumptions implicit in the step heating approach that might impact its ability to represent environmental temperature time series. The basis of that approach seems to be in the laboratory testing of thermal properties of materials where a step change in heat input can be applied, but this may be different from the type of input from environmental radiative forcings. I recall here the comparison of PDE and Force Restore approaches for surface temperature modeling in Johnson et al. 1991 Boundary Layer Meteorology 56, 275-294 that dealt with some similar issues.
Specific Comments
Line 23 Does the choice of walls with albedo = 0.5 make the test somewhat conservative for Ts, given Ts would be expected to be lower in this case (compared to a wall with lower albedo?) Perhaps this choice is related to expected MRT?
Line 55-57 The Fourier's law of diffusion is just for the conduction of heat - so a surface energy balance approach could potentially replace Fourier's law with an alternative -e.g. a force restore approach could be used to resolve conduction while still using a surface energy balance approach.
Line 62 And possibly also to the season or time period of that calibration.
Line 67 Thermal effussivity is equivalent to thermal admittance (see e.g. Oke et al. 2017 Urban Climates) - which may be familiar to readers with an urban climate background so it could be useful to point out this equivalence.
Line 76-79. I wonder if a bit more context is needed - e.g. SOLWEIG is not quite a full energy balance model - and success at modeling Ts will also be related in part to the performance of the convective and conductive parameterizations (as well as radiation).
Line 84 Are there any issues given that the environmental forcing is different from a Dirac heat pulse as assumed in the reference?
Equation 1: Equation 1) has surface temperature on both sides of the equation, since outgoing longwave that is part of omega is determined by Ts – this may be worth noting.
Line 102 The matching is probably best for walls a few hours after sunset (e.g. to allow west-facing walls that may have been directly irradiated to cool closer to Tair)
Equation 13 I guess fsun could also depend on the canyon geometry? and some of the empiricism is related to expected obstructions within a canyon and/or microscale facet structure such as awnings that creates self shading?
Equation 15 I'm not sure why the "fraction" terminology is used when the other values are view factors? Essentially you are calculating the view factor of the opposite wall of the street canyon?
Equation 16 If psi_g is SVF at ground level then psi_v is SVF at...? Is this for ground-level vegetation (grass?) or tree canopy? or... This is not to clear to me, and for tree canopy things get more complicated? Maybe some more description would help the reader here and/or given an example for vegetation.
Line 181 I was initially surprised at the very small distance of the IRR (infrared radiometer) from the surface, however Figure 1 is helpful for putting this small distance in context - so perhaps adding some words about the need for a simple lightweight mount attached to the surface of interest might be useful.
As a more general comment: The IFOV of the wall surface for the IRR would be pretty small with such a small distance from the surface - and presumably the angle of observation is set to try and avoid any shadow issues within the IFOV? I guess this helps with not having a long period of sunlit/shadow transition on the wall if that is an issue for this position; on the other hand a larger distance would get a larger spatial sample and reduce any impacts of the sensor on the surface itself. The wooden wall also shows some microscale structure from the edge of the boards – possibly these are on a slight angle and have some microshadows beneath them? Do you think this has any impact on the model vs observed relation? Is the angle of view of the IRR set to try to avoid any such effects?
Line 182 Have you tried any sensitivity testing on the choice of these thermal parameters?
Line 185 I guess this is average shortwave radiation incident on the plane of the wall (not horizontal incident shortwave at the top of the canyon?) Probably should make this explicit.
Line 187 Re emissivity setting of the IRR - it may be worth checking on what assumptions it is making regarding the incident longwave radiation on the target , which may not always be true in urban environments.
Line 305-307 Some older models might also be included? e.g. the wall temperature comparison done for the TUF3d model (Krayenhoff & Voogt 2007) and Henon et al. 2012 with SOLENE? And there may be others too.
Line 325 I might describe Ts as 'influencing' the sensible heat flux, since there are probably larger scale processes occurring within the street canyon that may influence the extent and direction of sensible heat flux.
Line 327-328 I agree and the success of relatively simple models for estimating 3d urban temperature in applications of urban thermal anisotropy also point to the dominance of radiation as a control on surface temperature.
Line 368 We may know the range of parameters but isn't a challenge knowing how to assign such parameters to the variety of buildings in an urban area and not necessarily knowing the cross-sectional structure of the walls for all buildings? Presumably this leads to a level of uncertainty that would be greater than for a test case with a building where more precise information can be determined.
Wording/typographical suggestions
Line 20 suggest “and quantify how it influences”
Line 41 Odd phrasing? The spatial differences in urban air temperature are typically small in most environments by day? (Key element here is daytime - since at night time the spatial differences in Tair are larger and better linked to surface characteristics).
Line 105 Suggest "The omega term (in eq. 1) ..." (just to avoid starting the sentence with the symbol)
Line 108 Suggest: "The absorbed shortwave radiation" (the equations are representing the absorbed, not just the incident shortwave radiation)
Line 121 Suggest lower case w for where
Line 134 I think it would be useful to be explicit on what Lup and Ldown are.
Line 135 Suggest: "of sunlit non-sky"
Line 147 Here an upper case version of the symbol is used compared to equations 5 and 10; suggest it be changed to match.
Lines 171-174 Perhaps this sentence should be moved to section 2.5?
Line 178 Specify the wall aspect explicitly here? Looks to be south-east based on Figure 1.
Line 222 suggest “and underestimates by 2-3°C in nighttime.”
Line 223 I'm not sure i understand the wording - do you mean they overestimate the surface temperature during times of intense solar radiation?
Line 233 I would move this sentence to the previous sentence that begins on line 230.
Line 239 suggest “are subject to partly cloudy weather”
Fig 4 The label of the regression line omits the + sign of the equation; a reader would presumably infer this but might be useful to give all the signs given it is written as an equation.
Line 316-322 Move this to the methods?