SEARCH · Engineering Papers

Results for “surface complexation modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

NASA NTRS2011

Microwave Processing of Planetary Surfaces for Volatile Extraction

In-Situ Resource Utilization will be necessary for sustained exploration of space. Volatiles are present in planetary soils, but water by far has the strongest potential for effective utilization. The presence of water at the lunar poles, Mars, and possibly on Phobos opens the possibility of producing LOX for propellant. Water is also a useful radiation shielding material and water (and oxygen) are expendables that are also required for habitation in space. Because of the strong function of water vapor pressure with temperature, heating soil effectively liberates water vapor by sublimation. Microwave energy will penetrate soil and heat from within much more efficiently than heating from the surface with radiant heat. This is especially true under vacuum conditions since the heat transfer rate is very low. The depth of microwave penetration is a strong function of the microwave frequency and to a lesser extent on soil dielectric properties. Methods for measuring the complex electric permittivity and magnetic permeability are being developed and have been measured for some lunar soil simulants at 0.5, 2.45, and 10 GHz from room temperature down to liquid nitrogen temperature. A new method for delivery of microwaves deep into a planetary surface is being prototyped with laboratory experiments and modeled with COMSOL MultiPhysics. We have plans to set up a planetary testbed in a large vacuum chamber in the coming year. Recent results will be presented.

Ethridge, Edwin C.↗

NASA NTRS2021

The Role of Complex Ionized Absorber in the Soft X-ray Spectra of Intermediate Polars

In magnetic Cataclysmic Variables (mCVs), X-ray radiation originates from the shock heated multi-temperature plasma in the post-shock region near the white dwarf surface. These X-rays are modified by a complex distribution of absorbers in the pre-shock region. The presence of photo-ionized lines and warm absorber features in the soft X-ray spectra of these mCVs suggests that these absorbers are ionized. We developed the ionized complex absorber model zxipab, which is represented by a power-law distribution of ionized absorbers in the pre-shock ow. Using the ionized absorber model zxipab along with a cooling ow model and a reflection component, we model the broadband Chandra/HETG and NuSTAR spectra of two IPs: NY Lup and V1223 Sgr. We nd that this model describes well many of the H and He like emission lines from medium Z elements, which arises from the collisionally excited plasma. However the model fails to account for some of the He like triplets from medium Z elements, which points towards its photo-ionization origin. We do not find a compelling evidence for a blackbody component to model the soft excess seen in the residuals of the Chandra/HETG spectra, which could be due to the uncertainties in estimation of the interstellar absorption of these sources using Chandra/HETG data and/or excess fluxes seen in some photo-ionized emission lines which are not accounted by the cooling ow model. We describe the implications of this model with respect to the geometry of the pre-shock region in these two IPs

Nazma Islam Syeda↗

NASA NTRS2013

Experimental Investigation of Charging Properties of Interstellar Type Silica Dust Grains by Secondary Electron Emissions

The dust charging by electron impact is an important dust charging processes in astrophysical and planetary environments. Incident low energy electrons are reflected or stick to the grains charging the dust grains negatively. At sufficiently high energies electrons penetrate the grains, leading to excitation and emission of electrons referred to as secondary electron emission (SEE). Available classical theoretical models for calculations of SEE yields are generally applicable for neutral, planar, or bulk surfaces. These models, however, are not valid for calculations of the electron impact charging properties of electrostatically charged micron/submicron-size dust grains in astrophysical environments. Rigorous quantum mechanical models are not yet available, and the SEE yields have to be determined experimentally for development of more accurate models for charging of individual dust grains. At the present time, very limited experimental data are available for charging of individual micron-size dust grains, particularly for low energy electron impact. The experimental results on individual, positively charged, micron-size lunar dust grains levitated carried out by us in a unique facility at NASA-MSFC, based on an electrodynamic balance, indicate that the SEE by electron impact is a complex process. The electron impact may lead to charging or discharging of dust grains depending upon the grain size, surface potential, electron energy, electron flux, grain composition, and configuration (Abbas et al, 2010, 2012). In this paper, we discuss SEE charging properties of individual micron-size silica microspheres that are believed to be analogs of a class of interstellar dust grains. The measurements indicate charging of the 0.2m silica particles when exposed to 25 eV electron beams and discharging when exposed to higher energy electron beams. Relatively large size silica particles (5.2-6.82m) generally discharge to lower equilibrium potentials at both electron energies. These measurements conducted on silica microspheres are qualitatively similar in nature to our previous SEE measurements on lunar Apollo missions dust samples.

Tankosic, D.↗

NASA NTRS1997

Inconsistencies in the Understanding of Solder Joint Reliability Physics

Over the years, many analytical and experimental research studies have aimed to improve the state-of-the-art assessment of solder joint integrity from a physics-of-failure perspective. Although much progress has been made, there still exist many inconsistent and even contradictory correlations and conclusions. Before discussing some of the prominent inconsistencies found in the literature, this paper reviews the fundamental physics underlying the nature of solder failure...Using the complex constitutive properties of solder, fundamental mechanical and thermomechanical proccesses can be modeled to demonstrate some of the inconsistencies in the literature.

solder joints solder solder joint failure joint fa↗

NASA NTRS2001

Assimilation of GOES Land Surface Data into a Mesoscale Models

A technique has been developed for assimilating Geostationary Operational Environmental Satellite (GOES)-derived skin temperature tendencies and insolation into the surface energy budget equation of a mesoscale model so that the simulated rate of temperature change closely agrees with the satellite observations. A critical assumption of the technique is that the availability of moisture (either from the soil or vegetation) is the least known term in the model's surface energy budget. Therefore, the simulated latent heat flux, which is a function of surface moisture availability, is adjusted based upon differences between the modeled and satellite-observed skin temperature tendencies. An advantage of this technique is that satellite temperature tendencies are assimilated in an energetically consistent manner that avoids energy imbalances and surface stability problems that arise from direct assimilation of surface shelter temperatures. The fact that the rate of change of the satellite skin temperature is used rather than the absolute temperature means that sensor calibration is not as critical. The assimilation technique has been applied to the Oklahoma-Kansas region during the spring-summer 2000 time period when dynamic changes in vegetation cover occur. In April, central Oklahoma is characterized by large NDVI associated with winter wheat while surrounding areas are primarily rangeland with lower NDVI. In July the vegetation pattern reverses as the central wheat area changes to low NDVI due to harvesting and the surrounding rangeland is greener than it was in April. The goal of this study is to determine if assimilating satellite land surface data can improve simulation of the complex spatial distribution of surface energy and water fluxes across this region. The PSU/NCAR NM5 V3 system is used in this study. The grid configuration consists of a 36-km CONUS domain and a 12-km nest over the area of interest. Bulk verification statistics (BIAS and RMSE) of surface air temperature and dewpoint indicates that assimilation of the satellite data results reduces both the bias and RMSE for both state variables. In addition, comparison of model data with ARM/CART EBBR flux observations reveals that the assimilation technique adjusts the bowen ratio in a realistic fashion.

Lapenta, William M.↗

NASA NTRS1979

Solar p-mode oscillations as a tracer of radial differential rotation

Photoelectric observations of solar p-modes obtained with improved wavenumber and frequency resolution are presented. The observations are compared with model calculations of the p-modes, and the degree of spatial and temporal coherence of the observed wave pattern is investigated. It is found that the p-mode oscillations pervade the visible surface of the sun with a high degree of coherence in space and time, so that the whole complex pattern of standing waves with its nodes and antinodes can be regarded as a fixed pattern corotating with the solar surface layers. The p-modes are introduced as a tracer of solar rotational flow velocities. The equatorial differential rotation is estimated as a function of effective depth on the basis of the theoretical contribution functions for the p-modes recently derived by Ulrich et al. (1978). The results strongly indicate that the angular speed of rotation is not uniform even in the relatively shallow layer extending about 20,000 km below the photosphere.

Deubner, F.-L.↗

NASA NTRS1991

Wing tip vortex calculations with an unstructured adaptive-grid Euler solver

A solution-adaptive grid method has been developed for computing tip-vortex flowfields around rectangular wings. This method uses subdivision in order to locally refine the grid in regions with high vorticity. Two different flow solvers are used. Each solves the three-dimensional Euler equations on unstructured grids. Computed results are compared to experimentally measured surface pressures and vortex velocities on a NACA 0015 rectangular wing. Predicted results for surface pressures and integrated lift agree well with the experimental data. The predicted size of the rotational vortex core is larger than the experimentally-measured value and the peak velocities are less. This discrepancy appears to be caused by deficiencies in the inviscid Euler-equation model. This model cannot capture the complex viscous effects at the tip that determine the detailed structure of the resulting vortex. In spite of this limitation, the present Euler unstructured adaptive-grid method demonstrates the ability to convert vortical flows with low numerical diffusion. Applications for modeling helicopter rotor wake systems are discussed.

Strawn, Roger C.↗

NASA NTRS1993

Constraints on the Martian cratering rate imposed by the SNC meteorites and Vallis Marineris layered deposits

Following two independent lines of evidence -- estimates of the age and formation time of a portion of the Martian geologic column exposed in the layered deposits and the crystallization and ejection ages of the SNC meteorites -- it appears that the Martian cratering rate must be double the lunar rate or even higher. This means models such as NHII or NHIII (Neukum and Hiller models II and III), which estimate the Martian cratering rate as being several times lunar are probably far closer to reality on Mars than lunar rates. The effect of such a shift is profound: Mars is transformed from a rather Moon-like place into a planet with vigorous dynamics, multiple large impacts, erosion, floods, and volcanism throughout its history. A strong shift upward in cratering rates on Mars apparently solves some glaring problems; however, it creates others. The period of time during which Earth-like atmospheric conditions existed, the liquid water era on Mars, persists in NHIII up to only 0.5 b.y. ago. Scenarios of extended Earth-like conditions on Mars have been discounted in the past because they would have removed many of the craters from the early bombardment era found in the south. It does appear that some process of crater removal was quite vigorous in the north during Mars' past. Evidence exists that the northern plains may have been the home of long-lived seas or perhaps even a paleo-ocean, so models exist for highly localized destruction of craters in the north. However, the question of how the ancient crater population could be preserved in the south under a long liquid-water era found in any high-cratering-rate models is a serious question that must be addressed. It does appear to be a higher-order problem because it involves low-energy dynamics acting in localized areas, i.e., erosion of craters in the south of Mars, whereas the two problems with the low-cratering-rate models involve high-energy events acting over large areas: the formation of the Vallis Marineris, the SNC ejecting impacts, and the global atmospheric pressure and temperature conditions that allow liquid water to exist as a robust entity anywhere on the Martian surface. In any case, it appears Mars is a more complex and dynamic planet than previously supposed. It has canyons dating from the middle to late period of its history that contain apparent lake sediments bedded deeper than most sediments on Earth. Recent multiple, violent impacts on Mars have apparently provided us with multiple random samples of its surface that all crystallized less than 1.5 b.y. ago. These things cannot be accommodated in our present cratering chronologies of Mars, based on 1x lunar cratering rates, without great difficulties. These difficulties suggest that a new chronology, probably based on NHII or even NHIII, should be adopted; this new chronology will provide us with a new view of Mars as a dynamic planet of rich history.

Brandenburg, J. E.↗

NASA NTRS2002

Expression and affinity purification of recombinant proteins from plants

With recent advances in plant biotechnology, transgenic plants have been targeted as an inexpensive means for the mass production of proteins for biopharmaceutical and industrial uses. However, the current plant purification techniques lack a generally applicable, economic, large-scale strategy. In this study, we demonstrate the purification of a model protein, beta-glucuronidase (GUS), by employing the protein calmodulin (CaM) as an affinity tag. In the proposed system, CaM is fused to GUS. In the presence of calcium, the calmodulin fusion protein binds specifically to a phenothiazine-modified surface of an affinity column. When calcium is removed with a complexing agent, e.g., EDTA, calmodulin undergoes a conformational change allowing the dissociation of the calmodulin-phenothiazine complex and, therefore, permitting the elution of the GUS-CaM fusion protein. The advantages of this approach are the fast, efficient, and economical isolation of the target protein under mild elution conditions, thus preserving the activity of the target protein. Two types of transformation methods were used in this study, namely, the Agrobacterium-mediated system and the viral-vector-mediated transformation system. Copyright 2002 Elsevier Science (USA).

Non-NASA Center↗

NASA NTRS1999

Advanced Electromagnetic Methods for Aerospace Vehicles

The imminent destructive threats of Lightning on helicopters and other airborne systems has always been a topic of great interest to this research grant. Previously, the lightning induced currents on the surface of the fuselage and its interior were predicted using the finite-difference time-domain (FDTD) method as well as the NEC code. The limitations of both methods, as applied to lightning, were identified and extensively discussed in the last meeting. After a thorough investigation of the capabilities of the FDTD, it was decided to incorporate into the numerical method a subcell model to accurately represent current diffusion through conducting materials of high conductivity and finite thickness. Because of the complexity of the model, its validity will be first tested for a one-dimensional FDTD problem. Although results are not available yet, the theory and formulation of the subcell model are presented and discussed here to a certain degree. Besides lightning induced currents in the interior of an aircraft, penetration of electromagnetic fields through apertures (e.g., windows and cracks) could also be devastating for the navigation equipment, electronics, and communications systems in general. The main focus of this study is understanding and quantifying field penetration through apertures. The simulation is done using the FDTD method and the predictions are compared with measurements and moment method solutions obtained from the NASA Langley Research Center. Cavity-backed slot (CBS) antennas or slot antennas in general have many applications in aircraft-satellite type of communications. These can be flushmounted on the surface of the fuselage and, therefore, they retain the aerodynamic shape of the aircraft. In the past, input impedance and radiation patterns of CBS antennas were computed using a hybrid FEM/MoM code. The analysis is now extended to coupling between two identical slot antennas mounted on the same structure. The predictions are performed using both the hybrid FEM/MoM and the FDTD NEWS code. The results are compared with each other as well as with measurements performed in the ElectroMagnetic Anechoic Chamber (EMAC) of ASU. In addition to self and mutual impedances versus frequency, the comparisons include mutual coupling S(sub 12) as a function of distance for various slot orientations.

Balanis, Constantine A.↗

NASA NTRS2000

Applications of Fractal Analytical Techniques in the Estimation of Operational Scale

The observational scale and the resolution of remotely sensed imagery are essential considerations in the interpretation process. Many atmospheric, hydrologic, and other natural and human-influenced spatial phenomena are inherently scale dependent and are governed by different physical processes at different spatial domains. This spatial and operational heterogeneity constrains the ability to compare interpretations of phenomena and processes observed in higher spatial resolution imagery to similar interpretations obtained from lower resolution imagery. This is a particularly acute problem, since longterm global change investigations will require high spatial resolution Earth Observing System (EOS), Landsat 7, or commercial satellite data to be combined with lower resolution imagery from older sensors such as Landsat TM and MSS. Fractal analysis is a useful technique for identifying the effects of scale changes on remotely sensed imagery. The fractal dimension of an image is a non-integer value between two and three which indicates the degree of complexity in the texture and shapes depicted in the image. A true fractal surface exhibits self-similarity, a property of curves or surfaces where each part is indistinguishable from the whole, or where the form of the curve or surface is invariant with respect to scale. Theoretically, if the digital numbers of a remotely sensed image resemble an ideal fractal surface, then due to the self-similarity property, the fractal dimension of the image will not vary with scale and resolution, and the slope of the fractal dimension-resolution relationship would be zero. Most geographical phenomena, however, are not self-similar at all scales, but they can be modeled by a stochastic fractal in which the scaling properties of the image exhibit patterns that can be described by statistics such as area-perimeter ratios and autocovariances. Stochastic fractal sets relax the self-similarity assumption and measure many scales and resolutions to represent the varying form of a phenomenon as the pixel size is increased in a convolution process. We have observed that for images of homogeneous land covers, the fractal dimension varies linearly with changes in resolution or pixel size over the range of past, current, and planned space-borne sensors. This relationship differs significantly in images of agricultural, urban, and forest land covers, with urban areas retaining the same level of complexity, forested areas growing smoother, and agricultural areas growing more complex as small pixels are aggregated into larger, mixed pixels. Images of scenes having a mixture of land covers have fractal dimensions that exhibit a non-linear, complex relationship to pixel size. Measuring the fractal dimension of a difference image derived from two images of the same area obtained on different dates showed that the fractal dimension increased steadily, then exhibited a sharp decrease at increasing levels of pixel aggregation. This breakpoint of the fractal dimension/resolution plot is related to the spatial domain or operational scale of the phenomenon exhibiting the predominant visible difference between the two images (in this case, mountain snow cover). The degree to which an image departs from a theoretical ideal fractal surface provides clues as to how much information is altered or lost in the processes of rescaling and rectification. The measured fractal dimension of complex, composite land covers such as urban areas also provides a useful textural index that can assist image classification of complex scenes.

Emerson, Charles W.↗

NASA NTRS2024

PRODEM: An Annual Series of Summer DEMs (2019 through 2022) of the Marginal Areas of the Greenland Ice Sheet

Surface topography across the marginal zone of the Greenland Ice Sheet is constantly evolving in response to changing weather, season, climate, and ice dynamics. However, current digital elevation models (DEMs) for the ice sheet are usually based on data from a multi-year period, thus obscuring these changes over time. Here we present four 500 m resolution summer DEMs (PRODEMs) of the Greenland Ice Sheet marginal zone for 2019 through 2022. The PRODEMs cover the marginal zone from the ice edge to 50 km inland, hence capturing all Greenland outlet glaciers. Each PRODEM is based on data fusion of CryoSat-2 radar altimetry and ICESat-2 laser altimetry using regionally varying kriging of elevation anomalies relative to ArcticDEM. The PRODEMs are validated using leave-one-out cross-validation, and PRODEM19 is further validated against an external data set, showcasing their ability to correctly represent surface elevations within the associated spatially varying prediction uncertainties. We observe a general lowering of surface elevations during the 4-year PRODEM period, but the spatial pattern of change is highly complex and with annual changes superimposed. The PRODEMs enable detailed studies of the marginal ice sheet elevation changes. With their high spatio-temporal resolution, the PRODEMs will be of value to a wide range of researchers and users studying ice sheet dynamics and monitoring how the ice sheet responds to changing environmental conditions. PRODEMs from summer 2019 through 2022 are available at https://doi.org/10.22008/FK2/52WWHG (Winstrup, 2024), and we plan to annually update the product henceforth.

Digital Elevation Models↗

NASA NTRS1992

Meso-beta scale numerical simulation studies of terrain-induced jet streak mass/momentum perturbations

Work performed during the report period is summarized. The first numerical experiment which was performed on the North Carolina Supercomputer Center's CRAY-YMP machine during the second half of FY92 involved a 36 hour simulation of the CCOPE case study. This first coarse-mesh simulation employed the GMASS model with a 178 x 108 x 32 matrix of grid points spaced approximately 24 km apart. The initial data was comprised of the global 2.5 x 2.5 degree analyses as well as all available North American rawinsonde data valid at 0000 UTC 11 July 1981. Highly-smoothed LFM-derived terrain data were utilized so as to determine the mesoscale response of the three-dimensional atmosphere to weak terrain forcing prior to including the observed highly complex terrain of the northern Rocky Mountain region. It was felt that the model should be run with a spectrum of terrain geometries, ranging from observed complex terrain to no terrain at all, to determine how crucial the terrain was in forcing the mesoscale phenomena. Both convection and stratiform (stable) precipitation were not allowed in this simulation so that their relative importance could be determined by inclusion in forth-coming simulations. A full suite of planetary boundary layer forcing was allowed in the simulation, including surface sensible and latent heat fluxes employing the Blakadar PBL formulation. The details of this simulation, which in many ways could be considered the control simulation, including the important synoptic-scale, meso-alpha scale, and meso-beta scale circulations is described. These results are compared to the observations diagnosed by Koch and his colleagues as well as hypotheses set forth in the project proposal for terrain-influences upon the jet stream and their role in the generation of mesoscale wave phenomenon. The fundamental goal of the analyses being the discrimination among background geostrophic adjustment, terrain influences, and shearing instability in the initiation and maintainance of mesoscale internal wave phenomena. Based upon these findings, FY93 plans are discussed. A review of linear theory and theoretical modeling of a geostrophic zonal wind anomaly is included.

Lin, Yuh-Lang↗

NASA NTRS1983

A theoretical study on the mechanism of electronic to vibrational energy transfer in Hg/3P/ + CO

The mechanism of electronic-to-vibrational (E-V) energy transfer in Hg(3P) + CO collisions has been studied theoretically. The configuration interaction (CI) method was employed to calculate potential energy surfaces of the collision system. A simplified theoretical model, based on the reaction coordinate concept and the calculated potential energy characteristics, was used to discuss the mechanism of the singlet-triplet transition and the energy disposal in the collision. The results obtained were that: (a) the quenching process processed via a collision complex mechanism; and that (b) the triplet-singlet transition occurs near the collinear geometry. A model classical trajectory calculation gives a product CO vibrational distribution in good agreement with the experimental result.

Kato, S.↗

NASA NTRS1988

Life prediction modeling based on cyclic damage accumulation

A high temperature, low cycle fatigue life prediction method was developed. This method, Cyclic Damage Accumulation (CDA), was developed for use in predicting the crack initiation lifetime of gas turbine engine materials, where initiation was defined as a 0.030 inch surface length crack. A principal engineering feature of the CDA method is the minimum data base required for implementation. Model constants can be evaluated through a few simple specimen tests such as monotonic loading and rapic cycle fatigue. The method was expanded to account for the effects on creep-fatigue life of complex loadings such as thermomechanical fatigue, hold periods, waveshapes, mean stresses, multiaxiality, cumulative damage, coatings, and environmental attack. A significant data base was generated on the behavior of the cast nickel-base superalloy B1900+Hf, including hundreds of specimen tests under such loading conditions. This information is being used to refine and extend the CDA life prediction model, which is now nearing completion. The model is also being verified using additional specimen tests on wrought INCO 718, and the final version of the model is expected to be adaptable to most any high-temperature alloy. The model is currently available in the form of equations and related constants. A proposed contract addition will make the model available in the near future in the form of a computer code to potential users.

Nelson, Richard S.↗

NASA NTRS1993

Development of a large scale Chimera grid system for the Space Shuttle Launch Vehicle

The application of CFD techniques to large problems has dictated the need for large team efforts. This paper offers an opportunity to examine the motivations, goals, needs, problems, as well as the methods, tools, and constraints that defined NASA's development of a 111 grid/16 million point grid system model for the Space Shuttle Launch Vehicle. The Chimera approach used for domain decomposition encouraged separation of the complex geometry into several major components each of which was modeled by an autonomous team. ICEM-CFD, a CAD based grid generation package, simplified the geometry and grid topology definition by provoding mature CAD tools and patch independent meshing. The resulting grid system has, on average, a four inch resolution along the surface.

Pearce, Daniel G.↗

NASA NTRS2012

The Performance of Current Atmospheric Radiation Codes in Phase I of CIRC

The Continual Intercomparison of Radiation Codes (CIRC) is intended as an evolving and regularly updated reference source for evaluation of radiative transfer (RT) codes used in Global Climate Models and other atmospheric applications. In our presentation we will discuss our evaluation of the performance of 13 shortwave and 11 longwave RT codes that participated in Phase I of CIRC. CIRC differs from previous intercomparisons in that it relies on an observationally validated catalogue of cases. The seven CIRC Phase I baseline cases, five cloud-free, and two with overcast liquid clouds, are built around observations by the Atmospheric Radiation Measurements (ARM) program that satisfy the goals .of Phase I, namely to examine RT model performance in realistic, yet not overly complex, atmospheric conditions. Besides the seven baseline cases, additional idealized "subcases" are also examined to facilitate interpretation of model errors. We will quantify individual model performance with respect to reference line-by-line calculations, and will also highlight RT code behavior for conditions of doubled CO2 , aspects of utilizing a spectral specification of surface albedo, and the impact of the inclusion of scattering in the thermal infrared. Our analysis suggests that RT codes should work towards improving their calculation of diffuse shortwave flux, shortwave absorption, treatment of spectral surface albedo, and shortwave CO2 forcing. Despite practical difficulties in comparing our results to previous results by the Intercomparison of Radiation Codes in Climate Models (ICRCCM) conducted about 20 years ago, it appears that the current generation of RT codes do indeed perform better than the codes of the ICRCCM era. By enhancing the range of conditions under which participating codes are tested, future CIRC phases will hopefully allow even more rigorous examination of RT code performance.

Oreopoulos, L.↗

NASA NTRS2019

Prebiotic Chemistry of Pluto

We present the case for the presence of complex organic molecules, such as amino acids and nucleobases, formed by abiotic processes on the surface and in near-subsurface regions of Pluto. Pluto's surface is tinted with a range of non-ice substances with colors ranging from light yellow to red to dark brown; the colors match those of laboratory organic residues called tholins. Tholins are broadly characterized as complex, macromolecular organic solids consisting of a network of aromatic structures connected by aliphatic bridging units (e.g., Imanaka et al.,2004; Materese et al.,2014, 2015). The synthesis of tholins in planetary atmospheres and in surface ices has been explored in numerous laboratory experiments, and both gas- and solid-phase varieties are found on Pluto. A third variety of tholins, exposed at a site of tectonic surface fracturing called Virgil Fossae, appears to have come from a reservoir in the subsurface. Eruptions of tholin-laden liquid H2O from a subsurface aqueous repository appear to have covered portions of Virgil Fossae and its surroundings with a uniquely colored deposit (D.P. Cruikshank, personal communication) that is geographically correlated with an exposure of H2O ice that includes spectroscopically detected NH3 (C.M. Dalle Ore, personal communication). The subsurface organic material could have been derived from presolar or solar nebula processes, or might have formed in situ. Photolysis and radiolysis of a mixture of ices relevant to Pluto's surface composition (N2, CH4, CO) have produced strongly colored, complex organics with a significant aromatic content having a high degree of nitrogen substitution similar to the aromatic heterocycles pyrimidine and purine (Materese et al.,2014, 2015; Cruikshank et al.,2016). Experiments with pyrimidines and purines frozen in H2O-NH3 ice resulted in the formation of numerous nucleobases, including the biologically relevant guanine, cytosine, adenine, uracil, and thymine (Materese et al.,2017). The red material associated with the H2O ice may contain nucleobases resulting from energetic processing on Pluto's surface or in the interior. Some other Kuiper Belt objects also exhibit red colors similar to those found on Pluto and may therefore carry similar inventories of complex organic materials. The widespread and ubiquitous nature of similarly complex organic materials observed in a variety of astronomical settings drives the need for additional laboratory and modeling efforts to explain the origin and evolution of organic molecules. Pluto observations reveal complex organics on a small body that remains close to its place of origin in the outermost regions of the Solar System.

Cruikshank, D. P.↗