Search NASA⌕ Search

SEARCH · Search NASA

Results for “difference map”

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 1,459 records · Page 81

Asymmetries in Core-Collapse Supernovae from Maps of Radioactiver 44Ti in Cassiopeia A

Asymmetry is required by most numerical simulations of stellar core-collapse explosions, but the form it takes differs significantly among models. The spatial distribution of radioactive 44Ti, synthesized in an exploding star near the boundary between material falling back onto the collapsing core and that ejected into the surroundingmedium1, directly probes the explosion asymmetries. Cassiopeia A is a young2, nearby3, core-collapse4 remnant from which 44Ti emission has previously been detected5-8 but not imaged. Asymmetries in the explosion have been indirectly inferred from a high ratio of observed 44Ti emission to estimated 56Ni emission9, from optical light echoes10, and from jet-like features seen in the X-ray11 and optical12 ejecta. Here we report spatial maps and spectral properties of the 44Ti in Cassiopeia A. This may explain the unexpected lack of correlation between the 44Ti and iron X-ray emission, the latter being visible only in shock-heated material. The observed spatial distribution rules out symmetric explosions even with a high level of convective mixing, as well as highly asymmetric bipolar explosions resulting from a fast-rotating progenitor. Instead, these observations provide strong evidence for the development of low-mode convective instabilities in core-collapse supernovae.

stellarcore-collapse explosions↗

AN INDUCTIVE MAPPING WITH CONVOLUTIONAL REPRESENTATIONS FOR HUMAN SETTLEMENT DETECTION: PRELIMINARY RESULTS

Undoubtedly, deep convolutional learning methods continue to improve performance in image-level classification for computer vision and remote sensing applications. However, the spatio-temporal nature of remote sensing imagery offers other interesting challenges in multiscale image understanding. Emerging opportunities include seeking fine-grained and neighborhood mapping with overhead imagery. Limitations due to lack of relevant scale ground-truth often mandates that these challenges are pursued as disjoint. We test this premise and explore representation spaces from a single deep convolutional network and their visualization to argue for a unified feature extraction framework. The objective is to utilize and re-purpose trained feature extractors without the need for network retraining on three remote sensing tasks i.e. superpixel mapping, pixel-level segmentation and semantic based image visualization. By leveraging the same convolutional feature extractors and viewing them as visual information extractors that encode different image representation spaces, we demonstrate surprising potential on multiscale experiments that incorporate edge-level details up to semantic-level informationWe test this premise and explore representation spaces from a single deep convolutional network and their visualization to argue for a novel unified feature extraction framework. The objective is to utilize and re-purpose trained feature extractors without the need for network retraining on three remote sensing tasks i.e. superpixel mapping, pixel-level segmentation and semantic based image visualization. By leveraging the same convolutional feature extractors and viewing them as visual information extractors that encode different image representation spaces, we demonstrate a preliminary inductive transfer learning potential on multiscale experiments that incorporate edge-level details up to semantic-level informationWe test this premise and explore representation spaces from a single deep convolutional network and their visualization to argue for a novel unified feature extraction framework. The objective is to utilize and re-purpose trained feature extractors without the need for network retraining on three remote sensing tasks i.e. superpixel mapping, pixel-level segmentation and semantic based image visualization. By leveraging the same convolutional feature extractors and viewing them as visual information extractors that encode different image representation spaces, we demonstrate a preliminary inductive transfer learning potential on multiscale experiments that incorporate edge-level details up to semantic-level information.

Lunga, Dalton↗

Ecosystems and Networks Integrated with Genes and Molecular Assemblies (ENIGMA). Final report

The goal of the Stahl lab is to use a combination of isotopic, genomic, and activity-based analyses in the field, of selected field isolates, and defined synthetic communities to identify sources and sinks of different nitrogen species in the field, focusing on the FRC at the ORNL. The most general research objective is to characterize variants in the pathway for denitrification among field isolates and synthetic microbial assemblies as a foundation for predicting the biotic and abiotic controls of nitrogen and carbon transformations in the field. Specifically, we have focused on characterizing environmental controls of nitrous oxide emissions, an important greenhouse gas and ozone depleting molecule. Studies have mapped the flux of nitrous oxide and controlling genetic elements within the FRC subsurface, indicating the importance of different enzyme variants in nitrous oxide production and consumption in controlling surface emissions. We have developed methods to quantify and catalog the isotopic signatures imparted on nitrous oxide by different nitrous oxide reductase variants in isolated organisms and model synthetic communities. We have developed methods to quantify and catalog the nitrous oxide affinity of field relevant isolates and their biomass yield when respiring nitrous oxide on different carbon substrates. We continue to expand the collection of field-relevant microbial isolates within ENIGMA by focused enrichment and isolation from field material. These objectives have been facilitated by the development of new technologies, including: 1) the development of a multiplexed automated optical density measurement system for quantifying growth kinetics of microorganisms dependent on gaseous substrates in monoculture or gas exchange in synthetic communities and 2) the development of a microcalorimetric analysis format to quantify the impact of environmental stress on microbial maintenance energy. Since field relevant processes are mediated by assemblies of organisms interacting with each other and their local environment, we have also collaborated with other ENIGMA investigators on the assembly and characterization of synthetic communities. These studies have evaluated the role of organisms lacking a complete pathway for denitrification in the control of field processes, measuring process rates of assemblies of organisms that individually lack the full denitrification pathway but mediate the complete denitrification when growing together. These efforts have so far resulted in the publication of 18 peer reviewed papers ranging from ecological and physiological drivers of pathway optimization to field-scale process characterization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

ALMA Spectral Line and Imaging Survey of a Low and a High Mass-Loss Rate AGB Star Between 335 and 362 GHz

Context. Low and intermediate mass stars are known to power strong stellar winds when evolving through the asymptotic giant branch (AGB) phase. Initial mass, luminosity, temperature, and composition determine the pulsation characteristics of the star and the dust species formed in the pulsating photospheric layers. Radiation pressure on these grains triggers the onset of a stellar wind. However, as of today, we still cannot predict the wind mass-loss rates and wind velocities from first principles neither do we know which species are the first to condense in the upper atmospheric regions. Aims. We aim to characterise the dominant physical, dynamical, and chemical processes in the inner wind region of two archetypical oxygen-rich (C/O<1) AGB stars, that is, the low mass-loss rate AGB star R Dor (M~ 1x10(exp 7) M☉/yr) and the high mass-loss rate AGB star IK Tau (M~5x10(exp 6) M☉/yr). The purpose of this study is to observe the key molecular species contributing to the formation of dust grains and to cross-link the observed line brightnesses of several species to the global and local properties of the star and its wind. Methods. A spectral line and imaging survey of IK Tau and R Dor was made with ALMA between 335 and 362 GHz (band 7) at a spatial resolution of ~150 mas, which corresponds to the locus of the main dust formation region of both targets. Results. Some two hundred spectral features from 15 molecules (and their isotopologues) were observed, including rotational lines in both the ground and vibrationally excited states (up to v=5 for SiO). Detected species include the gaseous precursors of dust grains such as SiO, AlO, AlOH, TiO, and TiO2.We present a spectral atlas for both stars and the parameters of all detected spectral features. A clear dichotomy for the sulphur chemistry is seen: while CS, SiS, SO, and SO2 are abundantly present in IK Tau, only SO and SO2 are detected in R Dor. Also other species such as NaCl, NS, AlO, and AlOH display a completely different behaviour. From some selected species, the minor isotopologues can be used to assess the isotopic ratios. The channel maps of many species prove that both large and small-scale inhomogeneities persist in the inner wind of both stars in the form of blobs, arcs, and/or a disk. The high sensitivity of ALMA allows us to spot the impact of these correlated density structures in the spectral line profiles. The spectral lines often display a half width at zero intensity much larger than expected from the terminal velocity, v1, previously derived for both objects (36 km/s versus v(∞)~17.7 km/s for IK Tau and 23 km/s versus v(∞)~5.5 km/s for R Dor). Both a more complex 3D morphology and a more forceful wind acceleration of the (underlying) isotropic wind can explain this trend. The formation of fractal grains in the region beyond _400 mas can potentially account for the latter scenario. From the continuum map, we deduce a dust mass of ~3.7 x 10(exp 7) M☉ and ~2 x 10(ex0 8) M☉ for IK Tau and R Dor, respectively. Conclusions. The observations presented here provide important constraints on the properties of these two oxygen-dominated AGB stellar winds. In particular, the ALMA data prove that both the dynamical and chemical properties are vastly different for this high mass-loss rate (IK Tau) and low mass-loss rate (R Dor) star.

Stars: AGB and post-AGB, Stars: mass loss, Stars: ↗

Simulating water dynamics related to pedogenesis across space and time: Implications for four-dimensional digital soil mapping

Digital soil mapping (DSM) relies on machine-learning and geostatistics to represent soil property observations across space. DSM techniques are powerful but often empirical, being limited to the quality and density of point samples. Water dynamics are closely related to soil variability, and the physics that govern water movement are well known. Hydrological properties can hence be simulated by physical models through space and time, unveiling key characteristics about soils. We propose the use of hydrologic models to map soils across the surface (2D), depth (1D), and time (1D)–which provides a 4D approach to digital soil mapping (4DSM). The Distributed Hydrology Soil Vegetation Model (DHSVM) was applied to a watershed currently under pasture. Moisture sensors and wells were installed at different depths in the watershed on summit, sideslope and toeslope positions to validate the model. DHSVM simulations of soil moisture distribution and depth to saturation were performed during the hydrological year (October 2008-September 2009). Clusters of similar pixels based on soil moisture values were determined using Dynamic Time Warping (DTW) to align temporal data and K-means. Clustering was performed both seasonally and for the entire year. Temporal patterns simulated by DHSVM matched measurements given by moisture sensors and wells. Seasonal clusters differed from the annual cluster. Distinct clusters were observed for each season and with depth, showing that spatiotemporal soil variability is lost when statically assessing soils. Spatiotemporal clusters corroborated field observations of fragipan occurrence not explicitly spatially mapped by Soil Survey Geographic Database (SSURGO). If a connection can be made between water and soils, static and dynamic soil variability can be predicted using physically based hydrologic models. Hydrologic models can benefit soil mapping by enabling reliable 4D simulation of water dynamics, which are fundamental to soil variability and soil classification and directly relate to biological, physical and chemical soil processes not captured by typical soil sampling protocols.

54 ENVIRONMENTAL SCIENCES↗

Ice Thicknesses and Driving Stresses of Greenland Ice Sheet

Surface elevations from satellite radar altimetry (Geosat, Seasat, and ERS-1) and bedrock topography from airborne radar sounding (Simon Ekholm's Danish compilation) are combined to derive maps of the driving stresses in the Greenland ice sheet. The stress vector, tau = rho g h sin(alpha), is calculated using surface slope vectors, alpha, from surface elevations and ice thicknesses, h, from the difference between surface and basal elevations. Since the 5-km scale of the surface slope is only about 2 times the ice thickness, the stress maps show spatial variations indicative of longitudinal stress gradients associated with topographic undulations. Values of alpha generally vary from near zero at the ice divides to maxima values around 120 kpa, returning to near zero in a narrow band at the edges. The distribution of alpha's peaks at 60 kpa with an approximate sigma of +/- 20 kpa. Areas of very low alpha near the origin of the northeast ice stream may indicate small sub-glacial lakes. The profile of alpha, down the ice stream from near the ice divide, increases to a maximum of about 120 kpa near the margin, which is characteristic of East Antarctic outlet glaciers and in contrast to West Antarctic ice streams where alpha has maximum values 400 to 500 km inland from the grounding lines. Overall distributions of alpha values are compared with those for the Antarctic ice sheet and the Mars Northern ice cap.

Zwally, H. J.↗

Leveraging gradient weighted class activation mapping to improve classification effectiveness: Case study in transportation infrastructure characterization

Roadway “corners†are common for pedestrian use, whether designated with markings or not. Different types of markings have been deployed, ranging from simple parallel lines to more complex designs. Understanding the impact of different types of crosswalks is important for public safety. In this work we explore methods to improve the logging of marked crosswalk types. We used the Roadway Information Database from the Second Strategic Highway Research Project and used active learning methods with transfer learning to identify the crosswalk types (marked or unmarked). Upon completion we found our classifiers were unable to perform above roughly 94% correct classifications. To improve their efficacy, we separated the crosswalks into their “fine grained†types and used Gradient-Weighted Class Activation Mapping to isolate and study the features that classified the crosswalks. We compared this with sampled manually marked crosswalks and present findings. We believe this use case can represent a process to improve the active learning method for some visual machine learning applications.

Karnowski, Thomas↗

Investigation and Expression Analysis of R2R3-MYBs and Anthocyanin Biosynthesis-Related Genes during Seed Color Development of Common Bean ( Phaseolus vulgaris )

Anthocyanins are responsible for the coloration of common bean seeds, and their accumulation is positively correlated with the expression level of anthocyanin biosynthetic genes. The MBW (MYB-bHLH-WD40) complex is thought to regulate the expression of these genes, and MYB proteins, which are a key factor in activating anthocyanin pathway genes, have been identified in several plants. This study demonstrated gene structures, chromosomal placements, gene duplications of R2R3-MYBs, miRNAs associated with R2R3-MYBs, and the interaction of these genes with other flavonoid regulatory genes. qRT-PCR was used to investigate the role of specific R2R3-MYBs and flavonoid genes in common bean seed color development. As a result of a comprehensive analysis with the help of in silico tools, we identified 160 R2R3-MYB genes in the common bean genome. We divided these genes into 16 classes on the basis of their intron-exon and motif structures. Except for three, the rest of the common bean R2R3-MYB members were distributed to all chromosomes with different densities, primarily located on chromosomes 3 and 8. We identified a total of 44 duplicated gene pairs dispersed across 11 chromosomes and evolved under purifying selection (Ka/Ks < 1), 19 of which were derived from a whole-genome duplication. Our research uncovered 25 putative repressor PvMYB proteins that contain the EAR motif. Additionally, fifty different cis-regulatory elements regulated by light, stress, and hormone were identified. Within the genome of the common bean, we discovered a total of 36 microRNAs that target a total of 72 R2R3-MYB transcripts. The effect of 16 R2R3-MYB genes and 16 phenylpropanoid pathway genes, selected on the basis of their interaction in the protein-protein interaction map, playing role in the regulation of seed coat color development was evaluated using qRT-PCR in 5 different tissues at different developmental stages. The results revealed that these specific genes have different expression levels during different developmental periods, with higher levels in the pod filling and early pod stages than in the rest of the developmental periods. Furthermore, it was shown that PvTT8 (bHLH), PvTT2 (PvMYB42), PvMYB113, PvTTG1, and PvWD68 genes have effects on the regulation of seed coat color. The findings of this study, which is the first to use whole-genome analysis to identify and characterize the R2R3-MYB genes in common bean, may serve as a reference for future functional research in the legume.

59 BASIC BIOLOGICAL SCIENCES↗

Bias correcting regional scale Earth system model projections: novel approach using empirical mode decomposition

Bias correction is a crucial step in using Earth system model outputs for assessments, as it adjusts systematic errors by comparing the model to observations. However, standard methods – ranging from mean-based linear scaling to distribution-based quantile mapping typically treat bias correction as a single-scale process, overlooking the fact that biases can manifest differently across daily, seasonal, and annual timescales. In this study, we propose a novel, timescale-aware bias-correction approach built on Empirical Mode Decomposition. By decomposing the meteorological signal into multiple oscillatory components and aggregating them to represent distinct timescales, we apply targeted corrections to each component, thereby preserving both short- and long-term structure in the data. Experimental illustrations show that the timescale-aware EMDBC framework matches the performance of conventional quantile-delta mapping (QDM) at the native daily scale and achieves progressively larger bias reductions at bi-weekly, seasonal, and annual scales. As a result, the proposed approach offers a more robust path to accurate and reliable Earth system projections, strengthening their utility for resilience and adaptation planning.

Ganguli, Arkaprabha [Argonne National Laboratory (↗

Dissolution of a cylindrical disk in Hele-Shaw flow: a conformal-mapping approach

Here we apply conformal mapping to find the evolving shapes of a dissolving cylinder in a potential flow. Similar equations can be used to describe melting in a flowing liquid phase. Results are compared with microfluidic experiments and numerical simulations. Shapes predicted by conformal mapping agree almost perfectly with experimental observations, after a modest (20 %) rescaling of the time. Finite-volume simulations show that the differences with experiment are connected to the underlying assumptions of the analytical model: potential flow and diffusion-limited dissolution. Approximate solutions of the equations describing the evolution of the shape of the undissolved solid can be derived from a Laurent expansion of the mapping function from the unit circle. Asymptotic expressions for the evolution of the area of the disk and the shift in its centre of mass have been derived at low and high Péclet number. Analytic approximations to the leading-order Laurent coefficients provide additional insight into the mechanisms underlying pore-scale dissolution.

58 GEOSCIENCES↗

Application of ERTS-1 imagery to detecting and mapping modern erosion features, and to monitoring erosional changes, in southern Arizona

The author has identified the following significant results. The red MSS band 5 gives the sharpest definition of modern arroyos. On the best images, modern arroy0s can be distinguished as narrow as 150 to 200 feet in reaches where their contrast with adjacent areas is only moderate, and as narrow as 60 to 75 feet where their contrast is high. Both the red and infrared bands show differences is soils and vegetation. In the late fall and winter imagery, band 7 generally is the most useful for mapping the areas of the more erodible soils. A map at 1:1,000,000 scale has been prepared that shows all the arroyos within the 17,000 square mile study area that have been identified from ERTS-1 images. Also, from U-2 color infrared airphotos, a 1:125,000 scale map has been made of a 50 mile reach along San Simon Wash, in southeastern Arizona. This map shows not only the arroyo channels and narrow flood plains that have developed since 1890, but also areas within a few miles of the wash that are severely guilled, severely sheet-eroded, and moderately sheet-eroded. Two important effects of the third largest recorded flood of the upper Gila River also have been determined from the ERTS-1 images. The inundated area is best displayed on band 7, and the areas of severe sand/gravel erosion/deposition show best on band 5.

Morrison, R. B.↗

High resolution sub-millimetre mapping of starburst galaxies: Comparison with CO emission

Researchers present first results from a program of submillimeter continuum mapping of starburst galaxies, and comparison of their dust and CO emission. This project was prompted by surprising results from the first target, the nearby starburst M82, which shows in the dust continuum a morphology quite unlike that of its CO emission, in contrast to what might be expected if both CO and dust are accurately tracing the molecular hydrogen. Possible explanations for this striking difference are discussed. In the light of these results, the program has been extended to include sub-mm mapping of the nearby, vigorously star forming spirals, M83 and Maffei 2. The latter were also observed extensively in CO, in order to study excitation conditions in its central regions. The James Clerk Maxwell Telescope was used in these studies.

Smith, P. A.↗

Effect of different crystallographic properties on the electrical conductivity of two polymorphs of a spin crossover complex

In this study, the structure and transport properties of two polymorphs, nanoparticles and nanorods, of the iron(II) triazole [Fe(Htrz) 2 (trz)](BF 4 ) spin crossover complex were compared. Conductive atomic force microscopy was used to map the electrical conductivity of individual nanoparticles and nanorods. The [Fe(Htrz) 2 (trz)](BF 4 ) nanorods showed significantly higher conductivity compared to nanoparticles. This difference in electrical conductivity is partially associated to the different Fe–N bond lengths in each of the polymorphs, with an inverse relationship between Fe–N bond length and conductivity. Transport measurements were done on the nanorods for both high spin (at 380 K) and low spin (at 320 K) states under dark and illuminated conditions. The conductance is highest for the low spin state under dark conditions. In illumination, the conductance change is much diminished.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modeling the effects of molecular disorder on the properties of Frenkel excitons in organic molecular semiconductors

The Frenkel exciton model offers a convenient framework for simulating electronic excitations in organic conjugated systems that are too large to address with atomistic level electronic structure methods. Parameterization of this model is typically based on analytical expressions that incompletely describe the spatial and temporal correlations that are inherent to many condensed phase molecular systems. In this manuscript, we introduce a general procedure for including these correlations in the Frenkel exciton model, by mapping them directly from all-atom molecular configurations, for instance from classical molecular dynamics. Regardless of system morphology, this mapping automatically captures the spatial and temporal molecular correlations that are otherwise difficult or impossible to represent in terms of low-dimensional correlation functions. We apply this procedure to study the excited state properties of condensed phase materials made up of thiophene oligomers. We demonstrate that Frenkel model parameters can be mapped from a series of single molecule electronic structure calculations, and that for these materials efficient semi-empirical methods are sufficient to accurately reproduce experimental spectral measurements. By analyzing the statistics of model parameters derived from materials with different characteristic morphologies, we highlight failures in some assumptions that are commonly applied when generating model parameters. Finally, by simulating exciton dynamics on a mapped Frenekel exciton model, we demonstrate the ability to quantify the effect of material morphology on the dynamic properties of excitons.

74 ATOMIC AND MOLECULAR PHYSICS↗

Concentric Spherical GNN for 3D Representation Learning

Learning 3D representations that generalize well to arbitrarily oriented inputs is a challenge of practical importance in applications varying from computer vision to physics and chemistry. We propose a novel multi-resolution convolutional architecture for learning over concentric spherical feature maps, of which the single sphere representation is a special case. Our hierarchical architecture is based on alternatively learning to incorporate both intra-sphere and inter-sphere information. We show the applicability of our method for two different types of 3D inputs, mesh objects, which can be regularly sampled, and point clouds, which are irregularly distributed. We also propose an efficient mapping of point clouds to concentric spherical images, thereby bridging spherical convolutions on grids with general point clouds. We demonstrate the effectiveness of our approach in improving state-of-the-art performance on 3D classification tasks with rotated data.

97 MATHEMATICS AND COMPUTING↗

Kepler Planet Detection Metrics: Per-Target Flux-Level Transit Injection Tests of TPS for Data Release 25

Quantifying the ability of a transiting planet survey to recover transit signals has commonly been accomplished through Monte-Carlo injection of transit signals into the observed data and subsequent running of the signal search algorithm (Gilliland et al., 2000; Weldrake et al., 2005; Burke et al., 2006). In order to characterize the performance of the Kepler pipeline (Twicken et al., 2016; Jenkins et al., 2017) on a sample of over 200,000 stars, two complementary injection and recovery tests are utilized:1. Injection of a single transit signal per target into the image or pixel-level data, hereafter referred to as pixel-level transit injection (PLTI), with subsequent processing through the Photometric Analysis (PA), Presearch Data Conditioning (PDC), Transiting Planet Search (TPS), and Data Validation (DV) modules of the Kepler pipeline. The PLTI quantification of the Kepler pipeline's completeness has been described previously by Christiansen et al. (2015, 2016); the completeness of the final SOC 9.3 Kepler pipeline acting on the Data Release 25 (DR25) light curves is described by Christiansen (2017).2. Injection of multiple transit signals per target into the normalized flux time series data with a subsequent transit search using a stream-lined version of the Transiting Planet Search (TPS) module. This test, hereafter referred to as flux-level transit injection (FLTI), is the subject of this document. By running a heavily modified version of TPS, FLTI is able to perform many injections on selected targets and determine in some detail which injected signals are recoverable. Significant numerical efficiency gains are enabled by precomputing the data conditioning steps at the onset of TPS and limiting the search parameter space (i.e., orbital period, transit duration, and ephemeris zero-point) to a small region around each injected transit signal.The PLTI test has the advantage that it follows transit signals through all processing steps of the Kepler pipeline, and the recovered signals can be further classified as planet candidates or false positives in the exact same manner as detections from the nominal (i.e., observed) pipeline run (Twicken et al., 2016, Thompson et al., in preparation). To date, the PLTI test has been the standard means of measuring pipeline completeness averaged over large samples of targets (Christiansen et al., 2015, 2016; Christiansen, 2017). However, since the PLTI test uses only one injection per target, it does not elucidate individual-target variations in pipeline completeness due to differences in stellar properties or astrophysical variability. Thus, we developed the FLTI test to provide a numerically efficient way to fully map individual targets and explore the performance of the pipeline in greater detail. The FLTI tests thereby allow a thorough validation of the pipeline completeness models (such as window function (Burke and Catanzarite, 2017a), detection efficiency (Burke Catanzarite, 2017b), etc.) across the spectrum of Kepler targets (i.e., various astrophysical phenomena and differences in instrumental noise). Tests during development of the FLTI capability revealed that there are significant target-to-target variations in the detection efficiency.

DR25↗

Nonlinear dimension reduction for surrogate modeling using gradient information

We introduce a method for the nonlinear dimension reduction of a high-dimensional function $u:{\mathbb{R}}^d\rightarrow{\mathbb{R}}$, $d\gg 1$. Our objective is to identify a nonlinear feature map $g:{\mathbb{R}}^d\rightarrow{\mathbb{R}}^m$, with a prescribed intermediate dimension $m\ll d$, so that $u$ can be well approximated by $f\circ g$ for some profile function $f:{\mathbb{R}}^m\rightarrow{\mathbb{R}}$. We propose to build the feature map by aligning the Jacobian $\nabla g$ with the gradient $\nabla u$, and we theoretically analyze the properties of the resulting $g$. Once $g$ is built, we construct $f$ by solving a gradient-enhanced least squares problem. Our practical algorithm uses a sample $\{{\textbf{x}}^{(i)},u({\textbf{x}}^{(i)}),\nabla u({\textbf{x}}^{(i)})\}_{i=1}^N$ and builds both $g$ and $f$ on adaptive downward-closed polynomial spaces, using cross validation to avoid overfitting. We numerically evaluate the performance of our algorithm across different benchmarks, and explore the impact of the intermediate dimension $m$. We show that building a nonlinear feature map $g$ can permit more accurate approximation of $u$ than a linear $g$, for the same input data set.

97 MATHEMATICS AND COMPUTING↗

Forecast of Wildfire Potential Across California USA Using a Transformer

Wildfires are a major issue facing the United States, a matter further exacerbated by an ever-changing climate. In California alone, wildfires are responsible for billions of dollars in damages and take lives each year. Accurately predicting fire danger conditions allows preparation awareness before wildfires start. Transformers are a class of deep learning models designed to identify patterns in sequential datasets. In recent years, transformers have gained popularity through their impressive performance in natural language processing and other applications of signal recognition. This analysis demonstrates the ability of a transformer with a residual connection to forecast fire danger potential over the state of California. Wildland fire potential index (WFPI) maps collected from the US Geological Survey database from January 1st 2020 to December 31st 2023 were used to tune, train and evaluate the transformer. Meteorological inputs (provided by Daymet daily weather and climatological summaries), the normalized difference vegetation index (NDVI) (calculated from the Moderate Resolution Imaging Spectroradiometer (MODIS)), and outputs from the Scott and Burgman fire behavior fuel models (to characterize maps of fuel types), were used as inputs. Our results show that a transformer can effectively emulate the US Forest Service modeled WFPI maps of California USA for four week long forecasts over the month of July, 2023, with correlations ranging from 0.85 – 0.98.

Limber, Russell [ORNL]↗