The Impact of Landscape Fragmentation on Atmospheric Flow: A Wind-Tunnel Study

Landscape discontinuities such as forest edges play an important role in determining the characteristics of the atmospheric flow by generating increased turbulence and triggering the formation of coherent tree-scale structures. In a fragmented landscape, consisting of surfaces of different heights and roughness, the multiplicity of edges may lead to complex patterns of flow and turbulence that are potentially difficult to predict. Here, we investigate the effects of different levels of forest fragmentation on the airflow. Five gap spacings (of length approximately 5h, 10h, 15h, 20h, 30h, where h is the canopy height) between forest blocks of length 8.7h, as well as a reference case consisting of a continuous forest after a single edge, were investigated in a wind tunnel. The results reveal a consistent pattern downstream from the first edge of each simulated case, with the streamwise velocity component at tree top increasing and turbulent kinetic energy decreasing as gap size increases, but with overshoots in shear stress and turbulent kinetic energy observed at the forest edges. As the gap spacing increases, the flow appears to change monotonically from a flow over a single edge to a flow over isolated forest blocks. The apparent roughness of the different fragmented configurations also decreases with increasing gap size. No overall enhancement of turbulence is observed at any particular level of fragmentation.     

Development and Validation of an Evaporation Duct Model.Part Ⅰ:Model Establishment and Sensitivity Experiments

Based on the Coupled Ocean-Atmospheric Response Experiment（COARE）bulk algorithm and the Naval Postgraduate School（NPS）model,a universal evaporation duct（UED）model that can flexibly accommodate the latest improvements in component（such as stability function,velocity roughness,and scalar roughness）schemes for different stratification and wind conditions,is proposed in this paper.With the UED model,the sensitivity of the model-derived evaporation duct height（EDH）to stability function（Ψ）,ocean wave effect under moderate to high wind speeds,and scalar roughness length parameterization,is investigated,and relative contributions of these factors are compared.The results show that the stability function is a key factor influencing the simulated EDH values.Under unstable conditions,the EDH values from stability functions of Fairall et al.（1996）and Hu and Zhang（1992）are generally higher than those from others;while under stable conditions,unreasonably high EDHs can be avoided by use of the stability functions of Hu and Zhang（1992）and Grachev et al.（2007）.Under moderate to high wind speeds,the increase in velocity roughness length z₀ due to consideration of the true ocean wave effect acts to reduce modeled EDH values;this trend is more pronounced under stable conditions.Although the scalar roughness length parameterization has a minor effect on the model-derived EDH,a positive correlation is found between the scalar roughness length z_0qand the model-derived EDH.     

Pollutant Dispersion in Boundary Layers Exposed to Rural-to-Urban Transitions: Varying the Spanwise Length Scale of the Roughness

Both large-eddy simulations (LES) and water-tunnel experiments, using simultaneous stereoscopic particle image velocimetry and laser-induced fluorescence, have been used to investigate pollutant dispersion mechanisms in regions where the surface changes from rural to urban roughness. The urban roughness was characterized by an array of rectangular obstacles in an in-line arrangement. The streamwise length scale of the roughness was kept constant, while the spanwise length scale was varied by varying the obstacle aspect ratio l / h between 1 and 8, where l is the spanwise dimension of the obstacles and h is the height of the obstacles. Additionally, the case of two-dimensional roughness (riblets) was considered in LES. A smooth-wall turbulent boundary layer of depth 10h was used as the approaching flow, and a line source of passive tracer was placed 2h upstream of the urban canopy. The experimental and numerical results show good agreement, while minor discrepancies are readily explained. It is found that for \(l/h=2\) the drag induced by the urban canopy is largest of all considered cases, and is caused by a large-scale secondary flow. In addition, due to the roughness transition the vertical advective pollutant flux is the main ventilation mechanism in the first three streets. Furthermore, by means of linear stochastic estimation the mean flow structure is identified that is responsible for street-canyon ventilation for the sixth street and onwards. Moreover, it is shown that the vertical length scale of this structure increases with increasing aspect ratio of the obstacles in the canopy, while the streamwise length scale does not show a similar trend.     

Experimental Investigation on Shear-Stress Partitioning for Flexible Plants with Approximately Zero Basal-to-Frontal Area Ratio in a Wind Tunnel

Shear-stress partitioning is investigated for one type of flexible plant for very small values of the basal-to-frontal area ratio σ (0.001–0.007). The plant model is made of plastic with irregular structures, which are different from previously investigated rigid regular or flexible roughness elements with larger σ values. The distribution of the surface shear stress and the total shear stress at four plant densities with five plant heights are measured in a wind tunnel using Irwin-type sensors and a load cell, respectively. The wind-tunnel experiments prove that, for these flexible plants, the plant height and lateral cover usually decrease with increasing friction velocity, especially for taller plants, while the plant coverage generally increases. However, these characteristics may be inconsistent with flexible roughness elements with very large σ values (and usually very low aspect ratios) because these elements are less flexible. The present flexible plants generally result in lower shear-stress ratios compared with other roughness elements, which is also proven by the higher values of β (the ratio of the drag coefficient of an isolated roughness element to that of the bare surface) and a constant m (accounting for the difference between the average and peak surface shear stresses) from the present experiments (β?=?184–210 and m?=?0.68–0.79). The peak mean stress ratio of the present flexible plants is not a constant (1.07–1.54) because it is affected by the lateral cover, which is different from previous studies that consider the ratio to be constant without regard for the lateral cover.     

Classification of Vertical Wind Speed Profiles Observed Above a Sloping Forest at Nighttime Using the Bulk Richardson Number

Wind speed profiles above a forest canopy relate to scalar exchange between the forest canopy and the atmosphere. Many studies have reported that vertical wind speed profiles above a relatively flat forest can be classified by a stability index developed assuming wind flow above a flat plane. However, can such a stability index be used to classify vertical wind speed profiles observed above a sloping forest at nighttime, where drainage flow often occurs? This paper examines the use of the bulk Richardson number to classify wind speed profiles observed above a sloping forest at nighttime. Wind speed profiles above a sloping forest were classified by the bulk Richardson number Ri _B . Use of Ri _B distinguished between drainage flow, shear flow, and transitional flow from drainage flow to shear flow. These results suggest that Ri _B is useful to interpret nighttime CO₂ and energy fluxes above a sloping forest. Through clear observational evidence, we also show that the reference height should be high enough to avoid drainage-flow effects when calculating Ri _B .     

Wind-Speed Undulations Over Swell: Field Experiment and Interpretation

Results of field measurements of the swell-induced undulation of the wind speed taken from a Black Sea platform are presented. The wind speed and its fluctuations were measured at several heights between 1.3 and 21 m above the mean sea level under various wind and swell conditions. Parameters of the swell-induced undulations were derived from cross spectra of the wind-speed fluctuations and the sea-surface displacement. As found, the phase and the amplitude of the wind speed undulation in the layer from k _p z = 0.1 to k _p z = 3 (k _p is the swell wavenumber) are in good agreement with the theory of inviscid shear flow over a wavy surface. The main feature of the vertical profile of the swell-induced undulation is the exponential attenuation of its amplitude with height typical for the potential flow over the fast running waves. At the lowest levels the potential undulations are significantly distorted by the wind-speed variations caused by the vertical displacements of the shear airflow relative to a fixed sensor. No direct impact of swell on the mean properties of the turbulent boundary layer at k _p z > 0.1 is revealed. In particular, the mean wind-speed profile and spectra of the horizontal velocity in the inertial subrange obey Monin-Obukhov similarity theory.     

On the extension of the wind profile over homogeneous terrain beyond the surface boundary layer

Analysis of profiles of meteorological measurements from a 160 m high mast at the National Test Site for wind turbines at Høvsøre (Denmark) and at a 250 m high TV tower at Hamburg (Germany) shows that the wind profile based on surface-layer theory and Monin-Obukhov scaling is valid up to a height of 50–80 m. At higher levels deviations from the measurements progressively occur. For applied use an extension to the wind profile in the surface layer is formulated for the entire boundary layer, with emphasis on the lowest 200–300 m and considering only wind speeds above 3 m s^?1 at 10 m height. The friction velocity is taken to decrease linearly through the boundary layer. The wind profile length scale is composed of three component length scales. In the surface layer the first length scale is taken to increase linearly with height with a stability correction following Monin-Obukhov similarity. Above the surface layer the second length scale (L _MBL ) becomes independent of height but not of stability, and at the top of the boundary layer the third length scale is assumed to be negligible. A simple model for the combined length scale that controls the wind profile and its stability dependence is formulated by inverse summation. Based on these assumptions the wind profile for the entire boundary layer is derived. A parameterization of L _MBL is formulated using the geostrophic drag law, which relates friction velocity and geostrophic wind. The empirical parameterization of the resistance law functions A and B in the geostrophic drag law is uncertain, making it impractical. Therefore an expression for the length scale, L _MBL , for applied use is suggested, based on measurements from the two sites.     

Projection methods for error suppression in the meteorological fields before calculation of derivatives

In finite-difference approximation of derivatives on a grid with step of h, the error decreases with decreasing h. However, if the original function is given on a grid with an error of order of ?, then the approximate differentiation error grows like C?h ^k , where k is the derivative order and C is a constant depending on a chosen finite-difference approximation. Projection methods are discussed that are designed to decrease the errors of this kind. On the latitude-longitude grids, the step near the poles along the latitude circles is very small. Methods are proposed for projecting scalar and vector meteorological fields on the subspaces of smooth fields. The projected fields obtained take into account asymptotics of the smooth fields in the vicinity of the pole in the spherical coordinates, separately for each Fourier harmonic in longitude. Results are presented of such smoothing in the polar regions on the fields of objective analysis by Hydrometeorological Center of Russia and on the first-guess fields at standard isobaric heights.     

Velocity Adjustment and Passive Scalar Diffusion in and Above an Urban Canopy in Response to Various Approach Flows

We used wind-tunnel experiments to investigate velocity-field adjustment and scalar diffusion behaviour in and above urban canopies located downwind of various roughness elements. Staggered arrays of rectangular blocks of various heights H and plan area ratios λ_p were used to model the urban canopies. The velocity field in the roughness sublayer (height \({z \lesssim 2H}\)) reached equilibrium at distances proportional to \({\sqrt{L_{\rm c}H}}\) where L _c is the canopy-drag length scale determined as a function of λ_p and the block side length L. A distance of about \({20\sqrt{L_{\rm c}H}}\) was required for adjustment at z = H/2 (in the canopy), and a distance of about \({10\sqrt{L_{\rm c}H}}\) was required at z = 2H (near the top of the roughness sublayer). Diffusion experiments from a ground emission source revealed that differences in upwind roughness conditions had negligible effects on the plume growth near the source (up to a few multiples of L from the source) if the source was located at a fetch F larger than about \({10\sqrt{L_{\rm c}H}}\) from the upwind edge of the canopy. However, at locations farther downwind (more than several multiples of L from the source), upwind conditions had considerable effects on the plume growth. For a representative urban canopy, it was shown that a much larger fetch than required for velocity-field adjustment in the roughness sublayer was necessary to eliminate the effects of upwind conditions on plume widths at 24L downwind from the source.     

Characteristics of Turbulent Transfer during Episodes of Heavy Haze Pollution in Beijing in Winter 2016/17

We analyzed the structure and evolution of turbulent transfer and the wind profile in the atmospheric boundary layer in relation to aerosol concentrations during an episode of heavy haze pollution from 6 December 2016 to 9 January 2017. The turbulence data were recorded at Peking University’s atmospheric science and environment observation station. The results showed a negative correlation between the wind speed and the PM_2.5 concentration. The turbulence kinetic energy was large and showed obvious diurnal variations during unpolluted (clean) weather, but was small during episodes of heavy haze pollution. Under both clean and heavy haze conditions, the relation between the non-dimensional wind components and the stability parameter z/L followed a 1/3 power law, but the normalized standard deviations of the wind speed were smaller during heavy pollution events than during clean periods under near-neutral conditions. Under unstable conditions, the normalized standard deviation of the potential temperature σ _θ /|θ_*| was related to z/L, roughly following a –1/3 power law, and the ratio during pollution days was greater than that during clean days. The three-dimensional turbulence energy spectra satisfied a –2/3 power exponent rate in the high-frequency band. In the low-frequency band, the wind velocity spectrum curve was related to the stability parameters under clear conditions, but was not related to atmospheric stratification under polluted conditions. In the dissipation stage of the heavy pollution episode, the horizontal wind speed first started to increase at high altitudes and then gradually decreased at lower altitudes. The strong upward motion during this stage was an important dynamic factor in the dissipation of the heavy haze.     

A scattering-based over-land rainfall retrieval algorithm for South Korea using GCOM-W1/AMSR-2 data

Heavy summer rainfall is a primary natural disaster affecting lives and properties in the Korean Peninsula. This study presents a satellite-based rainfall rate retrieval algorithm for the South Korea combining polarization-corrected temperature (PCT) and scattering index (SI) data from the 36.5 and 89.0 GHz channels of the Advanced microwave Scanning Radiometer 2 (AMSR-2) onboard the Global Change Observation Mission (GCOM)-W1 satellite. The coefficients for the algorithm were obtained from spatial and temporal collocation data from the AMSR-2 and groundbased automatic weather station rain gauges from 1 July - 30 August during the years, 2012-2015. There were time delays of about 25 minutes between the AMSR-2 observations and the ground raingauge measurements. A new linearly-combined rainfall retrieval algorithm focused on heavy rain for the PCT and SI was validated using ground-based rainfall observations for the South Korea from 1 July - 30 August, 2016. The validation presented PCT and SI methods showed slightly improved results for rainfall > 5 mm h^-1 compared to the current ASMR-2 level 2 data. The best bias and root mean square error (RMSE) for the PCT method at AMSR-2 36.5 GHz were 2.09 mm h^-1 and 7.29 mm h^-1, respectively, while the current official AMSR-2 rainfall rates show a larger bias and RMSE (4.80 mm h^-1 and 9.35 mm h^-1, respectively). This study provides a scatteringbased over-land rainfall retrieval algorithm for South Korea affected by stationary front rain and typhoons with the advantages of the previous PCT and SI methods to be applied to a variety of spaceborne passive microwave radiometers.     

Field Experiments and Numerical Modeling of Wind Speed and Surface Waves in Medium-size Inland Reservoirs

An attempt is made to apply the modern methods of surface wave simulation developed for oceanic conditions to the modeling of waves in medium-size inland reservoirs (10–100 km). The results of field measurements of wind speed and waves are described, and on their basis the parameterization C _D (U ₁₀) is proposed. WAVEWATCH III spectral wave model was adapted to the conditions of a medium-size inl and reservoir. The simulated data are compared with the field data. The use of the new parameterization C _D (U ₁₀) allowed reducing the values of the wind wave growth rate that improved consistency in data from the field experiment and numerical modeling concerning the height of significant waves. Further steps towards improving the quality of prediction of the adapted WAVEWATCH III model are discussed.     

Effects of a Fence on Pollutant Dispersion in a Boundary Layer Exposed to a Rural-to-Urban Transition

Simultaneous particle-image velocimetry and laser-induced fluorescence combined with large-eddy simulations are used to investigate the flow and pollutant dispersion behaviour in a rural-to-urban roughness transition. The urban roughness is characterized by an array of cubical obstacles in an aligned arrangement. A plane fence is added one obstacle height h upstream of the urban roughness elements, with three different fence heights considered. A smooth-wall turbulent boundary layer with a depth of 10h is used as the approaching flow, and a passive tracer is released from a uniform line source 1h upstream of the fence. A shear layer is formed at the top of the fence, which increases in strength for the higher fence cases, resulting in a deeper internal boundary layer (IBL). It is found that the mean flow for the rural-to-urban transition can be described by means of a mixing-length model provided that the transitional effects are accounted for. The mixing-length formulation for sparse urban canopies, as found in the literature, is extended to take into account the blockage effect in dense canopies. Additionally, the average mean concentration field is found to scale with the IBL depth and the bulk velocity in the IBL.     

Aerodynamic Roughness Length of Fresh Snow

This study presents the results from a series of wind-tunnel experiments designed to investigate the aerodynamic roughness length z ₀ of fresh snow under no-drift conditions. A two-component hot-film anemometer was employed to obtain vertical profiles of velocity statistics in a zero pressure gradient turbulent boundary layer for flow over naturally deposited snow surfaces. The roughness of these snow surfaces was measured by means of digital photography to capture characteristic length scales that can be related to z ₀. Our results show that, under aerodynamically rough conditions, the mean value of the roughness length for fresh snow is \({\langle{z}_{0}\rangle= 0.24}\) mm with a standard deviation σ(z ₀) = 0.05 mm. In this study, we show that variations in z ₀ are associated with variations in the roughness geometry. The roughness measurements suggest that the estimated values of z ₀ are consistent with the presence of irregular roughness structures that develop during snowfalls that mimic ballistic deposition processes.     

Influence of the atmospheric conditions on PM<Subscript>10</Subscript> concentrations in Poznań, Poland

This study investigates atmospheric conditions’ influence on the mean and extreme characteristics of PM₁₀ concentrations in Poznań during the period 2006–2013. A correlation analysis was carried out to identify the most important meteorological variables influencing the seasonal dynamics of PM₁₀ concentrations. The highest absolute correlation values were obtained for planetary boundary layer height (r = ?0.57), thermal (daily minimum air temperature: r = ?0.51), anemological (average daily wind speed: r = ?0.37), and pluvial (precipitation occurrence: r = ?0.36) conditions, however the highest correlations were observed for temporal autocorrelations (1 day lag: r = 0.70). As regulated by law, extreme events were identified on the basis of daily threshold value i.e. 50 μg m^?3. On average, annually there are approximately 71.3 days anywhere in the city when the threshold value is exceeded, 46.6 % of those occur in winter. Additionally, 83.7 % of these cases have been found to be continuous episodes of a few days, with the longest one persisting for 22 days. The analysis of the macro-scale circulation patterns led to the identification of an easy-to-perceive seasonal relations between atmospheric fields that favour the occurrence of high PM₁₀ concentration, as well as synoptic situations contributing to the rapid air quality improvement. The highest PM₁₀ concentrations are a clear reaction to a decrease in air temperature by over 3 °C, with simultaneous lowering of PBL height, mean wind speed (by around 1 m s^?1) and changing dominant wind directions from western to eastern sectors. In most cases, such a situation is related to the expansion of a high pressure system over eastern Europe and weakening of the Icelandic Low. Usually, air quality conditions improve along with an intensification of westerlies associated with the occurrence of low pressure systems over western and central Europe. Opposite relations are distinguishable in summer, when air quality deterioration is related to the inflow of tropical air masses originating over the Sahara desert.     

Analysis of photosynthetically active radiation and applied parameterization model for estimating of PAR in the North China Plain

As photosynthetically active radiation (PAR) variability and PAR estimating methods play an important role in climate change and ecological process research, PAR variation trends and broadband global solar radiation (R _s ) ratios (PAR/R _s ) in the North China Plain (NCP) are examined using in situ PAR and R _s observed data for 2005 to 2011. The annual average PAR value found in the NCP is 22.9 mol m^?2 d^?1. The highest and lowest values were recorded at Changwu and Luancheng sites, respectively. The highest PAR/R _s value was found in Jiaozhouwan due to large water vapor volumes present in this area. PAR/R _s levels have increased in the NCP due to a decrease in fine aerosols and increase in water vapor concentration. From these analysis results, a parameterization model that can be applied to all sky conditions was checked. Empirical estimation model comparisons for obtaining PAR values indicate that model was least accurate when R _s was used independently. When the model included R _s, the clearness index (K _s) and the solar zenith angle, the model estimated PAR values with acceptable accuracy. A parameterization model was constructed by considering K _s and attenuation factors of PAR under clear weather conditions (ρ _clear). The improved parameterization model more accurately predicts values for local sites and for various observation sites.     

Retrieving air humidity,global solar radiation,and reference evapotranspiration from daily temperatures: development and validation of new methods for Mexico. Part I: humidity

The occurrence of flood and drought frequency is highly correlated with the temporal fluctuations of streamflow series; understanding of these fluctuations is essential for the improved modeling and statistical prediction of extreme changes in river basins. In this study, the complexity of daily streamflow fluctuations was investigated by using multifractal detrended fluctuation analysis (MF-DFA) in a large heterogeneous lake basin, the Poyang Lake basin in China, and the potential impacts of human activities were also explored. Major results indicate that the multifractality of streamflow fluctuations shows significant regional characteristics. In the study catchment, all the daily streamflow series present a strong long-range correlation with Hurst exponents bigger than 0.8. The q-order Hurst exponent h(q) of all the hydrostations can be characterized well by only two parameters: a (0.354 ≤ a ≤ 0.384) and b (0.627 ≤ b ≤ 0.677), with no pronounced differences. Singularity spectrum analysis pointed out that small fluctuations play a dominant role in all daily streamflow series. Our research also revealed that both the correlation properties and the broad probability density function (PDF) of hydrological series can be responsible for the multifractality of streamflow series that depends on watershed areas. In addition, we emphasized the relationship between watershed area and the estimated multifractal parameters, such as the Hurst exponent and fitted parameters a and b from the q-order Hurst exponent h(q). However, the relationship between the width of the singularity spectrum (Δα) and watershed area is not clear. Further investigation revealed that increasing forest coverage and reservoir storage can effectively enhance the persistence of daily streamflow, decrease the hydrological complexity of large fluctuations, and increase the small fluctuations.     

Formulation of stability-dependent empirical relations for turbulent intensities from surface layer turbulence measurements for dispersion parameterization in a lagrangian particle dispersion model

Season- and stability-dependent turbulence intensity (σ _u /u _*, σ _v /u _*, σ _w /u _*) relationships are derived from experimental turbulence measurements following surface layer scaling and local stability at the tropical coastal site Kalpakkam, India for atmospheric dispersion parameterization. Turbulence wind components (u′, v′, w′) measured with fast response UltraSonic Anemometers during an intense observation campaign for wind field modeling called Round Robin Exercise are used to formulate the flux–profile relationships using surface layer similarity theory and Fast Fourier Transform technique. The new relationships (modified Hanna scheme) are incorporated in a Lagrangian Particle Dispersion model FLEXPART-WRF and tested by conducting simulations for a field tracer dispersion experiment at Kalpakkam. Plume dispersion analysis of a ground level hypothetical release indicated that the new turbulent intensity formulations provide slightly higher diffusivity across the plume relative to the original Hanna scheme. The new formulations for σ _u , σ _v , σ _w are found to give better agreement with observed turbulent intensities during both stable and unstable conditions under various seasonal meteorological conditions. The simulated concentrations using the two methods are compared with those obtained from a classical Gaussian model and the observed SF₆ concentration. It has been found that the new relationships provide comparatively higher diffusion across the plume relative to the model default Hanna scheme and provide downwind concentration results in better agreement with observations.     

Multivariate drought frequency estimation using copula method in Southwest China

This study reveals the impacts of climatic variable trends on drought severity in Xinjiang, China. Four drought indices, including the self-calibrating Palmer drought severity index (sc-PDSI), Erinç’s index (I _m), Sahin’s index (I _sh), and UNEP aridity index (AI), were used to compare drought severity. The ensemble empirical mode decomposition and the modified Mann-Kendall trend test were applied to analyze the nonlinear components and trends of the climatic variable and drought indices. Four and six climatic scenarios were generated in sc-PDSI, I _m, I _sh, and AI with different combinations of the observed and detrended climatic variables, respectively. In Xinjiang, generally increasing trends in minimal, average, and maximal air temperature (T _min, T _ave, T _max) and precipitation (P) were found, whereas a decreasing trend in wind speed at 2 m height (U ₂) was observed. There were significantly increasing trends in all of the four studied drought indices. Drought relief was more obvious in northern Xinjiang than in southern Xinjiang. The strong influences of increased P on drought relief and the weak influences of increased T _min, T _ave, and T _max on drought aggravation were shown by comparing four drought indices under different climate scenarios. Decreased U ₂ had a weak influence on drought, as shown by the AI in different climate scenarios. The weak influences of T and U ₂ were considered to be masked by the strong influences of P on droughts. Droughts were expected to be more severe if P did not increase, but were likely milder without an increase in air temperature and with a decrease in U ₂.     

Uncertainty of global summer precipitation in the CMIP5 models: a comparison between high-resolution and low-resolution models

Worldwide, the majority of rapidly growing neighborhoods are found in the Global South. They often exhibit different building construction and development patterns than the Global North, and urban climate research in many such neighborhoods has to date been sparse. This study presents local-scale observations of net radiation (Q _* ) and sensible heat flux (Q _H ) from a lightweight low-rise neighborhood in the desert climate of Andacollo, Chile, and compares observations with results from a process-based urban energy-balance model (TUF3D) and a local-scale empirical model (LUMPS) for a 14-day period in autumn 2009. This is a unique neighborhood-climate combination in the urban energy-balance literature, and results show good agreement between observations and models for Q _* and Q _H . The unmeasured latent heat flux (Q _E ) is modeled with an updated version of TUF3D and two versions of LUMPS (a forward and inverse application). Both LUMPS implementations predict slightly higher Q _E than TUF3D, which may indicate a bias in LUMPS parameters towards mid-latitude, non-desert climates. Overall, the energy balance is dominated by sensible and storage heat fluxes with mean daytime Bowen ratios of 2.57 (observed Q _H /LUMPS Q _E )–3.46 (TUF3D). Storage heat flux (ΔQ _S ) is modeled with TUF3D, the empirical objective hysteresis model (OHM), and the inverse LUMPS implementation. Agreement between models is generally good; the OHM-predicted diurnal cycle deviates somewhat relative to the other two models, likely because OHM coefficients are not specified for the roof and wall construction materials found in this neighborhood. New facet-scale and local-scale OHM coefficients are developed based on modeled ΔQ _S and observed Q _* . Coefficients in the empirical models OHM and LUMPS are derived from observations in primarily non-desert climates in European/North American neighborhoods and must be updated as measurements in lightweight low-rise (and other) neighborhoods in various climates become available.