Revising the definition of anthropogenic heat flux from buildings: role of human activities and building storage heat flux

. Buildings are a major source of anthropogenic heat emissions, impacting energy use and human health in cities. The difference in magnitude and time lag between building energy consumption and building anthropogenic heat emission is poorly quantiﬁed. Energy consumption ( Q EC ) is a widely used proxy for the anthropogenic heat ﬂux from buildings ( Q F , B ). Here we revisit the latter’s deﬁnition. If Q F , B is the heat emission to the outdoor environment from human activities within buildings, we can derive it from the changes in energy balance ﬂuxes between occupied and unoccupied buildings. Our derivation shows that the difference be-tween Q EC and Q F , B is attributable to a change in the storage heat ﬂux induced by human activities ( (cid:49)S o-uo ) (i.e. Q F , B = Q EC − (cid:49)S o-uo ). Using building energy simulations (EnergyPlus) we calculate the energy balance ﬂuxes for a simpliﬁed isolated building (obtaining Q F , B , Q EC , (cid:49)S o-uo ) with different occupancy states. The non-negligible differences in diurnal patterns between Q F , B and Q EC are caused by thermal storage (e.g. hourly Q F , B to Q EC ratios vary between − 2.72 and 5.13 within a year in Beijing, China). Negative Q F , B can occur as human activities can reduce heat emission from a building but this is associated with a large storage heat ﬂux. Building operations (e.g. opening windows, use of space heating and cooling system) modify the Q F , B by affecting not only Q EC but also the (cid:49)S o-uo diurnal proﬁle. Air temperature and solar radiation are critical meteorological factors explaining day-to-day variability of Q F , B . Our new approach could be used to provide data for future parameterisations of both anthropogenic heat ﬂux and storage heat ﬂuxes from buildings. It is evident that storage heat ﬂuxes in cities could also be impacted by occupant behaviour


Introduction
Human activities that influence energy exchanges are critical to a wide variety of disciplines (e.g.meteorology, building design, geography, climatology, hydrology, engineering).As disciplines often have interests in different scales, purposes and/or boundary conditions, the terminology and acceptable assumptions differ.However, disciplines may provide data to each other or help improve assumptions used.In this study we are concerned with the interface between meteorology, climatology and building design in urban areas.
To model the weather and climate in urban areas, an important additional source of energy to the environment is the anthropogenic heat flux (Q F ).This is defined as the heat converted from consumption of biological, chemical and electrical energy and released to the atmosphere due to human activities (Oke et al., 2017).The Q F has three major sources, including metabolic (people and animals) activities (Q F,M ), transport (Q F,T ) and buildings (Q F,B ) (Grimmond, 1992).It can be large relative to incoming solar radiation in summer (e.g.43 % in an area of Beijing, Nie et al. (2014)) and increases air temperature in cities (e.g.Ichinose et al., 1999;Fan and Sailor, 2005), subsequently contributing to higher cooling demand for buildings (Santamouris et al., 2001;Takane et al., 2019).In winter Q F can contribute to the intensity of the urban heat island (Biggart et al., 2021).
Published by Copernicus Publications on behalf of the European Geosciences Union.
Y. Liu et al.: Revising the definition of anthropogenic heat flux from buildings Not all heat generated within the building volume is directly ejected into the outdoor environment immediately but subject to change in magnitude and time lag.For example, the heat generated from human activities inside buildings is released initially indoors (via heating or cooling application), then transported through the building fabric by conduction, allowing it to be transported into the atmosphere by turbulent sensible heat flux and outgoing longwave radiation.In this process the net storage heat flux ( Q S ) of a building is modified since the building fabric temperature is changed by absorbing more heat from the internal heat generation.
In urban areas, Q S is the net uptake or release of energy from urban volume.This term is an important determinant of urban climate and is regarded as a key process in the genesis of urban heat island (Goward, 1981).The change in building Q S is modified when heat is released by human activities but the timing of the external emissions are impacted by the building fabric characteristics and the conduction process.As prior studies often used energy consumption (Q EC ) as a proxy for Q F,B , derived from inventory-related approaches (e.g.Sailor and Lu, 2004;Iamarino et al., 2012) and building energy modelling (e.g.Heiple and Sailor, 2008;Nie et al., 2014), the impact on Q S is not addressed.To qualify the "real" Q F,B and change of Q S , we revisit the definition of Q F,B and attempt to understand how human activities affect the energy balance fluxes of buildings.
If Q F,B is the heat released from buildings into the atmosphere as a result of human activities inside the building (including human metabolism), when the building is completely unoccupied (e.g.no operational appliances, no people: such as "ghost cities" in China (Shepard, 2015) or vacant in Dublin, (Kelly and Scott, 2018)), then Q F,B is zero.However, heat released from the unoccupied building is non-zero as there is still heat exchange between the building and the ambient environment (see Eq. 1 and 2), as occurs in other environments with large mass, such as forests (e.g.Oliphant et al., 2004), and rocks (e.g.Wang et al., 2018).The Q F,B differs from building heat emission (BHE) (e.g.Hong et al., 2020;Ferrando et al., 2021) as the latter is the total heat flux released from buildings to the ambient air ) not due to human activities alone.Shortwave and longwave radiation can enter the unoccupied internal building space through windows and conduction through walls.It modifies the heat stored within the building volume and the temperature of the building envelope and indoor air, subsequently influencing the emission of heat via sensible heat flux, outgoing longwave radiation and air exchange.However, this energy leaving the unoccupied building is not anthropogenic heat flux.For an occupied building, the internal heat gain arises from 1. the equivalent sources and sinks as the unoccupied buildings; but also 2. the energy linked to the indoor human activities (metabolism, powered appliances and energy inputs to heating or cooling).
These will modify each of the energy balance fluxes.Some of this additional energy is transported out of buildings through indoor-outdoor ventilation exchange and/or heating ventilation and air conditioning (HVAC) systems, immediately contributes to Q F,B , while some is stored in the building fabric, and later is released outdoors through various pathways (convection, radiation, conduction) to become Q F,B with a time lag.Here, we derive Q F,B by looking at the difference of heat fluxes between occupied and unoccupied buildings.
If the energy balance for the building system (including the indoor air and building envelope) for an unoccupied dry building (assuming latent heat is not important in this case) is The radiation balance for an isolated unoccupied (uo) building can be expressed as: where Q * is the net all-wave radiation, K is the shortwave radiation incoming (↓) and outgoing (↑) to the external surfaces.The longwave (L) radiation exchanges depend on the view factors (F ) between the building of interest (boi), the surrounding facets of other surfaces or buildings (other b) and the sky: In Eq. ( 1), Q H is the turbulent sensible heat flux (convection) from external surfaces to the external ambient air and Q BAE is the net energy exchange from the buildings through air exchange (e.g.ventilation).When the building is sealed Q BAE is 0 W m −2 , otherwise (e.g.open windows, cracks) it can be a source or sink of energy (environment ← building, or inverse).The Q S is the net storage heat flux of the building volume (i.e.fabric, contents, including the air).The left-hand side (LHS) of Eq. ( 1) is the inputs or source of energy to the building, whereas the right-hand side (RHS) is the sink or energy dissipation outputs.With no human activities within the building, the internal heat generation from human and infrastructure activities is zero.When the building is occupied (o) (e.g.appliances operating, people present), additional terms are needed in Eq. (1) to account for the supply of energy into the building for these activities and the release of energy: The two additional sources of energy Here Q H only represents the convection heat transfer at building external surface (i.e.wall, roof and windows).Both Q Waste , and Q BAE will be incorporated into the turbulent sensible heat flux by the time they reach the inertial sublayer (ISL) or constant flux layer (CFL).Hence, sensors (e.g.eddy covariance or large aperture scintillometry) located in the ISL would observe this as Q H .The separation of these three terms is necessary for a better understanding of how human activities (e.g.open/closed windows, HVAC operation) influence each heat flux.Urban canopy parameterisation (UCP) can use this information about the separate sources and their roles in the urban energy balance to account for the modified fluxes by the time they reach the ISL.Additionally, it is clearer for multi-layer UCP where the energy should enter vertically.
To determine the impact of the occupancy (i.e.not just the physical building form) we can consider the difference between Eqs. ( 5) and (1).If the radiation balance for the occupied case is we assume that the incoming and outgoing shortwave radiation remain unchanged because the reflectivity, transmissivity and absorptivity do not change by occupancy activities then: The incoming longwave radiation is dependent on the surroundings, which are independent of the building state, so Thus, the difference in radiative fluxes between occupied and unoccupied buildings ( L ↑o−uo ) is Similarly, the difference of the heat transfer through air exchange is With the additional terms in Eq. ( 5) and the air exchange rate differences from the activities within the buildings, gives As the change in surface temperature influences the sensible heat fluxes and storage heat fluxes: By combining the Eqs.( 1) and ( 5), we obtain where the LHS accounts for the net available energy as a result of human activities in indoor environments and the RHS shows that these impact the longwave radiation, turbulent sensible and storage heat fluxes (in this dry case).With rearrangement: The additional energy generation associated with human activities to the whole building system (LHS) is apparent, as traditionally defined as Q F,B previously (Heiple and Sailor, 2008).Here, because the heat release from human metabolism indoors is considerably smaller than other sources, for simplicity of analysis we assume metabolic heat is also part of energy consumption ( ).In addition, some of additional energy is associated with the extra gain or release of stored heat within the building volume ( S o−uo ).The rest is the heat released to the outdoor environment from the building due to human activities, which is the Q F,B based on its definition: Eq. ( 14) demonstrates that the Q F,B is the relative heat emission at the exterior building boundary between unoccupied and occupied buildings through longwave radiation, convection, air exchange and waste heat from any mechanical heating/cooling system.The source of Q F,B within the building volume gives (by combining Eqs. ( 13) and ( 14) The sources of Q F,B are from both energy consumption (Q EC ) and differences of storage heat flux ( S o-uo ) between unoccupied and occupied buildings (Q F,B ) in this https://doi.org/10.5194/acp-22-4721-2022Atmos.Chem.Phys., 22, 4721-4735, 2022 study includes part of Q F,M from human metabolism).In most prior studies, the second term of Eq. ( 15) is ignored.Although the storage heat flux over a year should tend to zero, over short periods (e.g.sub-daily) S o-uo is not zero causing time lag and magnitude difference between Q F,B and Q EC .Therefore, estimation of Q F,B by differences in heat emission between occupied and unoccupied buildings can capture the impact of dynamic changes in the building storage heat flux especially in a sub-annual temporal cycle.
In this study, the objective is to understand the temporal profile of Q F,B and how and why it differs from Q EC at diurnal and seasonal time scales, by examining differences in energy balance fluxes between the same occupied and unoccupied building.A building energy simulation tool (Energy-Plus) is used to obtain the various energy balance fluxes from the building system.

Unoccupied (uo) and occupied (o) building energy simulation (BES)
Building energy simulation (BES) is widely used to estimate energy consumption, heat emission and heat storage within a building, while allowing changes in heat fluxes due to human activities to be estimated.Here we use Energy-Plus version 9.4 (DOE, 2020a) to study an isolated building (i.e.without a surrounding neighbourhood).The AS-NI/ASHRAE standard 140 Case 900 test model (ASHRAE, 2017) is used, which is developed in software-to-software comparative tests for validating building thermal load.It is a 48 m 2 one-story heavyweight rectangular prism with high mass fabrics (Fig. S1 in the Supplement), whose simple geometry is ideal to understand the process of how human activities change the building energy balance fluxes in a theoretical study.
Modifications of the original building model for this study, include: windows are reduced to one (6 m 2 south-facing) for more appropriate EnergyPlus single-sided ventilation calculations (Daish et al., 2016) and internal heat gain, ventilation control strategy and HVAC system operation are varied with different scenarios considered (Table 1).For the simulations, the building is assumed to be located in Beijing as the climate has both hot summer and cold winter conditions.Chinese standard weather data (CSWD) selected to create a typical meteorological year (TMY) (China Meteorological Bureau et al., 2005) are used as the meteorological forcing, as these data were developed for simulating building thermal load and energy use.
The modelling scenarios (Table 1) vary with building occupation state.Two types of unoccupied (uo) buildings are considered.Neither have internal heat gains or HVAC systems, but they differ based on air exchange between (1) unoccupied sealed (us) with no infiltration nor ventilation, and (2) unoccupied ventilated (uv) with 50 % of window area kept open.The single-sided natural ventilation rate is estimated by including both the wind-driven ventilation rate (V W , m 3 s −1 ) (Warren 1977): and the stack buoyancy-driven ventilation rate (V , m 3 s −1 ) (Warren 1977): where A eff is the effective opening area (m 2 ), U W is the reference wind speed at the height of opening (m s −1 ).C d is the discharge coefficient (usually taken as 0.6, Wang and Chen ( 2012)), T is the indoor and outdoor air temperature difference ( • C), H is the height of opening (m), g the gravitational acceleration (m s −2 ) and T ave is the average indoor and outdoor air temperature ( • C).The combined ventilation rate is (Fan et al., 2021) The three occupied (o) building simulations assume that occupant behaviour modifies internal heat generation, natural ventilation and HVAC systems (ov).First, ov 1 has internal heat gains (Q Internal, o ) from human metabolism, lighting and other appliances based on local building code (MOHURD, 2018), with window always open (50 %, as uv).The internal heat gains are held constant allowing the fraction of heat in Q F,B and Q S to be impacted by building and climate conditions but not the diurnal variability of human heat generation.
Second, ov 2 considers natural ventilation based on passive cooling and thermal comfort.The window opening is controlled automatically.It is opened (50 % of window area) when the indoor air temperature is higher than both outdoor air temperature and the ventilation setting point (23 • C for "warm limit" in the bedroom (Oikonomou et al., 2012)).Otherwise, it is closed to reduce heat loss and keep the building warm.Third, since natural ventilation alone may not satisfy indoor thermal comfort, mixed mode ventilation with an auxiliary HVAC system (e.g.Wang and Chen, 2013;Wang and Greenberg, 2015;Chen et al., 2017) is considered in ov 3 .The mechanical heating and cooling systems are active when the indoor temperature reaches the threshold (18 • C for heating and 26 • C for cooling, MOHURD, 2018).The ventilation control strategy in ov 3 is the same as ov 2 , but the EnergyPlus hybrid ventilation manager (DOE, 2020b) turns the HVAC off when natural ventilation is active to prevent simultaneous operation.

Determination of anthropogenic heat flux
The simulated hourly heat fluxes by radiation, convection, air exchange and waste heat generated from the HVAC sys- tem between the isolated building and the atmosphere (Table S3 in the Supplement) are analysed for each case (Table S2 in the Supplement).If cooling occurs, the waste heat consists of the cooling load and electrical energy consumed by the air conditioner (Q HVAC ).The Q HVAC is predicted using a static coefficient of performance (COP) for the air conditioner, and the heat removed by an air conditioner (Q AC ) to the total amount of electricity consumed: With a centralised heating system (as Beijing has), for simplicity we assume all energy associated with the heating system is released indoors, and waste heat due to boiler efficiency and pipe heat loss are not considered: Combining mechanical heating and cooling, the energy consumption and corresponding waste heat from HVAC system gives Each term in Eq. ( 14) is determined using an occupied (o) and unoccupied (uo) building result to determine Q F,B and the other fluxes.The results are analysed by season as spring (March, April and May, MAM), summer (June, July, August, JJA), autumn (September, October, November, SON) and winter (December, January, February, DJF)) using the median (50 %) and interquartile range (IQR) between the 25th and 75th percentiles to assess the diurnal patterns.

Ratio of anthropogenic heat flux to energy consumption
If the energy consumed within the building is ejected immediately into the atmosphere (Heiple and Sailor, 2008), the change in Q S is not accounted for, and therefore Q F,B is assumed to be only from energy consumption (Q EC ).The variation of Q S associated with human activities is considered when using the relative heat emissions in Eqs. ( 14) and ( 15).We use the ratio R = Q EC to determine the relative importance of building operation modes and choice of baselines on the discrepancy between Q F,B and Q EC .

Results and discussion
Building energy balance fluxes vary through each day and season (Fig. 1) associated with when a building is occupied and people's activities inside the building.First, we consider one case in detail -an occupied building with both natural ventilation and HVAC (ov 3 , Table 1) relative to an unoccupied sealed building (us, Table 1) -their difference (ov 3 -us) allows us to obtain the fluxes needed (Sect.1).
As noted (Sect.1), the shortwave and incoming longwave radiation fluxes for all cases (Table 1) are assumed to be identical, but all other terms of the building energy balance differ.Hence, the change in outgoing longwave radiation ( L ↑o-uo , Fig. 1c) is equivalent to the net allwave radiation difference (Q * o-uo , Fig. 1a-b) for the occupied and unoccupied buildings.The positive sensible heat flux difference (Eq. 10, H o-uo , Fig. 1c) and L ↑o-uo indicate the building is warmed up by internal heat gains (Q Internal, o ) with higher exterior surface temperatures.Their small magnitudes and flat patterns indicate small relative importance compared to the heat exchange from ventilation differences (Eq.8, BAE o-uo , Fig. 1c).The latter not only contributes the largest fraction of anthropogenic heat flux (Q F,B , Fig. 1c), but also has a diurnal pattern consistent with Q F,B , especially during spring and autumn (Fig. 1c, i).Rarely, heat (Q Waste, o , Fig. 1i) is emitted by the air conditioner in the mid-afternoon (shading) at this time of year, but more importantly in summer (Fig. 1f) when cooling demand increases.
The Q F,B (Eq.14, Fig. 1c) has four components of emitted heat, whereas energy consumption (Q EC , Fig. 1c) only has (in this case, constant) internal heat gains (Q Internal, o = 11.8W m −2 , Fig. 1b  systems (Q HVAC , Fig. 1b).Their difference is the storage heat flux difference (Eq.15 S o-uo in Fig. 1c).If S o-uo is positive, the building acts as a heat sink and stores the extra heat generated by human activities, or stored heat is released when S o-uo is negative.Hence, we can identify the impacts of seasonal varying human activities and building operations on the diurnal variability in S o-uo , Q EC and Q F,B .

Impact of human activities on seasonal and diurnal variations of the fluxes
For the same ov 3 -us case (Table 1, Fig. 1), we consider the diurnal and seasonal variability of the fluxes.In spring and autumn (Fig. 1a-c, g-i), natural ventilation is the dominant factor contributing to diurnal variation in S o-uo and Q F,B , while Q EC has minimal variability.The Q EC is slightly larger than Q Internal, o because of some short periods of HVAC use in the mid-afternoon (IQR shading in Fig. 1i).There is a clear diurnal cycle of Q F,B (Fig. 1c) with the median varying between 8 W m −2 (07:00) and 15 W m −2 (15:00) relative to the constant internal heat gain (11.8 W m −2 ).The difference between Q F,B and Q EC ( S o-uo ) is largely impacted by natural ventilation.During the night and early morning with closed windows, only part of the consumed energy is transferred externally to the atmosphere.The rest of the heat is stored in the building fabric (positive S o-uo ), hence Q F,B is lower than Q EC .However, when overheating may occur during the middle of the day, occupants keep the window opened (air conditioner is less frequently used) to cool the building down, with stored heat released (negative S o-uo ).This is consistent with the diurnal variability of BAE o-uo , which has a minimum at night (window closed) and maximum in the mid-afternoon (window open).
In summer, the daytime natural ventilation is replaced by air conditioning as natural ventilation alone could not maintain thermal comfort indoors.Natural ventilation and waste heat from the air conditioner (Q Waste, o ) contribute to one peak Q F,B at nighttime and daytime, respectively (Fig. 1f).Q F,B is higher than Q EC around these two peak periods (05:00-07:00 and 13:00-21:00).The peak Q F,B at night reaches 14 W m −2 (median) at 05:00, which is mainly attributed to natural ventilation when outdoor air temperature is cooler than indoors.Conversely, in the afternoon when outdoor temperature is warmer, occupants "choose" mechanical cooling to achieve thermal comfort.The peak Q F,B is 22 W m −2 at 16:00, approximately 22 % higher than Q EC .It indicates that using Q EC for the anthropogenic heat flux from buildings (e.g.Heiple and Sailor, 2008) may underestimate the effect of Q F,B on urban atmospheric processes especially during the late afternoon/early evening.In addition, Q F,B is always smaller than Q Waste, o because of the negative L ↑o-uo and H o-uo causing a cooler exterior surface.This suggests using Q Waste, o as Q F,B (e.g.Chow et al., 2014) may overestimate Q F,B in summer.
However, in winter mechanical heating and thermal mass effect shape the temporal pattern of Q F,B (Fig. 1i).The cool outdoor air temperature before sunrise results in a substantial heating supply and peak Q EC (16.43 W m −2 for median line) at 08:00.This heat is stored in building fabric (positive S o-uo ) and has a relatively stable release through convection and longwave radiation.Therefore the diurnal profile Q F,B is rather flatter and S o-uo has a highly consistent temporal pattern to Q EC .
Overall, this analysis recognises the crucial role of S o-uo in distinguishing Q F,B from Q EC , which is highly dependent on HVAC operation and natural ventilation (i.e.human activity of window opening).These two factors can rapidly increase or decrease Q F,B while convection and longwave radiation cannot.Whereas in winter, the larger IQR (shading) of Q F,B than Q EC indicates more day-to-day variation in Q F,B diurnal profile than Q EC .Estimates of Q F,B using satellite remote sensing found heat storage plays an important role in moderating energy use within buildings (Yu et al., 2021).As the storage heat flux change modifies the diurnal sensible heat flux pattern it modifies the surface temperature increment (Q F,B in remote sensing approach) and hence the apparent energy consumption.
The diurnal profiles of S o-uo are not identical between seasons as people use different actions to achieve thermal comfort in different weather conditions.This suggests that the Q F,B and Q EC differences may vary between climates and with cultural practices.In inventory methods the diurnal profiles may be limited, e.g.large scale urban consumption of energy (LUCY, Allen et al. (2011), weekday/weekend by country, and ignore seasonal variations.However, S o-uo behaviour type classes may benefit from distinguishing diurnal variation for different climates.

Impact of different building operation modes on seasonal and diurnal variations
Figure 2 illustrates the impact of different building operation modes (Table 1: ov 1 , ov 2 , ov 3 ; cf.us) on the Q F,B diurnal profiles.It suggests that the different ventilation strategies and HVAC systems do change Q F,B in both temporal pattern and magnitude, but their impacts vary among seasons.
In spring and autumn, different natural ventilation control strategies completely modify the Q F,B diurnal profile, whereas a HVAC system only increases the peak Q F,B slightly in autumn (Fig. 2i).The distinctly different (opposite) trend in diurnal Q F,B pattern for ov 1 cf.ov 2 or ov 3 (Fig. 2a-c, g-i) is largely explained by the diurnal change of BAE o-uo in the three cases.In ov 1 (window open, no control) the minimum outdoor air temperature before sunrise creates the maximum indoor and outdoor air temperature difference, therefore the highest BAE o-uo and peak Q F,B at 06:00 (30 W m −2 for the median in Fig. 2a).Whereas ov 2 and ov 3 have the window closed at night and early morning to avoid overcooling, therefore, the minimum Q F,B is in the https://doi.org/10.5194/acp-22-4721-2022Atmos.Chem.Phys., 22, 4721-4735, 2022 early morning (07:00).As outdoor air temperature increases through the day, Q F,B follows the reduced BAE o-uo in ov 1 , whereas natural ventilation is active in ov 2 and ov 3 , leading to an increase in BAE o-uo and Q F,B .Unlike ov 2 , ov 3 has a clear peak (16 W m −2 median, Fig. 2i) at 15:00, because natural ventilation alone cannot satisfy thermal comfort and ov 3 air conditioning is activated.But their overall patterns (IQR) are very consistent, indicating afternoon use of air conditioning could increase Q F,B magnitude but have a limited impact on other parts of the diurnal pattern.Surprisingly, negative Q F,B occurs around 17:00 in spring (Fig. 2a), suggesting the occupied building has less heat emission than the unoccupied building, because the natural ventilation at night and morning cools the building down and reduced fabric exterior surface temperature leads to a larger reduction in longwave radiation and convection ( L ↑o-uo and H o-uo ) than the increase in heat emission through natural ventilation ( BAE o-uo ) in the afternoon.Also, the reduced overall emissions are converted into increase in storage heat flux ( S o-uo ).Negative Q F,B also occurs when the unoccupied building is always ventilated (uv) and the occupied building is ventilated with control (ov 2 and ov 3 ) in spring (e.g.Fig. S7b-c in the Supplement).The window is closed to avoid excessive cooling at night in ov 2 .With BAE o-uo negative in this case, its magnitude is much larger than the increase in longwave radiation and convection ( L ↑o−uo and H o-uo ).The minimum Q F,B frequently corresponds to the peak S o-uo .
In summer, in ov 2 the window is open most of the time (as in ov 1 ) for thermal comfort, therefore, the Q F,B has no apparent difference to ov 1 .However for ov 3 , as the air conditioning runs from morning to late night there is a very different diurnal profile (cf.ov 2 and ov 1 ).Air conditioner use contributes to a much larger Q F,B (cf.ov 2 ) from 12:00 to 21:00.Not only is extra energy consumed, but it also removes heat from the building to the atmosphere in this period.In contrast, using natural ventilation as a cooling strategy (ov 1 and ov 2 ) contributes to a high Q F,B at night and early morning but very low even negative extra heat emission in the afternoon.This implies that natural ventilation as a passive cooling strategy could not only improve the thermal conditions indoors but could also contribute to the improvement of outdoor climate by modifying the diurnal pattern of anthropogenic heat emissions (Duan et al., 2019).
Consistent with results in the other seasons, different ventilation control strategies in winter cause a large change in the Q F,B profile between ov 1 and ov 2 .However, the temporal pattern of Q F,B (IQR) in ov 2 is quite similar to ov 3 because the supplied heat from the mechanical heating system does not immediately enhance Q F,B with a closed window.Ov 2 is the only scenario that has similar Q F,B and Q EC through the whole day.Comparison using an unoccupied ventilated (uv) baseline (Fig. S7) (cf.us Fig. 2) show that although Q F,B profiles differ, the impacts of different building operation modes are consistent when the same occupied buildings are used.The impact of baselines with different air exchange on Q F,B are analysed in Sect.3.3.

Impact of unoccupied baseline chosen
Here two unoccupied baselines (us -unoccupied sealed building, uv -unoccupied ventilated building with uncontrolled open window) are used to assess the impact.A ratio between Q F,B and Q EC (R) is used (Fig. 3) to normalise the impact of baselines on their difference with different building operation modes.The largest difference in R occurs on 23 December at 11:00, with values of 5.13 (ov 3 -uv) and −2.72 (ov 1 -us), reflecting the considerable difference between Q F,B and Q EC .
Two diurnal patterns of the R ratio are distinguished.When the window is always open (ov 1 in all seasons, ov 2 in summer), R>1 (Q F,B >Q EC ) at night/early morning (22:00-08:00), reaching its maximum around 05:00-07:00 (near sunrise in all seasons).For the remaining periods, which are relatively warm R<1.Whereas, when window opening/closing is controlled and HVAC is used for thermal comfort an almost inverse temporal pattern of R occurs, with R>1 during the afternoon when either the window is open or the air conditioner is activated.The peak R occurs at 15:00 when both outdoor temperature and solar radiation are high.
When different unoccupied baselines are used, the temporal patterns of R are similar for all cases, but their magnitudes differ significantly.R is close to 1 when window states between unoccupied and occupied buildings are similar (e.g.ov 1 -uv in all seasons, ov 2 -uv in summer).Hence, a greater difference occurs in heat transfer from ventilation or mechanical heating/cooling between occupied and unoccupied buildings (i.e.larger R).Thus, the baseline chosen impacts the results and requires appropriate consideration for incorporating Q F,B into atmospheric modelling.

Comparison between Q F,B and building heat emission (BHE)
Comparison of building heat emissions (BHE), determined using the Hong et al. (2020) approach, to Q F,B (this study) for one case (ov 3 -us) shows that the former is much larger than Q F,B during the day but smaller at night and has different diurnal patterns (Fig. 4).Convection from the exterior envelope (Q H , Fig. 1b, e, h, k) is the main contributor to BHE and therefore influences the BHE diurnal profile in each season.During the day, solar radiation is high and a major control whereas Q F,B is relatively small and consistent but modified by building-human interactions (e.g.opening windows, activation of mechanical heating and cooling systems).In this scenario shown, natural ventilation and mechanical cooling dominate Q F,B in summer and shoulder season (i.e.spring and autumn); while in winter in their absence, convection and longwave radiation are more important.

Daily variation of fluxes in relation to meteorological conditions
Ambient air temperature is one of the most crucial factors controlling building energy consumption (Sailor and Vasireddy, 2006).Hence, it is often used to determine daily variability of Q EC (e.g.Lindberg et al., 2013) and the resulting monthly variations (e.g.Allen et al., 2011).By accounting for S o-uo in this study, the response of Q F,B to ambient air temperature may differ to previous studies.To examine this we used the ov 3 -us case to consider the relations of daily mean (unless indicated) variables of air tempera-  ture (mean) , solar radiation (daily total) and simulated available energy to the building from human activities ( B o-uo ) with anthropogenic heat flux (Q F,B in Fig. 5a), energy consumption (Q EC in Fig. 5b) and their difference ( S o-uo in Fig. 5c).The overall trends between Q F,B and Q EC to ambient air temperature are consistent, with Q F,B and Q EC small-est when temperatures are between 10-15 • C.This coincides with the Nicol and Humphreys (2002) monthly balancepoint temperature of 12 • C, which has been regarded as the equivalent ambient air temperature with the minimum energy use within the building (e.g.Allen et al., 2011, Koralegedara et al., 2016).As the temperature increases (de- creases), Q EC increases proportionally with temperature due to mechanical cooling (heating).However, in contrast to Q EC , Q F,B has a much larger variability at the same temperature caused by a large range of S o-uo (−7.7 to 9.0 W m −2 ), which is highly dependent on human activities on the diurnal scale (Sect.3.1).
To understand the large daily variability of S o-uo , we use B o-uo (net available energy from human activities in buildings in Eq. ( 9)) to indicate the effect of human activities (heat addition or removal) in one day.Higher B o-uo (larger circles) are associated with higher S o-uo at the same ambient air temperature, especially in winter (Fig. 5c).This is not unexpected as buildings will absorb more heat when extra internal energy is added into the building.Inversely, negative B o-uo (small circles) contributes to much more heat release from heat storage (lower S o-uo through either natural ventilation or mechanical cooling.The sign and magnitude of B o-uo are linked to daily cumulative solar radiation.At the same ambient air temperature, more solar radiation enhances the need for larger heat removal or less heat addition to the building for thermal comfort, therefore leading to a smaller B o-uo and lower S o-uo .Consequently, we can conclude that both ambient air temperature and cumulative solar radiation are important meteorological factors to determine S o-uo and Q F,B .

Conclusions
Anthropogenic heat flux from buildings (Q F,B ) is defined as the additional heat released from the building into the atmosphere due to human activities.It is qualitatively and quantitatively different to building energy consumption (Q EC ) in temporal pattern and magnitude as a result of thermal inertia of the building (Iamarino et al., 2012).However, as there is no standard to quantify "real" Q F,B most studies use Q EC as a proxy via inventory and building energy modelling approaches.This paper proposes a new method to quantify a more appropriate Q F,B by utilising the difference in heat fluxes between an occupied and unoccupied building (i.e. the built structure with absolutely no energy use and no human metabolism).We show that the difference between Q EC and Q F,B is attributable to a change in the storage heat flux induced by human activities ( S o-uo ).Q F,B has four components based on its dissipation pathways, including outgoing longwave radiation, turbulent sensible heat flux (convection), heat release due to air exchange and waste heat from HVAC systems.We use one simplified case study in Beijing to demonstrate the analysis using building energy simulations to quantify the temporal difference between Q EC and Q F,B and to understand the relative importance of building operations for thermal comfort and meteorological conditions on Q F,B .The key conclusions are 1.Hourly ratios between Q F,B and Q EC can differ between −2.72 and 5.13 because of differences in occupancy use of the building (within a year, in Beijing's climate).Individual ratios frequently exceed 3 between 14:00 and 16:00 when controlled natural ventilation or mechanical cooling is activated in a shoulder season (i.e.spring and autumn).Thus, the differences in the definitions are large.
2. Natural ventilation ( BAE o-uo ) or HVAC operation (Q Waste, o for cooling and Q HVAC for heating) are two predominant contributors to the storage heat flux.Hence, different building operations to control thermal comfort determine the diurnal profile of Q F,B by affecting not only Q EC but also S o-uo .
3. The day-to-day variation of the Q F,B diurnal profile is broader than that of Q EC .
4. The diurnal profile of S o-uo varies with season as occupants modify their behaviour and the interaction with buildings to achieve thermal comfort (e.g.cooling in summer and heating in winter), indicating that differences between Q F,B and Q EC will vary with both climate and cultural norms.
6. Daily mean temperature only accounts for the day-today variability in Q EC rather than S o-uo .Both ambient air temperature and cumulative solar radiation are important meteorological factors to determine S o-uo and Q F,B .
Our new approach should be used to provide data for future parameterisations of both anthropogenic heat flux from buildings and storage heat fluxes for urban weather and climate modelling.We conclude that storage heat fluxes in cities are also being modified by occupant behaviour, particularly by natural ventilation and mechanical cooling.It is expected that the diurnal variation of S o-uo will vary with operation schedules for different building uses (e.g.residential vs. commercial buildings).Given that the release of stored heat has a critical influence on the nocturnal canopy layer urban heat island (CL-UHI), the impact of different HVAC operations on nocturnal UHI should be explored further.This is an important factor to determine diurnal pattern of Q F,B in the shoulder season and can be expressed more accurately.However, in different climates and with different social cul-tural practices the periods most influenced will change.Further studies are being conducted to explore the impacts of these, while also addressing feedback at the neighbourhood scale.
For developers of urban canopy parameterisations (UCP) there are several considerations because of computational efficiencies which are essential for undertaking weather and climate modelling: (1) human activities within a building are modifying both the storage heat flux and the anthropogenic heat flux; (2) assuming within an UCP that a "simple" building energy model (BEM) (cf. a full building energy simulation scheme such as EnergyPlus) will require some human activities to be simplified, such as using a fixed ventilation rate, instead of dynamic natural ventilation depending on both outdoor weather conditions and thermal comfort requirements and (3) with a multi-layer UCP the appropriate levels for the impact of these energy exchanges can be accounted for.Our current research is extending this analysis to consider moisture and exploring the role of building materials, construction, other aspects of building design and external meteorology.The outcome of this work will also have implications for UCP development, as this can help identify what can be simplified and what are critical controls in different climates and urban settings.

Figure 2 .
Figure 2. As Fig. 1c, f, i, j , but comparing three different building operation types (a, d, g, j) ov 1 : window is always open without control, no HVAC; (b, e, h, k) ov 2 : controlled natural ventilation for indoor thermal comfort, no HVAC; (c, f, i, l) ov 3 : mixed mode ventilation.

Figure 3 .
Figure 3. Q F,B to Q EC ratio I median (line) and IQR (shading) for (a-b) spring, (c-d) summer, (e-f) autumn and (g-h) winter, using two unoccupied baselines: (a, c , e, g) sealed (us), and (b, d, f, h) ventilation (uv); each with three occupancy types (colour): ov 1 : only internal heat gains are applied and window is fully open; ov 2 : internal heat gains and natural ventilation control are applied.Ov3: internal heat gains, natural ventilation control and HVAC system are applied.Ratio R = 1 (Black dotted line).

Figure 5 .
Figure 5. Daily results for the ov 3 -us case stratified by daily cumulative solar radiation (colour) and daily mean available energy to the building (size) Eq. (9) associated with human activities, with mean external air (ambient) temperature and (a) mean anthropogenic heat flux, (b) energy consumption and (c) difference in storage heat flux.

Table 1 .
Cases simulated differ based on building occupation state, internal heat gain (Q Internal, o ) and presence of natural ventilation and HVAC.Notation is defined in the text and nomenclature.