Spacecraft Data Research Articles

The Bioburden and Ionic Composition of Hypersaline Lake Ices: Novel Habitats on Earth and Their Astrobiological Implications.

We present thermophysical, biological, and chemical observations of ice and brine samples from five compositionally diverse hypersaline lakes in British Columbia's interior plateau. Possessing a spectrum of magnesium, sodium, sulfate, carbonate, and chloride salts, these low-temperature high-salinity lakes are analogs for planetary ice-brine environments, including the ice shells of Europa and Enceladus and ice-brine systems on Mars. As such, understanding the thermodynamics and biogeochemistry of these systems can provide insights into the evolution, habitability, and detectability of high-priority astrobiology targets. We show that biomass is typically concentrated in a layer near the base of the ice cover, but that chemical and biological impurities are present throughout the ice. Coupling bioburden, ionic concentration, and seasonal temperature measurements, we demonstrate that impurity entrainment in the ice is directly correlated to ice formation rate and parent fluid composition. We highlight unique phenomena, including brine supercooling, salt hydrate precipitation, and internal brine layers in the ice cover, important processes to be considered for planetary ice-brine environments. These systems can be leveraged to constrain the distribution, longevity, and habitability of low-temperature solar system brines-relevant to interpreting spacecraft data and planning future missions in the lens of both planetary exploration and planetary protection.

Determination of coronal mass ejection orientation and consequences for their propagation

Context. The configuration of the interplanetary magnetic field and features of the related ambient solar wind in the ecliptic and meridional plane are different. Therefore, one can expect that the orientation of the flux-rope axis of a coronal mass ejection (CME) influences the propagation of the CME itself. However, the determination of the CME orientation, especially from image data, remains a challenging task to perform. Aim. This study aims to provide a reference to different CME orientation determination methods in the near-Sun environment. Also, it aims to investigate the non-radial flow in the sheath region of the interplanetary CME (ICME) in order to provide the first proxy to relate the ICME orientation with its propagation. Methods. We investigated 22 isolated CME-ICME events in the period 2008–2015. We determined the CME orientation in the near-Sun environment using the following: (1) a 3D reconstruction of the CME with the graduated cylindrical shell (GCS) model applied to coronagraphic images provided by the STEREO and SOHO missions; and (2) an ellipse fitting applied to single spacecraft data from SOHO/LASCO C2 and C3 coronagraphs. In the near-Earth environment, we obtained the orientation of the corresponding ICME using in situ plasma and field data and also investigated the non-radial flow in its sheath region. Results. The ability of GCS and ellipse fitting to determine the CME orientation is found to be limited to reliably distinguish only between the high or low inclination of the events. Most of the CME-ICME pairs under investigation were found to be characterized by a low inclination. For the majority of CME-ICME pairs, we obtain consistent estimations of the tilt from remote and in situ data. The observed non-radial flows in the sheath region show a greater y direction to z direction flow ratio for high-inclination events, indicating that the CME orientation could have an impact on the CME propagation.

ULF Wave‐Induced Ion Pitch Angle Evolution in the Dayside Outer Magnetosphere

AbstractDrift‐bounce resonance between ultralow frequency (ULF) waves and ions is essential for ion energization in the magnetosphere. Here, we present the first comprehensive study of drift‐bounce resonance in the dayside outer magnetosphere, where off‐equatorial magnetic field minima would strongly distort ions' bounce and drift motion. A generalized theory is proposed, in which the effects of off‐equatorial minima, time‐evolving fields and ion bounce motion are taken into account. In consequence of these effects, ion pitch angle distributions undergo dramatic changes. In the presence of off‐equatorial minima, the time‐of‐flight effect of ion bounce motion forms latitude‐dependent dispersions besides “paw‐track shaped” structures, while evolving wave fields cause time‐dependent phase shifts in “paw‐tracks.” All the predicted signatures have been confirmed by 5 years of Magnetospheric Multiscale spacecraft data and numerical simulations. These results allow us to better understand the interactions between ULF waves and thermal ion species in global magnetospheric dynamics.

Arecibo S-band Radar Characterization of Local-scale Heterogeneities within Mercury’s North Polar Deposits

Abstract Ground-based planetary radar observations first revealed deposits of potentially nearly pure water ice in some permanently shadowed regions (PSRs) on Mercury’s poles. Later, the MESSENGER spacecraft confirmed the icy nature of the deposits, as well as their location within PSRs. Considering the geologic context provided by MESSENGER, we further characterized the north polar deposits by pairing spacecraft data with new Arecibo S-band radar observations. Here we show that some ice deposits within PSRs have a gradational pattern in their radar properties that is likely associated with differences in ice purity. Radar-bright features with a circular polarization ratio μ c > 1 can be characterized by water ice with ≳3% impurities by volume while those with μ c < 1 by ≳20% impurities. Furthermore, areas in PSRs with μ c < 1 typically surround locations of stronger radar backscatter with μ c > 1. Therefore, deposits of nearly pure water ice are likely surrounded by lower-purity material, such as water-ice-rich regolith, which could be the result of impact gardening or the crater’s thermal environment. However, such deposits are not always colocated within large polar craters where ice should be the most stable, even at the surface. In fact, we found that there is no significant difference between the radar backscattering properties of deposits thought to have surficial ice and those with buried ice. Our results also help improve the identification of icy reservoirs elsewhere, such as the Moon. Indeed, we found that μ c is not an adequate diagnostic, but rather the radar backscatter in each circular polarization independently provides information to identify water-ice deposits.

Geodetic investigations of the mission concept MAGIC to reveal Callisto's internal structure

Geodetic and geophysical investigations of the Galilean moon Callisto address fundamental questions regarding the formation and evolution of the Jovian system. Callisto's evolution and internal structure appear to significantly differ from the other Jovian satellites. Similarly-sized Ganymede is a highly evolved ice-rock moon with a differentiated interior, intrinsic magnetic field, and abundant surface evidence of internal activity. In contrast, Callisto's surface is ancient, and Galileo spacecraft data suggest its interior is only incompletely differentiated, despite the presumed presence of a sub-surface ocean. These properties make Callisto uniquely able to constrain the timing and nature of the Jovian system formation. The Magnetics, Altimetry, Gravity, and Imaging of Callisto (MAGIC) mission concept is conceived to fully characterize the properties of this enigmatic moon from its deep interior to the icy shell. Three main instruments are included as a scientific payload. Highly accurate measurements of Callisto's topography, magnetic field, and morphology are obtained by the onboard laser altimeter, magnetometer, and camera, respectively. The telecommunication system supports an additional gravity and radio science investigation. Long- and short-wavelength gravity anomalies afford powerful constraints on internal differentiation and the properties of the hydrosphere (water and ice). Comprehensive numerical simulations and covariance analyses of MAGIC mission scenarios presented in this paper show that the gravitational degree-2 normalized coefficients and the pole obliquity enable the determination of the moment of inertia with an accuracy better than 0.015%. The combination of gravity and altimetry measurements acquired by MAGIC are essential to the characterization of Callisto's interior if – as is likely – the degree-2 gravity includes non-hydrostatic terms. MAGIC's radio science data yield the estimation of Callisto's gravity field with spatial resolutions of <100 km. The combination of gravitational and deformation tides that are retrieved by the radio science and altimetry investigations, respectively, leads to the recovery of the rigid ice shell thickness to within ∼3 km. Together these datasets would resolve ambiguities inherent in Galileo flyby data, revealing Callisto's interior structure as well as the existence and properties of its postulated internal ocean.

Relativistic Electron Enhancements Through Successive Dipolarizations During a CIR‐Driven Storm

AbstractRelativistic electrons in the Earth's radiation belts are highly dynamic on a variety of timescales during the geomagnetic storm. Using Van Allen Probe spacecraft data, we investigate rapid enhancements of relativistic electrons in the outer radiation belt during a corotating interaction region (CIR) driven storm. Successive dipolarizations associated with 100keV‐MeV electron injections are identified. The evolution of energetic electrons is analyzed in the space of adiabatic invariants (μ, K and L*). Within less than a few hours, the phase space density (PSD) of the relativistic electrons promptly increases corresponding to injections of MeV electrons. The PSD of MeV electrons cumulatively increases by a factor of 4–10 at L* = 4.5–5.8 which is likely due to successive groups of dipolarizations and injections. Both near‐equatorial (small K) and off‐equatorial (large K) energetic electrons increase significantly. The increases in the near‐equatorial electrons are still dominant, suggesting the operation of betatron acceleration. The event study shows that successive dipolarizations associated with the CIR‐driven storm may rapidly affect relativistic electrons of the outer radiation belt over a wide range in the phase space.

A Remote Sensing Perspective on Mass Wasting in Contrasting Planetary Environments: Cases of the Moon and Ceres

Mass wasting, as one of the most significant geomorphological processes, contributes immensely to planetary landscape evolution. The frequency and diversity of mass wasting features on any planetary body also put engineering constraints on its robotic exploration. Mass wasting on other Solar System bodies shares similar, although not identical, morphological characteristics with its terrestrial counterpart, indicating a possible common nature for their formation. Thus, planetary bodies with contrasting environmental conditions might help reveal the effects of the atmosphere, subsurface fluids, mass accumulation/precipitation, and seismicity on mass wasting, and vice versa. Their relative positions within our Solar System and the environmental and geophysical conditions on the Moon and the dwarf planet Ceres are not only extremely different from Earth’s but from each other too. Their smaller sizes coupled with the availability of global-scale remote sensing datasets make them ideal candidates to understand mass wasting processes in widely contrasting planetary environments. Through this concept article, we highlight several recent advances in and prospects of using remote sensing datasets to reveal unprecedented details on lunar and Cerean mass wasting processes. We start with briefly discussing several recent studies on mass wasting using Lunar Reconnaissance Orbiter Camera (LROC) data for the Moon and Dawn spacecraft data for Ceres. We further identify the prospects of available remote sensing data in advancing our understanding of mass wasting processes under reduced gravity and in a scant (or absent) atmosphere, and we conclude the article by suggesting future research directions.

The Need for a System Science Approach to Global Magnetospheric Models

This perspective advocates for the need of a combined system science approach to global magnetospheric models and to spacecraft magnetospheric data to answer the question “Do simulations behave in the same manner as the magnetosphere does?” (instead of the standard validation question “How well do simulations reproduce spacecraft data?”). This approach will 1) validate global magnetospheric models statistically, without the need for a direct comparison against spacecraft data, 2) expose the deficiencies of the models, and 3) provide physics support to the system analysis performed on the magnetospheric system.

Reconstruction of Mercury's internal magnetic field beyond the octupole

Abstract. The reconstruction of Mercury's internal magnetic field enables us to take a look into the inner heart of Mercury. In view of the BepiColombo mission, Mercury's magnetosphere is simulated using a hybrid plasma code, and the multipoles of the internal magnetic field are estimated from the virtual spacecraft data using three distinct reconstruction methods: the truncated singular value decomposition, the Tikhonov regularization and Capon's minimum variance projection. The study shows that a precise determination of Mercury's internal field beyond the octupole up to the dotriacontapole is possible and that Capon's method provides the same high performance as the Tikhonov regularization, which is superior to the performance of the truncated singular value decomposition.

The Wave Telescope Technique

AbstractUnambiguous determination of wavevectors (wavelengths and propagation directions) is one of the major challenges of space plasma observations, since conventional wave analysis methods for the single spacecraft data are limited to the time series data analysis and the wave properties are determined as a function of the frequencies in the spacecraft frame. The wave telescope technique overcomes the problem of wavevector determination by making extensive use of the multi‐spacecraft interferometric data. The wave telescope techniques do not assume any wave mode or plasma mode, and is suited to testing for the wave and turbulence theories against the multi‐spacecraft data. This article presents an updated report on the wave telescope technique on the basis of adaptive filtering (Capon's minimum variance projection), the construction of an estimator for the multi‐point magnetic field data, the performance and limits in the wavevector analysis (aliasing and low‐wavenumber behavior), the applications to the Cluster fluxgate magnetometer data (dispersion relations, phase velocity diagrams, propagation pattern, energy spectra, helicities, and higher‐order statistics), and possible extensions of the technique.

Alkali magmatism on Mars: an unexpected diversity

Despite an apparent north/south topographic dichotomy that formed >4.0 Ga, the young Martian meteorites (<2.4 Ga) and first-order remote sensing observations revealed a surface of Mars that is uniformly basaltic. This simplistic vision has been challenged by the discovery of a brecciated meteorite and additional spacecraft data that all point to the presence of alkaline igneous rocks, thereby demonstrating an unexpected igneous diversity on Mars. In the present paper, we review a variety of effusive alkaline rocks (basalts to trachytes) recognized so far in the southern hemisphere of Mars as observed from a unique 4.47 Ga Martian meteorite, as well as ground, and orbital data. The complementary of effusive alkaline rocks and plutonic orthopyroxene-rich rocks in early Mars is discussed. We propose that mantle-derived magmas at high extent of melting at rather low pressure either erupted forming orthopyroxene-rich lavas, or crystallized at shallow crustal depths, fractionating orthopyroxene which sank to the bottom of the chamber and residual alkaline magmas which erupted at the surface of Mars. Widespread low pressure fractionation processes could also be related to heavy bombardment on the early Martian crust generating melt sheets that ultimately differentiated. The Noachian crust is more diverse than being merely basaltic.

On-orbit performance of the Spitzer Space Telescope: science meets engineering

The Spitzer Space Telescope operated for over 16 years in an Earth-trailing solar orbit, returning not only a wealth of scientific data but, as a by-product, spacecraft and instrument engineering data which will be of interest to future mission planners. These data will be particularly useful because Spitzer operated in an environment essentially identical to that at the L2 LaGrange point where many future astrophysics missions will operate. In particular, the radiative cooling demonstrated by Spitzer has been adopted by other infrared space missions, from JWST to SPHEREx. This paper aims to facilitate the utility of the Spitzer engineering data by collecting the more unique and potentially useful portions into a single, readily-accessible publication. We avoid discussion of less unique systems, such as the telecom, flight software, and electronics systems and do not address the innovations in mission and science operations which the Spitzer team initiated. These and other items of potential interest are addressed in references supplied in an appendix to this paper.

A Fast Bow Shock Location Predictor-Estimator From 2D and 3D Analytical Models: Application to Mars and the MAVEN Mission.

We present fast algorithms to automatically estimate the statistical position of the bow shock from spacecraft data, using existing analytical two‐dimensional (2D) and three‐dimensional (3D) models of the shock surface. We derive expressions of the standoff distances in 2D and 3D and of the normal to the bow shock at any given point on it. Two simple bow shock detection algorithms are constructed, one solely based on a geometrical predictor from existing models, the other using this predicted position to further refine it with the help of magnetometer data, an instrument flown on many planetary missions. Both empirical techniques are applicable to any planetary environment with a defined shock structure. Applied to the Martian environment and the NASA/MAVEN mission, the predicted shock position is on average within 0.15 planetary radius R p of the bow shock crossing. Using the predictor‐corrector algorithm, this estimate is further refined to within a few minutes of the true crossing (≈0.05R p). Between 2014 and 2021, we detect 14,929 clear bow shock crossings, predominantly quasi‐perpendicular. Thanks to 2D conic and 3D quadratic fits, we investigate the variability of the shock surface with respect to Mars Years (MY), solar longitude (Ls), and solar EUV flux levels. Although asymmetry in Y and Z Mars Solar Orbital coordinates is on average small, we show that for MY32 and MY35, Ls = [135°−225°] and high solar flux, it can become particularly noticeable, and is superimposed to the usual North‐South asymmetry due in part to the presence of crustal magnetic fields.

The Source of Auroral Omegas

AbstractThe auroral wave‐like structures called “omega bands” appear within the post‐midnight sector auroral oval with shapes resembling the Greek letter omega, and are typically associated with the recovery phase of substorms. Prior work and MHD simulations suggest both high speed earthward flows and post‐midnight flow shears are possible omega band source mechanisms. However, what produces omega bands is not well understood. It is most likely that the paucity of concurrent magnetospheric data has limited our ability to understand fully the mechanism responsible for the generation of the omega bands. We have identified about 263 auroral omegas in seven different THEMIS all sky cameras (ASCs) over a 10 year period. A fraction appear to form from north‐south streamers, but some appear to form from east‐west auroral arcs. Fifty‐one of the 263 Ω also have conjugate or near conjugate THEMIS and GOES spacecraft data. There is evidence in about 80% of these events that high speed earthward flows have occurred prior to or at about the same time as the auroral omega observation. This evidence consists of high speed earthward flows, magnetic field dipolarizations, particle injections into the inner magnetosphere, and auroral streamers. 11 events also show plasma flow data azimuthally along the inner magnetosphere, but for six of these events the Vy component is inconsistent with the auroral omega direction of motion. Of the remaining five events, four have Vy flow shear events but also show evidence of magnetic field dipolarizations, particle injections into the inner magnetosphere, and auroral streamers. Our observations suggest high speed earthward flows are more likely to be the source of auroral omega bands.

Numerical and experimental evidence for a new interpretation of residence times in space

Aims. We investigate the energy dependence of Jovian electron residence times, which allows for a deeper understanding of adiabatic energy changes that occur during charged particle transport, as well as of their significance for simulation approaches. Thereby we seek to further validate an improved approach to estimate residence times numerically by investigating the implications on previous analytical approaches and possible effects detectable by spacecraft data. Methods. Utilising a propagation model based on a stochastic differential equation (SDE) solver written in CUDA, residence times for Jovian electrons were calculated over the whole energy range dominated by the Jovian electron source spectrum. We analysed the interdependences both with the magnetic connection between the observer and the source as well as between the distribution of the exit (simulation) times and the resulting residence times. Results. We point out a linear relation between the residence time for different kinetic energies and the longitudinal shift of the 13 month periodicity typically observed for Jovian electrons and discuss the applicability of these findings to data. Furthermore, we utilise our finding that the simulated residence times are approximately linearly related to the energy loss for Jovian and Galactic electrons, and we develop an improved analytical estimation in agreement with the numerical residence time and the longitudinal shift observed by measurements.

Variations of flow direction in solar wind streams of different types

Studying the direction of the solar wind flow is a topical problem of space weather forecasting. As a rule, the quiet and uniform solar wind propagates radially, but significant changes in the solar wind flow direction can be observed, for example, in compression regions before the interplanetary coronal mass ejections (Sheath) and Corotating Interaction Regions (CIR) that precede high-speed streams from coronal holes. In this study, we perform a statistical analysis of the longitude (φ) and latitude (θ) flow direction angles and their variations on different time scales (30 s and 3600 s) in solar wind large-scale streams of different types, using WIND spacecraft data. We also examine the relationships of the value and standard deviations SD of the flow direction angles with various solar wind parameters, regardless of the solar wind type.&#x0D; We have established that maximum values of longitude and latitude angle modulus, as well as their variations, are observed for Sheath, CIR, and Rare, with the probability of large deviations from the radial direction (&gt;5°) increasing. The dependence on the solar wind type is shown to decrease with scale. We have also found that the probability of large values of SD(θ) and SD(φ) increases with increasing proton temperature (Tp) in the range 5–10 eV and with increasing proton velocity (Vp) in the range 400–500 km/s.

Вариации направления потока в различного типа течениях солнечного ветра

Studying the direction of the solar wind flow is a topical problem of space weather forecasting. As a rule, the quiet and uniform solar wind propagates radially, but significant changes in the solar wind flow direction can be observed, for example, in compression regions before the interplanetary coronal mass ejections (Sheath) and Corotating Interaction Regions (CIR) that precede high-speed streams from coronal holes. In this study, we perform a statistical analysis of the longitude (φ) and latitude (θ) flow direction angles and their variations on different time scales (30 s and 3600 s) in solar wind large-scale streams of different types, using WIND spacecraft data. We also examine the relationships of the value and standard deviations SD of the flow direction angles with various solar wind parameters, regardless of the solar wind type.&#x0D; &#x0D; We have established that maximum values of longitude and latitude angle modulus, as well as their variations, are observed for Sheath, CIR, and Rare, with the probability of large deviations from the radial direction (&gt;5°) increasing. The dependence on the solar wind type is shown to decrease with scale. We have also found that the probability of large values of SD(θ) and SD(φ) increases with increasing proton temperature (Tp) in the range 5–10 eV and with increasing proton velocity (Vp) in the range 400–500 km/s.

The study of the strongest solar event on a minimum of the 24th solar cycle

The strongest solar flares of the 24th solar cycle erupted on September 6, 2017, and it was the 8th strongest solar flare recorded since 1996. This extreme solar flare occurred at the minimum of the 24th solar cycle. The active region is located in the Western Hemisphere and produced the violent explosion of class X9.3 and X2.2 on September 6, X1.3 on September 7, and X8.2 on September 10, 2017. The injection duration of the solar energetic particles of the solar event was 17 minutes. All data for this solar event was collected from the Advanced Composition Explorer and simulated for particles’ motion using the transport equation and solved by the numerical technique. We obtained the injection time of the solar energetic particle propagation by comparing fitting between the simulation results and the spacecraft data. Injection time taken by high-energy particles to travel from the Sun to the Earth was found to be in the range of 39 to 743 minutes. At the peak of this solar flare, the coronal mass ejection was detected, which increased the injection time. The Kp-index of this solar flare was 4; thus, there was no effect on the Earth. The Kp-index value increased to 8 on September 7-8, 2017, due to another solar event from the same sunspot region, indicating the effect of solar flare and CME, which resulted in the appearance of aurora.

Analysis of Magnetohydrodynamic Perturbations in the Radial-field Solar Wind from Parker Solar Probe Observations

Abstract We report analysis of sub-Alfvénic magnetohydrodynamic (MHD) perturbations in the low-β radial-field solar wind employing the Parker Solar Probe spacecraft data from 2018 October 31 to November 12. We calculate wavevectors using the singular value decomposition method and separate MHD perturbations into three eigenmodes (Alfvén, fast, and slow modes) to explore the properties of sub-Alfvénic perturbations and the role of compressible perturbations in solar wind heating. The MHD perturbations show a high degree of Alfvénicity in the radial-field solar wind, with the energy fraction of Alfvén modes dominating (∼45%–83%) over those of fast modes (∼16%–43%) and slow modes (∼1%–19%). We present a detailed analysis of a representative event on 2018 November 10. Observations show that fast modes dominate magnetic compressibility, whereas slow modes dominate density compressibility. The energy damping rate of compressible modes is comparable to the heating rate, suggesting the collisionless damping of compressible modes could be significant for solar wind heating. These results are valuable for further studies of the imbalanced turbulence near the Sun and possible heating effects of compressible modes at MHD scales in low-β plasma.

Origin and structure of electromagnetic generator regions at the edge of the electron diffusion region

During the process of magnetic reconnection, converging magnetic fields transfer some of the stored magnetic energy into plasma energy. Simulations and spacecraft data suggest that at the edges of the electron diffusion region, “generators” can form where the local electromagnetic field gains energy from the plasma. We use a 2D Particle-in-cell (PIC) simulation and Magnetospheric Multiscale (MMS) spacecraft data to investigate the electron dynamics and energy budget of these generators. PIC simulation results suggest that these generators form where the electron population begins to remagnetize, while still not being fully frozen in to the magnetic fields. Super-Alfvénic electrons locally lose bulk kinetic energy to Hall electric fields along the outflow axis and to the out-of-plane reconnection electric field as they become remagnetized. These “partial” remagnetization regions also exhibit signatures of irreversible dissipation in MMS data, such as electron enthalpy flux and enhanced electron thermal energy density as remagnetized electrons scatter away from the neutral line within the Hall field structures of the exhaust.