Search NASA⌕ Search

SEARCH · Search NASA

Results for “data harmonization”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

A Critical Examination of Wind-Wave Spectral Functional Form

Traditionally, data from random ocean waves are presented in spectral functions. The spectra are the result of Fourier analysis. Fourier spectral analysis has dominated data analysis for, at least, the last hundred years. It has been the standard method for is examining the global amplitude-frequency distributions. Although Fourier transform valid under extremely general conditions, there are some crucial restrictions for the Fourier spectral analysis. The system must be linear, and the data must be stationary- otherwise, the resulting spectrum will make little physical sense. The stationarity requirement is also a common required criterion for most of other available data analysis methods. Nevertheless, few, if any, natural phenomena are linear and stationary. To compound these complications is the imperfection of our probes or numerical schemes the interactions of the imperfect probes even with a perfect linear system can make the final data nonlinear. Furthermore, all the available data are usually of finite duration. Under these conditions, Fourier analysis is of limited use, For lack of alternatives, however, Fourier analysis is still used to process such data. The loose application of Fourier analysis and the insouciant adoption of the stationary and linear assumptions may lead to misleading conclusions. Ocean waves are know to be nonlinear, and the wind system generating the wave field are seldom stationary- As a result, the traditional examination of the spectral form hardly made physical sense. A new method for analyzing nonlinear and nonstationary data has been developed. The key part is the Empirical Mode Decomposition (EMD) method with which any complicated data set can be decomposed into a finite and often small number of Intrinsic Mode Functions (IMF) that serve as the basis of the representation of the data, This decomposition method is adaptive, and, therefore, highly efficient. The IMFs admit well-behaved Hilbert transforms, and yield instantaneous energy and frequency as functions of time that give sharp identifications of imbedded structures. The final presentation of the results is an energy-frequency-time distribution, designated as the Hilbert Spectrum. Among the main conceptual innovations is the introduction of the instantaneous frequencies for complicated data sets, which eliminate the need of spurious harmonics to represent nonlinear and nonstationary signals. The spectral forms of the ocean waves are very different. This new method offers the first physical meaningful view of the wave spectrum. Data from laboratory and field will be presented to illustrate the differences.

Huang, Norden E.↗
DOE OSTI2025

Assessing methods in fusion and fitting for time series construction in remote sensing-based earth observations

This study evaluates the comparative performance of spatiotemporal fusion and time-series fitting methods for constructing high-spatiotemporal-resolution remote sensing time-series data. Due to in-class similarity of fusion methods and fitting methods, we employ the Fit-FC (Fitting, spatial Filtering, and residual Compensation) model as a representative fusion method and the linear harmonic fitting model as a representative fitting method. Both Fit-FC and the linear harmonic fitting are widely used for high-spatiotemporal-resolution time-series data construction, and we modify the original Fit-FC model to enable automatic time-series fusion. To ensure data representativeness, we use 3 years (2019–2021) of Harmonized Landsat and Sentinel-2 surface reflectance datasets and Terra MCD43A4 products. Eight experimental regions are selected worldwide to guarantee generalization of the comparative performance between fusion and fitting methods, covering diverse land-use types (cropland, developed land, forest, and grassland) and varying climatological conditions. Time-series of NDVI and surface reflectance are analyzed under both actual observations and simulated data-missing scenarios. The constructed time-series data reveals that (1) the modified Fit-FC and linear harmonic fitting model achieve excellent performance in constructing high-resolution time-series images; (2) the fusion method outperforms the fitting method in constructing time-series of NDVI and surface reflectance images in cropland-, forest-, and grassland-dominated regions; (3) both methods achieve comparable performance in developed-dominated regions; (4) the fusion method is more robust to missing data, and better captures abrupt phenological transitions under conditions of continuous missing data; (5) the fitting method is computationally more efficient, making it suitable for large-scale time-series image reconstruction. This study provides valuable insights for selecting optimal strategies to generate high-resolution time-series images across diverse application scenarios and lays a foundation for extensions to other vegetation indices or land surface variables.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1979

M2 ocean tide parameters and the deceleration of the moon's mean longitude from satellite orbit data

An estimation is made of the principal long-period spherical harmonic parameters in the representation of the M2 ocean tide from the orbital histories of the three satellites 1967-92A, Starlette, and GEOS 3. The data used are primarily the evolution of the orbital inclinations of the satellites in conjunction with the longitude of the ascending node from GEOS 3. Analysis procedure and analytic formulation, as well as ocean tidal parameter estimation and deceleration of the lunar mean longitude are outlined. The credibility of the M2 ocean tide solution is further enhanced by the close accord between the computed value for the deceleration of the lunar mean longitude and other recently reported estimates. It is evident from the results presented that studies of close earth satellite orbits are able to provide important information about the tidal forces acting on the earth.

Felsentreger, T. L.↗
NASA NTRS1993

Empirical wind model for the middle and lower atmosphere. Part 1: Local time average

The HWM90 thermospheric wind model was revised in the lower thermosphere and extended into the mesosphere and lower atmosphere to provide a single analytic model for calculating zonal and meridional wind profiles representative of the climatological average for various geophysical conditions. Gradient winds from CIRA-86 plus rocket soundings, incoherent scatter radar, MF radar, and meteor radar provide the data base and are supplemented by previous data driven model summaries. Low-order spherical harmonics and Fourier series are used to describe the major variations throughout the atmosphere including latitude, annual, semiannual, and longitude (stationary wave 1). The model represents a smoothed compromise between the data sources. Although agreement between various data sources is generally good, some systematic differences are noted, particularly near the mesopause. Root mean square differences between data and model are on the order of 15 m/s in the mesosphere and 10 m/s in the stratosphere for zonal wind, and 10 m/s and 4 m/s, respectively, for meridional wind.

Hedin, A. E.↗
NASA NTRS1999

Probing Ganymede's Magnetosphere with Field Line Resonances

We report on the spectrum of field line resonances identified in data acquired by the Galileo spacecraft within Ganymede's magnetosphere on a relatively low latitude pass. We infer properties of the plasma distribution and its transport from the observed spectrum. The harmonic structure in the spectrum of the magnetometer data agrees very well with the frequencies predicted for resonances of a dipole field. The spectrum implies a density of 2 amu per cubic centimeter near the equator on closed field lines of Ganymede's magnetic field inside of 2 RG (Ganymede radii). This density is significantly reduced relative to the local density of the Jovian plasma sheet near Ganymede (less than or equal to 8 electrons per cubic centimeter, or approximately equal to 100 amu per cubic centimeter for an average ion mass per charge of 20 amu and average charge of 1.5 electron charges). A shadowing effect of Ganymede for the flow of particles injected at a reconnection layer on the side of the moon downstream relative to the direction of torus plasma flow accounts for the marked density depletion. Implications for conducting paths near Ganymede's surface are considered.

Volwerk, M.↗
NASA NTRS1992

Annual and interannual variations of earth-emitted radiation based on a 10-year data set

The method of empirical orthogonal functions (EOFs) was applied to a 10-yr data set of outgoing longwave radiation. Spherical harmonic functions are used as a basis set for producing equal area map results. The first EOF accounts for 66 percent of the variance. After that, each EOF accounts for only a small variance, forming a slowly converging series. The first two EOFs describe mainly the annual cycle. The third EOF is primarily the semiannual cycle although many other EOFs also contain significant semiannual parts. These results reaffirm those based on a shorter data set. In addition, a much stronger spring/fall mode was found in the central equatorial Pacific Ocean for the second EOF than was found earlier. The fourth EOF describes much of the 1976 to 1977 and 1982 to 1983 ENSO phenomena. There is typically a gap in the spectrum between a semiannual peak and the annual cycle for all but the first EOF. A semiannual OLR dipole straddles the Asian-Australian monsoon track.

Bess, T. D.↗
NASA NTRS1986

DE 1 observations of harmonic auroral kilometric radiation

The plasma wave instrument on board the DE 1 spacecraft has observed several intervals of auroral kilometric radiation during which harmonic structure is clearly present. The paper presents evidence, some of which is based on unique capabilities of the DE instrument, which argues strongly that the harmonic structures are natural rather than instrumental in origin. The harmonic emissions occur infrequently, but when present may persist for intervals of up to an hour. The emissions are relatively narrow band and consist of a relatively weak fundamental accompanied by an even weaker second harmonic. The ratio of power in the fundamental band to the power in the harmonic ranges from 10 to 100. In all cases, polarization data indicate that the fundamental is a left-hand ordinary (L-O) mode emission while the harmonic is right-hand extraordinary (R-X) mode. These observations are consistent with predictions based on the cyclotron maser mechanism.

Mellott, M. M.↗
NASA NTRS2001

On the Application of Hilbert Spectral Analysis for Climate Studies

The Hilbert spectral analysis (Huang et al, 1998, Proceedings of the Royal Society of London, A 454, pp 903-995) consisted of two steps: First, the data has to be reduced into a finite number of Intrinsic Mode Function by the Empirical Mode Decomposition method, then the resulting Intrinsic Mode Functions are converted to time-frequency-energy distribution through Hilbert transform. In this approach, the Empirical Mode Functions served as the basis functions with which the data is expanded. This basis function is adaptive, and the decomposition is nonlinear. Furthermore, as the Hilbert transform is a singular transform, it retains a high degree of local information. The instantaneous frequency is determined by differentiation of the phase function; therefore, there is no restriction of the 'uncertainty principle' for all the time-frequency analysis resulting from a priori basis approach. With the adaptive basis and the instantaneous frequency, the Hilbert Spectral analysis can represent data from nonlinear and nonstationary processes without resorting to the harmonics. Another advantage of using instantaneous frequency is the ability to find out frequency from limited length of data, which is a critical problem in climate studies. As the processes driving the climate changes could be both nonlinear and nonstationary, the Hilbert Spectral Analysis could be of great use in examining the underlying mechanisms. A preliminary study based on the length of day data will be presented as example for the application of the Hilbert Spectral Analysis for climate study.

Huang, Norden E.↗
NASA NTRS1989

Annual and interannual variations of Earth-emitted radiation based on a 10-year data set

The method of empirical orthogonal functions (EOF) was applied to a 10-year data set of outgoing longwave radiation. Spherical harmonic functions are used as a basis set for producing equal area map results. The following findings are noted. The first EOF accounts for 66 percent of the variance. After that, each EOF accounts for only a small variance, forming a slowly converging series. The first two EOF's describe mainly the annual cycle. The third EOF is primarily the semiannual cycle although many other EOF's also contain significant semiannual parts. These results reaffirm those based on a shorter data set. In addition, a much stronger spring/fall mode was found in the central equatorial Pacific Ocean for the second EOF than was found earlier. This difference is attributed to the use of broadband radiometer data which were available for the present study. The earlier study used data from a window channel instrument which is not as sensitive to water vapor variations. The fourth EOF describes much of the 1976 to 1977 and 1982 to 1983 ENSO phenomena. There is typically a gap in the spectrum between a semiannual peak and the annual cycle for all but the first EOF. A semiannual OLR dipole straddles the Asian-Australian monsoon track.

Bess, T. Dale↗
NASA NTRS1983

Status of the geopotential

Satellite laser ranging, satellite altimetry, and improved measurements of surface gravitational anomalies have broadened the data base on intermediate and short wavelength regions of the earth gravity field. The global data set served to develop new geopotential models with a resolution in spherical harmonics out to degree 180. The resolution was made possible using Seasat altimetry data containing 56,761 values of 1 x 1 deg gravity anomalies. Satellite-to-satellite tracking techniques involving the Geos-3 and Apollo spacecraft data for the sea surface temperature have yielded accurate intermediate wavelength gravity variations which correlate well with residual depth anomalies. Oceanic gravity anomalies have been computed directly from satellite altimetry or through statistical estimation using oceanic geoid heights. The data sets for gravimetric geoids have been compared with altimetric surfaces to identify areas which were of interest for geophysical investigation. Future data sets could become available from a proposed satellite-to-satellite Doppler tracking system (Gravsat) launched by NASA.

Lerch, F. J.↗
NASA NTRS1982

Math modeling for helicopter simulation of low speed, low altitude and steeply descending flight

A math model was formulated to represent some of the aerodynamic effects of low speed, low altitude, and steeply descending flight. The formulation is intended to be consistent with the single rotor real time simulation model at NASA Ames Research Center. The effect of low speed, low altitude flight on main rotor downwash was obtained by assuming a uniform plus first harmonic inflow model and then by using wind tunnel data in the form of hub loads to solve for the inflow coefficients. The result was a set of tables for steady and first harmonic inflow coefficients as functions of ground proximity, angle of attack, and airspeed. The aerodynamics associated with steep descending flight in the vortex ring state were modeled by replacing the steady induced downwash derived from momentum theory with an experimentally derived value and by including a thrust fluctuations effect due to vortex shedding. Tables of the induced downwash and the magnitude of the thrust fluctuations were created as functions of angle of attack and airspeed.

Sheridan, P. F.↗
NASA NTRS1994

TOPEX/POSEIDON tides estimated using a global inverse model

Altimetric data from the TOPEX/POSEIDON mission will be used for studies of global ocean circulation and marine geophysics. However, it is first necessary to remove the ocean tides, which are aliased in the raw data. The tides are constrained by the two distinct types of information: the hydrodynamic equations which the tidal fields of elevations and velocities must satisfy, and direct observational data from tide gauges and satellite altimetry. Here we develop and apply a generalized inverse method, which allows us to combine rationally all of this information into global tidal fields best fitting both the data and the dynamics, in a least squares sense. The resulting inverse solution is a sum of the direct solution to the astronomically forced Laplace tidal equations and a linear combination of the representers for the data functionals. The representer functions (one for each datum) are determined by the dynamical equations, and by our prior estimates of the statistics or errors in these equations. Our major task is a direct numerical calculation of these representers. This task is computationally intensive, but well suited to massively parallel processing. By calculating the representers we reduce the full (infinite dimensional) problem to a relatively low-dimensional problem at the outset, allowing full control over the conditioning and hence the stability of the inverse solution. With the representers calculated we can easily update our model as additional TOPEX/POSEIDON data become available. As an initial illustration we invert harmonic constants from a set of 80 open-ocean tide gauges. We then present a practical scheme for direct inversion of TOPEX/POSEIDON crossover data. We apply this method to 38 cycles of geophysical data records (GDR) data, computing preliminary global estimates of the four principal tidal constituents, M(sub 2), S(sub 2), K(sub 1) and O(sub 1). The inverse solution yields tidal fields which are simultaneously smoother, and in better agreement with altimetric and ground truth data, than previously proposed tidal models. Relative to the 'default' tidal corrections provided with the TOPEX/POSEIDON GDR, the inverse solution reduces crossover difference variances significantly (approximately 20-30%), even though only a small number of free parameters (approximately equal to 1000) are actually fit to the crossover data.

Egbert, Gary D.↗
DOE OSTI2024

HarDWR - Harmonized Water Rights Records

For a detailed description of the database of which this record is only one part, please see the HarDWR meta-record. Here we present a new dataset of western U.S. water rights records. This dataset provides consistent unique identifiers for each spatial unit of water management across the domain, unique identifiers for each water right record, and a consistent categorization scheme that puts each water right record into one of 7 broad use categories. These data were instrumental in conducting a study of the multi-sector dynamics of intersectoral water allocation changes through water markets (Grogan et al., in review). Specifically, the data were formatted for use as input to a process-based hydrologic model, WBM, with a water rights module (Grogan et al., in review). While this specific study motivated the development of the database presented here, U.S. west water management is a rich area of study (e.g., Anderson and Woosly, 2005; Tidwell, 2014; Null and Prudencio, 2016; Carney et al, 2021) so releasing this database publicly with documentation and usage notes will enable other researchers to do further work on water management in the U.S. west. The raw downloaded data for each state is described in Lisk et al. (in review), as well as here. The dataset is a series of various files organized by state sub-directories. The first two characters of each file name is the abbreviation for the state the in which the file contains data for. After the abbreviation is the text which describes the contents of the file. Here is each file type described in detail: XXFullHarmonizedRights.csv: A file of the combined groundwater and surface water records for each state. Essentially, this file is the merging of XXGroundwaterHarmonizedRights.csv and XXSurfaceWaterHarmonizedRights.csv by state. The column headers for each of this type of file are: state - The name of the state the data comes from. FIPS - The two-digit numeric state ID code. waterRightID - The unique identifying ID of the water right, the same identifier as its state uses. priorityDate - The priority date associated with the right. origWaterUse - The original stated water use(s) from the state. waterUse - The water use category under the unified use categories established here. source - Whether the right is for surface water or groundwater. basinNum - The alpha-numeric identifier of the WMA the record belongs to. CFS - The maximum flow of the allocation in cubic feet per second (ft3s-1). Arizona is unique among the states, as its surface and groundwater resources are managed with two different sets of boundaries. So, for Arizona, the basinNum column is missing and instead there are two columns: surBasinNum - The alpha-numeric identifier of the surface water WMA the record belongs to. grdBasinNum - The alpha-numeric identifier of the groundwater WMA the record belongs to. XXStatePOD.shp: A shapefile which identifies the location of the Points of Diversion for the state's water rights. It should be noted that not all water right records in XXFullHarmonizedRights.csv have coordinates, and therefore may be missing from this file. XXStatePOU.shp: A shapefile which contains the area(s) in which each water right is claimed to be used. Currently, only Idaho and Washington provided valid data to include within this file. XXGroundwaterHarmonizedRights.csv: A file which contains only harmonized groundwater rights collected from each state. See XXFullHarmonizedRights.csv for more details on how the data is formatted. XXSurfaceWaterHarmonizedRights.csv: A file which contains only harmonized surface water rights collected from each state. See XXFullHarmonizedRights.csv for more details on how the data is formatted. Additionally, one file, stateWMALabels.csv, is not stored within a sub-directory. While we have referred to the spatial boundaries that each state uses to manage its water resources as WMAs, this term is not shared across all states. This file lists the proper name for each boundary set, by state. For those whom may be interested in exploring our code more in depth, we are also making available an internal data file for convenience. The file is in .RData format and contains everything described above as well as some minor additional objects used within the code calculating the cumulative curves. For completeness, here is a detailed description of the various objects which can be found within the .RData file: states: A character vector containing the state names for those states in which data was collected for. More importantly, the index of the state name is also the index in which that state's data can be found in the various following list objects. For example, if California is the third index in this object, the data for California will also be in the third index for each accompanying list. rightsByState_ground: A list of data frames with the cleaned ground water rights collected from each state. This object holds the the data that is exported to created the xxGroundwaterHarmonizedRights.csv files. rightsByState_surface: A list of data frames with the cleaned surface water rights collected from each state. This object holds the the data that is exported to created the xxSurfaceWaterHarmonizedRights.csv files. fullRightsRecs: A list of the combined groundwater and surface water records for each state. This object holds the the data that is exported to created the xxFullHarmonizedRights.csv files. projProj: The spatial projection used for map creation in the beginning of the project. Specifically, the World Geodetic System (WGS84) as a coordinate reference system (CRS) string in PROJ.4 format. wmaStateLabel: The name and/or abbreviation for what each state legally calls their WMAs. h2oUseByState: A list of spatial polygon data frames which contain the area(s) in which each water right is claimed to be used. It should be noted that not all water right records have a listed area(s) of use in this object. Currently, only Idaho and Washington provided valid data to be included in this object. h2oDivByState: A list of spatial points data frames which identifies the location of the Point of Diversion for the state's water rights. It should be noted that not all water right records have a listed Point of Diversion in this object. spatialWMAByState: A list of spatial polygon data frames which contain the spatial WMA boundaries for each state. The only data contained within the table are identifiers for each polygon. It is worth reiterating that Arizona is the only state in which the surface and groundwater WMA boundaries are not the same. wmaIDByState: A list which contains the unique ID values of the WMAs for each state. plottingDim: A character vector used to inform mapping functions for internal map making. Each state is classified as either "tall" or "wide", to maximize space on a typical 8x11 page. The code related to the creation of this dataset can be viewed within HarDWR GitHub Repository/dataHarmonization.

Economics↗
NASA NTRS1975

Error enhancement in geomagnetic models derived from scalar data

An investigation conducted by Backus (1970) regarding the possible existence of two harmonic functions of certain characteristics in three-dimensional space is considered. The derivation of a model of the main geomagnetic field from scalar data is discussed along with a numerical simulation study. It is found that experimental discrepancies between vector field observations and the predictions of the model may have a mathematical origin, related to the work of Backus.

Stern, D. P.↗
NASA NTRS1972

Gravitational field models for the earth (GEM 1 and 2)

Two models of the earth's gravitational field have been computed at Goddard Space Flight Center. The first, Goddard Earth Model 1 (GEM 1), has been derived from satellite tracking data. The second, Goddard Earth Model 2 (GEM 2), has been derived from a combination of satellite tracking and surface gravimetric data. The geopotential models are represented in spherical harmonics complete to degree and order 16 for the combined solution and complete to degree and order 12 for the satellite solution. Both solutions include zonal terms to degree 21 and related satellite resonant coefficients to degree 22. The satellite data consisted primarily of optical data processed on 300 weekly orbital arcs for 25 close earth satellites. Surface gravity data were employed in the form of 5 deg x 5 deg mean free-air gravity anomalies providing about 70% world coverage. Station locations were obtained for 46 tracking sites by combining electronic, laser, and additional optical tracking data with the above satellite data. Analysis of the radial positions of these stations and a value of mean gravity on the geoid indicated a mean equatorial radius for the earth of about 6378145 meters. Results of geopotential tests on satellite data not used in the solution show that better agreement was obtained with the GEM 1 and GEM 2 models than with the 1969 Smithsonian Standard Earth 2 model.

Lerch, F. J.↗
NASA NTRS1981

Deconvolution of wide field-of-view radiometer measurements of earth-emitted radiation. II - Analysis of first year of Nimbus 6 ERB data

The theory of deconvolution considered by Smith and Green (1981) is applied to Nimbus 6 data in order to interpret the data with the maximum obtainable resolution. Attention is given to the data sample, sun-contaminated measurements, the measurement model, the deconvolution, the degree variance, the spherical harmonic coefficients, the geographical distribution of longwave radiation, time histories of zonal coefficients, and the effect of a grid system. Degree variance plots for 12 months of longwave radiation data show that the limit for a spherical harmonic representation of the Nimbus 6 wide field-of-view longwave data is 12th degree. The degree variance plots also reveal that most of the power is in the lower degree terms. Contour maps of the radiation field show the geographical distribution of earth emitted radiant exitance for each month.

Bess, T. D.↗
NASA NTRS1990

A Solar-luminosity Model and Climate

Although the mechanisms of climatic change are not completely understood, the potential causes include changes in the Sun's luminosity. Solar activity in the form of sunspots, flares, proton events, and radiation fluctuations has displayed periodic tendencies. Two types of proxy climatic data that can be related to periodic solar activity are varved geologic formations and freshwater diatom deposits. A model for solar luminosity was developed by using the geometric progression of harmonic cycles that is evident in solar and geophysical data. The model assumes that variation in global energy input is a result of many periods of individual solar-luminosity variations. The 0.1-percent variation of the solar constant measured during the last sunspot cycle provided the basis for determining the amplitude of each luminosity cycle. Model output is a summation of the amplitudes of each cycle of a geometric progression of harmonic sine waves that are referenced to the 11-year average solar cycle. When the last eight cycles in Emiliani's oxygen-18 variations from deep-sea cores were standardized to the average length of glaciations during the Pleistocene (88,000 years), correlation coefficients with the model output ranged from 0.48 to 0.76. In order to calibrate the model to real time, model output was graphically compared to indirect records of glacial advances and retreats during the last 24,000 years and with sea-level rises during the Holocene. Carbon-14 production during the last millenium and elevations of the Great Salt Lake for the last 140 years demonstrate significant correlations with modeled luminosity. Major solar flares during the last 90 years match well with the time-calibrated model.

Perry, Charles A.↗
DOE OSTI2026

From models to reality: a systematic review on simulated and measured residential heat pump energy savings

High-performance HVAC solutions are central to residential energy management. A substantial share of these are electric, reversible-cycle systems, with heat pumps representing the largest portion of current and near-term adoption. This review synthesizes peer-reviewed and grey literature on residential space heating and cooling heat pumps. The academic literature is dominated by modeling (73.8%), with limited field measurement (13.1%). Grey literature from United States serve as a supplemental resource providing measured savings. Conversions from electric-resistance heating consistently show the largest site energy reductions, while oil/propane baselines yield moderate savings, and gas baseline scenario often deliver small and region-dependent savings. This study cross-checks the grey literature measured data with simulation data filtered from the ResStock dataset. The comparison indicates a discrepancy between simulations and measured data: simulated site EUIs are typically lower than measured EUIs, but percentage energy savings fall in similar ranges, implying simulations capture directional effects while underestimating energy use. Factors associated with variability and model–measurement differences include system characterization and control representation (e.g., backup heat engagement, thermostat/setpoint strategies, commissioning/installation quality), occupant behavior, weather normalization, metering scope, and envelope characterization. This paper also outlines the proposed methodology for comparing simulation and measured data for heat pumps. It emphasizes the metrics used for comparison and units harmonization, building characteristics matching, and compact metadata are needed for simulations to match measured data. The proposed methodology is expected to improve the credibility of simulated savings as measured evidence grows.

Yu, Lili↗