Perspectives on the Future of Land Surface Models and the Challenges of Representing Complex Terrestrial Systems
Not Available
SEARCH · Search NASA
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.
Not Available
In this work we employed machine learning-augmented density functional theory (DFT) thermodynamic calculations to assess the stability of different AgO x structures under catalytic ethylene epoxidation reaction conditions. We found that there are multiple AgO x surface motifs that could co-exist under the relevant conditions. These included Ag surface oxides (e.g., AgO_p(4 × 4) and Ag 1.83 O) and atomic oxygen-covered Ag(111) surfaces. Furthermore, we employed DFT calculations to evaluate the energetics of different reaction mechanisms by which ethylene and oxygen can react on these surfaces. These studies revealed several energetically viable reaction pathways for ethylene epoxidation. Microkinetic modeling analyses, based on the DFT-calculated reaction pathways, showed that ethylene epoxidation can proceed on all surfaces and that multiple pathways, including those involving Langmuir–Hinshelwood and Eley–Rideal mechanisms, could be involved in selective and unselective reactions. The diversity of mechanisms that we discovered in the context of the relatively simple ethylene epoxidation reaction on Ag suggests that the richness and complexity of surface chemistry are most likely a rule rather than an exception in heterogeneous catalytic chemical transformations on metal surfaces and that the concept of a single or even a dominant mechanism and reaction intermediates might need to be revisited for many reactions.
Dynamic nuclear polarization surface-enhanced nuclear magnetic resonance (NMR) spectroscopy has enabled the determination of the three-dimensional configuration of surface sites, in particular supported metal complexes of relevance to single-site heterogeneous catalysis. These approaches have chiefly leveraged the application of NMR double-resonance experiments that either reveal the complex conformation via point-to-point intramolecular distances between spin-labeled atoms or the complex-surface orientation via distances between the spins and the surface plane. Either method typically requires expensive isotope labeling and each reports on different structural features. The application of an experiment that simultaneously reveals both types of distances with chemical resolution would be ideal. Here, in this article, we describe an 17 O{ 1 H} pseudo-3D correlation experiment that achieves this goal. Specifically, Si–O–Si and Si–O–M oxygens are well-resolved by 17 O NMR; therefore, distances can be simultaneously measured radially, between Si– 17 O–M and the 1 H’s of the ligands, and vertically to the Si– 17 O–Si linkages of the silica support. We demonstrate the experiment using supported yttrium and zirconium complexes. Good agreement is obtained when comparing the experimental results to theoretical predictions from density functional theory calculations, highlighting the reliability of this relatively simple experiment.
Surface crack detection and dimensional measurement at active mining sites present significant safety and operational challenges. Manual inspection methods are labor-intensive, spatially incomplete, and expose personnel to hazardous environments, while existing automated approaches have been developed primarily for concrete civil infrastructure and have not been validated on the complex, variable surfaces characteristic of mining environments. This dissertation presents an automated pipeline that integrates deep learning semantic segmentation with Structure-from-Motion photogrammetry to detect surface cracks and measure their aperture, length, and vertical displacement from standard RGB imagery acquired during routine Uncrewed Aerial Vehicle (UAV) survey operations, without requiring additional sensor hardware or manual measurement. The pipeline combines a U-Net architecture with an EfficientNet-B0 encoder, pretrained on the SDNET2018 concrete crack dataset and fine-tuned on a mining-specific dataset spanning laboratory concrete specimens, coal refuse impoundment embankments, and post-blast limestone quarry benches. Photogrammetric reconstruction is performed using COLMAP Structure-from-Motion and Multi-View Stereo, with crack segmentation masks projected into the reconstructed point cloud to enable three-dimensional vertical displacement measurement through local plane fitting and bimodal surface detection. The pipeline was validated across 36 controlled laboratory specimens at three imaging distances and four vertical displacement levels, achieving aperture measurement RMSE of 0.047 cm and R² of 0.954, and vertical displacement RMSE of 0.140 cm and R² of 0.966, against independent caliper measurements. Field application at a coal refuse impoundment in southwestern Pennsylvania detected 71 crack components across the embankment crest, with a dominant longitudinal crack exhibiting aperture values reaching 28 cm and a 95th percentile vertical displacement of 35.53 cm, consistent in magnitude and spatial distribution with simultaneously acquired LiDAR-derived estimates. Application across four post-blast limestone quarry bench datasets in California successfully characterized blast-induced fracture networks at ground sampling distances ranging from 0.59 to 1.23 cm/pixel, with detected crack geometries physically consistent with observable surface conditions at each site. The results demonstrate that deep learning-based crack detection and photogrammetric measurement can be integrated into routine UAV inspection workflows at mining sites, providing repeatable, scalable, and quantitative crack characterization across surface types, crack scales, and displacement magnitudes not previously addressed in the literature. The pipeline requires no dedicated surveying equipment beyond the UAV platforms already deployed at mine sites for survey and monitoring purposes, supporting practical adoption within existing operational workflows.
Lunar floor-fractured craters are primarily located near the maria and frequently contain ponded mare units and dark mantling deposits. Fracturing is confined to the crater interior, often producing a moat-like feature near the floor edge, and crater depth is commonly reduced by uplift of the crater floor. Although viscous relaxation of crater topography can produce such uplift, the close association of modification with surface volcanism supports a model linking floor fracture to crater-centered igneous intrusions. The consequences of two intrusion models for the lunar interior are quantitatively explored. The first model is based on terrestrial laccoliths and describes a shallow intrusion beneath the crater. The second model is based on cone sheet complexes where surface deformation results from a deeper magma chamber. Both models, their fit to observed crater modifications and possible implications for local volcanism are described.
We thank G.P. Mikitik and Yu.V. Sharlai for contributing this note and the cordial exchange about it. First and foremost, we note that the aim of our paper is to report a methodology to diagnose topological (semi)metals using magnetic quantum oscillations. Thus far, such diagnosis has been based on the phase offset of quantum oscillations, which is extracted from a “Landau fan plot”. A thorough analysis of the Onsager–Lifshitz–Roth quantization rules has shown that the famous π-phase shift can equally well arise from orbital or spin magnetic moments in topologically trivial systems with strong spin-orbit coupling or small effective masses. Therefore, the “Landau fan plot” does not by itself constitute a proof of a topologically nontrivial Fermi surface. In the paper at hand, we report an improved analysis method that exploits the strong energy dependence of the effective mass in linearly dispersing bands. This leads to a characteristic temperature dependence of the oscillation frequency which is a strong indicator of nontrivial topology, even for multi-band metals with complex Fermi surfaces. Three materials, Cd 3 As 2 , Bi 2 O 2 Se and LaRhIn5 served as test cases for this method. Linear band dispersions were detected for Cd 3 As 2 , as well as the F ≈ 7 T pocket in LaRhIn 5 .
We uncover a surface-mediated mechanism for visible-light-driven photocatalytic desulfurization of thiols on titanium dioxide (TiO 2 ) nanoparticles under LED illumination (>420 nm). Combined diffuse reflectance, infrared, and cryogenic electron paramagnetic resonance spectroscopy verify the formation of a surface thiol–Ti charge-transfer complex that introduces visible-light absorption into an otherwise UV-active reaction system. Selective photoexcitation of this surface complex drives interfacial charge transfer, directly activating the adsorbed thiol and inducing C–S bond scission without the band gap excitation of TiO 2 . The resulting carbon-centered radical intermediates are intercepted by styrene to form C–C coupling products. These results uncover a distinct surface-mediated photocatalytic mechanism in which charge-transfer complex formation governs both light absorption and chemical activation, providing a mechanistic framework for extending visible-light reactivity to wide-bandgap metal oxide photocatalysts.
High-resolution space-based spectral imaging of the Earth's surface delivers critical information for monitoring changes in the Earth system as well as resource management and utilization. Orbiting spectrometers are built according to multiple design parameters, including ground sampling distance (GSD), spectral resolution, temporal resolution, and signal-to-noise ratio. Different applications drive divergent instrument designs, so optimization for wide-reaching missions is complex. The Surface Biology and Geology component of NASA's Earth System Observatory addresses science questions and meets applications needs across diverse fields, including terrestrial and aquatic ecosystems, natural disasters, and the cryosphere. The algorithms required to generate the geophysical variables from the observed spectral imagery each have their own inherent dependencies and sensitivities, and weighting these objectively is challenging. Here, we introduce intrinsic dimensionality (ID), a measure of information content, as an applications-agnostic, data-driven metric to quantify performance sensitivity to various design parameters. ID is computed through the analysis of the eigenvalues of the image covariance matrix, and can be thought of as the number of significant principal components. This metric is extremely powerful for quantifying the information content in high-dimensional data, such as spectrally resolved radiances and their changes over space and time. We find that the ID decreases for coarser GSD, decreased spectral resolution and range, less frequent acquisitions, and lower signal-to-noise levels. This decrease in information content has implications for all derived products. ID is simple to compute, providing a single quantitative standard to evaluate combinations of design parameters, irrespective of higher-level algorithms, products, applications, or disciplines.
The UTD-based NEWAIR3 code is modified such that it can model modern aircraft by composite plates. One good model of conductor-backed composites is the impedance boundary condition where the composites are replaced by surfaces with complex impedances. This impedance-plate model is then used to model the composite plates in the NEWAIR3 code. In most applications, the aircraft distorts the desired radiation pattern of the antenna. However, test examples conducted in this report have shown that the undesired scattered fields are minimized if the right impedance values are chosen for the surface impedance plates.
Airless planetary bodies are subject to space weathering processes such as micrometeoroid bombardment and solar wind irradiation. These processes alter the microstructural, chemical, and optical properties of planetary regoliths and, in turn, complicate interpretations of surface composition from remote sensing data and matching of meteorites to their parent bodies. Unlike lunar and ordinary chondrite materials, space weathering studies of primitive carbonaceous chondrites are still evolving [1]. Investigations that examine the spectral, microstructural, and chemical effects of space weathering on carbonaceous materials are critical for maximizing the scientific return of ongoing and future missions to C-complex asteroids, like the NASA OSIRIS-REx (Bennu) and JAXA Hayabusa2 (Ryugu) missions[2,3]. Here, we present results from our experimental simulation of solar wind weathering on organic-rich, C-complex asteroid surfaces in which we irradiated a slab of the Murchison CM2 chondrite with 1 keV/amu H+ and He+ ions.
Surfaces of normal, cultured, and mitogen-stimulated mouse lymphoid cells were examined by scanning electron microscopy (SEM). Lymphocytes with smooth, highly villous and intermediate surfaces were observed in cell suspensions from both spleens and thymuses of normal mice and from spleens of congenitally athymic (nude) mice. Several strain-specific surface features were noted, including the spine-like appearance of microvilli on C57B1/6 lymphocytes. Although thymus cell suspensions contained somewhat more smooth cells than did spleen cell preparations, lymphocyte derivation could not be inferred from SEM examination. Studies of cells stimulated with mitogenic agents for thymus-derived lymphocytes (concanavalin A) or for bone marrow-derived lymphocytes (lipopolysaccharide) suggested that, in the mouse, development of a complex villous surface is a general concomitant of lymphocyte activation and transformation.
The initial condition problem for a binary neutron star system requires a Poisson equation solver for the velocity potential with a Neumann-like boundary condition on the surface of the star. Difficulties that arise in this boundary value problem are: (a) the boundary is not known a priori, but constitutes part of the solution of the problem; (b) various terms become singular at the boundary. In this work, we present a new method to solve the fluid Poisson equation for irrotational/spinning binary neutron stars. The advantage of the new method is that it does not require complex fluid surface fitted coordinates and it can be implemented in a Cartesian grid, which is a standard choice in numerical relativity calculations. This is accomplished by employing the source term method proposed by Towers, where the boundary condition is treated as a jump condition and is incorporated as additional source terms in the Poisson equation, which is then solved iteratively. The issue of singular terms caused by vanishing density on the surface is resolved with an additional separation that shifts the computation boundary to the interior of the star. We present two-dimensional tests to show the convergence of the source term method, and we further apply this solver to a realistic three-dimensional binary neutron star problem. By comparing our solution with the one coming from the initial data solver cocal, we demonstrate agreement to approximately 1%. We report our method can be used in other problems with non-smooth solutions like in magnetized neutron stars.
A large-scale, fully resolved direct simulation Monte Carlo (DSMC) computation of a non-equilibrium, reactive flow of pure oxygen over a double cone is presented. Under the simulated near-continuum conditions, the computational demands are shown to be significant because of the wide range of length scales that must be resolved. Therefore, robust grid adaption capabilities and efficient parallelization of the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code that is utilized in this work are essential. The thermochemical and transport collision models were selected for efficiency and simplicity. First-principles data, obtained from the highly accurate direct molecular simulation method, were used to inform the collision models’ parameters. Importantly, because SPARTA implements molecular collision models using collision-specific energies, the resulting macroscopic relaxation rates were evaluated a posteriori via zero-dimensional heat bath simulations. The comparisons of surface properties, namely heat flux and pressure, show very close agreement with previous computational fluid dynamics (CFD) results. Differences with the measurements were found to be similar to the CFD simulations. The unresolved discrepancy with the measurements could be due to inconsistent free stream conditions with the actual experimental data or missing physical phenomena altogether, for example atomic and molecular oxygen electronically excited states, three-dimensional effects, or more complex gas–surface interactions. As shown in this work, the advantages of obtaining a DSMC particle solution for these flows reside in the method's ability to be directly informed from first principles and to seamlessly describe internal energy non-equilibrium for all modes. With the advent of exascale computing and beyond, particle methods will be an increasingly important tool to verify the validity of physical assumptions in reduced-order models via fully resolved, experimental-scale simulations, down to the level of molecular-level distributions.
The immobilization of molecularly precise metal complexes to substrates, such as silica, provides an attractive platform for the design of active sites in heterogeneous catalysts. Specific steric and electronic variations of the ligand environment enable the development of structure–activity relationships and the knowledge-driven design of catalysts. At present, however, the three-dimensional environment of the precatalyst, much less the active site, is generally not known for heterogeneous single-site catalysts. We explored the degree to which NMR-based surface-to-complex interatomic distances could be used to solve the three-dimensional structures of three silica-supported metal complexes. Furthermore, the structure solution revealed unexpected features related to the environment around the metal that would be difficult to discern otherwise. This approach appears to be highly robust and, due to its simplicity, is readily applied to most single-site catalysts with little extra effort.
We present Keck/NIRC2 Ks band high-contrast coronagraphic imaging of the luminous debris disk around the nearby, young A star HD 32297 resolved at a projected separation of r = 0.3-2.5 arcse (approx 35-280 AU). The disk is highly warped to the north and exhibits a complex, "wavy" surface brightness profile interior to r approx 110 AU, where the peaks/plateaus in the profiles are shifted between the NE and SW disk lobes. The SW side of the disk is 50 - 100% brighter at r = 35 - 80 AU, and the location of its peak brightness roughly coincides with the disk's mm emission peak. Spectral energy distribution modeling suggests that HD 32297 has at least two dust populations that may originate from two separate belts likely at different locations, possibly at distances coinciding with the surface brightness peaks. A disk model for a single dust belt including a phase function with two components and a 5-10 AU pericenter offset explains the disk's warped structure and reproduces some of the surface brightness profile's shape (e.g. the overall "wavy" profile, the SB peak/plateau shifts) but more poorly reproduces the disk's brightness asymmetry. Although there may be alternate explanations, agreement between the SW disk brightness peak and disk's peak mm emission is consistent with an overdensity of very small, sub-blowout-sized dust and large, 0.1-1 mm-sized grains at approx 45 AU tracing the same parent population of planetesimals. New near-IR and submm observations may be able to clarify whether even more complex grain scattering properties or dynamical sculpting by an unseen planet are required to explain HD 32297's disk structure.
We present Keck/NIRC2 K(sub s) band high-contrast coronagraphic imaging of the luminous debris disk around the nearby, young A star HD 32297 resolved at a projected separation of r = 0.3 - 2.5" (approx equals 35 - 280 AU). The disk is highly warped to the north and exhibits a complex, "wavy" surface brightness profile interior to r approx equals 110 AU, where the peaks/plateaus in the profiles are shifted between the NE and SW disk lobes. The SW side of the disk is 50 - 100% brighter at r = 35 - 80 AU, and the location of its peak brightness roughly coincides with the disk's mm emission peak. Spectral energy distribution modeling suggests that HD 32297 has at least two dust populations that may originate from two separate belts likely at different locations, possibly at distances coinciding with the surface brightness peaks. A disk model fur a single dust belt including a phase function with two components and a 5 - 10 AU pericenter offset explains the disk's warped structure and reproduces some of the surface brightness profile's shape (e.g. the overall "wavy" profile, the SB peak/plateau shifts) but more poorly reproduces the disk's brightness asymmetry and the profile at wider separations (r > 110 AU). Although there may be a1ternate explanations, agreement between the SW disk brightness peak and disk's peak rom emission is consistent with an overdensity of very small, sub-blowout-sized dust and large, 0.1 - 1 mm-sized grains at approx equal 45 AU tracing the same parent population of planetesimals. New near-IR and submm observations may be able to clarify whether even more complex grain scattering properties or dynamical sculpting by an unseen planet are required to explain HD 32297's disk structure.
An understanding of the physical and chemical process occurring in a nuclear explosion enables predictions of the effects of nuclear weapons, including characteristics of radioactive fallout resulting from the explosion. Near-surface nuclear explosions are of particular interest due to the potential for significant amounts of environmental material to interact with and alter the physical and chemical behavior of the fireball. Such interactions have the potential to affect the distribution of radioactive species in the fireball and subsequently become incorporated into fallout through a process known as radiochemical fractionation. Studying variations in fallout formed in different historical testing environments allows us to understand the influence of local environments on fallout formation processes. In particular, constraining variations in thermal evolution and redox conditions during the evolution of the fireball can be useful to understanding how sensitive fallout radiochemical fractionation may be to the local explosion environment. However, untangling these conditions in complex, multicomponent fallout is a challenge. Here we present one method of constraining and interpreting fallout formation conditions by relating computationally derived phase stability predictions to observations in historic fallout. Development of such approaches will help improve physics-based models of fallout formation and radiochemical fractionation in complex, near surface nuclear detonations.
During the Medieval Climate Anomaly (MCA), Western North America experienced episodes of intense aridity that persisted for multiple decades or longer. These megadroughts are well documented in many proxy records, but the causal mechanisms are poorly understood. General circulation models (GCMs) simulate megadroughts, but do not reproduce the temporal clustering of events during the MCA, suggesting they are not caused by the time history of volcanic or solar forcing. Instead, GCMs generate megadroughts through (1) internal atmospheric variability, (2) sea-surface temperatures, and (3) land surface and dust aerosol feedbacks. While no hypothesis has been definitively rejected, and no GCM has accurately reproduced all features (e.g., timing, duration, and extent) of any specific megadrought, their persistence suggests a role for processes that impart memory to the climate system (land surface and ocean dynamics). Over the 21st century, GCMs project an increase in the risk of megadrought occurrence through greenhouse gas forced reductions in precipitation and increases in evaporative demand. This drying is robust across models and multiple drought indicators, but major uncertainties still need to be resolved. These include the potential moderation of vegetation evaporative losses at higher atmospheric [CO2], variations in land surface model complexity, and decadal to multidecadal modes of natural climate variability that could delay or advance onset of aridification over the the next several decades. Because future droughts will arise from both natural variability and greenhouse gas forced trends in hydroclimate, improving our understanding of the natural drivers of persistent multidecadal megadroughts should be a major research priority.