Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “horizon 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 19 records

Horizon Based Orientation Estimation for Planetary Surface Navigation

Planetary rovers navigate in extreme environments for which a Global Positioning System (GPS) is unavailable, maps are restricted to relatively low resolution provided by orbital imagery, and compass information is often lacking due to weak or not existent magnetic fields. However, an accurate rover localization is particularly important to achieve the mission success by reaching the science targets, avoiding negative obstacles visible only in orbital maps, and maintaining good communication connections with ground. This paper describes a horizon solution for precise rover orientation estimation. The detected horizon in imagery provided by the on board navigation cameras is matched with the horizon rendered over the existing terrain model. The set of rotation parameters (roll, pitch yaw) that minimize the cost function between the two horizon curves corresponds to the rover estimated pose.

localization↗

What Can Geolocated Sferics Tell Us About Terrestrial Gamma-Ray Flashes?

The Fermi Gamma-ray Burst Monitor (GBM) has been detecting TGFs with increasing sensitivity over the past two years, owing to changes in flight software that have lowered its threshold for triggering and, recently, allowed a search for TGFs weaker than those which would cause an onboard trigger. Associations between TGFs detected in the first 18 months of operation and sferics detected using the World Wide Lightning Location Network (WWLLN) show that TGF peaks and lightning discharges are simultaneous to within tens of microseconds, and that GBM triggered on TGFs that occurred up to a distance of 300 km from the sub-spacecraft position. In the work presented here, we look for associations between TGFs detected by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) and WWLLN sferics over the same 18 months, and we compare the match rate and detection horizon of the two instruments. We also look for associations between WWLLN sferics and more recent GBM TGFs, both triggered events and weaker TGFs uncovered in our untriggered search. We discuss whether in this new mode, GBM is detecting TGFs that are more distant from the sub-spacecraft point than 300 km, or whether the weaker TGFs are instead indicative of a luminosity distribution, either because the weaker ones originate deeper in the atmosphere or because they are intrinsically dimmer.

Connaughton, V.↗

Flight Demonstration of the Tailored Arrival Manager

A flight demonstration of arrival traffic management automation was conducted in partnership between NASA, FAA and Boeing as an element of the latter’s ecoDemonstrator2020 flight program. For the demonstration, NASA’s prototype Tailored Arrival Manager(TAM) was used to compute efficient trajectory-based solutions to traffic management problems representing those encountered during time-based metering operations today. TAM solutions involving route modifications were uplinked to a Boeing 787 airplane using Controller Pilot Data Link Communications and seamlessly loaded into the airplane’s Flight Management System (FMS). The paper describes the concept and technology behind TAM along with the data exchanges and procedures involved with the demonstration. All TAM solutions were delivered in a timely manner and successfully integrated with the airplane’s FMS, thereby demonstrating the basic feasibility of trajectory-based arrival management using currently available data communications and avionics equipage. Although TAM solutions were not executed during this initial demonstration, a limited study of trajectory prediction accuracy was possible given that the airplane flew uninterrupted, automated descents with known route and speed intent. Analysis revealed that TAM meter-fix arrival time predictions were accurate to within ± 30 seconds for time horizons of 30 minutes or less, which matched closely with FMS predictions acquired through real-time data exchange. Top-of-descent predictions over similar time horizons were found accurate to within ± 5 nautical miles.

air traffic management↗

Flight Demonstration of the Tailored Arrival Manager

A flight demonstration of arrival traffic management automation was conducted in partnership between NASA, FAA and Boeing as an element of the latter’s ecoDemonstrator2020 flight program. For the demonstration, NASA’s prototype Tailored Arrival Manager(TAM) was used to compute efficient trajectory-based solutions to traffic management problems representing those encountered during time-based metering operations today. TAM solutions involving route modifications were uplinked to a Boeing 787 airplane using Controller Pilot Data Link Communications and seamlessly loaded into the airplane’s Flight Management System (FMS). The paper describes the concept and technology behind TAM along with the data exchanges and procedures involved with the demonstration. All TAM solutions were delivered in a timely manner and successfully integrated with the airplane’s FMS, thereby demonstrating the basic feasibility of trajectory-based arrival management using currently available data communications and avionics equipage. Although TAM solutions were not executed during this initial demonstration, a limited study of trajectory prediction accuracy was possible given that the airplane flew uninterrupted, automated descents with known route and speed intent. Analysis revealed that TAM meter-fix arrival time predictions were accurate to within ± 30 seconds for time horizons of 30 minutes or less, which matched closely with FMS predictions acquired through real-time data exchange. Top-of-descent predictions over similar time horizons were found accurate to within ± 5 nautical miles.

Air traffic managment↗

Flexible Wing Base Micro Aerial Vehicles: Towards Flight Autonomy: Vision-Based Horizon Detection for Micro Air Vehicles

Recently substantial progress has been made towards design building and testifying remotely piloted Micro Air Vehicles (MAVs). This progress in overcoming the aerodynamic obstacles to flight at very small scales has, unfortunately, not been matched by similar progress in autonomous MAV flight. Thus, we propose a robust, vision-based horizon detection algorithm as the first step towards autonomous MAVs. In this paper, we first motivate the use of computer vision for the horizon detection task by examining the flight of birds (biological MAVs) and considering other practical factors. We then describe our vision-based horizon detection algorithm, which has been demonstrated at 30 Hz with over 99.9% correct horizon identification, over terrain that includes roads, buildings large and small, meadows, wooded areas, and a lake. We conclude with some sample horizon detection results and preview a companion paper, where the work discussed here forms the core of a complete autonomous flight stability system.

Nechyba, Michael C.↗

The psychophysical function for perceived gravitational-inertial force does not depend on the orientation of the otolith organs

It has generally been believed that the perceived intensity of a gravitational-inertial force depends on both the magnitude and orientation of the force with respect to the otolith organs, as does the elevator illusion. In this study, we examine the perceived intensity of Gz force and the elevator illusion as a function of the applied force and the orientation of the subject's head. Methods: Each of 7 male subjects was seated upright in a swinging chair mounted in the Ames 20-G Human Centrifuge while he set a visual target to his apparent horizon and judged the perceived intensity of Gz forces by cross-modal matches on a hand dynamometer. Plateau Gz levels were 1.00 1.25, 1.50, 2.00, 2.25, and 2.50; a 30 second ramp to plateau was used in all cases, and the duration of exposure at each plateau was 120 seconds. All measures were obtained both with the subject's head erect and pitched forward 30 degrees. Results: Although the elevator illusion changed with head orientation (F(6,60) = 7.56; p less than 0.001) the perceived intensity of Gz was essentially the same for both orientations of the head (F (6,60) = 0.61; p greater than 50). Conclusions: The results of this experiment suggest that the perceived intensity of gravitational-inertial force does not depend on otolith mechanisms in the same way as does the elevator illusion and that somesthetic, tactile, and other proprioceptive inputs are important for the psychophysical function.

Cohen, M.↗

Remote Sensing of the Ocean Surface Refractive Index via Short-Wave Infrared Polarimetry

The linear polarization of sunlight reflected by ocean surfaces in the shortwave infrared (SWIR), at geometries where specular reflection dominates the signal, is a direct function of the refractive index of the surface microlayer(SML). This simple physical concept is at the base of a novel technique presented in this study. We invert observations obtained by the airborne Research Scanning Polarimeter (RSP) in the sun glint region, where each pixel's radiance is dominated by the signal originating from the wave slopes oriented precisely to cause specular reflection. The SWIR wavelength ensures minimization of aerosol interference when radiance travels through the atmosphere; strong absorption by the water body then limits the penetration depth to the first micrometer or so, effectively probing the SML. The resulting Degree of Linear Polarization (DoLP) is then governed by the refractive index via the Fresnel law for the specfic pixel geometry, independently of the wind speed. The selected dataset concerns several field deployments from both low- and high-altitude aircraft, including total reflectance measurements with the sole purpose of accounting for the residual aerosol effect.Stable retrievals from transects above pure seawater yielded values of refractive index that match the values published in the literature within an accuracy of 5 x 10 [superscript 4]. Flying over the oil spill caused by the explosion of the Deepwater Horizon platform, detected variations were found compatible with the presence of an oil slick. The robustness of the results, guaranteedby the high RSP polarimetric accuracy (less than or equal to 0:2%), opens the possibility for remote-sensing detection of other entities that similarly affect the refractive index including whitecaps, microplastics, biological gels, seaweed and grass mats.

Surface Microlayer↗

An Updated Secondary Lunar Meteoroid Ejecta Model for Engineering Design

Introduction: The surface of the Moon is constantly being bombarded by a flux of meteoroids of various sizes. Impacts due to these meteoroids produce secondary ejecta material at much lower speeds but with a total mass larger than the original impactor. Details about the secondary ejecta are important for planning missions on the lunar surface. In this work, an updated ejecta model is presented called the Meteoroid Model of Secondary Ejecta (MeMoSeE), to re-place the Apollo-era ejecta model, NASA SP-8013 [1], in the SLS-SPEC 159 Design Specification for Natural Environments (DSNE) [2]. The model pro-duces secondary ejecta flux environments for a user-specified location on the lunar surface, and sorts the incoming secondary flux by angular direction and speed. Methods: MeMoSeE is separated into three parts: the inputs, the conversion step, and the integration of fluxes. Inputs to the model include the primary meteoroid fluxes and the primary near-Earth object (NEO) fluxes. Meteoroid fluxes, both asteroidal and cometary, are calculated using the Meteoroid Engineering Model (MEM3) [3] for different locations on the Moon. For each surface location, an ephemeris is generated using the JPL HORIZONS System [4] that feeds into MEM3. The NEO fluxes are approximated by the high-density population of MEM3 (i.e., only the directionality), where the speed distribution of the NEO fluxes is renormalized to match observations [5]. The regolith properties are used as defined in the DSNE [2]. The conversion step utilizes scaling laws given by Housen & Holsapple 2011 [6] to convert the primary impactor flux to the total mass of secondary ejecta. The ejecta distribution, at the point-of-impact (POI), is separated into a zenith angle and azimuthal angle distribution. The zenith angle distribution follows a beta distribution where the peak depends on the impact altitude angle and the impact azimuth [7]. We employ an ejecta azimuth distribution that is based on Rival & Mandeville 1999 [8] which focuses ejecta in the downstream direction for more oblique impacts. Finally, during the integration step, we sum secondary ejecta number fluxes at a particular region-of-interest (ROI) that originated from many POI locations over the entire surface of the Moon. We keep track of both altitude and azimuth angle bins as well as a range of speed bins, following the igloo gridding as done in MEM3 [9]. The ejecta particle size distribution and density is assumed to be the same as the lunar regolith [10, 2]. Results: The primary fluxes are computed for one Metonic cycle (19 years) for various locations over the lunar surface with a fixed orientation. Both the angular and speed distributions of the primary fluxes are dependent on the latitude and longitude. The over-all primary fluxes show a roughly 13% increase from the eastern limb to the western limb. In general, the speed distribution of the secondary fluxes span from a user-defined minimum speed to the escape speed of the Moon (2.38 km/s), roughly following a power-law relation [6]. Different parts of the speed distribution come from different primary impact locations on the Moon. The secondary ejecta is dominated by the slowest speeds, where these particles originate nearby the ROI. For speeds around 71% of the escape speed, the secondary ejecta originates from locations near the antipodal point. On the other hand, for secondary ejecta speeds that exceed roughly 90% the escape speed, the ejecta particles come from all over the lunar surface to the ROI. Comparing the secondary ejecta fluxes from MeMoSeE with NASA SP-8013 [1], there is a reduction by about 2-3 orders of magnitude for secondary ejecta particles greater than 1 µg. These estimates agree with recent findings from Bjorkman & Christiansen 2019 [11]. The secondary ejecta fluxes are also compared with the primary fluxes, where the ejecta fluxes are roughly an order of magnitude greater than the primary fluxes. References: [1] Cour-Palais, B. G., (1969) NASA SP-8013. [2] NASA SLS-SPEC-159 Rev. H (2020). [3] Moorhead, A. V., et al. (2019) JS&R, 1-17. [4] Giorgini, J. D., (2015) IAUGA, 29, 2256293. [5] Moorhead, A. V., (2020) Memo OSMA/MEO/Lunar-001. [6] Housen, K. R., and Holsapple, K. A., (2011) Icarus, 211(1), 856-875. [7] Gault, D. E., and Wedekind, J. A., (1978) L&PSCP, 9, 3843-3875. [8] Rival, M., and Mandeville, J., (1999) Space Debris, 1(1), 45-57. [9] Moorhead, A. V., (2019) MEM3 User Guide. [10] Carrier, W. D., (1973) The Moon, 6(3-4), 250-263. [11] Bjorkman, M. D. & Christiansen, E. L., (2019) ODC

impacts↗

Ultraviolet spectroscopy of the zodiacal light at 20-deg elongation

The zodiacal light at 20-deg elongation and 10-deg inclination was observed by rocket ultraviolet spectrometers at 10-15-A resolution in the spectral range 1200-3200 A during an experiment designed to observe comet Kohoutek (1973 XII). The data were obtained above 180 km when scattered horizon light in the startracker caused a loss of tracking on the comet. Airglow emission due to NO and O(+), identified spectroscopically and by its variation with altitude, is significant between 1900 and 2500 A. Longward of 2600 A, the spectrum matches that of the sun, and the derived value of the color ratio, relative to the visible, is 0.90 + or - 0.20. At 1600 A, an upper limit on the zodiacal-light emission of 0.07 R per A or 7 hundred-millionths erg/s per sq cm/sterad per A is obtained.

Feldman, P. D.↗

Electromagnetic signatures from supermassive binary black holes approaching merger

We present fully relativistic predictions for the electromagnetic emission produced by accretion disks surrounding spinning and nonspinning supermassive binary black holes on the verge of merging. We use the code Bothros to post-process data from 3D General Relativistic Magnetohydrodynamic (GRMHD) simulations via ray-tracing calculations. These simulations model the dynamics of a circumbinary disk and the mini-disks that form around two equal-mass black holes orbiting each other at an initial separation of 20 gravitational radii, and evolve the system for more than 10 orbits in the inspiral regime. We model the emission as the sum of thermal blackbody radiation emitted by an optically thick accretion disk and a power-law spectrum extending to hard X-rays emitted by a hot optically thin corona. We generate time-dependent spectra, images, and light curves at various frequencies to investigate intrinsic periodic signals in the emission, as well as the effects of the black hole spin. We find that prograde black hole spin makes mini-disks brighter since the smaller ISCO angular momentum demands more dissipation before matter plunges to the horizon. However, compared to mini-disks in larger separation binaries with spinning black holes, our mini-disks are less luminous: unlike those systems, their mass accretion rate is lower than in the circumbinary disk, and they radiate with lower efficiency because their inflow times are shorter. Compared to a single black hole system matched in mass and accretion rate, these binaries have spectra noticeably weaker and softer in the UV. Finally, we discuss the implications of our findings for the potential observability of these systems.

Eduardo M. Gutiérrez↗

Comparison of MPEG digital video with super VHS tape for diagnostic echocardiographic readings

BACKGROUND: Digital recording of echocardiographic studies is on the clinical horizon. However, full digital capture of complete echocardiographic studies in traditional video format is impractical, given current storage capacity and network bandwidth. To overcome these constraints, we evaluated the diagnostic image quality of digital video by using MPEG (Motion Picture Experts Group) compression. METHODS AND RESULTS: Fifty-eight complete, consecutive studies were recorded simultaneously with the use of MPEG-1 and sVHS videotape. Each matched MPEG and sVHS study pair was reviewed by two from a total of six readers, and findings were recorded with the use of a detailed, computerized reporting tool. Intrareader and interreader discrepancies were characterized as major or minor and analyzed in total and for specific subgroups of findings (left and right ventricular parameters, valvular insufficiency, and left ventricular regional wall motion). Intrareader discrepancies were reviewed by a consensus panel for agreement with either MPEG or sVHS findings. There was an exact concordance between MPEG and sVHS readings in 83% of findings. The majority of discrepancies were minor, with major discrepancies in only 2.7% of findings. There was no difference in the rate of consensus panel agreement with MPEG or sVHS for instances of intrareader discrepancy, either in total or for any subgroup of findings. Interreader discrepancy rates were nearly identical for both MPEG and sVHS. CONCLUSIONS: MPEG-1 digital video is equivalent to sVHS videotape for diagnostic echocardiography. MPEG increases the range of practical options for digital echocardiography and offers, for the first time, the advantages of digital recording in a familiar video format.

Non-NASA Center↗

Integrated Demand Management: Concepts and Procedures

This report provides a comprehensive description of the Integrated Demand Management concept. Motivation: NASA’s Integrated Demand Management (IDM) research explores the idea that, under certain conditions, time-based flow management (TBFM) arrival operations can benefit from the coordinated use of a strategic traffic management initiative (TMI) to “precondition” the inbound demand. The research was motivated by the observation that TBFM was usually turned off during convective weather, even in facilities where it was routinely used. Our hypothesis was that strategic adjustments to the inbound traffic so that it provided a better match to the off-nominal changes in capacity observed during these conditions could enable TBFM scheduling to continue to provide effective support for arrival traffic management. Concept: IDM proposes that a TMI (e.g., a Collaborative Trajectory Options Program, or CTOP) be used to adjust the rate and/or geographic distribution across flows of the traffic inbound to a high-demand, TBFM-managed airport before that traffic reaches the TBFM planning horizon. After this strategic preconditioning, TBFM can then tactically fine-tune the demand to deliver a well-managed, orderly feed to the destination airport. Coordinated use of these two flow management capabilities is intended to improve system performance in terms of: • Equity of ground delay assignment, avoiding excessive ground delay for TBFM-scheduled departures, without penalizing longer flights; • Throughput, by distributing traffic to maximize use of available capacity; • Predictability for operators, providing advance notice about the impact on individual flights; • Increased flexibility, supporting operator mitigation strategies such as slot swapping or trajectory options; • Efficiency of flight operations, using ground delay more effectively and reducing airborne delay. The operational description in this document highlights how the IDM concept builds upon already existing tools and procedures, and also indicates where tool enhancements could facilitate conduct of IDM operations. However, enhanced tools are not a requirement for concept introduction. In fact, initial deployment that focused on training procedures and rationale for coordinated use of TFMS and TBFM, without changes to existing tools, might be a simpler way to introduce and to familiarize traffic managers with the idea of preconditioning. The concept and procedures described in this document can hopefully provide useful guidance for introduction of IDM into field operations..

IDM↗

Environmental modeling and recognition for an autonomous land vehicle

An architecture for object modeling and recognition for an autonomous land vehicle is presented. Examples of objects of interest include terrain features, fields, roads, horizon features, trees, etc. The architecture is organized around a set of data bases for generic object models and perceptual structures, temporary memory for the instantiation of object and relational hypotheses, and a long term memory for storing stable hypotheses that are affixed to the terrain representation. Multiple inference processes operate over these databases. Researchers describe these particular components: the perceptual structure database, the grouping processes that operate over this, schemas, and the long term terrain database. A processing example that matches predictions from the long term terrain model to imagery, extracts significant perceptual structures for consideration as potential landmarks, and extracts a relational structure to update the long term terrain database is given.

Lawton, D. T.↗

The Nustar Spectrum of Mrk 335: Extreme Relativistic Effects Within Two Gravitational Radii of the Event Horizon?

We present 3-50 keV NuSTAR observations of the active galactic nuclei Mrk 335 in a very low flux state. The spectrum is dominated by very strong features at the energies of the iron line at 5-7 keV and Compton hump from 10-30 keV. The source is variable during the observation, with the variability concentrated at low energies, which suggesting either a relativistic reflection or a variable absorption scenario. In this work, we focus on the reflection interpretation, making use of new relativistic reflection models that self consistently calculate the reflection fraction, relativistic blurring and angle-dependent reflection spectrum for different coronal heights to model the spectra. We find that the spectra can be well fitted with relativistic reflection, and that the lowest flux state spectrum is described by reflection alone, suggesting the effects of extreme light-bending occurring within approx. 2 gravitational radii (RG) of the event horizon. The reflection fraction decreases sharply with increasing flux, consistent with a point source moving up to above 10 RG as the source brightens. We constrain the spin parameter to greater than 0.9 at the 3(sigma) confidence level. By adding a spin-dependent upper limit on the reflection fraction to our models, we demonstrate that this can be a powerful way of constraining the spin parameter, particularly in reflection dominated states. We also calculate a detailed emissivity profile for the iron line, and find that it closely matches theoretical predictions for a compact source within a few RG of the black hole.

galaxies: active aEuro" galaxies: individual: Mrk ↗

Lunar science strategy: Exploring the Moon with humans and machines

Important scientific questions that can be addressed from the lunar surface are reviewed for a number of scientific disciplines. A successful strategy for human exploration of the Moon is outlined. It consists of several elements: thorough preparation; a means of extending the human reach; measurement of the mobility of both human and robotic components; and flexible technologies so as to be able to take the most effective path as successive decision points occur. Part of thorough preparation involves concurrent development of a set of science goals and objectives as well as a supporting information base; neither can evolve independently of the other. This matched set will drive the definition of missions and technologies used to satisfy the requirements of various science disciplines. No single site on the Moon will satisfy all requirements. Thus, global accessibility is a goal of the current Lunar and Mars Exploration Program science strategy. Human mobility on the surface is limited to a few kilometers without the use of vehicles. Unpressurized crew carrying rovers could take explorers to distances tens of kilometers from an outpost; the distance is primarily limited by health and safety concerns. Pressurized rovers could extend the range to hundreds of kilometers, but size, mass, and costs limit such vehicles to missions beyond current planning horizons. The establishment of several outposts instead of one would provide opportunities for effective use of the unique capabilities of humans. Extending the human reach to global dimensions may be accomplished through teleoperation or telepresence. The most effective mix of these techniques is a decision that will evolve as experience is gained on the surface. Planning and technology must be flexible enough to allow a variety of options to be selected.

Morrison, Donald A.↗

Estimates of Minimum Energy Requirements for Range-Controlled Return of a Nonlifting Satellite from a Circular Orbit

Existing expressions are used to obtain the minimum propellant fraction required for return from a circular orbit as a function of vacuum trajectory range. trajectory are matched to those of the atmospheric trajectory to obtain a complete return from orbit to earth. The results are restricted by the assumptions of (1) impulsive velocity change, (2) nearly circular transfer trajectory, ( 3) spherical earth, atmosphere, and gravitational field, (4) exponential atmospheric density variation with attitude, and (5) a nonrotating atmosphere. The solutions for the parameters of the vacuum Calculations are made t o determine the effects of longitudinal and lateral range on required propeUant fraction and reentry loading for a nonrotating earth and for several orbital altitudes. the single- and two-impulse method of return is made and the results indicate a "trade off" between propellant fraction required and landing- position accuracy. A comparison of An example of a return mission from a polar orbit is discussed where the initial deorbit point is the intersection of the North Pole horizon with the satellite orbit. Some effects of a rotating earth are also considered. It is found that, for each target-orbital-plane longitudinal difference, there exists a target latitude for which the required propellant fraction is a minimum.

Jackson, Charlie M., Jr.↗

Investigating Low-Altitude Constellations of Ad-Hoc Lunar PNT System for Distributed Spacecraft Autonomy

In this study, we examine a low-altitude Lunar Position, Navigation, and Timing (LPNT) constellations and the localization performance of Centralized Extended Kalman Filter (CEKF) and Decentralized Extended Kalman Filter (DEKF) algorithms. The primary investigation involves a 100-node swarm operating at a 100 km altitude, in contrast to previous studies that examined a 21-node asset in a frozen-orbit at 5,500 km. The autonomous operation of large-scale swarm is based on two-way Inter-Satellite Link (ISL) measurements, which involve pseudoranges and relative velocities among swarm nodes. We perform a numerical assessment of the two filtering approaches, utilizing ‘fully sampled’ measurements from all available assets as well as ‘two ISL’ measurements where each spacecraft is restricted to only two antennas. This research includes an analysis of CEKF under 2-ISL constraints and evaluates the performance of DEKF in a 100-node swarm, which has not been explored in previous studies. In addition, we examine the impact of increasing the sampling frequency for DEKF, showing that the update cycle can be shortened from a 10-minute interval. A novel approach for ‘2-ISL limited’ DEKF will also be introduced, using a matching formulation that exhaustively enumerates all potential matches. This study provides valuable insights into large-scale distributed swarm operations, considering various filter configurations, sampling frequencies, matching strategies, and scalability of CEKF and DEKF for low-altitude LPNT applications. The Lunar PNT technology plays a key role in providing reliable and robust navigation services on the Moon's surface and the South pole, where the primary Lunar missions are planned. To support upcoming Lunar missions, including small satellites from NASA's Commercial Lunar Payload Services program, the Lunar PNT system must be adaptable to smaller platforms like CubeSats. Driven by the growing involvement of public and private exploration partnerships, the traditional low Earth orbit missions are shifting to beyond geosynchronous orbit [1]. These upcoming missions aim to foster a sustainable and innovative exploration program, in collaboration with commercial and international partners, to facilitate human expansion throughout the solar system and return new knowledge and opportunities to Earth [2]. As part of this trend, there are increasing efforts to utilize science missions in Lunar orbit to develop a non-dedicated and ad-hoc PNT network system. Two traditional approaches, the Deep Space Network (DSN) and the weak signal Global Positioning System (GPS), are established deep-space navigation technologies for missions beyond the geosynchronous orbit. Beginning in 1958, the DSN was developed to communicate with the Explorer 1 spacecraft based on the use of radiometric tracking in spacecraft navigation [3]. The DSN is capable of providing nearly unfettered coverage to spacecraft beyond low-Earth orbit (LEO), however, increased space mission volume has created concerns about future expectations of DSN usage for spacecraft navigation [4]. For cislunar mission applications, the position accuracy using DSN achieves 100 m (3σ) with at least three geometrically diverse ground stations when using radiometric tracking alone [5]. The DSN's dependence on Earth-based ground stations restricts its operational capabilities to periods of Earth visibility. This limitation, coupled with its poor localization performance, renders the DSN unsuitable for future lunar missions that demand continuous tracking and precise positioning. To satisfy the increasing requirements of DSN in Lunar applications, spacecrafts are also required to improve their onboard antenna power and efficiency of the transmission. However, there is an important aggregate cost trade between adding capabilities to every spacecraft and adding to a capacity on the ground that serves multiple spacecraft [6]. A weak GPS system can provide PNT service while the user spacecraft is bound to the Moon, leveraging a single, steerable high gain antenna with the relatively narrow beam which includes all the sources in its field of view [7]. However, the higher the altitude the receiver is above the GPS constellations, the poorer and the weaker are the relative geometry and the received signal powers, respectively, leading to a significant navigation accuracy reduction [8]. The transmitted power becomes weaker with increasing distance from the Earth as well as signals tracked from one of the side lobes of the GPS antenna pattern. As a results, the number of visible satellites and relative geometric condition of the GPS satellites at very high altitude drops dramatically and reduces the navigation solution accuracy. Therefore, the weak GPS system is also not an ideal way to provide PNT service to upcoming Lunar missions when considering its limited geometric condition and the recued navigation accuracy. Another navigation approach on the Moon is being developed, similar to the Global Navigation Satellite System (GNSS) on Earth, aiming to offer navigation service with continuous 24/7 coverage across the entire Lunar surface. For example, lunar communications relay and navigation systems (LCRNS) by NASA and Lunar navigation satellite systems (LNSS) by JAXA are designed to serve as dedicated Position, Navigation, and Timing (PNT) systems for the Moon. However, designing a dedicated LNSS and PNT service involves additional challenges, which are unique to the lunar environment, including limited payload capacity for the CubeSat platform, i.e., the size, weight, and power (SWaP) of the onboard clock, limited lunar ground monitoring stations, and limited financial investment as compared to the legacy Earth-GPS [9]. NASA’s focus on utilizing CubeSat platforms on the Moon leads to an alternative Lunar navigation platform that leverages the existing Lunar science and exploration assets. The small satellites used in Lunar missions can be used to create a low-cost, autonomous, ad-hoc, and on-demand mission-centric Lunar PNT swarm capable of providing PNT services to these low-cost lunar missions [10]. As upcoming Lunar missions will often operate at low-altitude about 30 km to 100 km for scientific observations and mapping purposes, the low-altitude orbital constellations could be employed to create an ad-hoc Lunar PNT system. However, several issues must be addressed, such as the instability of these orbits, which often require maintenance or are only suitable for short-duration missions, operating for fewer than 90 days. Additionally, at an altitude of 100 km, the satellites have a limited period during which they are above the horizon and capable of providing PNT service to users. The implementation of a non-dedicated, ad-hoc Lunar navigation constellation facilitates on-demand PNT services. A preliminary study of ad-hoc Lunar PNT system was conducted using 21 spacecraft in 5,5000 km altitude frozen orbits to test its feasibility and a basic performance of orbital asset localization among ad-hoc Lunar constellations in small satellites format [10]. These swarm assets are designed for autonomous localization with minimal Earth interaction, reducing dependency on bandwidth and ground resources. The design in [10] demonstrated the feasibility of a decentralized PNT approach, specifically employing a DEKF approach for state estimation, which helps minimize onboard operating costs. The DEKF method distributes computation across individual satellites, which lightens the computational load while maintaining accuracy in orbit ephemeris and clock offsets, similar to centralized systems [11]. In a follow-on study [12], each spacecraft was limited to 2 communications antennae, forcing the selection of measurements and scheduling spacecraft activities to perform the measurements. A matching algorithm is implemented to select the best measurements and schedule position estimation updates. The decentralized localization performance is also investigated with increasing levels of network degradation for swarm assets considering the impact of intermittent and permanent communication failure, to demonstrate the robustness and fidelity of the decentralized Lunar PNT service [13]. This study confirmed that the ad-hoc PNT constellations in frozen orbit are highly robust and resilient to communication failures. However, unlike frozen orbit swarm assets, the low-altitude satellites have a limited ground view at an altitude of 100 km, where the ad-hoc Lunar constellation consists of 98 low-altitude satellites, evenly distributed across seven circular polar orbital planes, alongside two satellites in a frozen orbit at an altitude of 5,500 km (Figure 1). Therefore, the number of satellites visible to ground users is significantly limited in low-altitude orbit constellations. As each visibility of a spacecraft remains intact for only a few ticks before it moves out of the field of view, the ground user encounters challenges in maintaining continuous navigation service, resulting in sparse availability and provision of Lunar PNT system. Consequently, service availability is primarily restricted to the Lunar South Pole region (Figure 2). Given these limitations and concerns, the localization performance of low-altitude swarm assets will be assessed in this study. We focus on the investigation of the localization performance of low-altitude swarm assets and ground users near the Lunar South Pole. The overall flow of the Lunar PNT simulation incorporates the DEKF approach of asset localization and the weighted least-squares approach in user localization (Figure 3). The autonomous Lunar PNT simulation is primarily implemented in MATLAB, where the DEKF based on the matching scheduler is implemented with Google’s OR-tools as a model builder and Gurobi optimization tool as a backend solver. The General Mission Analysis Tool (GMAT) is utilized to generate ephemeris data for swarm assets, and accounts for satellite orbital details, mass, and perturbations like solar radiation pressure and drag coefficients. Each ephemeris dataset is produced in the Moon International Celestial Reference Frame (ICRF) inertial coordinate system. For state estimation, the distributed swarm assets rely on two-way Inter-Satellite Link (ISL) measurements, which involve tracking pseudoranges and relative velocities between visible satellites and anchor nodes during each observation. Numerical evaluations of the decentralized localization process are conducted to demonstrate the feasibility of the low-altitude PNT system in providing reliable navigation services. The main approach involves using DEKF and CEKF to localize 100 satellites in low-altitude constellations, where the CEKF is implemented to serve as a baseline for comparing the performance of distributed algorithms. In both cases, we evaluate ‘fully sampled’ measurements from all available assets, and ‘two ISL’ measurements when spacecraft are constrained to have only two antennas. We test four estimation techniques: CEKF fully sampled, CEKF two ISL, DEKF fully sampled, and DEKF two ISL filters. As the DEKF update cycle is comprised of network setup, communication, and computations, a global broadcast network and 2-way ISL network setup will take from 4 to 6 minutes as maximum [12]. In this simulation, the DEKF update cycle is set to 10 minutes, including a 4-minute latency for obtaining and computing the actual measurement updates. We experiment an increased update cycle to demonstrate the feasibility and evaluate the impact on localization performance using various tuning values for measurement noise covariances (Figures 4 and 5). By comparing centralized and decentralized approaches using a matching algorithm, we analyze the influence of cross-correlation factors in the covariance matrix, assuming 100% reliability of all assets and measurements. The increased frequency and the adjustments of tuning parameters reveal distinct error patterns between the two scenarios. The localization accuracy of the swarm assets and ground users is assessed by taking the median error across 100 assets and one ground user (84.9°S, 137.5°E) over 7-day simulation period (Table 1). Since the user localization accuracy is significantly affected by the performance of the swarm assets, it is crucial to maintain high localization accuracy within the swarm. This study will continue to explore decentralized filtering for autonomous LPNT operations, with further investigation of an 'iterative' matching approach which enumerates every valid matching pair, planned for the following month.

Yeji Kim↗