Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed estimation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Distribution Function Instead of Steady-State Assumption in Time-Series Simulation

The quasi-steady-state assumption in time-series simulation is inadequate to model phenomena of interest to energy system integration such as inverter clipping, sell-back of excess electricity to a utility and utility hosting capacity, and battery throughput. Researchers are working on stochastic modeling, higher-resolution time series data, and machine learning approaches to address this need. This poster describes a distribution function that allows a maximum value, minimum value, and shape to the curve that can vary within each time-step based only on data inputs at the resolution of that time step (eg. Hourly). The distribution function is not used to synthesize high-resolution data, rather scalar integrals are used to calculate the quantities of interest within each time step. The form of the distribution function shows significant reduction in error when compared to 1 minute data and commercial software employing the distribution function (HOMER Pro version 14.3 and higher) shows a much improved estimate of inverter clipping in an example.

battery↗
DOE OSTI2020

Model-Free State Estimation Using Low-Rank Canonical Polyadic Decomposition

As electric grids experience high penetration levels of renewable generation, fundamental changes are required to address real-time situational awareness. Here, we utilize unique traits of tensors to devise a model-free situational awareness and energy forecasting framework for distribution networks. This work formulates the state of the network at multiple time instants as a three-way tensor; hence, recovering full state information of the network is tantamount to estimating all the values of the tensor. Given measurements received from µphasor measurement units and/or smart meters, the recovery of unobserved quantities is carried out using the low-rank canonical polyadic decomposition of the state tensor—that is, the state estimation task is posed as a tensor imputation problem utilizing observed patterns in measured (sampled) quantities. Two structured sampling schemes are considered, namely, asynchronous slab and fiber sampling. For both schemes, we present sufficient conditions on the number of sampled slabs and fibers that guarantee identifiability of the factors of the state tensor. Numerical results demonstrate the ability of the proposed framework to achieve high estimation accuracy in multiple sampling scenarios.

42 ENGINEERING↗
DOE OSTI2021

Towards robust autonomous impedance spectroscopy analysis: A calibrated hierarchical Bayesian approach for electrochemical impedance spectroscopy (EIS) inversion

Distribution-based analyses, such as the distribution of relaxation times (DRT) and the distribution of diffusion times (DDT), present model-free alternatives to equivalent circuit modeling for analysis of electrochemical impedance spectroscopy (EIS) data. However, reconstructing such distributions from noisy impedance data is an ill-posed problem that must be solved with specialized inversion algorithms, requiring careful control and tuning. Furthermore, most inversion algorithms developed to date can only solve problems of limited complexity. Herein, we present a new hierarchical Bayesian method for EIS inversion, leveraging efficient algorithms for optimization and Hamiltonian Monte Carlo (HMC) sampling to solve models of arbitrary complexity. We overcome the challenge of ad-hoc parameter tuning by encoding intrinsic characteristics of the DRT and DDT into flexible prior distributions and “pre-calibrating” the model to simulated data. This approach is versatile, highly robust to noise, and provides quantitative estimates of both the error structure of the data and the uncertainty in the recovered distributions. The model is validated with simulated data to demonstrate accurate recovery of the DRT and the DDT. The method also shows promise for simultaneous recovery of multiple distributions, raising the intriguing possibility of semi-autonomous EIS analysis and ad-hoc model construction. Finally, the practical utility of the method is illustrated with experimental data. Throughout, we draw comparisons to several recently published EIS inversion methodologies.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Baselining the Indirect Effect by Improving Quantification of Sea Spray and Marine Sources at Ascension Island (Final Report)

Oceans cover two-thirds of the Earth and understanding the interactions of aerosols with clouds in these large marine regions requires quantifying the man-made contributions to the budget of cloud-drop forming particles (known as cloud condensation nuclei, or CCN) relative to the non-manmade “baseline” conditions and understanding the meteorology of boundary layer clouds. Modeling studies have shown substantial uncertainties and sensitivities to natural marine CCN sources, meaning that to reduce uncertainties in indirect effects we must be able to better quantify the CCN budget in ocean regions. While models provide important constraints on these uncertainties, actually reducing uncertainties requires substantial observations in open-ocean and coastal regions in order to establish the baseline on which manmade emissions are added. The tropical South Atlantic Ocean is one of the least-sampled regions of the planet, making the comprehensive measurements of the Department of Energy Atmospheric Radiation Measurement Layered Atlantic Smoke Interactions with Clouds (LASIC) campaign provides the longest record of aerosol size distribution measurements from a differential mobility analyzer in a cloud-influenced marine location in the ARM database, with 17 months of ground-based aerosol size distribution measurements. This project addressed the research topic of Aerosol-Cloud Interactions (ACI) by supporting three publications from the LASIC measurements: 1. The ARM measurements were combined with a new value-added technique that was pioneered by the Russell group for quantifying sea salt. This quantification of sea salt uses supermicron scattering measurements for retrieving reasonable sea spray mass concentrations, providing the best-available, observationally-constrained estimate of the sea spray mode properties when supermicron size distribution measurements are not available. 2. The fitted modes of the distribution were used to investigate the signatures of cloud processing for very clean to very smoky aerosol conditions, revealing not only differences in the particles that activate in clouds but also in the mechanisms that control that droplet formation process. In clean air, the size required to form a cloud droplet is influenced by the number of particles, as well as how quickly particles take up water during growth in cloud. 3. Building on this aerosol characterization, cloud and meteorological conditions were used to evaluate aerosol-related changes in cloud albedo and optical depth. We introduced a new method of decomposing the impact of aerosols on clouds known as the Twomey effect by incorporating retrieved supersaturation from two independent sets of observations to constrain the feedback of aerosol particles on cloud properties. The method quantifies the reduction of the Twomey effect at high aerosol concentrations, which has never been explained quantitatively by observations. In addition, the results provide the first direct validation for the conditions observed in the tropical South Atlantic of a parcel-based approach that is embedded in many climate models. Together these findings illustrate how ARM extended field campaigns in stratocumulus-covered regions can be used to constrain ACI processes with direct observations, providing specific radiative effects without models. By making such process-specific constraints available to improve ACI in global climate models, ARM observations play a key role in supporting model development.

58 GEOSCIENCES↗
DOE OSTI2026

Quantitative Evaluation of Potassium Iodide Implementation Strategies for Emergency Preparedness and Response

This study evaluates the effectiveness of potassium iodide (KI) distribution strategies in mitigating exposure to radioiodine during severe nuclear power plant accidents. This analysis quantitatively compares various KI distribution methods (pre-distributed versus stockpiles), including scenarios with and without KI administration. The results indicate that differences in distribution strategies impact the projected thyroid dose by at least an order of magnitude. The results also indicate that the timing of KI administration is critical, as expected. For scenarios involving delayed releases of significant quantities of radionuclides, evacuation is the most effective protection strategy regardless of KI distribution method. For scenarios involving rapid releases, retrieving KI from stockpiles can have a detrimental effect. Pre-distributed KI is potentially the most effective approach when used as a supplement to evacuation and sheltering. However, these model results are based on idealized conditions for KI distribution and administration; the actual benefits of KI prophylaxis are likely to be less than estimated in this report due to many variables. The results highlight the importance of considering the cost and rigor of different distribution programs and public compliance with emergency instructions. This report includes suggested research to explore KI distribution plans for advanced reactors.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Disc dichotomy signature in the vertical distribution of [Mg/Fe] and the delayed gas infall scenario

Context. Analysis of the Apache Point Observatory Galactic Evolution Experiment project (APOGEE) data suggests the existence of a clear distinction between two sequences of disc stars in the [α/Fe] versus [Fe/H] abundance ratio space, known as the high- and low-α sequence, respectively. This dichotomy also emerges from an analysis of the vertical distribution of the [α/Fe] abundance ratio. Aims. We aim to test whether the revised two-infall chemical evolution models designed to reproduce the low- and high-α sequences in the [α/Fe] versus [Fe/H] ratios in the solar neighbourhood are also capable of predicting the disc bimodality observed in the vertical distribution of [Mg/Fe] in APOGEE DR16 data. Methods. Along with the chemical composition of the simple stellar populations born at different Galactic times predicted by our reference chemical evolution models in the solar vicinity, we provide their maximum vertical height above the Galactic plane |zmax| computed assuming the relation between the vertical action and stellar age in APOGEE thin-disc stars. Result. The vertical distribution of the [Mg/Fe] abundance ratio predicted by the reference chemical evolution models is in agreement with that observed when combining the APOGEE DR16 data (chemical abundances) with the astroNN catalogue (stellar ages, orbital parameters) for stars younger than 8 Gyr (only low-α sequence stars). Including the high-α disc component, the dichotomy in the vertical [Mg/Fe] abundance distribution is reproduced considering the observational cut in the Galactic height of |z|< 2 kpc. However, our model predicts an overly flat (almost constant) growth of the maximum vertical height |zmax| quantity as a function of [Mg/Fe] for high-α objects in contrast with the median values from APOGEE data. Possible explanations for such a tension are that: (i) the APOGEE sample with |z|< 2 kpc is more likely than ours to be contaminated by halo stars, causing the median values to be kinematically hotter, and (ii) external perturbations – such as minor mergers – that the Milky Way experienced in the past could have heated up the disc, and the heating of the orbits cannot be modeled by only scattering processes. Assuming a disc dissection based on chemistry for APOGEE-DR16 stars (|z|< 2 kpc), the observed |zmax| distributions for high-α and low-α sequences are in good agreement with our model predictions if we consider the errors in the vertical action estimates in the calculation. Moreover, a better agreement between predicted and observed stellar distributions at different Galactic vertical heights is achieved if asteroseismic ages are included as a constraint in the best-fit model calculations. Conclusions. The signature of a delayed gas infall episode, which gives rise to a hiatus in the star formation history of the Galaxy, are imprinted both in the [Mg/Fe] versus [Fe/H] relation and in vertical distribution of [Mg/Fe] abundances in the solar vicinity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Distribution System Model Calibration for GMLC 3.3.3 "Incipient Failure Identification for Common Grid Asset Classes" - Project Summary

Distribution system model calibration is a key enabling task for incipient failure identification within the distribution system. This report summarizes the work and publications by Sandia National Laboratories on the GMLC project titled “Incipient Failure Identification for Common Grid Asset Classes”. This project was a joint effort between Sandia National Laboratories, Lawrence Livermore National Laboratory, National Energy Technology Laboratory, and Oak Ridge National Laboratory. The included work covers distribution system topology identification, transformer groupings, phase identification, regulator and tap position estimation, and the open-source release and implementation of the developed algorithms.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Real Time Effective Impedance Estimation for Power System State Estimation

Here, we introduce two tools-the Effective Impedance State Estimator (EISE), and the Real Time Effective Impedance Estimator (RTEIE). The EISE is novel in that it is an entirely data-driven state estimator, and does not require a network model. Thus, the EISE is well suited for applications in which an accurate network model is not available. Instead of using a network impedance model, the EISE uses an effective impedance network model that is learned online by the RTEIE. The RTEIE recursively estimates the effective impedance network model using phasor measurement unit (PMU) measurements. The mathematical reasoning for using an effective impedance model, rather than a full network impedance model, is given. The recursive estimation framework allows the RTEIE to take the entire history of phasor measurements into account to produce a model estimate, while still adjusting to changes in the network impedance due to external factors such as topology changes and temperature. The merits of the RTEIE and EISE are demonstrated by comparing the accuracy of the EISE with a conventional state estimator on a toy network.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Characterization and Valuation of the Uncertainty of Calibrated Parameters in Microsimulation Decision Models

We evaluated the implications of different approaches to characterize the uncertainty of calibrated parameters of microsimulation decision models (DMs) and quantified the value of such uncertainty in decision making. We calibrated the natural history model of CRC to simulated epidemiological data with different degrees of uncertainty and obtained the joint posterior distribution of the parameters using a Bayesian approach. We conducted a probabilistic sensitivity analysis (PSA) on all the model parameters with different characterizations of the uncertainty of the calibrated parameters. We estimated the value of uncertainty of the various characterizations with a value of information analysis. We conducted all analyses using high-performance computing resources running the Extreme-scale Model Exploration with Swift (EMEWS) framework. The posterior distribution had a high correlation among some parameters. The parameters of the Weibull hazard function for the age of onset of adenomas had the highest posterior correlation of -0.958. When comparing full posterior distributions and the maximum-a-posteriori estimate of the calibrated parameters, there is little difference in the spread of the distribution of the CEA outcomes with a similar expected value of perfect information (EVPI) of $\$$653 and $\$$685, respectively, at a willingness-to-pay (WTP) threshold of $\$$66,000 per quality-adjusted life year (QALY). Ignoring correlation on the calibrated parameters’ posterior distribution produced the broadest distribution of CEA outcomes and the highest EVPI of $\$$809 at the same WTP threshold. Different characterizations of the uncertainty of calibrated parameters affect the expected value of eliminating parametric uncertainty on the CEA. Ignoring inherent correlation among calibrated parameters on a PSA overestimates the value of uncertainty.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Statistical estimates of the binary properties of rotational variables

ABSTRACT We present a model to estimate the average primary masses, companion mass ranges, the inclination limit for recognizing a rotational variable, and the primary mass spreads for populations of binary stars. The model fits a population’s binary mass function distribution and allows for a probability that some mass functions are incorrectly estimated. Using tests with synthetic data, we assess the model’s sensitivity to each parameter, finding that we are most sensitive to the average primary mass and the minimum companion mass, with less sensitivity to the inclination limit and little to no sensitivity to the primary mass spread. We apply the model to five populations of binary spotted rotational variables identified in ASAS-SN, computing their binary mass functions using RV data from APOGEE. Their average primary mass estimates are consistent with our expectations based on their CMD locations ($\sim 0.75 \, {\rm M}_{\odot }$ for lower main sequence primaries and $\sim 0.9$–$1.2 \, {\rm M}_{\odot }$ for RS CVn and sub-subgiants). Their companion mass range estimates allow companion masses down to $M_2/M_1\simeq 0.1$, although the main sequence population may have a higher minimum mass fraction ($\sim 0.4$). We see weak evidence of an inclination limit $\gtrsim 50^{\circ }$ for the main sequence and sub-subgiant groups and no evidence of an inclination limit in the other groups. No groups show strong evidence for a preferred primary mass spread. We conclude by demonstrating that the approach will provide significantly better estimates of the primary mass and the minimum mass ratio and reasonable sensitivity to the inclination limit with 10 times as many systems.

Phillips, Anya (ORCID:000900051914974X)↗
DOE OSTI2023

Winter connectivity and leapfrog migration in a migratory passerine

Abstract Technological advances in migratory tracking tools have revealed a remarkable diversity in migratory patterns. One such pattern is leapfrog migration, where individuals that breed further north migrate to locations further south. Here, we analyzed migration patterns in the Painted Bunting ( Passerina ciris ) using a genetic‐based approach. We started by mapping patterns of genetic variation across geographic space (called a genoscape) using 386 individuals from 25 populations across the breeding range. We then genotyped an additional 230 samples from 31 migration stopover locations and 178 samples from 16 wintering locations to map patterns of migratory connectivity. Our analyses of genetic variation across the breeding range show the existence of four genetically distinct groups within the species: Eastern, Southwestern, Louisiana, and Central groups. Subsequent assignment of migrating and wintering birds to genetic groups illustrated that birds from the Central group migrated during the fall via western Mexico or southern Texas, spent the winter from northeastern Mexico to Panama, and migrated north via the Gulf Coast of Mexico. While Louisiana birds overlapped with Central birds on their spring migratory routes along the Gulf Coast, we found that Louisiana birds had a more restricted wintering distribution in the Yucatan Peninsula and Central America. Further estimation of the straight‐line distance from the predicted breeding location to the wintering location revealed that individuals sampled at lower winter latitudes traveled to greater distances (i.e., the predicted breeding area was further north; p > .001), confirming that these species exhibit a leapfrog migration pattern. Overall, these results demonstrate the utility of a genoscape‐based approach for identifying range‐wide patterns of migratory connectivity such as leapfrog migration with a high degree of clarity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Scalable statistical inference of photometric redshift via data subsampling

Handling big data has largely been a major bottleneck in traditional statistical models. Consequently, when accurate point prediction is the primary target, machine learning models are often preferred over their statistical counterparts for bigger problems. But full probabilistic statistical models often outperform other models in quantifying uncertainties associated with model predictions. We develop a data-driven statistical modeling framework that combines the uncertainties from an ensemble of statistical models learned on smaller subsets of data carefully chosen to account for imbalances in the input space. We demonstrate this method on a photometric redshift estimation problem in cosmology, which seeks to infer a distribution of the redshift—the stretching effect in observing the light of far-away galaxies—given multivariate color information observed for an object in the sky. Our proposed method performs balanced partitioning, graph-based data subsampling across the partitions, and training of an ensemble of Gaussian process models.

data subsampling↗
DOE OSTI2022

Micromechanics of Fracture Propagation During Multistage Stress Relaxation and Creep in Brittle Rocks

Time-dependent rock deformation caused by the initiation and growth of fractures leads to the weakening of the rock mass. Understanding the fracturing mechanisms involved in the time-dependent behavior in brittle rocks is very important and to achieve this goal, a systematic series of three types of experiments was performed on double-flawed prismatic Barre granite specimens under unconfined compression. The first series aimed to identify the failure mechanism in the short-term failure mode under monotonic loading, the, second series involved multistage relaxation (constant strain) experiments to analyze the damage at different strain levels, and the third series explored the fracture propagation under multistage creep (constant load) experiments. The spatial and temporal evolution of cracking mechanisms were evaluated using the acoustic emission (AE) and two-dimensional digital image correlation (2D-DIC) techniques to observe the whole crack growth process as well as the accumulated inelastic strain at the specified region of interest. Results suggest that in the case of multistage creep experiments, the time to failure was less compared to the multistage relaxation, when loaded above the crack damage threshold (CD) estimated from the monotonic testing. Furthermore, the frequency magnitude distribution of the AE events generated in the three loading conditions followed the Gutenberg Richter model. A relatively lower b-value was obtained for the creep experiments, indicative of high energy AE events and faster crack growth. In addition, the AE and DIC results also revealed high evolution of tensile cracks at-different stages of creep and relaxation compared to shear and mixed-mode cracks.

58 GEOSCIENCES↗
DOE OSTI2023

Empirical evidence that mean soil contaminant concentration is not a conservative indicator of external exposure to wildlife

Adequate information on external exposure is often a limiting factor in ecological risk assessments because of a lack of exposure data for free-ranging animals as they interact with a heterogeneous distribution of contaminants. Instead, most wildlife exposure data are estimated from model simulations. However, empirical exposure data are rarely available to verify model results, or to test fundamental assumptions inherent in exposure assessments. We used novel, GPS-coupled radiation dosimeters on reindeer (Rangifer tarandus) to quantify their external exposure for five months in a Chernobyl fallout affected mountainous area of Norway. The exposure data were then used to test a fundamental assumption in screening-level risk assessments: Mean soil contaminant concentrations can be used to conservatively estimate chronic external exposures to representative animals. External exposure and GPS coordinates were obtained hourly and compared to model simulations of exposure based on airborne surveys of 137 Cs soil contamination. Our objectives were to determine 1) if modelled dose estimates based on mean 137 Cs soil concentrations were conservative, and 2) if modelled dose estimates could be improved by adjusting for spatial–temporal heterogeneity of contaminants using 137 Cs concentrations within reindeers’ GPS-defined home ranges and core areas. We found that external exposure modelled from mean soil contaminant concentrations underestimated exposure by 70% compared to empirical measurements from GPS-dosimeter collars worn on individual animals. Using area-weighted mean contaminant levels within the animals’ home range improved the model estimates, but still underestimated the field dosimeter data by 53%. The best estimate of external exposure was obtained when high-use areas within their home range, along with spatial heterogeneity of the contaminant, was accounted for by using area-weighted mean soil contaminant levels within GPS-derived core areas. Accounting for temporal use of core areas within animals’ home range provided a 24% improvement. Here, a comparison of our analyses to similar work performed on Chernobyl wolves found the two studies concur that risk assessors should be cautious in relying on the paradigm of using averaged soil contaminant levels to provide conservative estimates of external exposure to representative individuals. Results of modelled external exposure based on mean soil contaminant concentrations could misguide assessors and cause inappropriate management decisions.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2023

Shedding light on the economic costs of long-duration power outages: A review of resilience assessment methods and strategies

Here this paper provides a literature review of methods and modeling techniques to estimate the cost of power system outages, along with the value of outage mitigation or system resilience. Regulators, policymakers, and infrastructure owners have a growing need to understand the methods for estimating the benefits of resilience improvements of electric infrastructure against natural and man-made disasters. There is a broad literature that estimates the cost of short-duration outages and a small but developing literature on estimating the cost of long-duration outages. This article reviews the models used to estimate the cost of outages and discusses their relative strengths. Additionally, this paper identifies key questions from stakeholders regarding resilience investment and maps them to the relevant models that would help answer them. We include recommendations for future work to include recent advances in regional economic modeling that can estimate region and demographic-specific costs and the distributional consequences of potential resilience projects.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Range-wide population assessments for subalpine fir indicate widespread disturbance-driven decline

Subalpine forests in western North America are threatened by rapid climate change, increased activity by endemic and exotic insects and diseases, and changing wildfire regimes. The interactive effects of these stressors have resulted in pronounced population declines in many subalpine tree species; however, a systematic assessment of the status and trends of subalpine forests is lacking. Subalpine fir (Abies lasiocarpa) is a widespread species across the western United States, with documented population declines in many parts of its distribution. Here we use subalpine fir as an initial leverage point to build a more complete understanding of subalpine forest baseline conditions and responses to environmental change. Specifically, we leverage the USDA Forest Service Forest Inventory and Analysis (FIA) database to (1) ask how subalpine fir populations are changing across the species’ distribution in the western US, (2) assess the drivers of recent subalpine fir population trends, and (3) explore whether those changes imply generalized species-wide and/or system-wide decline. We found that subalpine fir abundance and basal area are declining concurrently across ~ 62% of the species’ distribution, and increasing across ~ 19%. Range-wide, we estimated 25.02 ± 2.74 % subalpine fir mortality between 2000 and 2009 and 2010–2019 FIA inventory periods, with higher mortality concentrated in the eastern Oregon Cascades, central Idaho, and parts of southern Colorado. High regeneration density did not predict positive population trajectories, which were instead associated with higher rates of adult recruitment. While the importance of different mortality agents varied substantially between ecoregions, 83.4% of total range-wide mortality was related to fire or biological disturbance. Declining subalpine fir basal area coincided with declines in the basal area of other co-occurring tree species in 39% of subalpine forest area, and with increases in conspecific basal area in 22% of forest area. Fire disturbance was the single largest cause of subalpine fir mortality; however, even where subalpine fir fire mortality was high, mortality among other species was primarily caused by insects. In conclusion, our results suggest that subalpine fir declines across large portions of the western United States are driven by forest disturbance, and that declines in subalpine fir populations may be indicative of negative change in subalpine forest systems broadly.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Comparing land surface temperature and mean radiant temperature for urban heat mapping in Philadelphia

Many cities are experiencing more frequent extreme heat during hot summers. With the rise of global temperature, the thermal comfort in urban areas become even worse. Quantitative information of the spatial distributions of urban heat has become increasingly important for resilience and adaptation to climate change in cities. This study compares satellite-derived land surface temperature (LST) and urban microclimate modeling-based mean radiant temperature (Tmrt) for mapping the urban heat distributions in Philadelphia, Pennsylvania, USA. Here, the LST was estimated based on Landsat 8 thermal imagery with a spatial resolution of around 100 m, while the Tmrt was simulated based on high resolution LiDAR and national aerial imagery program multispectral aerial imageries with a spatial resolution of 1 m. Result shows that both LST and Tmrt show a similar general pattern of the urban heat across the study area, while the Tmrt presents much more details of the heat variations street by street and neighborhood by neighborhood. The LST tends to have a stronger relationship with the Tmrt on building roofs, which are usually not the place for human activities. This studyprovides evidence for choosing more appropriate metrics in urban heat-related studies.

54 ENVIRONMENTAL SCIENCES↗