Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spectral simulations”

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,279 records · Page 71

Spatial frequency spectrum of the x-ray scatter distribution in CBCT projections

Purpose: X-ray scatter is a source of significant image quality loss in cone-beam computed tomography (CBCT). The use of Monte Carlo (MC) simulations separating primary and scattered photons has allowed the structure and nature of the scatter distribution in CBCT to become better elucidated. This work seeks to quantify the structure and determine a suitable basis function for the scatter distribution by examining its spectral components using Fourier analysis.Methods: The scatter distribution projection data were simulated using a CBCT MC model based on the EGSnrc code. CBCT projection data, with separated primary and scatter signal, were generated for a 30.6 cm diameter water cylinder [single angle projection with varying axis-to-detector distance (ADD) and bowtie filters] and two anthropomorphic phantoms (head and pelvis, 360 projections sampled every 1°, with and without a compensator). The Fourier transform of the resulting scatter distributions was computed and analyzed both qualitatively and quantitatively. A novel metric called the scatter frequency width (SFW) is introduced to determine the scatter distribution's frequency content. The frequency content results are used to determine a set basis functions, consisting of low-frequency sine and cosine functions, to fit and denoise the scatter distribution generated from MC simulations using a reduced number of photons and projections. The signal recovery is implemented using Fourier filtering (low-pass Butterworth filter) and interpolation. Estimates of the scatter distribution are used to correct and reconstruct simulated projections.Results: The spatial and angular frequencies are contained within a maximum frequency of 0.1 cm{sup −1} and 7/(2π) rad{sup −1} for the imaging scenarios examined, with these values varying depending on the object and imaging setup (e.g., ADD and compensator). These data indicate spatial and angular sampling every 5 cm and π/7 rad (∼25°) can be used to properly capture the scatter distribution, with reduced sampling possible depending on the imaging scenario. Using a low-pass Butterworth filter, tuned with the SFW values, to denoise the scatter projection data generated from MC simulations using 10{sup 6} photons resulted in an error reduction of greater than 85% for the estimating scatter in single and multiple projections. Analysis showed that the use of a compensator helped reduce the error in estimating the scatter distribution from limited photon simulations by more than 37% when compared to the case without a compensator for the head and pelvis phantoms. Reconstructions of simulated head phantom projections corrected by the filtered and interpolated scatter estimates showed improvements in overall image quality.Conclusions: The spatial frequency content of the scatter distribution in CBCT is found to be contained within the low frequency domain. The frequency content is modulated both by object and imaging parameters (ADD and compensator). The low-frequency nature of the scatter distribution allows for a limited set of sine and cosine basis functions to be used to accurately represent the scatter signal in the presence of noise and reduced data sampling decreasing MC based scatter estimation time. Compensator induced modulation of the scatter distribution reduces the frequency content and improves the fitting results.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Radiation transfer of models of massive star formation. III. The evolutionary sequence

We present radiation transfer simulations of evolutionary sequences of massive protostars forming from massive dense cores in environments of high mass surface densities, based on the Turbulent Core Model. The protostellar evolution is calculated with a multi-zone numerical model, with the accretion rate regulated by feedback from an evolving disk wind outflow cavity. The disk evolution is calculated assuming a fixed ratio of disk to protostellar mass, while the core envelope evolution assumes an inside-out collapse of the core with a fixed outer radius. In this framework, an evolutionary track is determined by three environmental initial conditions: the core mass M{sub c} , the mass surface density of the ambient clump Σ{sub cl}, and the ratio of the core's initial rotational to gravitational energy β {sub c}. Evolutionary sequences with various M{sub c} , Σ{sub cl}, and β {sub c} are constructed. We find that in a fiducial model with M{sub c} = 60 M {sub ☉}, Σ{sub cl} = 1 g cm{sup –2}, and β {sub c} = 0.02, the final mass of the protostar reaches at least ∼26 M {sub ☉}, making the final star formation efficiency ≳ 0.43. For each of the evolutionary tracks, radiation transfer simulations are performed at selected stages, with temperature profiles, spectral energy distributions (SEDs), and multiwavelength images produced. At a given stage, the envelope temperature depends strongly on Σ{sub cl}, with higher temperatures in a higher Σ{sub cl} core, but only weakly on M{sub c} . The SED and MIR images depend sensitively on the evolving outflow cavity, which gradually widens as the protostar grows. The fluxes at ≲ 100 μm increase dramatically, and the far-IR peaks move to shorter wavelengths. The influence of Σ{sub cl} and β {sub c} (which determines disk size) are discussed. We find that, despite scatter caused by different M{sub c} , Σ{sub cl}, β {sub c}, and inclinations, sources at a given evolutionary stage appear in similar regions of color-color diagrams, especially when using colors with fluxes at ≳ 70 μm, where scatter due to inclination is minimized, implying that such diagrams can be useful diagnostic tools for identifying the evolutionary stages of massive protostars. We discuss how intensity profiles along or perpendicular to the outflow axis are affected by environmental conditions and source evolution and can thus act as additional diagnostics of the massive star formation process.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulated JWST Data Sets for Multispectral and Hyperspectral Image Fusion

The James Webb Space Telescope (JWST) will provide multispectral and hyperspectral infrared images of a large number of astrophysical scenes. Multispectral images will have the highest angular resolution, while hyperspectral images (e.g., with integral field unit spectrometers) will provide the best spectral resolution. This paper aims at providing a comprehensive framework to generate an astrophysical scene and to simulate realistic hyperspectral and multispectral data acquired by two JWST instruments, namely, NIRCam Imager and NIRSpec IFU. We want to show that this simulation framework can be resorted to assess the benefits of fusing these images to recover an image of high spatial and spectral resolutions. To do so, we make a synthetic scene associated with a canonical infrared source, the Orion Bar. We develop forward models including corresponding noises for the two JWST instruments based on their physical features. JWST observations are then simulated by applying the forward models to the aforementioned synthetic scene. We test a dedicated fusion algorithm we developed on these simulated observations. We show that the fusion process reconstructs the high spatio-spectral resolution scene with a good accuracy on most areas, and we identify some limitations of the method to be tackled in future works. The synthetic scene and observations presented in the paper can be used, for instance, to evaluate instrument models, pipelines, or more sophisticated algorithms dedicated to JWST data analysis. Besides, fusion methods such as the one presented in this paper are shown to be promising tools to fully exploit the unprecedented capabilities of the JWST.

79 ASTRONOMY AND ASTROPHYSICS↗

Chemical Source Inversion using Assimilated Constituent Observations in an Idealized Two-dimensional System

We present a source inversion technique for chemical constituents that uses assimilated constituent observations rather than directly using the observations. The method is tested with a simple model problem, which is a two-dimensional Fourier-Galerkin transport model combined with a Kalman filter for data assimilation. Inversion is carried out using a Green's function method and observations are simulated from a true state with added Gaussian noise. The forecast state uses the same spectral spectral model, but differs by an unbiased Gaussian model error, and emissions models with constant errors. The numerical experiments employ both simulated in situ and satellite observation networks. Source inversion was carried out by either direct use of synthetically generated observations with added noise, or by first assimilating the observations and using the analyses to extract observations. We have conducted 20 identical twin experiments for each set of source and observation configurations, and find that in the limiting cases of a very few localized observations, or an extremely large observation network there is little advantage to carrying out assimilation first. However, in intermediate observation densities, there decreases in source inversion error standard deviation using the Kalman filter algorithm followed by Green's function inversion by 50% to 95%.

Tangborn, Andrew↗

System Identification and Uncertainty Quantification Using Orthogonal Excitations and the Semi-span Super Sonic Transport (S4T) Model

Orthogonal harmonic multisine excitations were utilized in a wind tunnel test and in simulation of the SemiSpan Supersonic Transport model to assess aeroservoelastic characteristics. Fundamental issues associated with analyzing sinusoidal signals were examined, including spectral leakage, excitation truncation, and uncertainties on frequency response functions and mean-square coherence. Simulation allowed for evaluation of these issues relative to a truth model, while wind tunnel data introduced real-world implementation issues.

Heeg, Jennifer↗

A Universal Power-law Prescription for Variability from Synthetic Images of Black Hole Accretion Flows

We present a framework for characterizing the spatiotemporal power spectrum of the variability expected from the horizon-scale emission structure around supermassive black holes, and we apply this framework to a library of general relativistic magnetohydrodynamic (GRMHD) simulations and associated general relativistic ray-traced images relevant for Event Horizon Telescope (EHT) observations of Sgr A*. We find that the variability power spectrum is generically a red-noise process in both the temporal and spatial dimensions, with the peak in power occurring on the longest timescales and largest spatial scales. When both the time-averaged source structure and the spatially integrated light-curve variability are removed, the residual power spectrum exhibits a universal broken power-law behavior. On small spatial frequencies, the residual power spectrum rises as the square of the spatial frequency and is proportional to the variance in the centroid of emission. Beyond some peak in variability power, the residual power spectrum falls as that of the time-averaged source structure, which is similar across simulations; this behavior can be naturally explained if the variability arises from a multiplicative random field that has a steeper high-frequency power-law index than that of the time-averaged source structure. We briefly explore the ability of power spectral variability studies to constrain physical parameters relevant for the GRMHD simulations, which can be scaled to provide predictions for black holes in a range of systems in the optically thin regime. We present specific expectations for the behavior of the M87* and Sgr A* accretion flows as observed by the EHT.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling Soft X‐Ray Emissions at the Dayside Magnetopause

In this study, we simulate the Solar Wind Charge Exchange (SWCX) soft X-ray emissions at dayside magnetosheath and cusps by using magnetohydrodynamic (MHD) and LAtmos TEst Particle (LaTeP) models. MHD models are unable to resolve the particle kinetic effects, such as the different behaviors of ions with different q/m, or distinguish the magnetospheric plasma from the solar wind plasma. We investigate these effects with the LaTeP model. As the LaTeP model does not self-compute magnetic and electric field, the magnetic and electric field data obtained from Open Geospace General Circulation Model (OpenGGCM) and Lagrangian version of the piecewise parabolic method (PPMLR) MHD model are used as the input to LaTeP model. The soft X-ray emissivity maps simulated from pure OpenGGCM and PPMLR MHD approaches and from LaTeP-OpenGGCM and LaTeP-PPMLR approaches are presented and compared. The results indicate that the LaTeP model can well resolve the kinetic effects and can be used to investigate the individual spectral characteristics. Therefore, the LaTeP model is a complementary approach for simulating the X-ray emissions near the dayside magnetopause. We also calculate the ratio of integrated OVII/OVIII line intensities, produced by charge exchange of O7+ ions and O8+ ions, respectively. We find a relatively higher ratio at the bow shock compared to the surrounding areas, suggesting that this ratio can be an effective parameter to identify the bow shock location.

Qiuyu Xu↗

Modeling hot gas flow in the low-luminosity active galactic nucleus of NGC 3115

Based on the dynamical black hole (BH) mass estimates, NGC 3115 hosts the closest billion solar mass BH. Deep studies of the center revealed a very underluminous active galactic nucleus (AGN) immersed in an old massive nuclear star cluster. Recent 1 Ms Chandra X-ray visionary project observations of the NGC 3115 nucleus resolved hot tenuous gas, which fuels the AGN. In this paper we connect the processes in the nuclear star cluster with the feeding of the supermassive BH. We model the hot gas flow sustained by the injection of matter and energy from the stars and supernova explosions. We incorporate electron heat conduction as the small-scale feedback mechanism, the gravitational pull of the stellar mass, cooling, and Coulomb collisions. Fitting simulated X-ray emission to the spatially and spectrally resolved observed data, we find the best-fitting solutions with χ{sup 2}/dof = 1.00 for dof = 236 both with and without conduction. The radial modeling favors a low BH mass <1.3 × 10{sup 9} M {sub ☉}. The best-fitting supernova rate and the best-fitting mass injection rate are consistent with their expected values. The stagnation point is at r {sub st} ≲ 1'', so that most of the gas, including the gas at a Bondi radius r{sub B} = 2''-4'', outflows from the region. We put an upper limit on the accretion rate at 2 × 10{sup –3} M {sub ☉} yr{sup –1}. We find a shallow density profile n∝r {sup –β} with β ≈ 1 over a large dynamic range. This density profile is determined in the feeding region 0.''5-10'' as an interplay of four processes and effects: (1) the radius-dependent mass injection, (2) the effect of the galactic gravitational potential, (3) the accretion flow onset at r ≲ 1'', and (4) the outflow at r ≳ 1''. The gas temperature is close to the virial temperature T{sub v} at any radius.

79 ASTRONOMY AND ASTROPHYSICS↗

A Characteristic-Based Spectral Element Method for Moving-Domain Problems

Here, we introduce a characteristic-based numerical procedure for simulating incompressible flows in domains with moving boundaries. Our approach utilizes an operator-integration-factor splitting technique to help produce an efficient and stable numerical scheme. Using the spectral element method and an arbitrary Lagrangian-Eulerian formulation, we explore flows where the convective acceleration effects are non-negligible. Several examples, ranging from laminar to turbulent flows, are considered. Comparisons with a standard, semi-implicit time-stepping procedure illustrate the improved performance of our scheme.

97 MATHEMATICS AND COMPUTING↗

Effects of pore size on water dynamics in mesoporous silica

Water confined in mesoporous silica plays a central role in its many uses ranging from gas sorption to nanoconfined chemical reactions. Here in this paper, the influence of pore diameter (2.5–5.4 nm) on water hydrogen bond (H-bond) dynamics in MCM41 and SBA15 mesoporous silicas is investigated using femtosecond infrared vibrational spectroscopy and molecular dynamics simulations on selenocyanate (SeCN-) anions dissolved in the pores. As shown recently, SeCN- spectral diffusion is a reliable probe of surrounding water H-bond structural motions. Additionally, the long CN stretch vibrational lifetime facilitates measurement of the full range of confined dynamics, which are much slower than in bulk water. The simulations shed light on quantitative details that are inaccessible from the spatially averaged observables. The dependence of SeCN- orientational relaxation and that of spectral diffusion on the distance from the silica interface are quantitatively described with an exponential decay and a smoothed step-function, respectively. The distance-dependence of both quantities is found to be independent of the diameter of the pores, and the spatial distribution of SeCN- is markedly non-uniform, reaching a maximum between the interface and the pore center. The results indicate that the commonly invoked two-state, or core–shell, model is a more appropriate description of spectral diffusion. Using these insights, we model the full time-dependence of the measured dynamics for all pore sizes and extract the “core” and “shell” dynamical correlation functions and SeCN- spatial probability distributions. The results are critically compared to those for water confined in reverse micelles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Convective and Turbulent Motions in Nonprecipitating Cu. Part II: LES Simulated Cloud Represented by a Starting Plume

The dynamic structure of a small trade wind cumulus (Cu) is analyzed using a novel approach. Cu developing in a shear-free environment is simulated by 10-m-resolution LES model with spectral bin microphysics. The aim is to clarify the dynamical nature of cloud updraft zone (CUZ) including entrainment and mixing in growing Cu. The validity of concept stating that a cloud at developing state can be represented by a parcel or a jet is tested. To investigate dynamical entrainment in CUZ performed by motions with scales larger than the turbulence scales, the modeled fields of air velocity were filtered by wavelet filter that separated convective motions from turbulent ones. Two types of objects in developing cloud were investigated: small volume ascending at maximal velocity (point parcel) and CUZ. It was found that the point parcel representing the upper part of cloud core is adiabatic. The motion of the air in this parcel ascending from cloud base determines cloud-top height. The top-hat (i.e., averaged) values of updraft velocity and adiabatic fraction in CUZ are substantially lower than those in the point parcel. Evaluation of the terms in the dynamical equation typically used in 1D cloud parcel models show that this equation can be applied for calculation of vertical velocities at the developing stage of small Cu, at least up to the heights of the inversion layer. Dynamically, the CUZ of developing cloud resembles the starting plume with the tail of nonstationary jet. Both the top-hat vertical velocity and buoyancy acceleration linearly increase with the height, at least up to the inversion layer. An important finding is that lateral entrainment of convective (nonturbulent) nature has a little effect on the top-hat CUZ velocity and cannot explain the vertical changes of conservative variables qt and θl. In contrast, entrained air lifting inside CUZ substantially decreases top-hat liquid water content and its adiabatic fraction. Possible reasons of these effects are discussed.

54 ENVIRONMENTAL SCIENCES↗

Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE): National Geospatial-Intelligence Agency Calibration Target Placements Field Campaign Report

The purpose of our research has been to develop a capability of simulating image collection from any sensor so as to account for sensor configuration, ground spectral radiance, and atmospheric contributions. Atmosphere is arguably the most difficult aspect of the image pipeline to accurately characterize and integrate into an image simulation. Some of the reasons for this include the highly variable spatial content of different absorbing and scattering species in the atmosphere. Additionally, the interaction of land and atmosphere can significantly impact the total radiation budget, especially in the southern Great Plains.

54 ENVIRONMENTAL SCIENCES↗

Process-oriented evaluation of quasi-stationary Rossby waves and their impact on surface air temperature extremes in dynamical downscaling over North America

Quasi-stationary Rossby waves are a crucial component of the general circulation and play a significant role in regional water and energy cycles, as well as in extreme events. However, process-oriented evaluation for Rossby waves is rarely performed for dynamical downscaling simulations. To close this gap, we evaluate three classes of dynamical downscaling approaches, with a focus on quasi-stationary Rossby waves and their impact on surface air temperature over North America during Northern Hemisphere summer. The three classes of models differ in the way large-scale forcing is provided: a limited-area model (LAM) constrained only by lateral boundary conditions, represented by RegCM4 from the North American branch of the Coordinated Regional Downscaling Experiment (NA-CORDEX), a LAM with spectral nudging to maintain consistency in large-scale dynamics with the forcing data, represented by the Weather Research and Forecasting (WRF) model simulation in NA-CORDEX, and a global variable-resolution model with smoothly varying grid spacings, represented by the Community Atmosphere Model version 5.4, with the Model for Prediction Across Scales (MPAS) as its dynamical core (CAM-MPAS). With no constraints on the atmospheric dynamics, CAM-MPAS exhibits several mean biases in the upper-level circulations over the Pacific Coast region: a weaker subtropical jet, a northward-shifted mid-latitude jet, and an overestimated southerly flow. With the lateral boundary constraint alone, RegCM4 also exhibits weaker jets and overestimated southerly winds off the West Coast. Rossby ray theory reveals that those wind biases direct incoming Rossby waves northward. The erroneously routed Rossby waves distort the relationship between the accumulation of wave activity over the US West Coast and surface temperature anomalies over the Southern Great Plains, which emerges approximately 4 d after the convergence of wave-activity flux in the ERA-Interim reanalysis. Furthermore, the response of heatwaves to the extreme wave activity flux is not reproduced by the two models, a serious drawback as a dynamical downscaling framework is expected to connect large-scale forcing to local-scale phenomena. The WRF model employing spectral nudging is largely free from the aforementioned problems. A pair of sensitivity simulations suggests that spectral nudging is the key to improving the dynamics of quasi-stationary Rossby waves and their impact on surface air temperature. Our results also demonstrate the effectiveness of Rossby wave diagnostics that allow for realistic background flows for assessing the credibility of dynamical downscaling over North America, where incoming Rossby waves propagate through complex circulation patterns before traveling across the continent.

54 ENVIRONMENTAL SCIENCES↗

FAST VARIABILITY AND MILLIMETER/IR FLARES IN GRMHD MODELS OF Sgr A* FROM STRONG-FIELD GRAVITATIONAL LENSING

We explore the variability properties of long, high-cadence general relativistic magnetohydrodynamic (GRMHD) simulations across the electromagnetic spectrum using an efficient, GPU-based radiative transfer algorithm. We focus on both standard and normal evolution (SANE) and magnetically arrested disk (MAD) simulations with parameters that successfully reproduce the time-averaged spectral properties of Sgr A* and the size of its image at 1.3 mm. We find that the SANE models produce short-timescale variability with amplitudes and power spectra that closely resemble those inferred observationally. In contrast, MAD models generate only slow variability at lower flux levels. Neither set of models shows any X-ray flares, which most likely indicates that additional physics, such as particle acceleration mechanisms, need to be incorporated into the GRMHD simulations to account for them. The SANE models show strong, short-lived millimeter/infrared (IR) flares, with short (≲1 hr) time lags between the millimeter and IR wavelengths, that arise from the combination of short-lived magnetic flux tubes and strong-field gravitational lensing near the horizon. Such events provide a natural explanation for the observed IR flares with no X-ray counterparts.

79 ASTRONOMY AND ASTROPHYSICS↗

MOJAVE: Monitoring of jets in active galactic nuclei with VLBA experiments. XI. Spectral distributions

We have obtained milliarcsecond-scale spectral index distributions for a sample of 190 extragalactic radio jets through the Monitoring of Jets in Active Galactic Nuclei with the VLBA Experiments (MOJAVE) project. The sources were observed in 2006 at 8.1, 8.4, 12.1, and 15.4 GHz, and we have determined spectral index maps between 8.1 and 15.4 GHz to study the four-frequency spectrum in individual jet features. We have performed detailed simulations to study the effects of image alignment and (u, v)-plane coverage on the spectral index maps to verify our results. We use the spectral index maps to study the spectral index evolution along the jet and determine the spectral distributions in different locations of the jets. The core spectral indices are on average flat with a mean value of +0.22 ± 0.03 for the sample, while the jet spectrum is in general steep with a mean index of –1.04 ± 0.03. A simple power-law fit is often inadequate for the core regions, as expected if the cores are partially self-absorbed. The overall jet spectrum steepens at a rate of about –0.001 to –0.004 per deprojected parsec when moving further out from the core with flat spectrum radio quasars having significantly steeper spectra (mean –1.09 ± 0.04) than the BL Lac objects (mean –0.80 ± 0.05). However, the spectrum in both types of objects flattens on average by ∼0.2 at the locations of the jet components indicating particle acceleration or density enhancements along the jet. The mean spectral index at the component locations of –0.81 ± 0.02 corresponds to a power-law index of ∼2.6 for the electron energy distribution. We find a significant trend that jet components with linear polarization parallel to the jet (magnetic field perpendicular to the jet) have flatter spectra, as expected for transverse shocks. Compared to quasars, BL Lacs have more jet components with perpendicular magnetic field alignment, which may explain their generally flatter spectra. The overall steepening of the spectra with distance can be explained with radiative losses if the jets are collimating or with the evolution of the high-energy cutoff in the electron spectrum if the jets are conical. This interpretation is supported by a significant correlation with the age of the component and the spectral index, with older components having steeper spectra.

79 ASTRONOMY AND ASTROPHYSICS↗

Understanding the Space Weathering of Mercury Through Laboratory Experiments

Introduction: Airless surfaces across the solar system are continually modified by energetic particles from solar wind and micrometeoroid bombardment [1,2]. This process is known as space weathering, and it alters the chemical, microstructural, and optical properties of surface regoliths on airless bodies, including Mercury. On the Moon and S-type asteroids, the reflectance spectral signatures of space weathering include reddening (increasing reflectance with increasing wavelength), darkening (lowering of reflectance), and the attenuation of characteristic absorption bands [2]. Such spectral changes are driven by the production of Fe-bearing nanoparticles(npFe) through both solar wind irradiation and micrometeoroid bombardment. While our understanding of space weathering for the Moon and near-Earth S-types asteroids is advanced, insight into how these processes operate on other planetary bodies is limited. In particular, Mercury experiences a uniquely intense space weathering environment than planetary counterparts at 1 AU, including a moreintense solar wind flux and higher velocity micrometeoroid impacts [4]. Additionally, Mercury has a surface composition unique in the inner solar system, including regions of the surface with very low albedo known as the low reflectance material (LRM), which is enriched in carbon, likely graphite, up to 4wt.% [5].In addition, the concentration of Fe across Mercury’s surface islow (<2 wt.%) compared to the Moonor S-type asteroids asteroids[6]. Our understanding of the effects of space weathering on C-rich and Fe-poor phases is limited. Since Fe plays a critical role inthe development of space weathering characteristicson other airless surfaces(e.g., npFe), its limited availability may significantly affect the development of space weathering features in Mercury surface materials. We can simulate space weathering processes in the laboratory to explore their effects on the microstructural, chemical, and spectral characteristics of Mercury surface materials[7]. Here we used pulsed laser irradiation to simulate the short duration, high-temperature events associated with micrometeoroid impacts. We performed coordinated analyses including reflectance spectroscopy and electron microscopy to investigate the spectral, chemical, and microstructural changes in these mercurian analog samples. Methods: For these experiments, we usedforsteritic olivine with varying FeOcontents, a mineral phase proposed to be abundant on the surface of Mercury. We mixed each sample with graphite to simulate LRM regions of the surface. Wesynthesized the olivinesamplesat 1-bar at NASA’s Johnson Space Centerand prepared pressed powder pellets for laser irradiation [8].We prepared three samples, each with a base layer of olivineto maintain structural integrity and topped witha surface layer containing the graphite-olivine mixture: 1) Sample SC-001 San Carlos olivine(Fo90.91),2)Sample F-S-002 with0.05 wt.% FeO olivine, and 3) F-T-004 with 0.53 wt.% FeO olivine. Each sample was mixed with 5 wt.% powdered graphite and had grain sizes ranging from 45to 125μm. We irradiated each sample usinga pulsed Nd-YAG laser, (l=1064 nm, ~6 ns pulse duration, energy of 48 mJ/pulse) while undervacuumat Northern Arizona University. The laser was rastered1x and then 5x over the surfaceof each sample to simulate progressive space weathering. We collected in situreflectance spectra from the samples after each laser pulse witha Nicolet IS50 Fourier-Transform Infrared spectrometer (lfrom 0.65-2.5 μm). We used an FEI Nova NanoSEM200scanning electron microscope (SEM) and a Hitachi TM4000 Plus benchtop SEM at Purdue University to image the surface morphology and topography of the samples. We extracted thin sections for analysis in the transmission electron microscope(TEM)using the FEI Helios NanoLab 660 focused ion beam (FIB) SEM at the University of Arizona.We performed analysis of the microstructural and chemical characteristics of the samples using the 200 keV JEOL 2500 scanning TEM at Johnson Space Center. Reflectance Spectroscopy Results:Reflectance spectra for each sample are shown in Fig. 1.SC-001:The spectrum of the unirradiated sample exhibits a weak 1.0 μm absorption feature, associated with Fe2+in the olivine,and low overall reflectance (Fig. 1a). Thereflectance and the depth of the absorption band increases after 1x laser raster but are at their lowest after 5x laser rasters.F-T-004:The unirradiated sample has a blue-sloped spectrum with low reflectance without identifiable absorption features (Fig. 1b). With progressive laser irradiation, the sample reflectance increases and becomes strongly red-sloped.F-S-002:The unirradiated sample exhibits a dark, blue-slopedspectrum. The brightness of the sample increases significantly from <0.2average reflectance over >0.8 reflectance in the most irradiated sampleand thespectral slope also becomesslightly reddened (Fig. 1c). Microstructural and Chemical Analysis: Two primary alteration textures were observed in the samples exposed to simulated space weathering: 1) fluffy C-rich,and 2) vesiculated melt. The fluffy C-rich texture is composed oflow-densitydeposits distributed across the surface of the sample(Fig. 2A). Analysis of a FIB section extracted from a low-density C-rich region in sample SC-001 reveals multiple globule-type deposits, discrete from stacked graphite, likely produced via melting from the laser irradiation [9].The vesiculated melt textureis smooth and uniformly distributed across isolated regions of the sample surface. The vesicles measure up to 100s of nmin diameter. Analysis of a FIB section from this texture was extracted from sample F-T-004 reveals a layer of amorphous melt material, close to 100 nm thick and uniform across the FIB section (Fig. 2B). Isolated regions of this melt layer contain small nanoparticles, <5 nm in diameter. Chemical analysis through energy dispersive X-ray spectroscopy reveals the composition of this layer is enriched in Si and depleted in Mg and O compared to the underlying sample. Implications for Space Weathering on Mercury: Previous experiments simulating space weathering of Mercury have showndarkening and reddening of spectra[7,10].However, our use of low-Fe materials and graphite to create a sample set more analogous to the mercurian surface. Our results indicate that sample composition plays a significant and important role in the space weathering of Mercury. In particular, our spectral data demonstrates a strong correlation between spectral slope, Fe content, and simulated space weathering. While the variation in FeO content between samples F-S-002 and F-T-004 is <0.6 wt.%, the spectra deviate from flat to strongly red-sloped(F-T-004). This reddeningmay be linked to the presence of very small nanoparticles observed in the melt textures extracted from sample F-T-004. For the SC-001 sample, the fluffy C-rich textures may be developed by the amalgamation of small graphite particles into these unique morphologies. Such observations indicate that space weathering on Mercury may result in both familiar and new microstructural and chemical characteristics. References: [1]Hapke B. (2001) J. Geophys. Res.-Planet.,106,10039–10073. [2]Pieters C.M. and Noble S.K. (2016) J. Geophys. Res-Planet., 121, 1865–1884. [3] Lucey P.G., and Riner, M.A. (2011) Icarus,212, 451-462.[4]CintalaM.J.(1992)J. Geophys. Res.-Planet.,97,947–973.[5]Klima R.L.et al.(2018)Geophys.Res.Letters, 45, 2945–2953. [6]Nittler L.R., et al. (2011) Science 333, 1847-1850.[7]Sasaki S. and Kurahashi E. (2004) Space weathering on Mercury, Adv.Space Res., 33, 2152-2155.[8] Vander KaadenK.E., et al. (2018) LPSCXLIX, Abstract 1230. [9] McGlaun M.L. et al. (2019) LPSCL, Abstract 2019. [10] TrangD.et al. (2018)LPSCXLIX,Abstract2083

M S Thompson↗

DEEP WIDEBAND SINGLE POINTINGS AND MOSAICS IN RADIO INTERFEROMETRY: HOW ACCURATELY DO WE RECONSTRUCT INTENSITIES AND SPECTRAL INDICES OF FAINT SOURCES?

Many deep wideband wide-field radio interferometric surveys are being designed to accurately measure intensities, spectral indices, and polarization properties of faint source populations. In this paper, we compare various wideband imaging methods to evaluate the accuracy to which intensities and spectral indices of sources close to the confusion limit can be reconstructed. We simulated a wideband single-pointing (C-array, L-Band (1–2 GHz)) and 46-pointing mosaic (D-array, C-Band (4–8 GHz)) JVLA observation using a realistic brightness distribution ranging from 1 μ Jy to 100 mJy and time-, frequency-, polarization-, and direction-dependent instrumental effects. The main results from these comparisons are (a) errors in the reconstructed intensities and spectral indices are larger for weaker sources even in the absence of simulated noise, (b) errors are systematically lower for joint reconstruction methods (such as Multi-Term Multi-Frequency-Synthesis (MT-MFS)) along with A-Projection for accurate primary beam correction, and (c) use of MT-MFS for image reconstruction eliminates Clean-bias (which is present otherwise). Auxiliary tests include solutions for deficiencies of data partitioning methods (e.g., the use of masks to remove clean bias and hybrid methods to remove sidelobes from sources left un-deconvolved), the effect of sources not at pixel centers, and the consequences of various other numerical approximations within software implementations. This paper also demonstrates the level of detail at which such simulations must be done in order to reflect reality, enable one to systematically identify specific reasons for every trend that is observed, and to estimate scientifically defensible imaging performance metrics and the associated computational complexity of the algorithms/analysis procedures.

79 ASTRONOMY AND ASTROPHYSICS↗

A theoretical and numerical study of polarimetric scattering and emission from random rough surfaces with anisotropic directional spectrum

In this paper, theoretical and numerical results of the polarimetric scattering and emission from random rough surfaces with anisotropic directional spectrum are presented for the remote sensing of ocean and soil surfaces. The polarimetric scattered field for rough dielectric surfaces is derived to the second order by the small perturbation method (SPM). It is found that the second-order scattered field is coherent in nature, and its coefficients for different polarizations present the lowest-order corrections to the Fresnel reflection coefficients of the surfaces. In addition, the cross-polarized (HV and VH) components of the coherent fields are reciprocal and not zero for surfaces with anisotropic directional spectrum when the azimuth angle of the incident direction is not aligned with the symmetry directions of surfaces. In order to verify the energy conservation condition of the theoretical results, which is important if the theory is to be applied to the passive polarimetry of rough surfaces, a Monte Carlo simulation is performed to numerically calculate the polarimetric reflectivities of one-dimensional random rough surfaces which are generated with a prescribed power-law spectrum in the spectral domain and transformed to the spatial domain by the FFT. The surfaces simulated by this approach are periodic with the period corresponding to the low-wavenumber cutoff. To calculate the scattering from periodic dielectric surfaces, the authors present a new numerical technique which applies the Floquet theorem to reduce the problem to one period and does not require the evaluation of one-dimensional periodic Green's function used in the conventional method of moment formulation. Once the scattering coefficients are obtained, the polarimetric Stokes vectors for the emission from the random surfaces are then calculated according to the Kirchhoff's law and are illustrated as functions of relative azimuth observation and row directions. The second-order SPM is also reduced to this one-dimensional case, and theoretical results for all four Stokes parameters are in excellent agreement with the numerical results obtained from the Monte Carlo simulation. Moreover, the second order coherent fields are found to be necessary for the theoretical evaluation of the third and fourth Stokes parameters. Without that, the reflectivities of random rough surfaces are significantly over-estimated, and the signs of the third and the fourth Stokes parameters are found to be incorrect, if calculated by using the first-order incoherent fields alone.

Yueh, S. H.↗