Search

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.

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.

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.

Small bodies record the chemical, physical, and dynamical processes that gave birth to and shaped the solar system. The great variety of small bodies reflects the diversity of both their genesis and their histories. The DARe mission conducts a critical test of how small body populations reflect a history of planetary migration and planetesimal scattering. This understanding is crucial for planning future NASA missions and placing current and past missions into context.

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 Deep Space One spacecraft flew by Comet 19P/Borrelly on September 22, 2001 and returned a rich array of imagery with resolutions of up to 48 m/pixel. These images provide a window into the surface structure, processes, and geological history of a comet. Additional information is contained in the original extended abstract.

The NASA-JPL Deep Space 1 Mission (DS1) encountered the short-period Jupiter-family comet 19/P Borrelly on September 22, 2001, about 8 days after perihelion. DS1's payload contained a remote-sensing package that included a 1024 square CCD and a near IR spectrometer with ~12nm resolution. Prior to its closest approach, the remote-sensing package on the spacecraft obtained 25 CCD images of the comet and 45 near-IR spectra.

The rate of detection of meteoroid impacts on the moon by the lunar seismic network shows a characteristic diurnal variation. Assuming that these meteoroids have a flux and a preimpact orbital distribution similar to that of fireballs observed by terrestrial camera networks, one can compute the expected diurnal variation for a given set of parameters that describe the seismic wave generation and transmission on the moon. An iterative process to match the theoretical variation with the observed one has led us to the following results: (1) the majority of the detected impact events occur within a closer range of the network than was believed earlier. This results in higher meteoroid flux estimates from lunar seismic data that agree with the terrestrially measured flux. (2) For meteoroid masses smaller than 1000 g, seismic amplitude is approximately proportional to the one-fifth power of the impact speed; for larger masses it is approximately proportional to the eight-fifth power, provided that the terrestrial meteor data used for analysis are not biased. (3) Seismic efficiency of meteoroids smaller than 1000 g is significantly less than that of large meteoroids. (4) Using orbits of fireballs that represent meteorites, we predict that the share of meteorites among the detected impacts is approximately 15 percent assuming that seismic efficiency of the high-density meteorites is the same as that of average meteoroids. A greatly increased seismic efficiency for these high-density objects is not likely.

Thermoelastic stress calculations show that if the moon was initially molten only in the outer few hundred kilometers, as in the magma ocean model of the moon, the highlands crust should be aseismic. In contrast, if the moon was initially totally molten, high stress (1 to more than about 3 kbar), shallow (0 to about 6 km deep), compressional moonquakes should be occurring in the highlands crust. Calculations of the minimum stress drops made for the 28 observed shallow moonquakes suggest that 3 of them probably have stress drops in the kbar range. Thus, these very limited seismic data are consistent with the model that the moon was initially totally molten.

Error probability upper and lower bounds determined for self synchronized binary PSK communications systems, presenting maximum-likelihood and Monte Carlo computer simulation

Error probability bounds for self synchronized binary PSK communication systems, simulating decision feedback, phase doubling and maximum likelihood systems

Double error probability in differentially coherent phase shift keyed communications systems, noting computerized simulation substantiation of theory

Single channel binary PSK communication systems performance influenced by degree of RF coherence between transmitter and receiver

Performance optimization for transmitted reference and signal channel phase shift keying systems