Search NASA⌕ Search

SEARCH · Search NASA

Results for “data and proxy modeling”

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 487 records · Page 27

Dynamics of Water Molecules and Ions in Concentrated Lithium Chloride Solutions Probed with Ultrafast 2D IR Spectroscopy

Water and ion dynamics in concentrated LiCl solutions were studied using ultrafast 2D IR spectroscopy with the methyl thiocyanate (MeSCN) CN stretch as the vibrational probe. The IR absorption spectrum of MeSCN has two peaks, one peak for water associated with the nitrogen lone pair of MeSCN (W) and the other peak corresponding to Li + associated with the lone pair (L). To extract the spectral diffusion (structural dynamics) of W and L species, we developed a method that isolates the peak of interest by subtracting the 2D Gaussian proxies of multiple interfering peaks. Center line slope data (normalized frequency–frequency correlation function) for 2D bands from the W and L are fit with triexponential functions. The fastest component (1.1–1.6 ps) is assigned to local hydrogen bond length fluctuations. The intermediate timescale (~4.0 ps) corresponds to the hydrogen bond network rearrangement. The slowest component decays in ~40 ps and corresponds to ion pair and ion cluster dynamics. The very similar W and L spectral diffusion indicates that the motions of the water and ions are strongly coupled. Orientational relaxations of the W and L species were extracted using a new method to eliminate the effects of overlapping peaks. The results show that MeSCN bound to water undergoes orientational relaxation significantly faster than MeSCN bound to Li + . The orientational and spectral diffusion results are compared. In conclusion, a Stark coupling model is used to extract the root mean square average electric field caused by the ion clouds along the CN moiety as a function of concentration.

14 SOLAR ENERGY↗

Machine Learning for Solar Accessibility: Implications for Low-Income Solar Expansion and Profitability

The solar industry in the US typically uses a credit score such as the FICO score as an indicator of consumer utility payment performance and credit worthiness to approve customers for new solar installations. Using data on over 800,000 utility payment performance and over 5,000 demographic variables, we compare machine learning and econometric models to predict the probability of default to credit-score cutoffs. Here, we compare these models across a variety of measures, including how they affect consumers of different socio-economic backgrounds and profitability. We find that a traditional regression analysis using a small number of variables specific to utility repayment performance greatly increases accuracy and LMI inclusivity relative to FICO score, and that using machine learning techniques further enhances model performance. Relative to FICO, the machine learning model increases the number of low-to-moderate income consumers approved for community solar by 1.1% to 4.2% depending on the stringency used for evaluating potential customers, while decreasing the default rate by 1.4 to 1.9 percentage points. Using electricity utility repayment as a proxy for solar installation repayment, shifting from a FICO score cutoff to the machine learning model increases profits by 34% to 1882% depending on the stringency used for evaluating potential customers.

Davuluri, Sruthi↗

Modeling precipitation <i>δ</i><sup>18</sup>O variability in East Asia since the Last Glacial Maximum: temperature and amount effects across different timescales

Water isotopes in precipitation have played a key role in the reconstruction of past climate on millennial timescales and longer. However, for midlatitude regions like East Asia with complex terrain, the reliability behind the basic assumptions of the temperature effect and amount effect is based on modern observational data and still remains unclear for past climate. In the present work, we reexamine the two basic effects on seasonal, interannual, and millennial timescales in a set of time slice experiments for the period 22–0 ka using an isotope-enabled atmospheric general circulation model (AGCM). Our study confirms the robustness of the temperature and amount effects on the seasonal cycle over China in the present climatic conditions, with the temperature effect dominating in northern China and the amount effect dominating in the far south of China but no distinct effect in the transition region of central China. However, our analysis shows that neither temperature nor amount effect is significantly dominant over China on millennial and interannual timescales, which is a challenge to those classic assumptions in past climate reconstruction. Our work helps shed light on the interpretation of the proxy record of δ 18 O from a modeling point of view.

54 ENVIRONMENTAL SCIENCES↗

Reduction of structural impacts and distinction of photosynthetic pathways in a global estimation of GPP from space-borne solar-induced chlorophyll fluorescence

Quantifying global photosynthesis remains a challenge due to a lack of accurate remote sensing proxies. Solar-induced chlorophyll fluorescence (SIF) has been shown to be a good indicator of photosynthetic activity across various spatial scales. However, a global and spatially challenging estimate of terrestrial gross primary production (GPP) based on satellite SIF remains unresolved due to the confounding effects of species-specific physical and physiological traits and external factors, such as canopy structure or photosynthetic pathway (C3orC4). Here we analyze an ensemble of far-red SIF data from OCO-2 satellite and ground observations at multiple sites, using the spectral invariant theory to reduce the effects of canopy structure and to retrieve a structure-corrected total canopy SIF emission (SIF total). We find that the relationships between observed canopy-leavingSIF and ecosystem GPP vary significantly among biomes. In contrast, the relationships between SIF total and GPP converge around two unique models, one for C3and one for C4plants. We show that the two single empirical models can be used to globally scale satellite SIF observations to terrestrial GPP. We obtain an independent estimate of global terrestrial GPP of 129.56 ± 6.54 PgC/year for the 2015–2017 period, which is consistent with the state-of-the-art data- and process-oriented models. The new GPP product shows improved sensitivity to previously undetected ‘hotspots’ of productivity, being able to resolve the double-peak in GPP due to rotational cropping systems. We suggest that the direct scheme to estimate GPP presented here, which is based on satelliteSIF, may open up new possibilities to resolve the dynamics of global terrestrial GPP across space and time.

scale satellite SIF observations to terrestrial GP↗

A First Approach to Global Runoff Simulation using Satellite Rainfall Estimation

Many hydrological models have been introduced in the hydrological literature to predict runoff but few of these have become common planning or decision-making tools, either because the data requirements are substantial or because the modeling processes are too complicated for operational application. On the other hand, progress in regional or global rainfall-runoff simulation has been constrained by the difficulty of measuring spatiotemporal variability of the primary causative factor, i.e. rainfall fluxes, continuously over space and time. Building on progress in remote sensing technology, researchers have improved the accuracy, coverage, and resolution of rainfall estimates by combining imagery from infrared, passive microwave, and space-borne radar sensors. Motivated by the recent increasing availability of global remote sensing data for estimating precipitation and describing land surface characteristics, this note reports a ballpark assessment of quasi-global runoff computed by incorporating satellite rainfall data and other remote sensing products in a relatively simple rainfall-runoff simulation approach: the Natural Resources Conservation Service (NRCS) runoff Curve Number (CN) method. Using an Antecedent Precipitation Index (API) as a proxy of antecedent moisture conditions, this note estimates time-varying NRCS-CN values determined by the 5-day normalized API. Driven by multi-year (1998-2006) Tropical Rainfall Measuring Mission (TRMM) Multi-satellite Precipitation Analysis, quasi-global runoff was retrospectively simulated with the NRCS-CN method and compared to Global Runoff Data Centre data at global and catchment scales. Results demonstrated the potential for using this simple method when diagnosing runoff values from satellite rainfall for the globe and for medium to large river basins. This work was done with the simple NRCS-CN method as a first-cut approach to understanding the challenges that lie ahead in advancing the satellite-based inference of global runoff. We expect that the successes and limitations revealed in this study will lay the basis for applying more advanced methods to capture the dynamic variability of the global hydrologic process for global runoff monlto~ngin real time. The essential ingredient in this work is the use of global satellite-based rainfall estimation.

Hong, Yang↗

Hole formation effect on surface morphological response of plasma-facing tungsten

We report simulation results on the effect of helium (He) bubble bursting-mediated surface hole formation on the surface morphological response of tungsten plasma-facing components (PFCs) in nuclear fusion devices. Our analysis is based on an atomistically informed, continuum-scale model, which is capable of accessing the spatiotemporal scales relevant to the fuzz nanostructure formation process on the surface of PFC tungsten. Our simulations account, in an empirical fashion, for two types of subsurface bubble dynamical phenomena in the nanobubble region of PFC tungsten during He plasma irradiation, involving bubble bursting and surface crater formation. We demonstrate that the hole formation effect on the PFC tungsten surface accelerates the growth rate of nanotendrils and the onset of fuzz formation. As a result, the predicted incubation time for surface nanotendril growth is shortened, in agreement with experimental data of incubation fluence at comparable plasma exposure conditions. We also explore systematically the dependence of the PFC surface morphological response on the areal density of holes introduced at regular time intervals onto the He-implanted tungsten surface, a parameter in our analysis that serves as a proxy for the rate of He bubble bursting. More importantly, our simulations capture fine surface features in the PFC tungsten surface morphology and predict that the average spacing between nanotendrils is on the order of 100 nm, consistent with the experimental findings.

Chen, Chao-Shou (ORCID:0000000307910042)↗

Carbon isotope fractionation by an ancestral rubisco suggests that biological proxies for CO 2 through geologic time should be reevaluated

The history of Earth’s carbon cycle reflects trends in atmospheric composition convolved with the evolution of photosynthesis. Fortunately, key parts of the carbon cycle have been recorded in the carbon isotope ratios of sedimentary rocks. The dominant model used to interpret this record as a proxy for ancient atmospheric CO 2 is based on carbon isotope fractionations of modern photoautotrophs, and longstanding questions remain about how their evolution might have impacted the record. Therefore, we measured both biomass (ε p ) and enzymatic (ε Rubisco ) carbon isotope fractionations of a cyanobacterial strain ( Synechococcus elongatus PCC 7942) solely expressing a putative ancestral Form 1B rubisco dating to ≫1 Ga. This strain, nicknamed ANC, grows in ambient pCO 2 and displays larger ε p values than WT, despite having a much smaller ε Rubisco (17.23 ± 0.61‰ vs. 25.18 ± 0.31‰, respectively). Surprisingly, ANC ε p exceeded ANC ε Rubisco in all conditions tested, contradicting prevailing models of cyanobacterial carbon isotope fractionation. Such models can be rectified by introducing additional isotopic fractionation associated with powered inorganic carbon uptake mechanisms present in Cyanobacteria, but this amendment hinders the ability to accurately estimate historical pCO 2 from geological data. Understanding the evolution of rubisco and the CO 2 concentrating mechanism is therefore critical for interpreting the carbon isotope record, and fluctuations in the record may reflect the evolving efficiency of carbon fixing metabolisms in addition to changes in atmospheric CO 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the predominance of unstable atmospheric conditions in the marine boundary layer offshore of the U.S. northeastern coast

Abstract The marine boundary layer of the northeastern U.S. is studied with focus on wind speed, atmospheric stability, and turbulent kinetic energy (TKE), the three most relevant properties in the context of offshore wind power development. Two long‐term observational data sets are analyzed. The first one consists of multilevel meteorological variables measured up to 60 m during 2003–2011 at the offshore Cape Wind tower, located near the center of the Nantucket Sound. The second data set comes from the 2013–2014 IMPOWR campaign (Improving the Modeling and Prediction of Offshore Wind Resources), in which wind and wave data were collected with new instruments on the Cape Wind platform, in addition to meteorological data measured during 19 flight missions offshore of New York, Connecticut, Rhode Island, and Massachusetts. It is found that, in this region: (1) the offshore wind resource is remarkable, with monthly average wind speeds at 60 m exceeding 7 m s −1 all year round, highest winds in winter (10.1 m s −1 ) and lowest in summer (7.1 m s −1 ), and a distinct diurnal modulation, especially in summer; (2) the marine boundary layer is predominantly unstable (61% unstable vs. 21% neutral vs. 18% stable), meaning that mixing is strong, heat fluxes are positive, and the wind speed profile is often nonlogarithmic (~40% of the time); and (3) the shape of the wind speed profile (log versus nonlog) is an effective qualitative proxy for atmospheric stability, whereas TKE alone is not.

54 ENVIRONMENTAL SCIENCES↗

Synergistic Use of Hyperspectral UV-Visible OMI and Broadband Meteorological Imager MODIS Data for a Merged Aerosol Product

The retrieval of optimal aerosol datasets by the synergistic use of hyperspectral ultraviolet(UV)–visible and broadband meteorological imager (MI) techniques was investigated. The Aura Ozone Monitoring Instrument (OMI) Level 1B (L1B) was used as a proxy for hyperspectral UV–visible instrument data to which the Geostationary Environment Monitoring Spectrometer (GEMS) aerosol algorithm was applied. Moderate-Resolution Imaging Spectroradiometer (MODIS) L1B and dark target aerosol Level 2 (L2) data were used with a broadband MI to take advantage of the consistent time gap between the MODIS and the OMI. First, the use of cloud mask information from the MI infrared (IR) channel was tested for synergy. High-spatial-resolution and IR channels of the MI helped mask cirrus and sub-pixel cloud contamination of GEMS aerosol, as clearly seen in aerosol optical depth (AOD) validation with Aerosol Robotic Network (AERONET) data. Second, dust aerosols were distinguished in the GEMS aerosol-type classification algorithm by calculating the total dust confidence index (TDCI) from MODIS L1B IR channels. Statistical analysis indicates that the Probability of Correct Detection (POCD) between the forward and inversion aerosol dust models (DS) was increased from 72% to 94% by use of the TDCI for GEMS aerosol-type classification, and updated aerosol types were then applied to the GEMS algorithm. Use of the TDCI for DS type classification in the GEMS retrieval procedure gave improved single-scattering albedo (SSA) values for absorbing fine pollution particles (BC) and DS aerosols. Aerosol layer height (ALH) retrieved from GEMS was compared with Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) data, which provides high-resolution vertical aerosol profile information. The CALIOP ALH was calculated from total attenuated backscatter data at 1064 nm, which is identical to the definition of GEMS ALH. Application of the TDCI value reduced the median bias of GEMS ALH data slightly. The GEMS ALH bias approximates zero, especially for GEMS AOD values of>~0.4 and GEMS SSA values of<~0.95.Finally, the AOD products from the GEMS algorithm and MI were used in aerosol merging with the maximum-likelihood estimation method, based on a weighting factor derived from the standard deviation of the original AOD products. With the advantage of the UV–visible channel in retrieving aerosol properties over bright surfaces, the combined AOD products demonstrated better spatial data availability than the original AOD products, with comparable accuracy. Furthermore, pixel-level error analysis of GEMS AOD data indicates improvement through MI synergy.

aerosol↗

Towards quantifying the impact of triaxiality on optical signatures of galaxy clusters: weak lensing and galaxy distributions

ABSTRACT We present observational evidence of the impact of triaxiality on radial profiles that extend to 40 Mpc from galaxy cluster centres in optical measurements. We perform a stacked profile analysis from a sample of thousands of nearly relaxed galaxy clusters from public data releases of the Dark Energy Survey and the Dark Energy Camera Legacy Survey. Using the central galaxy elliptical orientation angle as a proxy for galaxy cluster orientation, we measure cluster weak lensing and excess galaxy density axis-aligned profiles, extracted along the central galaxy’s major or minor axes on the plane of the sky. Our measurements show a ≳ 2σ–3σ difference per radial bin between the normalized axis-aligned profiles. The profile difference between each axis-aligned profile and the azimuthally averaged profile ($\sim \pm 10\,\rm per\ cent-20~{{\ \rm per\ cent}}$ along major/minor axis) appears inside the clusters (∼0.4 Mpc) and extends to the large-scale structure regime (∼10–20 Mpc). The magnitude of the difference appears to be relatively insensitive to cluster richness and redshift, and extends further out in the weak lensing surface mass density than in the galaxy overdensity. Looking forward, this measurement can easily be applied to other observational or simulation data sets and can inform the systematics in cluster mass modelling related to triaxiality. We expect imminent upcoming wide-area deep surveys, such as the Vera C. Rubin Observatory’s Legacy Survey of Space and Time, to improve our quantification of optical signatures of cluster triaxiality.

79 ASTRONOMY AND ASTROPHYSICS↗

Decomposing Fuel Economy and Greenhouse Gas Regulatory Standards in the Energy Conversion Efficiency and Tractive Energy Domain

The three foundational elements that determine mobile source energy use and tailpipe carbon dioxide (CO2) emissions are the tractive energy requirements of the vehicle, the on-cycle energy conversion efficiency of the propulsion system, and the energy source. The tractive energy requirements are determined by the vehicle's mass, aerodynamic drag, tire rolling resistance, and parasitic drag. Oncycle energy conversion of the propulsion system is dictated by the tractive efficiency, non-tractive energy use, kinetic energy recovery, and parasitic losses. The energy source determines the mobile source CO2 emissions. For current vehicles, tractive energy requirements and overall energy conversion efficiency are readily available from the decomposition of test data. For future applications, plausible levels of mass reduction, aerodynamic drag improvements, and tire rolling resistance can be transposed into the tractive energy domain. Similarly, by combining thermodynamic, mechanical efficiency, and kinetic energy recovery fundamentals with logical proxies, achievable levels of energy conversion efficiency can be established to allow for the evaluation of future powertrain requirements. Combining the plausible levels of tractive energy and on-cycle efficiency provides a means to compute sustainable vehicle and propulsion system scenarios that can achieve future regulations. Using these principles, the regulations established in the United States (U.S.) for fuel consumption and CO2 emissions are evaluated. Fleet-level scenarios are generated and compared to the technology deployment assumptions made during rule-making. When compared to the rule-making assumptions, the results indicate that a greater level of advanced vehicle and propulsion system technology deployment will be required to achieve the model year 2025 U.S. standards for fuel economy and CO2 emissions.

CAFE↗

Probing the limits of metal plasticity with molecular dynamics simulations

Ordinarily, the strength and plasticity properties of a metal are defined by dislocations—line defects in the crystal lattice whose motion results in material slippage along lattice planes. Dislocation dynamics models are usually used as mesoscale proxies for true atomistic dynamics, which are computationally expensive to perform routinely. However, atomistic simulations accurately capture every possible mechanism of material response, resolving every “jiggle and wiggle” of atomic motion, whereas dislocation dynamics models do not. We present fully dynamic atomistic simulations of bulk single-crystal plasticity in the body-centred-cubic metal tantalum. Our goal is to quantify the conditions under which the limits of dislocation-mediated plasticity are reached and to understand what happens to the metal beyond any such limit. In our simulations, the metal is compressed at ultrahigh strain rates along its [001] crystal axis under conditions of constant pressure, temperature and strain rate. To address the complexity of crystal plasticity processes on the length scales (85–340 nm) and timescales (1 ns–1μs) that we examine, we use recently developed methods of in situ computational microscopy to recast the enormous amount of transient trajectory data generated in our simulations into a form that can be analysed by a human. Our simulations predict that, on reaching certain limiting conditions of strain, dislocations alone can no longer relieve mechanical loads; instead, another mechanism, known as deformation twinning (the sudden re-orientation of the crystal lattice), takes over as the dominant mode of dynamic response. Below this limit, the metal assumes a strain-path-independent steady state of plastic flow in which the flow stress and the dislocation density remain constant as long as the conditions of straining thereafter remain unchanged. In this distinct state, tantalum flows like a viscous fluid while retaining its crystal lattice and remaining a strong and stiff metal.

36 MATERIALS SCIENCE↗

Global variations in critical drought thresholds that impact vegetation

ABSTRACT Identifying the thresholds of drought that, if crossed, suppress vegetation functioning is vital for accurate quantification of how land ecosystems respond to climate variability and change. We present a globally applicable framework to identify drought thresholds for vegetation responses to different levels of known soil-moisture deficits using four remotely sensed vegetation proxies spanning 2001–2018. The thresholds identified represent critical inflection points for changing vegetation responses from highly resistant to highly vulnerable in response to drought stress, and as a warning signal for substantial vegetation impacts. Drought thresholds varied geographically, with much lower percentiles of soil-moisture anomalies in vegetated areas covered by more forests, corresponding to a comparably stronger capacity to mitigate soil water deficit stress in forested ecosystems. Generally, those lower thresholds are detected in more humid climates. State-of-the-art land models, however, overestimated thresholds of soil moisture (i.e. overestimating drought impacts), especially in more humid areas with higher forest covers and arid areas with few forest covers. Based on climate model projections, we predict that the risk of vegetation damage will increase by the end of the twenty-first century in some hotspots like East Asia, Europe, Amazon, southern Australia and eastern and southern Africa. Our data-based results will inform projections on future drought impacts on terrestrial ecosystems and provide an effective tool for drought management.

54 ENVIRONMENTAL SCIENCES↗

Jason-1 and Jason-2 POD Using GPS

The Jason-2 satellite, launched in June 2008, is the latest follow-on to the successful Jason-1 altimetry satellite mission launched in December 7, 2001. Both, Jason-2 and Jason-1 are equipped with a GPS dual-frequency receiver, a laser retroreflector array, and a DORIS receiver for precise orbit determination (POD). A series of dynamic and reduced-dynamic Jason-2 orbits computed at NASA GSFC, based on GPS-only data and the std0905 standards, have been completed till cy74through cycle 74 using the IGS05 framework. These orbits, now publicly available, have been shown to agree radially at 1 cm RMS with the GSFC std0905 SLR/DORIS orbits and in comparison with orbits produced by JPL, ESA and CNES. In this paper, we describe the implementation of the IGS08 and repro1 framework for the Jason-2 and Jason-1 GPS POD processing with the NASA GSFC GEODYN software. . In our updated GPS POD, ambiguity fixing and updated time variable and static gravity fields. We also evaluate the implementation of non-tidal and degree-1 loading displacement as forward modeling to the tracking stations. Reduced-dynamic versus dynamic orbit differences are used to characterize the remaining force model errors and TRF instability. In particular, we assess their consistency radially and the stability of the altimeter satellite reference frame in the North/South direction as a proxy to assess the consistency of the reference frame.

Melachroinos, Stavros↗

Phase-Angle Dependence of Determinations of Diameter, Albedo, and Taxonomy: A Case Study of NEO 3691 Bede

Parameters important for NEO risk assessment and mitigation include Near-Earth Object diameter and taxonomic classification, which translates to surface composition. Diameters of NEOs are derived from the thermal fluxes measured by WISE, NEOWISE, Spitzer Warm Mission and ground-based telescopes including the IRTF and UKIRT. Diameter and its coupled parameters Albedo and IR beaming parameter (a proxy for thermal inertia and/or surface roughness) are dependent upon the phase angle, which is the Sun-target-observer angle. Orbit geometries of NEOs, however, typically provide for observations at phase angles greater than 20 degrees. At higher phase angles, the observed thermal emission is sampling both the day and night sides of the NEO. We compare thermal models for NEOs that exclude (NEATM) and include (NESTM) night-side emission. We present a case study of NEO 3691 Bede, which is a higher albedo object, X (Ec) or Cgh taxonomy, to highlight the range of H magnitudes for this object (depending on the albedo and phase function slope parameter G), and to examine at different phase angles the taxonomy and thermal model fits for this NEO. Observations of 3691 Bede include our observations with IRTF+SpeX and with the 10 micrometer UKIRT+Michelle instrument, as well as WISE and Spitzer Warm mission data. By examining 3691 Bede as a case study, we highlight the interplay between the derivation of basic physical parameters and observing geometry, and we discuss the uncertainties in H magnitude, taxonomy assignment amongst the X-class (P, M, E), and diameter determinations. Systematic dependencies in the derivation of basic characterization parameters of H-magnitude, diameter, albedo and taxonomy with observing geometry are important to understand. These basic characterization parameters affect the statistical assessments of the NEO population, which in turn, affects the assignment of statistically-assessed basic parameters to discovered but yet-to-be-fully-characterized NEOs.

Albedo and taxonomy↗

Simulation-based inference for parameter estimation of complex watershed simulators

High-resolution, spatially distributed process-based (PB) simulators are widely employed in the study of complex catchment processes and their responses to a changing climate. However, calibrating these PB simulators using observed data remains a significant challenge due to several persistent issues, including the following: (1) intractability stemming from the computational demands and complex responses of simulators, which renders infeasible calculation of the conditional probability of parameters and data, and (2) uncertainty stemming from the choice of simplified representations of complex natural hydrologic processes. Here, we demonstrate how simulation-based inference (SBI) can help address both of these challenges with respect to parameter estimation. SBI uses a learned mapping between the parameter space and observed data to estimate parameters for the generation of calibrated simulations. To demonstrate the potential of SBI in hydrologic modeling, we conduct a set of synthetic experiments to infer two common physical parameters – Manning's coefficient and hydraulic conductivity – using a representation of a snowmelt-dominated catchment in Colorado, USA. We introduce novel deep-learning (DL) components to the SBI approach, including an “emulator” as a surrogate for the PB simulator to rapidly explore parameter responses. We also employ a density-based neural network to represent the joint probability of parameters and data without strong assumptions about its functional form. While addressing intractability, we also show that, if the simulator does not represent the system under study well enough, SBI can yield unreliable parameter estimates. Approaches to adopting the SBI framework for cases in which multiple simulator(s) may be adequate are introduced using a performance-weighting approach. The synthetic experiments presented here test the performance of SBI, using the relationship between the surrogate and PB simulators as a proxy for the real case.

54 ENVIRONMENTAL SCIENCES↗

A Generative Model for Realistic Galaxy Cluster X-Ray Morphologies

Abstract The X-ray morphologies of clusters of galaxies display significant variations, reflecting their dynamical histories and the nonlinear dependence of X-ray emissivity on the density of the intracluster gas. Qualitative and quantitative assessments of X-ray morphology have long been considered a proxy for determining whether clusters are dynamically active or “relaxed.” Conversely, the use of circularly or elliptically symmetric models for cluster emission can be complicated by the variety of complex features realized in nature, spanning scales from megaparsecs down to the resolution limit of current X-ray observatories. In this work, we use mock X-ray images from simulated clusters from The Three Hundred project to define a basis set of cluster image features. We take advantage of the clusters’ approximate self-similarity to minimize the differences between images before encoding the remaining diversity through a distribution of high-order polynomial coefficients. Principal component analysis then provides an orthogonal basis for this distribution, corresponding to natural perturbations from an average model. This representation allows novel, realistically complex X-ray cluster images to be easily generated, and we provide code to do so. The approach provides a simple way to generate training data for cluster image analysis algorithms and could be straightforwardly adapted to generate clusters displaying specific types of features or selected by physical characteristics available in the original simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Estimation of the Relationship Between Remotely Sensed Anthropogenic Heat Discharge and Building Energy Use

This paper examined the relationship between remotely sensed anthropogenic heat discharge and energy use from residential and commercial buildings across multiple scales in the city of Indianapolis, Indiana, USA. The anthropogenic heat discharge was estimated with a remote sensing-based surface energy balance model, which was parameterized using land cover, land surface temperature, albedo, and meteorological data. The building energy use was estimated using a GIS-based building energy simulation model in conjunction with Department of Energy/Energy Information Administration survey data, the Assessor's parcel data, GIS floor areas data, and remote sensing-derived building height data. The spatial patterns of anthropogenic heat discharge and energy use from residential and commercial buildings were analyzed and compared. Quantitative relationships were evaluated across multiple scales from pixel aggregation to census block. The results indicate that anthropogenic heat discharge is consistent with building energy use in terms of the spatial pattern, and that building energy use accounts for a significant fraction of anthropogenic heat discharge. The research also implies that the relationship between anthropogenic heat discharge and building energy use is scale-dependent. The simultaneous estimation of anthropogenic heat discharge and building energy use via two independent methods improves the understanding of the surface energy balance in an urban landscape. The anthropogenic heat discharge derived from remote sensing and meteorological data may be able to serve as a spatial distribution proxy for spatially-resolved building energy use, and even for fossil-fuel CO2 emissions if additional factors are considered.

Multi-scale↗