Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sparse Data”

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,261 records · Page 70

Reproduction drives changes in space use and habitat selection in a highly adaptable invasive mammal

Abstract For ungulates, it is not well understood how the interaction between habitat and reproduction affects movement behavior, space use, and habitat selection. We used known parturition (farrowing) data to validate First Passage Time (FPT) movement analysis to identify wild pig (Sus scrofa) parturition events from data collected with GPS collars. We examined home range and habitat selection during three physiologically distinct time periods in the reproductive cycle: late-gestation; farrowing; and neonate-care period. Sows exhibited a behavioral change 1-2 days prior to farrowing, suggesting FPT analysis can reliably identify parturition events for wide-ranging species that have a focused birthing area. Home range was smallest during the farrowing period, likely reflective of nest building, parturition, and protection of neonates. Home range size during the neonate-care period was intermediate between the late-gestation and farrowing periods, indicative of offspring care that may restrict maternal movement. Across all periods, sows avoided developed areas that have sparse canopy and ground cover and are associated with human activities. During late-gestation, sows avoided pine forests that have an open understory and less vegetative cover. During late-gestation and neonate-care, sows selected bottomland hardwood forests, habitat associated with ample food, cover, and water. During farrowing and neonate-care periods, sows selected upland hardwood forests, habitat with high quality food and cover for neonates. The physiological requirements of reproduction drive female habitat selection and spatial scale of movement patterns associated with parturition. Our study contributes to delineation of the appropriate scale at which to analyze movement data to provide insight about where individuals chose to place their home range, how much space to use, and how they use resources on the landscape to maximize reproductive success and fitness.

Zoology↗

Development of a Global Reference Surface Reflectance and BRDF Datasets from Geostationary Satellite Observations and AERONET Measurements

Surface reflectances and their dependency on illumination-view geometries (i.e., BRDF) are the foundation of many high-level satellite products for land and water monitoring. Yet it is difficult to evaluate the quality of satellite-based surface reflectances with ground-based measurements due to the spatial scale differences. In order to fill the gap, here we develop a reference dataset of surface reflectance and BRDF at the global AERONET sites with data streams from operational geostationary sensors including Himawari 8/9 AHI, GK-2A AMI, and GOES 16/17 ABI. Taking the top-of-atmosphere (TOA) reflectance and the site measured atmospheric aerosol optical depth (AOD) as the main inputs, we apply the GeoNEX-AC algorithm to performance accurate atmospheric correction and derive 10-minute surface reflectance and daily Ross-Thick-Li-Sparse (RTLS) BRDF parameters at AERONET sites where coincident AOD measurements and TOA observations are available from 2016 (for Himawari) or 2018 (for GOES) onwards. The algorithm ensures that the retrieved surface BRDF parameters, along with the site-measured AOD, allow the atmospheric radiative transfer model, SHARM, accurately simulate the observed TOA reflectance at diurnal and longer time scales. They are our best estimates of the surface optical properties and thus can serve as the “reference” to evaluate the performance of operational atmospheric correction algorithms (where AOD is assumed unknown and needs to be retrieved). The reference BRDF also allow us to evaluate the spectral band ratios between the SWIR (e.g., 2200 nm) and the visible (e.g., 650 nm) regions, which are commonly used in operational atmospheric correction algorithms. Finally, we demonstrate that the reference dataset can be used to develop potential data synergies between different GEO satellites as well as GEO-LEO sensors.

Weile Wang↗

ChatMPI: LLM-Driven MPI Code Generation for HPC Workloads

The Message Passing Interface (MPI) standard plays a crucial role in enabling scientific applications for parallel computing and is an essential component in high-performance computing (HPC). However, implementing MPI code manually—especially applying a proper domain decomposition and communication pattern—is a challenging and error-prone task. We present ChatMPI, an AI assistant for MPI parallelization of sequential C codes. In our analysis, we focus on testing six essential HPC workloads, which are based on Basic Linear Algebra Subprograms levels 1, 2, and 3 as well as sparse, stencil, and iterative operations. We analyze the process of creating ChatMPI by using the ChatHPC library. This lightweight large language model (LLM)–based infrastructure enables HPC experts to efficiently create and supervise trustworthy AI capabilities for critical HPC software tasks. We study the data required for training (fine-tuning) ChatMPI to generate parallel codes that not only use MPI syntax correctly but also apply HPC techniques to reduce memory communication and maximize performance by using proper work decomposition. With a relatively small training dataset composed of a few dozen prompts and fewer than 15 minutes of fine-tuning on one node equipped with two NVIDIA H100 GPUs, ChatMPI elevates trustworthiness for MPI code generation of current LLMs (e.g., Code Llama, ChatGPT-4o and ChatGPT 5). Additionally, we evaluate the performance of the MPI codes generated by ChatMPI in comparison with the ones generated by ChatGPT-4o and ChatGPT-5. The codes generated by ChatMPI provide up to a 4 × boost in performance by using better problem decomposition, communication patterns, and HPC techniques (e.g., communication avoiding).

Valero Lara, Pedro [ORNL] (ORCID:0000000214794310)↗

Role of physics in physics-informed machine learning

Physical systems are characterized by inherent symmetries, one of which is encapsulated in the units of their parameters and system states. These symmetries enable a lossless order-reduction, e.g., via dimensional analysis based on the Buckingham theorem. Despite the latter's benefits, machine learning (ML) strategies for the discovery of constitutive laws seldom subject experimental and/or numerical data to dimensional analysis. We demonstrate the potential of dimensional analysis to significantly enhance the interpretability and generalizability of ML-discovered secondary laws. Our numerical experiments with creeping fluid flow past solid ellipsoids show how dimensional analysis enable both deep neural networks and sparse regression reproduce old results, e.g., Stokes law for a sphere, and generate new ones, e.g., an expression for an ellipsoid misaligned with the flow direction. Furthermore, our results suggest the need to incorporate other physics-based symmetries and invariances into ML-based techniques for equation discovery.

97 MATHEMATICS AND COMPUTING↗

Precision X-ray Timing of RX J0806.3+1527 with CHANDRA: Evidence for Gravitational Radiation from an Ultracompact Binary

RX J0806.3+1527 is a candidate double degenerate binary with possibly the shortest known orbital period. The source shows an approximately equal to 100% X-ray intensity modulation at the putative orbital frequency of 3.11 mHz (321.5 s). If the system is a detached, ultracompact binary gravitational radiation should drive spin-up with a magnitude of nu(sup dot) approximately 10(exp -16) Hz per second. Efforts to constrain the X-ray frequency evolution to date have met with mixed success, principally due to the sparseness of earlier observations. Here we describe the results of the first phase coherent X-ray monitoring campaign on RX J0806.3+1527 with Chandra. We obtained a total of 70 ksec of exposure in 6 epochs logarithmically spaced over 320 days. With these data we conclusively show that the X-ray frequency is increasing at a rate of 3.77 plus or minus 0.8 x 10(exp -16) Hz per second. Using the ephemeris derived from the new data we are able to phase up all the earlier Chandra and ROSAT data and show they are consistent with a constant nu(sup dot) = 3.63 plus or minus 0.06 x 10(exp -16) Hz per second over the past decade. This value appears consistent with that recently derived by Israel et al. largely from monitoring of the optical modulation, and is in rough agreement with the solutions reported initially by Hakala et al., based on ground-based optical observations. The large and stable nu(sup dot) over a decade is consistent with gravitational radiation losses driving the evolution. An intermediate polar (IP) scenario where the observed X-ray period is the spin period of an accreting white dwarf appears less tenable because the observed nu(sup dot) requires an m(sup dot) approximately equal to 4 x 10 (exp -8) solar mass yr(sup -l), that is much larger than that inferred from the observed X-ray luminosity (although this depends on the uncertain distance and bolometric corrections), and it is difficult to drive such a high m(sup dot) in a binary system with parameters consistent with all the multiwavelength data. If the ultracompact scenario is correct, then the X-ray flux cannot be powered by stable accretion which would drive the components apart, suggesting a new type of energy source (perhaps electromagnetic) may power the X-ray flux.

Strohymayer, Tod E.↗

Large Area Lunar Dust Flux Measurement Instrument

The instrument under development is designed to characterize the flux and size distribution of the lunar micrometeoroid and secondary ejecta environment. When deployed on the lunar surface, the data collected will benefit fundamental lunar science as well as enabling more reliable impact risk assessments for human lunar exploration activities. To perform this task, the instrument requirements are demanding. It must have as large a surface area as possible to sample the very sparse population of the larger potentially damage-inducing micrometeorites. It must also have very high sensitivity to enable it to measure the flux of small (<10 micron) micrometeorite and secondary ejecta dust particles. To be delivered to the lunar surface, it must also be very low mass, rugged and stow compactly. The instrument designed to meet these requirements is called FOMIS. It is a large-area thin film under tension (i.e. a drum) with multiple fiber optic displacement (FOD) sensors to monitor displacements of the film. This sensor was chosen since it can measure displacements over a wide dynamic range: 1 cm to sub-Angstrom. A prototype system was successfully demonstrated using the hypervelocity impact test facility at the University of Kent (Canterbury, UK). Based on these results, the prototype system can detect hypervelocity (approx.5 km/s) impacts by particles as small as 2 microns diameter. Additional tests using slow speeds find that it can detect secondary ejecta particles (which do not penetrate the film) with momentums as small as 15 pico-gram 100m/s, or nominally 5 microns diameter at 100 m/s.

Corsaro, R.↗

Smc-Last Point Source Catalog Spitzer'S Last Survey of the Small Magellanic Cloud

We present a processing for the point source catalog derived from Spitzer's final two surveys of the Small Magellanic Cloud (SMC). These 2017-2018 maps duplicate the footprint from the SAGE-SMC program from 2008 and were obtained at 3.6 and 4.5 μm with the Infrared Array Camera (IRAC). They provide a nearly 10-year baseline of all observations of the SMC made by Spitzer. We extracted point sources from both the basic calibrated data (BCD) as well as the final processed mosaics. We did not use any prior catalogs as inputs in order to search for moving objects (e.g., foreground brown dwarfs) and other transient phenomena. The extractions from two test fields are shown here: a field sparsely populated with stars and a second containing more complicated structures. Comparing our extractions between both of our epochs and with the SAGE epochs shows the completeness and range of our survey. The mosaics are currently available on the Infrared Science Archive (IRSA) and the Harvard Dataverse. The final point source catalog will be available at both archives, as well. This work is supported by NASA Award 80NSSC19K0585.

Thomas Kuchar↗

Geothermal Play-Fairway Analysis of Washington State Prospects: Final Report

The Washington State Geothermal Play-Fairway Analysis overcomes the exploration challenges posed by dense vegetation, glacial deposits, and extreme precipitation. The geothermal play-fairways we target are locations where heat, permeability, and saturated porosity are present in sufficient volume to provide adequate heat exchange at depths accessible by modern drilling technology. The three study areas lie along the Cascade Range magmatic arc and are near Mount Baker, Mount St. Helens, and the Wind River Valley. The seven-year project is divided into three phases. In Phase 1 we build on a previous statewide assessment of geothermal resources and develop an initial modeling approach. The results are a series of favorability, uncertainty, and risk maps for three targeted study areas. Based on these initial results, we collect new geologic and geophysical data to further refine our modeling and reduce exploration uncertainty in Phase 2. We improve the modeling method to handle the new data and update the favorability, uncertainty, and risk maps. We also update the conceptual geothermal resource models. In Phase 3 we validate our modeling approach by drilling two temperature-gradient holes and collecting and analyzing core, image logs, and new geochemistry. Our modeling approach improves on an earlier statewide method through a more-rigorous and detailed assessment of heat and permeability. Permeability potential is assessed through geomechanical modeling of the deformation that can generate and maintain reservoir porosity and permeability. Metrics to inform heat potential include temperature-gradient wells, which are sparse in Washington; proximity of Quaternary volcanic vents and young intrusive rock; spring temperature; and reservoir temperature inferred from geothermometry. We weight the individual components using an expert-guided approach known as the Analytical Hierarchy Process. During Phase 2 we also develop a fluid-filled fracture model, and an infrastructure model that helps to delineate areas which are more favorable for geothermal development based on proximity to transmission lines, elevation, land ownership and use restrictions, and availability of process water. New geologic and geophysical data is collected during Phase 2 in each of our three main study areas. At Mount Baker and north of Mount St. Helens we conduct 1:24,000-scale geologic mapping and lidar analysis to better constrain the location and character of surface faults; detailed mapping in the Wind River Valley was completed just prior to the start of this project. Ages of intrusive rocks are determined with 40 Ar/ 39 Ar geochronology, though all of our samples are Miocene or older. We collect ground based gravity observations (a total of 1,580 new stations) in all of our study areas and ground-based magnetic lines (a total of 93 km) at Mount Baker. These data are combined with existing gravity and aeromagnetic data and used to constrain fault locations and geometry. Two to three cross sections are constructed at each study area using the mapped surface geology and forward-modeling of the gravity and magnetic data; these cross sections form the basis for our updated conceptual models. We collect magnetotelluric surveys at Mount Baker and Mount St. Helens and these data are inverted to form a resistivity model from the surface to about 10 km depth; each model shows conductive zones that can be interpreted as upwelling geothermal fluids. At Mount St. Helens we deploy a passive seismic array and use the newly detected events to refine the location of the Saint Helens seismic zone. We also employ ambient-noise tomography to develop a detailed seismic-velocity model for the study area and use this model to help constrain our cross sections and conceptual model. Based on the new data collected during Phase 2—and our updated models—we develop a campaign of temperature-gradient holes and core analysis to validate our modeling in Phase 3. Drill hole MB76-31 is located near Little Park Creek, 11 km west-southwest of the summit of Mount Baker, and is 1,471 ft deep. About 410 ft of core from the lower portion of the hole—and image logs from ~175 ft below ground surface to the bottom—are collected and analyzed. Water samples are collected and processed for geothermometry. Drill hole MSH17-24 is located along upper Schultz Creek, 16 km north-northeast of Mount St. Helens and has core from 470 ft to the bottom at 1,053 ft. We did not collect image logs due to borehole stability concerns, but water samples are collected and analyzed for geothermometry. Repeat temperature-gradient measurements are made at both sites and thermal conductivity is measured from core samples. At MB76-31, the equilibrated temperature gradient of 64°C/km and calculated heat flow of 141–159 mW/m 2 is more than twice the regional average. Detailed mapping and analysis of the core, coupled with correlation to the image logs, indicates a history of permeability generation consistent with our predictions of high permeability. Because the site has high favorability in the Phase 2 model, we consider the results a positive validation of the modeling. At site MSH17-24, the equilibrated temperature gradient of ~15°C/km and calculated heat flow of 41–43 mW/m 2 are similar to regional. Geochemical analysis of the water samples indicates a meteoric source without any geothermal component. Detailed outcrop-based mapping of fault exposures near the drill site and analysis of image logs from nearby boreholes indicates a history of permeability generation consistent with our predictions. Because the site has low favorability in the Phase 2 model, we consider the results a positive validation of the modeling. Together, the two sites provide a reasonably positive validation of the Phase 2 modeling and should encourage future use of this modeling approach.

15 GEOTHERMAL ENERGY↗

Deep learning of free boundary and Stefan problems

Free boundary problems appear naturally in numerous areas of mathematics, science and engineering. These problems present a great computational challenge because they necessitate numerical methods that can yield an accurate approximation of free boundaries and complex dynamic interfaces. In this work, we propose a multi-network model based on physics-informed neural networks to tackle a general class of forward and inverse free boundary problems called Stefan problems. Specifically, we approximate the unknown solution as well as any moving boundaries by two deep neural networks. Besides, we formulate a new type of inverse Stefan problems that aim to reconstruct the solution and free boundaries directly from sparse and noisy measurements. We demonstrate the effectiveness of our approach in a series of benchmarks spanning different types of Stefan problems, and illustrate how the proposed framework can accurately recover solutions of partial differential equations with moving boundaries and dynamic interfaces. All code and data accompanying this manuscript are publicly available at https://github.com/PredictiveIntelligenceLab/DeepStefan.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Precision X-ray Timing of RX J0806.3+1527 with CHANDRA: Evidence for Gravitational Radiation from an Ultracompact Binary

RX J0806.3+1527 is a candidate double degenerate binary with possibly the shortest known orbital period. The source shows an $\approx 100 \%$ X-ray intensity modulation at the putative orbital frequency of 3.11 mHz (321.5 s). If the system is a detached, ultracompact binary gravitational radiation should drive spin-up with a magnitude of $\dot\nu \sim 10(exp -16)$ Hi s$(exp -l}$. Efforts to constrain the X-ray frequency evolution to date have met with mixed success, principally due to the sparseness of earlier observations. Here we describe the results of the first phase coherent X-ray monitoring campaign on RX J0806.3+1527 with {\it Chandra) . We obtained a total of 70 ksec of exposure in 6 epochs logarithmically spaced over 320 days. With these data we conclusively show that the X-ray frequency is increasing at a rate of $3.77 \pm\ 0.8 \times l0(exp -16)$ Hz s$(exp -l)$. Using the ephemeris derived from the new data we are able to phase up all the earlier {\it Chandra) and ROSAT data and show they are consistent with a constant $\dot\nu = 3.63 \pm 0.06 \times 10(exp -16)$ Hz s$(exp -1)s over the past decade. This value appears consistent with that recently derived by Israel et al. largely from monitoring of the optical modulation, and is in rough agreement with the solutions reported initially by Hakala et al., based on ground-based optical observations. The large and stable $\dot\nu$ over a decade is consistent with gravitational radiation losses driving the evolution. An intermediate polar (IP) scenario where the observed X-ray period is the spin period of an accreting white dwarf appears less tenable because the observed $\dot\nu$ requires an $\dot m \approx 2 \times 10(exp -8)$ $M-{\odot)$ yr$(exp -l)$, that is much larger than that inferred from the observed X-ray luminosity (although this depends on the uncertain distance and bolometric corrections), and it is difficult to drive such a high $\dot m$ in a binary system with parameters consistent with all the multiwavelength data. If the ultracompact scenario is

Strohmayer, Tod↗

Lunar Magnetic Field Models From Lunar Prospector and SELENE/Kaguya Along-Track Magnetic Field Gradients

We use L1-norm model regularization of |Br| component at the surface on magnetic monopoles bases and along-track magnetic field differences alone (without vector observations) to derive high quality global magnetic field models at the surface of the Moon. The practical advantages to this strategy are the following: monopoles are more stable at closer spacing in comparison to dipoles, improving spatial resolution; L1-norm model regularization leads to sparse models which may be appropriate for the Moon which has regions of localized magnetic field features; and along-track differences reduce the need for ad-hoc external field noise reduction strategies. We examine also the use of Lunar Prospector and SELENE/Kaguya magnetometer data, combined and separately, and find that the Lunar Prospector along-track vector field differences lead to surface field models that require weaker regularization and, hence, result in higher spatial resolution. Significantly higher spatial resolution (wavelengths of roughly 25–30 km) and higher amplitude surface magnetic fields can be derived over localized regions of high amplitude anomalies (due to their higher signal-to-noise ratio). These high-resolution field models are also compared with the results of Surface Vector Mapping approach of Tsunakawa et al. (2015, https://doi.org/10.1002/2014JE004785). Finally, the monopoles- as well as dipoles-based patterns of the Serenitatis high amplitude magnetic feature have characteristic textbook patterns of Br and Bθ component fields from a nearly vertically downwardly magnetized source region and it implies that the principal source of the anomaly was formed when the region was much closer to the north magnetic pole of the Moon.

D. Ravat↗

Midwest Water Resources: Developing an Evapotranspiration Climatology to Analyze Spatiotemporal Water Budget Patterns for Agriculture and Natural Resources Managers in the Midwest

Evapotranspiration (ET) is a climatic variable critical to the hydrologic cycle and is used to evaluate spatiotemporal trends in drought conditions. Although in-situ observations provide accurate ET information, these records are spatially sparse. The United States Department of Agriculture (USDA) Midwest Climate Hub, National Integrated Drought Information System, Minnesota Department of Agriculture, Michigan State University, and the USDA Foreign Agricultural Service have partnered with DEVELOP to gain new insights on spatiotemporal patterns of ET with NASA satellite data. This project evaluates the feasibility of using remotely sensed ET data products to understand trends from 2001 through 2020. Actual ET (aET) data were sourced from NASA’s Terra Moderate Resolution Imaging Spectroradiometer (MODIS), and reference ET (refET) and precipitation data were sourced from the Gridded Surface Meteorological Dataset (gridMET). Spatial and temporal differences in the MODIS and gridMET data sets were resolved such that gridded data could be compared. These data sets were then used to produce monthly Normals maps of mm/8-day of precipitation, aET, and refET. Using precipitation and ET Normals maps, hydrologic state maps were produced by subtracting either refET or aET from precipitation. These hydrologic state maps provide a proxy water balance by summarizing the difference between water entering and leaving the ground surface. Additionally, the ET products were compared in timeseries plots. The ET products from this project provided partners with comparable datasets to assess potential drought and flooding conditions to support Midwest agricultural and natural resource managers in decision-making.

Emma Myrick↗

Joint Channel Equalization and Symbol Detection for IoT Devices in Severe Multipath Channels

Internet of Things (IoT) devices often transmit very short message blocks without strong error correction coding or pilots for channel equalization. When the transmitted signals encounter a severe multipath channel, the receiver is often unable to equalize the resulting inter-symbol interference (ISI) via traditional methods, leading to a high retransmission rate. This paper proposes a joint channel equalization and symbol detection scheme for this scenario by utilizing the preamble sequence as a training pilot for equalization and the relatively weak cyclic redundancy check (CRC) 8-Dallas code for error correction. By assuming a sparse channel impulse response, the proposed joint channel equalization and symbol detection scheme formulates an l1 norm constrained optimization problem and solves it by a saddle-point algorithm. The proposed algorithm is applied to a set of real-world underwater animal tracking IoT data, and the results show improvement of correct message detection rate from 21:5% to 41:6%.

Jing, Shusen↗

The North Atlantic circulation in the early 1980s - An estimate from inversion of a finite-difference model

A finite-difference model of the North Atlantic is constructed for the purpose of making an estimate of the circulation through an inverse calculation. The data base is eclectic, and includes hydrography, oxygen, nutrients, current meter and float records, atmospheric momentum, heat and water vapor transfers, as well as estimates of certain integral fluxes. Owing to the available hydrographic database, the model resolution is restricted to 1 deg at best, and is much coarser in many aspects. This limited resolution is a major obstacle to accurate estimates of climatological fluxes. In its final form, there are about 9000 constraints in 29,000 formal unknowns plus 9000 noise unknowns. The system is solved as a tapered least-squares system by a sparse conjugate gradient algorithm. With the exception of a few float velocities, all constraints are found to be consistent within error estimates. The model produces estimates of large-scale fluxes and flux divergences for all conventional properties including heat and nutrients as well as carbon dioxide and alkalinity. Meridional fluxes of carbon are found to be indistinguishable from zero, whereas the North Atlantic tends to export nutrients to the south, but carry heat to the north. Traditional oceanographic depictions of the circulation through combination of nonsynoptic data into steady models may have reached their useful limit in the present calculation, as the conflicts between the data and physical requirements become quantitatively apparent.

Martel, Francoise↗

Investigating the Impact of UV Radiation on High-Altitude Shallow Lake Habitats, Life Diversity, and Life Survival Strategies: Clues for Mars' Past Habitability Potential?

We present data and results from an ongoing project of astrobiological high-altitude expeditions investigating the highest and least explored perennial lakes on Earth in the Bolivian and Chilean Andes, including several volcanic crater lakes nearing and beyond 6,000 m in elevation. In the next five years, they will provide the first integrated long-term astrobiological characterization and monitoring of lacustrine environments and their biology for such altitude. These extreme lakes are natural laboratories. They provide the field data missing beyond 4,000 m to complete our understanding of terrestrial lakes and biota. Research on the effects of UV has been performed in lower altitude lakes and models of UV flux over time are being developed. Lakes showing a high content of dissolved organic material (DOM) shield organisms from UV. DOM acts as a natural sunscreen as it influences the water transparency, therefore is a determinant of photic zone depth. In sparsely vegetated alpine areas, lakes are clearer and offer less protection from UV to organisms living in the water. Transparent water and high UV irradiance may maximize the penetration and effect of UV radiation. Shallow-water communities in these lakes are particularly sensitive to UV radiation. The periphyton can live on various susbtrates. While on rocks, it includes immobile species that cannot seek low UV refuges unlike sediment-dwelling periphyton or alpine phytoflagellates which undergo vertical migration. Inhibition of algal photosynthesis by UV radiation has been documented in laboratory and showed that phytoplankton production is reduced by formation of nucleic acid lesions or production of peroxides and free oxygen radicals. of peroxides and free oxygen radicals. Our project is providing the field data that is missing from natural laboratories beyond 4,000 m and will complement the vision of the effects of UV on life and its adaptation modes (or lack thereof).

Cabrol, A.↗

DS-GL: Advancing Graph Learning via Harnessing the Power of Nature within Dynamic Systems

With the rapid digitization of the world, an increasing number of real-world applications are turning to nonEuclidean data, modeled as graphs. Due to their intrinsic high complexity and irregularity, learning from graph data demands tremendous computational power. Recently, CMOS-compatible Ising machines, i.e., dynamic systems composed of CMOS components, have emerged as a new approach that harnesses the inherent power of natural annealing within dynamic systems to efficiently resolve binary optimization problems and have been adopted for traditional graph computation, such as max-cut. However, when performing complex Graph Learning (GL) tasks, Ising machines face significant hurdles: (i) they are inherently binary and thus ill-suited for real-valued problems; (ii) their expensive all-to-all coupling network that guarantees effective natural annealing poses daunting scalability concerns. To address these challenges, this paper proposes a nature-powered graph learning framework dubbed DS-GL, which is the first effort to transform the process of solving graph learning problems into the natural annealing process within a parameterized dynamic system embodied as a CMOS chip. To tackle the two major hurdles, DS-GL first augments the Ising machine architecture to modify the self-reaction term of its Hamiltonian function from linear to quadratic, effectively serving as an energy regulator. This adjustment maintains the system’s original physical interpretation while enabling it to process continuous, real-valued data. Second, to address the scaling issue, DS-GL further upgrades the real-valued dense Ising machine by decomposing it into a mesh-based multi-PE dynamic system that supports efficient distributed spatial-temporal co-annealing across different PEs through sparse interconnects. By exploiting the inherent sparsity and component structures in real-world graphs, DS-GL is able to map complex graph learning tasks onto the scalable dynamic system while maintaining high accuracy. Evaluations with three diverse GL applications across six real-world datasets, including traffic flow and COVID-19 prediction, show that DS-GL can deliver from 102× to 106× speedups and 500× energy reduction over Graph Neural Networks on GPUs, with 5% - 20% accuracy enhancement.

Song, Ruibing↗

Modeling the radiant transfers of sparse vegetation canopies

The scattering dynamics of sparse vegetation canopies were studied within the framework of the three-dimensional radiative transfer model of Kimes (1984). The model was upgraded by including an algorithm for the anisotropic scattering of a soil boundary. Validation of the model was carried out using measured directional reflectance data for two canopies exhibiting typical scattering behavior with low and intermediate vegetation density. The canopies were: an orchard grass canopy; and a hard wheat canopy. A number of factors were found contributing to the final reflectance distribution of the canopies, including: (1) the strong anisotropic scattering properties of the soil; (2) the geometric effect of the vegetation probability gap function on the soil anisotropy and solar irradiance; and (3) the anisotropic scattering of vegetation which is controlled by the phase function and the layering of leaves. The application of the theoretical results to the development of earth-observing sensor systems is discussed.

Kimes, D. S.↗

Extreme sparsification of physics-augmented neural networks for interpretable model discovery in mechanics

Data-driven constitutive modeling with neural networks has received increased interest in recent years due to its ability to easily incorporate physical and mechanistic constraints and to overcome the challenging and time-consuming task of formulating phenomenological constitutive laws that can accurately capture the observed material response. However, even though neural network-based constitutive laws have been shown to generalize proficiently, the generated representations are not easily interpretable due to their high number of trainable parameters. Sparse regression approaches exist that allow for obtaining interpretable expressions, but the user is tasked with creating a library of model forms which by construction limits their expressiveness to the functional forms provided in the libraries. Here, in this work, we propose to train regularized physics-augmented neural network-based constitutive models utilizing a smoothed version of $L^0$-regularization. This aims to maintain the trustworthiness inherited by the physical constraints, but also enables interpretability which has not been possible thus far on any type of machine learning-based constitutive model where model forms were not assumed a priori but were actually discovered. During the training process, the network simultaneously fits the training data and penalizes the number of active parameters, while also ensuring constitutive constraints such as thermodynamic consistency. We show that the method can reliably obtain interpretable and trustworthy constitutive models for compressible and incompressible hyperelasticity, yield functions, and hardening models for elastoplasticity, using synthetic and experimental data. This work aims to set a new paradigm for interpretable machine learning models in the broad area of solid mechanics where low and limited data is available along with prior knowledge of physical constraints that the learned maps need to obey. This paradigm can potentially be extended to a broader spectrum of scientific exploration.

Data-driven constitutive models↗