Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “feature matching”

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 883 records · Page 49

Detailed Infrared Study of Amorphous to Crystalline Propionitrile Ices Relevant to Observed Spectra of Titan's Stratospheric Ice Clouds

We have conducted a comprehensive study of propionitrile (C2H5CN) ice from the amorphous to crystalline phase in order to provide detailed information on this specific cyanide, which may potentially contribute to the chemical composition of the Haystack ice cloud observed in Titan's stratosphere by the Cassini Composite InfraRed Spectrometer (CIRS). Infrared transmission spectra of thin films of pure propionitrile ices deposited at low temperature (30-160 K) were collected from 50 cm1 to 11,700 cm1 (200-0.85 m). The far-infrared spectral region was specifically targeted to compare with CIRS far-infrared limb spectra. The temperature and time evolution of C2H5CN ice was thoroughly investigated to better understand discrepancies reported in pre- viously published laboratory studies on the crystalline phase of C2H5CN. Specifically, we observe peculiar temperature and time-driven ice phase transitions, revealed by significant spectral variations in the ice, which stabilizes once a complete crystalline phase is achieved. From these results, the crystalline phase of propionitrile ice was identified at deposition temperatures greater than or equal to 135 K and < 140 K. Our findings corro- borate previous studies that ruled out pure propionitrile ice as the sole chemical identity of Titan's observed Haystack emission feature. In order to understand and identify the Haystack cloud, we have initiated co-de- position experiments that incorporate mixtures of Titan-relevant organics, many of which have corresponding vapors that are abundantly present in Titan's stratosphere. In this paper, we present the result of one example of a co-deposited ternary ice mixture containing 16% hydrogen cyanide (HCN), 23% C2H5CN, and 61% benzene (C6H6). Although this co-condensed ice mixture is the best fit thus far obtained to match the broad width of the Haystack, it is still not the appropriate chemical candidate. However, it reveals an intriguing result: the strong lattice mode of pure C2H5CN ice is drastically altered by the surrounding molecules as a result of mixing in a co- condensed phase. The laboratory results reported here on propionitrile ice may help to further constrain the chemical identification of Titan's stratospheric Haystack ice cloud, as well as improve on the current state of knowledge of Titan's stratospheric ice cloud chemistry.

Nna-Mvondo, Delphine↗

A provably stable numerical method for the anisotropic diffusion equation in confined magnetic fields

We present a novel numerical method for solving the anisotropic diffusion equation in magnetic fields confined to a periodic box which is accurate and provably stable. We derive energy estimates of the solution of the continuous initial boundary value problem. A discrete formulation is presented using operator splitting in time with the summation by parts finite difference approximation of spatial derivatives for the perpendicular diffusion operator. Weak penalty procedures are derived for implementing both boundary conditions and parallel diffusion operator obtained by field line tracing. We prove that the fully-discrete approximation is unconditionally stable. Discrete energy estimates are shown to match the continuous energy estimate given the correct choice of penalty parameters. A nonlinear penalty parameter is shown to provide an effective method for tuning the parallel diffusion penalty and significantly minimises rounding errors. Several numerical experiments, using manufactured solutions, the “NIMROD benchmark” problem and a single island problem, are presented to verify numerical accuracy, convergence, and asymptotic preserving properties of the method. Finally, we present a magnetic field with chaotic regions and islands and show the contours of the anisotropic diffusion equation reproduce key features in the field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Improving Subsurface Stress Characterization for Carbon Dioxide Storage Projects by Incorporating Machine Learning Techniques

The overall objective of this project is to develop a framework for reliable characterization and prediction of the state of stress in the overburden and underburden (including the basement) in CO 2 storage reservoirs using machine learning and integrated geomechanics and geophysical methods. Specifically, we propose to develop workflow encompassing of technologies and/or methods to predict stress and pressure changes due to CO 2 injection in an active tertiary recovery site and their impacts on subtle fault activation, fractures and occurrence of microseismic events and compare responses to field observations. In this project, we anticipate using dataset from the Farnsworth field Unit (FWU) which is operated by Purdure Petroleum. A novel elastic-waveform VSP inversion technique will be used to estimate high-resolution spatial and temporal changes of elastic moduli in CO 2 storage reservoirs, which will be combined with velocity-stress relationship derived from laboratory tests to obtain subsurface pressure and stress. Clustered microseismic data will be jointly inverted for improved focal mechanisms. Least-squares reverse-time migration of microseismic waveform data will be performed to directly image fracture/fault zones. Additionally, a deep neural network machine learning technique with convolutional and recurrent layers will be used for learning the spectro-temporal structures in microseismic waveforms. The results of this geotechnical data analysis will be integrated to develop a high-resolution 3D mechanical earth model extending from the overburden sealing formations to the underburden including the basement. Mechanical properties will be derived through integration of mechanical logs, tests, available results from chemo-mechanical laboratory tests, and elastic inversion of seismic data using a combination of Bayesian and stochastic methods as well as machine learning technique. Failure features (faults/fractures) will be represented and/or modeled based on seismic and core data analysis. A transient hydrodynamic-geomechanical model will be developed through coupling with the calibrated FWU reservoir simulation model. The full physics coupled model will be used to train a reduced order proxy model using machine learning algorithm for estimating stress which will then be used with appropriate constitutive relationships and forward seismological models to simulate pressure changes and induced microseismicity. An advanced optimization framework will be developed to perform a history match to minimize error between field observations and simulated. The history matched proxy model will be verified against the full-physics equivalent. The field observations that will be used in the coupled model calibration process include pressure/stress inverted from VSP, moment magnitude from microseismic analysis, real time downhole pressure measurements, production and injection data. Parameter sensitivity and uncertainty analysis will be performed to characterize the impact of model parameter uncertainty on stress estimates. The proposed project will have significant impact on future field implementation of the proposed technology. Because the project field site is an ongoing CO 2 EOR development, the value of the new technology will be demonstrated in an operational context and evaluated as a viable risk mitigation strategy. Cost/benefit will be evaluated together with the various commercial incentives for CO 2 sequestration available to oil and gas operators. The extensive available dataset and ongoing data acquisition under the SWP Phase III work plan provides flexibility for investigation of multiple approaches and reduces technical risk.

58 GEOSCIENCES↗

Multi-spatio-temporal scales PIV in a turbulent buoyant jet discharging in a linearly stratified environment

Time-resolved particle image velocimetry is implemented with a camera array and several laser sheets; this results in a multi-spatio-temporal scale system that is modular and flexible. The setup is optimized to capture global flow features, while locally resolving in space and time near the Kolmogorov scale. The apparatus relies extensively on machine vision CMOS cameras; they are inexpensive and stream to computer hard drive with virtually continuous time-resolved records (up to one hour for the current system). This leads to statistically converged data and also helps in minimizing down time between experimental runs. Details of the implementation and design of experiment are reported. The system instruments a vertical buoyant jet discharging in a linearly stratified environment. Refractive index matched aqueous solutions enable precise optical deployment. The density difference is 3% and the fluids have similar dynamic viscosity. The jet Reynolds number is 2.00 x 10 4 and is above the mixing transition. Such flows are typically challenging to instrument and few velocity field data are available to date. The velocity statistics are reported as well as first insights gained from the campaign.

42 ENGINEERING↗

Combining Spectroscopic Techniques to Determine the Optical Constants of Powdered Lactose

A method for deriving the optical constants (n/k) of organic powdered materials using pressed pellets in the mid-infrared spectral range is introduced that combines variable angle spectroscopic ellipsometry and transmission spectroscopy. The approach is applied to lactose, in which three different forms of pellets were pressed and measured: A pure lactose pellet and a mixed lactose/KBr pellet with a large analyte percentage were used for ellipsometric measurements, and a KBr transmission pellet with only a small analyte percentage were used for transmission measurements. The transmission data provide an initial set of oscillators and improve the spectral fitting of weak absorption features (k<0.01). Ellipsometric data for the pure and mixed pellets are then fit simultaneously to derive the final n/k values for lactose from 6000-400 cm-1. An alternative method just using the ellipsometric data from the mixed pellet and the transmission data is also presented and shows good agreement with the multi-sample analysis, providing a simpler method for powders that do not press easily into pure pellets. Finally, the derived optical constants are used to model the reflectance data, demonstrating a good match with the measured reflectance spectra if non-idealities are included.

Yokosuk, Michael O.↗

Analysis and interpretation of MAGSAT anomalies over north Africa

Crustal anomaly detection with MAGSAT data is frustrated by inherent resolving power of the data and by contamination from external and core fields. Quality of the data might be tested by modeling specific tectonic features which produce anomalies that fall within proposed resolution and crustal amplitude capabilities of MAGSAT fields. To test this hypothesis, north African hotspots associated with Ahaggar, Tibesti and Darfur were modeled as magnetic induction anomalies. MAGSAT data were reduced by subtracting external and core fields to isolate scalar and vertical component crustal signals. Of the three volcanic areas, only the Ahaggar region had an associated anomaly of magnitude above error limits of the data. Hotspot hypothesis was tested for Ahaggar by seeing if predicted magnetic signal matched MAGSAT anomaly. Predicted model magnetic signal arising from surface topography of the uplift and the Curie isothermal surface was calculated at MAGSAT altitudes by Fourier transform technique modified to allow for variable magnetization. Curie isotherm surface was calculated using a method for temperature distribution in a moving plate above a fixed hotspot. Magnetic signal was calculated for a fixed plate as well as a number of plate velocities and directions.

Phillips, R. J.↗

Measuring Positions of Objects using Two or More Cameras

An improved method of computing positions of objects from digitized images acquired by two or more cameras (see figure) has been developed for use in tracking debris shed by a spacecraft during and shortly after launch. The method is also readily adaptable to such applications as (1) tracking moving and possibly interacting objects in other settings in order to determine causes of accidents and (2) measuring positions of stationary objects, as in surveying. Images acquired by cameras fixed to the ground and/or cameras mounted on tracking telescopes can be used in this method. In this method, processing of image data starts with creation of detailed computer- aided design (CAD) models of the objects to be tracked. By rotating, translating, resizing, and overlaying the models with digitized camera images, parameters that characterize the position and orientation of the camera can be determined. The final position error depends on how well the centroids of the objects in the images are measured; how accurately the centroids are interpolated for synchronization of cameras; and how effectively matches are made to determine rotation, scaling, and translation parameters. The method involves use of the perspective camera model (also denoted the point camera model), which is one of several mathematical models developed over the years to represent the relationships between external coordinates of objects and the coordinates of the objects as they appear on the image plane in a camera. The method also involves extensive use of the affine camera model, in which the distance from the camera to an object (or to a small feature on an object) is assumed to be much greater than the size of the object (or feature), resulting in a truly two-dimensional image. The affine camera model does not require advance knowledge of the positions and orientations of the cameras. This is because ultimately, positions and orientations of the cameras and of all objects are computed in a coordinate system attached to one object as defined in its CAD model.

Klinko, Steve↗

Crater Identification Algorithm for the Lost in Low Lunar Orbit Scenario

Recent emphasis by NASA on returning astronauts to the Moon has placed attention on the subject of lunar surface feature tracking. Although many algorithms have been proposed for lunar surface feature tracking navigation, much less attention has been paid to the issue of navigational state initialization from lunar craters in a lost in low lunar orbit (LLO) scenario. That is, a scenario in which lunar surface feature tracking must begin, but current navigation state knowledge is either unavailable or too poor to initiate a tracking algorithm. The situation is analogous to the lost in space scenario for star trackers. A new crater identification algorithm is developed herein that allows for navigation state initialization from as few as one image of the lunar surface with no a priori state knowledge. The algorithm takes as inputs the locations and diameters of craters that have been detected in an image, and uses the information to match the craters to entries in the USGS lunar crater catalog via non-dimensional crater triangle parameters. Due to the large number of uncataloged craters that exist on the lunar surface, a probability-based check was developed to reject false identifications. The algorithm was tested on craters detected in four revolutions of Apollo 16 LLO images, and shown to perform well.

Hanak, Chad↗

An In Situ Feed Monitoring System for Molten Salt Reactors with Fast Neutron Energy Spectrum Molten Salt Reactor Applications

Molten salt reactors (MSRs) are one of the six promising advanced reactor technologies selected for further research and development by the Generation IV International Forum. More than twenty MSR designs are actively being developed around the world. Several of these designs are liquid-fueled and intended for operation within the fast neutron energy spectrum.1 National regulations will require liquid-fueled MSRs to control and account for nuclear material within licensed facilities. Additionally, states with comprehensive safeguards agreements with the International Atomic Energy Agency (IAEA) are obligated to declare nuclear material quantities within facilities. In return, the IAEA Department of Safeguards independently verifies these quantities and provides assurance that the nuclear material and facility are being used only for peaceful purposes. One key distinction of liquid-fueled MSRs compared with other types of reactors is that in portions of the facility, the nuclear material is in bulk form rather than discrete items. Traditional nuclear material accounting techniques such as physical item counting and verification of serial numbers on fresh fuel assemblies do not translate directly to all process streams within liquid-fueled reactors. Liquid-fueled MSRs are typically designed with low excess reactivity. This feature provides safety benefits but also means that most MSRs require the addition of makeup fuel salt while a reactor is operational. The nuclear material in the initial fuel salt and in any makeup fuel salt must be quantified. Additionally, distinct nuclear material diversion and reactor misuse scenarios form the basis of the detection methods and monitoring systems developed for liquid-fueled MSRs. For example, the IAEA provides assurance that fuel salt containing nuclear material is not being diverted from the system, that the feed salt matches the reported actinide concentrations and uranium enrichment, and that no additional fertile material is being introduced into the system. Measurement systems currently used for nuclear material control and accounting are not directly applicable to achieving MSR safeguards goals. This paper concerns a system being designed to account for the nuclear material added to liquid-fueled MSRs and monitor for diversion and misuse scenarios related to MSR feed systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Operator-coached machine vision for space telerobotics

A prototype system for interactive object modeling has been developed and tested. The goal of this effort has been to create a system which would demonstrate the feasibility of high interactive operator-coached machine vision in a realistic task environment, and to provide a testbed for experimentation with various modes of operator interaction. The purpose for such a system is to use human perception where machine vision is difficult, i.e., to segment the scene into objects and to designate their features, and to use machine vision to overcome limitations of human perception, i.e., for accurate measurement of object geometry. The system captures and displays video images from a number of cameras, allows the operator to designate a polyhedral object one edge at a time by moving a 3-D cursor within these images, performs a least-squares fit of the designated edges to edge data detected with a modified Sobel operator, and combines the edges thus detected to form a wire-frame object model that matches the Sobel data.

Bon, Bruce↗

Computing the Properties of Matter with Leadership Computing Resources (Closeout Report for DE-SC0018121)

In order to add more capabilities to Halide, we have designed a new framework called Tiramisu and integrated this framework into Halide. Since Tiramisu enables Halide to target heterogeneous architectures, our development efforts have been refocused on Tiramisu. Most high-performance computer systems today are complex and increasingly heterogeneous; they may have CPUs, GPUs and FPGAs. Achieving best performance requires taking full advantage of all these different architectures. To address this issue, we have designed Tiramisu, an optimization framework that enables Halide (and other DSLs) to target heterogeneous architectures. Tiramisu is an optimization framework that takes as input a high level, architecture-independent representation of code and a set of scheduling and data mapping commands that guide code transformation. The input can either be generated by a domain-specific language (DSL) compiler such as Halide or directly written by a programmer. Tiramisu then applies the user-specified code and data-layout transformations and generates an architecture-specific, low-level intermediate representation (IR) that takes advantage of modern architectural features such as multicore parallelism, non-uniform memory (NUMA) hierarchies, clusters, and accelerators like GPUs and FPGAs. We integrated Tiramisu within Halide and implemented a representative set of benchmarks to evaluate this integration. Tiramisu is now open source and is available for public use (http://tiramisu-compiler.org/). A paper about Tiramisu was published, it shows that Tiramisu extends Halide with many new capabilities and that Tiramisu can generate efficient code for multicores, GPUs, FPGAs and distributed heterogeneous systems. The performance of code generated by the Tiramisu backends matches or exceeds hand optimized reference implementations. For example, the multicore backend matches the highly optimized Intel MKL library on many kernels and shows speedups reaching 4x over the original Halide. In addition to making Tiramisu more robust, we have used Tiramisu to implement a set of representative tensor operation for constructing baryon building blocks required for multi baryon contractions in LQCD. In order to implement this code, we needed to generalize Tiramisu in two ways: first we needed to support indirect array accesses, and second, we needed to add support for complex numbers to Tiramisu. The code generated by Tiramisu is 6x faster than the reference code. Our efforts towards an MPI based multi-node version of tiramisu have matured and the resulting code scales well on multiple nodes (tests up to 512 KNL nodes have been undertaken).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Binary pseudo-random array (BPRA) for inspection and calibration for cylindrical wavefront interferometry

High-accuracy metrology is vitally important in manufacturing ultra-high-quality free-form mirrors designed to manipulate X-ray light with nanometer-scale wavelengths. However, surface topography measurements are instrument dependent, and without the knowledge of how the instrument performs under the practical usage conditions, the measured data contain some degree of uncertainty. Binary Pseudo Random Array (BPRA) “white noise” artifact are effective and useful for characterizing the Instrument Transfer Function (ITF) of surface topography metrology tools and wavefront measurement instrument. BPRA artifact contains features with all spatial frequencies in the instrument bandpass with equal weight. As a result, power spectral density of the patterns has a deterministic white-noise-like character that allows direct determination of the ITF with uniform sensitivity over the entire spatial frequency range. The application examples include electron microscopes, x-ray microscopes, interferometric microscopes, and large field-of-view Fizeau Interferometers. Furthermore, we will introduce the application of BPRA method to characterizing the ITF of Cylindrical Wavefront Interferometry (CWI), by developing the BPRA artifact which matches the radius of curvature of the cylindrical wavefront. The data acquisition and analysis procedures for different applications of the ITF calibration technique developed are also discussed.

Munechika, K↗

A New CIERRA Gridded Product Climatology for LIS/OTD

The GLM-CIERRA processing generates two products: A level-2 cluster feature dataset that merges flashes that should have been a single flash following the clustering algorithm flash definition (for GLM, this is due to the LCFA thresholds); A set of level-3 gridded products including Flash Extent Density, Convective Probability, etc. The CIERRA Level-3 products are generated on a 0.1 degree grid. This is suboptimal given the spatial variations in GLM pixel size AND the size of gridpoints in the output grid. As a stepping-stone towards developing improved GLM-CIERRA grids, I am testing pixel matching techniques on the LIS/OTD data. The goal is to maintain accuracy while vastly improving computational efficiency. The CIERRA reclustering codes are also being applied to LIS / OTD to merge flashes split by the “first fit” clustering technique. The end result will be standardized grids at the nominal resolution of each instrument (5 km for LIS, 10 km for OTD) that take into account the pixel geometry for each event that comprises a given reclustered flash.

54 ENVIRONMENTAL SCIENCES↗

Star counts and galactic structure

Star counts are discussed as an aid in studies of galactic structure. Iterative computations using a method-of-color-bins to match numerical models with observational data are described. The star counts provide bounds for the luminosity function and the color-magnitude diagram and allow spheroid normalization in a copycat Galaxy model which posits a galaxy with a Population I disk and a Population II spheroid. Good agreement has been obtained between star counts predicted by the model and observed distributions, including color and magnitude ranges and direction. Modifications being considered for the model include consideration of the halo, a galactic bulge and a thick disk. Several galactic structural features which have yet to be modeled and which may be tractable to the copycat model corrected with star counts are identified.

Bahcall, John N.↗

Application of High-Resolution Thermal Infrared Remote Sensing and GIS to Assess the Urban Heat Island Effect

Day and night airborne thermal infrared image data at 5 m spatial resolution acquired with the 15-channel (0.45 micron - 12.2 micron) Advanced Thermal and Land Applications Sensor (ATLAS) over Alabama, Huntsville on 7 September, 1994 were used to study changes in the thermal signatures of urban land cover types between day and night. Thermal channel number 13 (9.6 micron - 10.2 micron) data with the best noise-equivalent temperature change (NEAT) of 0.25 C after atmospheric corrections and temperature calibration were selected for use in this analysis. This research also examined the relation between land cover irradiance and vegetation amount, using the Normalized Difference Vegetation Index (NDVI), obtained by ratioing the difference and the sum of the red (channel number 3: 0.60-0.63 micron) and reflected infrared (channel number 6: 0.76-0.90 micron) ATLAS data. Based on the mean radiance values, standard deviations, and NDVI extracted from 351 pairs of polygons of day and night channel number 13 images for the city of Huntsville, a spatial model of warming and cooling characteristics of commercial, residential, agricultural, vegetation, and water features was developed using a GIS approach. There is a strong negative correlation between NDVI and irradiance of residential, agricultural, and vacant/transitional land cover types, indicating that the irradiance of a land cover type is greatly influenced by the amount of vegetation present. The predominance of forests, agricultural, and residential uses associated with varying degrees of tree cover showed great contrasts with commercial and services land cover types in the center of the city, and favors the development of urban heat islands. The high-resolution thermal infrared images match the complexity of the urban environment, and are capable of characterizing accurately the urban land cover types for the spatial modeling of the urban heat island effect using a GIS approach.

Lo, C. P.↗

Lorentz symmetry and IR structure of the BFSS matrix model

The BFSS matrix model relates flat space M-theory to a large N limit of matrix quantum mechanics describing N non-relativistic D0-branes. M-theory, being a theory of gravity in flat space, has a rich infrared structure that includes various soft theorems and an infinite set of conserved charges associated to asymptotic symmetries. In this work, we ask: to what extent is this infrared structure present in BFSS? We find that all the salient features concerning the infrared structure of M-theory carry over naturally to the quantum mechanics dual. Moreover, we demonstrate that the dual statement of the soft graviton theorem in the matrix model implies that D0-brane scattering amplitudes in BFSS enjoy the full 11d Lorentz symmetry of M-theory, a claim which has been long anticipated. We also offer several first-principle consistency checks for our findings, including a computation of the soft theorem which does not presuppose the BFSS duality and a non-trivial match between several known symmetries of M-theory and BFSS that appear naturally in this formalism. These calculations give non-perturbative evidence in support of the BFSS duality as a model of flat space holography.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Modeling Transverse Space Charge effects in IOTA with pyORBIT

The role and mitigation of space charge effects are important aspects of the beam physics research to be performed at the Integrable Optics Test Accelerator (IOTA) at Fermilab. The impact of nonlinear integrability (partial and complete) on space charge driven incoherent and coherent resonances will be a new feature of this accelerator and has the potential to influence the design of future low energy proton synchrotrons. In this report we model transverse space charge effects using the PIC code pyORBIT. First we benchmark the single particle tracking results against MADX with checks on symplecticity, tune footprints, and dynamic aperture in a partially integrable lattice realized with a special octupole string insert. Our space charge calculations begin with an examination of the 4D symplecticity. Short term tracking is done first with the initial bare lattice and then with a matching of the rms values with space charge. Next, we explore slow initialization of charge as a technique to establish steady state and reduce emittance growth and beam loss following injection into a mismatched lattice. We establish values of space charge simulation parameters so as to ensure numerical convergence. Finally, we compare the simulated space charge induced tune shifts and footprints against theory.

43 PARTICLE ACCELERATORS↗

Investigations of the Mars Cratered Terrain: Smooth Plains Boundary in the Eastern Equatorial Region

As part of an ongoing study to place geologic constraints on, and to interpret the origin of the ancient cratered terrain - smooth plains boundary on Mars, structural features in the eastern hemisphere were mapped, looking at wide-scale variations in aereal percentage of crater vocer and distinctive landforms. Maps and a digitized data set for graben, scarps, ridges and elongate plateaus in the eastern hemisphere were completed for 30 deg N to 30 deg. In the equatorial region, results indicate a marked change from the N oriented ridges north of the boundary to the NW orientations of those in the inter and intracrater plains south of the boundary. The WNW orientation of the boundary itself is matched by similar orientations of graben and scraps in the southern terrain. In order to determine whether outlying plateaus and isolated knobby hills found north of the boundary are formed from material of the ancient cratered terrain, the Viking 2 approach color data were used to discriminate surface materials in the eastern hemisphere. In the southern half of the Amenthes quadrange (MC-14), several correlations exist among surficial spectral reflectance and observed geologic units. Old terrain on the eastern rim of the Isidis basin is much redder than the smooth volcanic plains on either side of the rim, even though the plains unit east of the rim is located several kilometers higher than the plains on the floor of Isidis. Based on statistics derived from the color data, the plains within Isidis were found to differ from those north of the boundary zone.

Maxwell, T. A.↗