New Analysis of Mercury Laser Altimeter Crossovers to Improve Geodetic Constraints by MESSENGER
No abstract available
Engineering topics
Publications and source records attributed to Neumann, G. A..
No abstract available
Images of asteroid (101955) Bennu acquired by the OSIRIS-REx mission reveal a rocky world covered in rubble; Shape deviates from hydrostatic surface; Internal friction and/or cohesion even if no tensile strength; Understanding the deviation of the surface from idealized shape may help constrain mechanical properties of the interior; Geologic evolution of Bennu is driven by downslope migration of surface material and rubble; May be caused by YORP-induced spin-up, re-accumulation, impact-induced seismic shaking, thermal stresses, or tidal disruption by close encounters to larger bodies.
The Lunar Orbiter Laser Altimeter (LOLA) aboard the Lunar Reconnaissance Orbiter (LRO) has collected nearly seven billion measurements of surface height on the Moon with an absolute accuracy of ∼1 m and a precision of ∼10 cm. Converting time-of-flight laser altimeter measurements to topographic elevations requires accurate knowledge of the laser pointing with respect to the spacecraft body-fixed coordinate system. To that end, we have utilized altimetric crossovers from LOLA, as well as bidirectional observations of the LOLA laser and receiver boresight via an Earth-based laser tracking ground station. Based on a sample of ∼780,000 globally distributed crossovers from the circular-orbit phase of LRO's mission (∼27 months), we derive corrections to the LOLA laser boresight. These corrections improve the cross-track and along-track agreement of the crossovers by 24% and 33%, respectively, yielding RMS residuals of ∼10 m. Since early in the LRO mission, the bidirectional laser tracking experiments have confirmed a pointing anomaly when the LOLA instrument is facing toward deep space or the night side of the Moon and have allowed the reconstruction of the laser far-field pattern and receiver telescope pointing. By conducting such experiments shortly after launch and nearly eight years later, we have directly measured changes in the laser characteristics and obtained critical data to understand the laser behavior and refine the instrument pointing model. The methods and results presented here are also relevant to the design, fabrication, and operation of future planetary laser altimeters and their long-term behavior in the space environment.
We are interested in illumination conditions and the temperature distribution within the upper two meters of regolith near the lunar poles. Here, areas exist receiving almost constant illumination near areas in permanent shadow, which were identified as potential exploration sites for future missions. For our study a numerical simulation of the illumination and thermal environment for lunar near-polar regions is needed. Our study is based on high-resolution, twenty meters per pixel and 400 x 400 km large polar Digital Terrain Models (DTMs), which were derived from Lunar Orbiter Laser Altimeter (LOLA) data. Illumination conditions were simulated by synthetically illuminating the LOLA DTMs using the horizon method considering the Sun as an extended source. We model polar illumination for the central 50 x 50 km subset and use it as an input at each time-step (2 h) to evaluate the heating of the lunar surface and subsequent conduction in the sub-surface. At surface level we balance the incoming insolation with the subsurface conduction and radiation into space, whereas in the sub-surface we consider conduction with an additional constant radiogenic heat source at the bottom of our two-meter layer. Density is modeled as depth-dependent, the specific heat parameter as temperature-dependent and the thermal conductivity as depth- and temperature-dependent. We implemented a fully implicit finite-volume method in space and backward Euler scheme in time to solve the one-dimensional heat equation at each pixel in our 50 x 50 km DTM. Due to the non-linear dependencies of the parameters mentioned above, Newton's method is employed as the non-linear solver together with the Gauss-Seidel method as the iterative linear solver in each Newton iteration. The software is written in OpenCL and runs in parallel on the GPU cores, which allows for fast computation of large areas and long time scales.
The small obliquity of Mercury causes topographic depressions located near its poles to cast persistent shadows. Many [1, 9, 15] have shown these permanently shadowed regions (PSRs) may trap water ice for geologic time periods inside cold-traps. More recently, direct evidence for the presence of water ice deposits inside craters was remotely sensed in RADAR [5] and visible imagery [3]. Albedo measurements (reflectence at 1064 nm) obtained by the MErcury Space ENviroment GEochemistry and Ranging Laser Altimeter (MLA) found unusually bright and dark areas next to Mercury's north pole [7]. Using a thermal illumination model, Paige et al. [8] found the bright deposits are correlated with surface cold-traps, and the dark deposits are correlated with subsurface cold-traps. They suggested these anomalous deposits were brought to the surface by comets and were processed by the magnetospheric radiation flux, removing hydrogen and mixing C-N-O-S atoms to form a variety of molecules which will darken with time. Here we use a thermal illumination model to find the link between the cold-trap area fraction of a rough surface and its albedo. Using this link and the measurements obtained by MESSENGER we derive a surface and a subsurface ice distribution map on Mercury's north pole below the MESSENGER spatial resolution, approximately 500 m. We find a large fraction of the polar ice on Mercury resides inside micro cold-traps (of scales 10 - 100 m) distributed along the inter-crater terrain.
We processed and analyzed one-way laser ranging data from International Laser Ranging Service ground stations to NASA's Lunar Reconnaissance Orbiter (LRO), obtained from June 13, 2009 until September 30, 2014. We pair and analyze the one-way range observables from station laser fire and spacecraft laser arrival times by using nominal LRO orbit models based on the GRAIL gravity field. We apply corrections for instrument range walk, as well as for atmospheric and relativistic effects. In total we derived a tracking data volume of approximately 3000 hours featuring 64 million Full Rate and 1.5 million Normal Point observations. From a statistical analysis of the dataset we evaluate the experiment and the ground station performance. We observe a laser ranging measurement precision of 12.3 centimeters in case of the Full Rate data which surpasses the LOLA (Lunar Orbiting Laser Altimeter) timestamp precision of 15 centimeters. The averaging to Normal Point data further reduces the measurement precision to 5.6 centimeters. We characterized the LRO clock with fits throughout the mission time and estimated the rate to 6.9 times10 (sup -8), the aging to 1.6 times 10 (sup -12) per day and the change of aging to 2.3 times 10 (sup -14) per day squared over all mission phases. The fits also provide referencing of onboard time to the TDB (Barycentric Dynamical Time) time scale at a precision of 166 nanoseconds over two and 256 nanoseconds over all mission phases, representing ground to space time transfer. Furthermore we measure ground station clock differences from the fits as well as from simultaneous passes which we use for ground to ground time transfer from common view observations. We observed relative offsets ranging from 33 to 560 nanoseconds and relative rates ranging from 2 times 10 (sup -13) to 6 times 10 (sup -12) between the ground station clocks during selected mission phases. We study the results from the different methods and discuss their applicability for time transfer.
Impact bombardment during the first billion years after the formation of the Moon produced at least several tens of basins. The Gravity Recovery and Interior Laboratory (GRAIL) mission mapped the gravity field of these impact structures at significantly higher spatial resolution than previous missions, allowing for detailed subsurface and morphological analyses to be made across the entire globe. GRAIL-derived crustal thickness maps were used to define the regions of crustal thinning observed in centers of lunar impact basins, which represents a less unambiguous measure of a basin size than those based on topographic features. The formation of lunar impact basins was modeled numerically by using the iSALE-2D hydrocode, with a large range of impact and target conditions typical for the first billion years of lunar evolution. In the investigated range of impactor and target conditions, the target temperature had the dominant effect on the basin subsurface morphology. Model results were also used to update current impact scaling relationships applicable to the lunar setting (based on assumed target temperature). Our new temperature-dependent impact-scaling relationships provide estimates of impact conditions and transient crater diameters for the majority of impact basins mapped by GRAIL. As the formation of lunar impact basins is associated with the first approximately 700 Myr of the solar system evolution when the impact flux was considerably larger than the present day, our revised impact scaling relationships can aid further analyses and understanding of the extent of impact bombardment on the Moon and terrestrial planets in the early solar system.
We used one-way laser ranging data from International Laser Ranging Service (ILRS) ground stations to NASA's Lunar Reconnaissance Orbiter (LRO) for a demonstration of orbit determination. In the one-way setup, the state of LRO and the parameters of the spacecraft and all involved ground station clocks must be estimated simultaneously. This setup introduces many correlated parameters that are resolved by using a priori constraints. More over the observation data coverage and errors accumulating from the dynamical and the clock modeling limit the maximum arc length. The objective of this paper is to investigate the effect of the arc length, the dynamical and modeling accuracy and the observation data coverage on the accuracy of the results. We analyzed multiple arcs using lengths of 2 and 7 days during a one-week period in Science Mission phase 02 (SM02,November2010) and compared the trajectories, the post-fit measurement residuals and the estimated clock parameters. We further incorporated simultaneous passes from multiple stations within the observation data to investigate the expected improvement in positioning. The estimated trajectories were compared to the nominal LRO trajectory and the clock parameters (offset, rate and aging) to the results found in the literature. Arcs estimated with one-way ranging data had differences of 5-30 m compared to the nominal LRO trajectory. While the estimated LRO clock rates agreed closely with the a priori constraints, the aging parameters absorbed clock modeling errors with increasing clock arc length. Because of high correlations between the different ground station clocks and due to limited clock modeling accuracy, their differences only agreed at the order of magnitude with the literature. We found that the incorporation of simultaneous passes requires improved modeling in particular to enable the expected improvement in positioning. We found that gaps in the observation data coverage over 12h (approximately equals 6 successive LRO orbits) prevented the successful estimation of arcs with lengths shorter or longer than 2 or 7 days with our given modeling.
Floor-fractured craters are a class of lunar crater hypothesized to form in response to the emplacement of a shallow magmatic intrusion beneath the crater floor. The emplacement of a shallow magmatic body should result in a positive Bouguer anomaly relative to unaltered complex craters, a signal which is observed for the average Bouguer anomaly interior to the crater walls. We observe the Bouguer anomaly of floor-fractured craters on an individual basis using the unfiltered Bouguer gravity solution from GRAIL and also a degree 100-600 band-filtered Bouguer gravity solution. The low-magnitude of anomalies arising from shallow magmatic intrusions makes identification using unfiltered Bouguer gravity solutions inconclusive. The observed anomalies in the degree 100-600 Bouguer gravity solution are spatially heterogeneous, although there is spatial correlation between volcanic surface morphologies and positive Bouguer anomalies. We interpret these observations to mean that the spatial heterogeneity observed in the Bouguer signal is the result of variable degrees of magmatic degassing within the intrusions.
We present initial calibration and results of passive radiometry collected by the Lunar Orbiter Laser Al- timeter onboard the Lunar Reconnaissance Orbiter over the course of 12 months. After correcting for time- and temperature-dependent dark noise and detector responsivity variations, the LOLA passive radiometry measurements are brought onto the absolute radiance scale of the SELENE Spectral Profiler. The resulting photometric precision is estimated to be approximately 5%. We leverage the unique ability of LOLA to measure normal albedo to explore the 1064 nm phase function's dependence on various geologic parameters. On a global scale, we find that iron abundance and optical maturity (quantified by FeO and OMAT) are the dominant controlling parameters. Titanium abundance (TiO2 ), surface roughness on decimeter to decameter scales, and soil thermophysical properties have a smaller effect, but the latter two are correlated with OMAT, indicating that exposure age is the driving force behind their effects in a globally-averaged sense. The phase function also exhibits a dependence on surface slope at approximately 300 m baselines, possibly the result of mass wasting exposing immature material and/or less space weathering due to reduced sky visibility. Modeling the photometric function in the Hapke framework, we find that, relative to the highlands, the maria exhibit decreased backscattering, a smaller opposition effect (OE) width, and a smaller OE amplitude. Immature highlands regolith has a higher backscattering fraction and a larger OE width compared to mature highlands regolith. Within the maria, the backscattering fraction and OE width show little dependence on TiO2 and OMAT. Variations in the phase function shape at large phase angles are observed in and around the Copernican-aged Jackson crater, including its dark halo, a putative impact melt deposit. Finally, the phase function of the Reiner Gamma Formation behaves more optically immature than is typical for its composition and OMAT, suggesting the visible-to-near-infrared spectrum and phase function respond differently to the unusual regolith evolution and properties at this location.
The Lunar Orbiter Laser Altimeter (LOLA) experiment on Lunar Reconnaissance Orbiter (LRO) is a laser altimeter that also measures the strength of the return pulse from the lunar surface. These data have been used to estimate the reflectance of the lunar surface, including regions lacking direct solar illumination. A new calibration of these data is presented that features lower uncertainties overall and more consistent results in the polar regions. We use these data, along with newly available maps of the distribution of lunar maria, also derived from LRO instrument data, to investigate a newly discovered dependence of the albedo of the lunar maria on latitude (Hemingway et al., [2015]). We confirm that there is an increase in albedo with latitude in the lunar maria, and confirm that this variation is not an artifact arising from the distribution of compositions within the lunar maria, using data from the Lunar Prospector Neutron Spectrometer. Radiative transfer modeling of the albedo dependence within the lunar maria is consistent with the very weak to absent dependence of albedo on latitude in the lunar highlands; the lower abundance of the iron source for space weathering products in the lunar highlands weakens the latitude dependence to the extent that it is only weakly detectable in current data. In addition, photometric mod- els and normalization may take into account the fact that the lunar albedo is latitude dependent, but this dependence can cause errors in normalized reflectance of at most 2% for the majority of near-nadir geometries. We also investigate whether the latitude dependent albedo may have obscured detection of small mare deposits at high latitudes. We find that small regions at high latitudes with low roughness similar to the lunar maria are not mare deposits that may have been misclassified owing to high albedos imposed by the latitude dependence. Finally, we suggest that the only modest correlations among space weathering indicators defined for the lunar samples may be due to mixing of soils from distinct latitudes.
We present initial calibration and results of passive radiometry collected by the Lunar Orbiter Laser Altimeter onboard the Lunar Reconnaissance Orbiter over the course of 12 months. After correcting for time- and temperature-dependent dark noise and detector responsivity variations, the LOLA passive radiometry measurements are brought onto the absolute radiance scale of the SELENE Spectral Profiler. The resulting photometric precision is estimated to be 5%. We leverage the unique ability of LOLA to measure normal albedo to explore the 1064 nm phase function's dependence on various geologic parameters. On a global scale, we find that iron abundance and optical maturity (quantified by FeO and OMAT) are the dominant controlling parameters. Titanium abundance (TiO2), surface roughness on decimeter to decameter scales, and soil thermo- physical properties have a smaller effect, but the latter two are correlated with OMAT, indicating that exposure age is the driving force behind their effects in a globally-averaged sense. The phase function also exhibits a dependence on surface slope at approximately 300 m baselines, possibly the result of mass wasting exposing immature material and/or less space weathering due to reduced sky visibility. Modeling the photometric function in the Hapke framework, we find that, relative to the highlands, the maria exhibit decreased backscattering, a smaller opposition effect (OE) width, and a smaller OE amplitude. Immature highlands regolith has a higher backscattering fraction and a larger OE width compared to mature highlands regolith. Within the maria, the backscattering fraction and OE width show little dependence on TiO2 and OMAT. Variations in the phase function shape at large phase angles are observed in and around the Copernican-aged Jackson crater, including its dark halo, a putative impact melt deposit. Finally, the phase function of the Reiner Gamma Formation behaves more optically immature than is typical for its composition and OMAT, suggesting the visible-to-near-infrared spectrum and phase function respond differently to the unusual regolith evolution and properties at this location. The phase difference map revealed additional geologically-influenced variations in the phase function's shape. In particular, variations were observed associated with the dark halo around Jackson crater, the impact ejecta of Copernicus and Giordano Bruno, and the Reiner Gamma Formation. For the latter, we 915 found that the phase function behaves more optically immature than the global phase function for its composition and OMAT, suggesting a difference in how the visible-to-near-IR spectrum and phase function respond to the unusual regolith evolution and properties at this location.
We present an improved lunar digital elevation model (DEM) covering latitudes within +/-60 deg, at a horizontal resolution of 512 pixels per degree ( approx.60 m at the equator) and a typical vertical accuracy approx.3 to 4 m. This DEM is constructed from approx.4.5 ×10(exp 9) geodetically-accurate topographic heights from the Lunar Orbiter Laser Altimeter (LOLA) onboard the Lunar Reconnaissance Orbiter, to which we co-registered 43,200 stereo-derived DEMs (each 1 deg×1 deg) from the SELENE Terrain Camera (TC) ( approx.10(exp 10) pixels total). After co-registration, approximately 90% of the TC DEMs show root-mean-square vertical residuals with the LOLA data of < 5 m compared to approx.50% prior to co-registration. We use the co-registered TC data to estimate and correct orbital and pointing geolocation errors from the LOLA altimetric profiles (typically amounting to < 10 m horizontally and < 1 m vertically). By combining both co-registered datasets, we obtain a near-global DEM with high geodetic accuracy, and without the need for surface interpolation. We evaluate the resulting LOLA + TC merged DEM (designated as "SLDEM2015") with particular attention to quantifying seams and crossover errors.
The Lunar Orbiter Laser Altimeter (LOLA) measures the backscattered energy of the returning altimetric laser pulse at its wavelength of 1064 nm, and these data are used to map the reflectivity of the Moon at zero-phase angle with a photometrically uniform data set. Global maps have been produced at 4 pixels per degree (about 8 kilometers at the equator) and 2 kilometers resolution within 20 deg latitude of each pole. The zero-phase geometry is insensitive to lunar topography, so these data enable characterization of subtle variations in lunar albedo, even at high latitudes where such measurements are not possible with the Sun as the illumination source. The geometric albedo of the Moon at 1064 nm was estimated from these data with absolute calibration derived from the Kaguya Multiband Imager and extrapolated to visual wavelengths. The LOLA estimates are within 2 sigma of historical measurements of geometric albedo. No consistent latitude-dependent variations in reflectance are observed, suggesting that solar wind does not dominate space weathering processes that modify lunar reflectance. The average normal albedo of the Moon is found to be much higher than that of Mercury consistent with prior measurements, but the normal albedo of the lunar maria is similar to that of Mercury suggesting a similar abundance of space weathering products. Regions within permanent shadow in the polar regions are found to be more reflective than polar surfaces that are sometimes illuminated. Limiting analysis to data with slopes less than 10 deg eliminates variations in reflectance due to mass wasting and shows a similar increased reflectivity within permanent polar shadow. Steep slopes within permanent shadow are also more reflective than similar slopes that experience at least some illumination. Water frost and a reduction in effectiveness of space weathering are offered as possible explanations for the increased reflectivity of permanent shadow; porosity is largely ruled out as the sole explanation. The south polar crater Shackleton is found to be among the most reflective craters in its size range globally but is not the most reflective, so mass wasting cannot be ruled out as a cause for the crater's anomalous reflectance. Models of the abundance of ice needed to account for the reflectance anomaly range from 3 to 14% by weight or area depending on assumptions regarding the effects of porosity on reflectance and whether ice is present as patches or is well mixed in the regolith. If differences in nanophase iron abundances are responsible for the anomaly, the permanently shadowed regions have between 50 and 80% the abundance of nanophase iron in mature lunar soil.
We have derived algorithms and techniques to precisely co-register laser altimeter profiles with gridded Digital Terrain Models (DTMs), typically derived from stereo images. The algorithm consists of an initial grid search followed by a least-squares matching and yields the translation parameters at sub-pixel level needed to align the DTM and the laser profiles in 3D space. This software tool was primarily developed and tested for co-registration of laser profiles from the Lunar Orbiter Laser Altimeter (LOLA) with DTMs derived from the Lunar Reconnaissance Orbiter (LRO) Narrow Angle Camera (NAC) stereo images. Data sets can be co-registered with positional accuracy between 0.13 m and several meters depending on the pixel resolution and amount of laser shots, where rough surfaces typically result in more accurate co-registrations. Residual heights of the data sets are as small as 0.18 m. The software can be used to identify instrument misalignment, orbit errors, pointing jitter, or problems associated with reference frames being used. Also, assessments of DTM effective resolutions can be obtained. From the correct position between the two data sets, comparisons of surface morphology and roughness can be made at laser footprint- or DTM pixel-level. The precise co-registration allows us to carry out joint analysis of the data sets and ultimately to derive merged high-quality data products. Examples of matching other planetary data sets, like LOLA with LRO Wide Angle Camera (WAC) DTMs or Mars Orbiter Laser Altimeter (MOLA) with stereo models from the High Resolution Stereo Camera (HRSC) as well as Mercury Laser Altimeter (MLA) with Mercury Dual Imaging System (MDIS) are shown to demonstrate the broad science applications of the software tool.
The Gravity Recovery and Interior Laboratory (GRAIL) spacecraft conducted the mapping of the gravity field of the Moon from March 1, 2012 to May 29, 2012. The twin spacecraft acquired highly precise K Band range-rate (KBRR) intersatellite ranging data and Deep Space Network (DSN) data during this prime mission phase from altitudes of 15 to 75 km above the lunar surface over three lunar months. We have processed these data using the NASA GSFC GEODYN orbit determination and geodetic parameter estimation program, and we have determined gravity fields up to degree and order 420 in spherical harmonics. The new gravity solutions show improved correlations with LOLA-derived topography to high degree and order and resolve many lunar features in the geopotential with a resolution of less than 30 km, including for example the central peak of the crater Tycho. We discuss the methodology used for the processing of the GRAIL data, the quality of the orbit determination on the GRAIL satellites and the derivation of the solutions, and their evaluation with independent data, including Lunar Prospector. We show that with these new GRAIL gravity solutions, we can now fit the low altitude, extended mission Lunar Prospector tracking data better than with any previous gravity model that included the LP data.
Shackleton crater, whose interior lies largely in permanent shadow, is of interest due to its potential to sequester volatiles. Observations from the Lunar Orbiter Laser Altimeter onboard the Lunar Reconnaissance Orbiter have enabled an unprecedented topographic characterization, revealing Shackleton to be an ancient, unusually well-preserved simple crater whose interior walls are fresher than its floor and rim. Shackleton floor deposits are nearly the same age as the rim, suggesting little floor deposition since crater formation over 3 billion years ago. At 1064 nm the floor of Shackleton is brighter than the surrounding terrain and the interiors of nearby craters, but not as bright as the interior walls. The combined observations are explainable primarily by downslope movement of regolith on the walls exposing fresher underlying material. The relatively brighter crater floor is most simply explained by decreased space weathering due to shadowing, but a 1-mm-thick layer containing approx 20% surficial ice is an alternative possibility.
Floor-fractured craters (FFCs) are a class of lunar craters defined by their distinctly shallow, often plate-like floors, and combinations of radial, con-centric, and polygonal floor-fractures; a variety of other interior features are often observed, such as moats, ridges, small dark-haloed pits, and patches of mare material. They were first classified by Schultz [1] , who recognized eight overall types of floor-fractured crater. These eight subtypes have widely differing appearances, a factor that could provide insight into formation mechanisms (different manifestations of the same mechanism, or indicators of varying formation mechanisms). Two formation mechanisms for FFCs were initially proposed: 1) magmatic intrusion [1], in which magma rising toward the surface in dikes encountered low-density breccia lenses beneath crater floors and spread laterally to form sills, raising and fracturing the crater floor. 2) viscous relaxation [2], in which the properties of the crust permitted viscous flow in the vicinity of the crater, causing long-wavelength relaxation of the topography and uplift and fracturing of the crater floor. Critical to distinguishing between these two end-member hypotheses and identifying others is a quantitative assessment of the topography of FFCs and knowledge of their regional and local settings. The purpose of this study is to use newly available Lunar Reconnaissance Orbiter (LRO) Lunar Orbiter Laser Altimeter (LOLA) altimeter and Lunar Reconnaissance Orbiter Camera (LROC) image data to provide an updated global catalog of the locations, classes, morphometric and morphologic characteristics of all lunar floor-fractured craters. We use the excellent 8-class system initially described in Schultz [1] as a starting point for classification and the enhanced LOLA/LROC data sets to examine and categorize all FFCs; we found evidence for a new FFC class, discernably different from the previously existing types. Our approach, and the global categorization of all FFCs, permits the spatial distribution of each FFC-subtype to be plotted and assessed allowing for further investigation into FFC formation mechanisms. Upon completion, the data set will be made available on our web site at http://www.planetary.brown.edu/html_pages/data.htm.