Performance of GPS-based accelerometry: A simulation experiment

Recent studies have shown that with the availability of high-quality CHAMP and GRACE gravity field models, it is feasible to determine accurate non-gravitational accelerations for low Earth orbiting satellites indirectly from precise GPS satellite-to-satellite observations. Possible applications of this so-called GPS-based accelerometry approach consist of accelerometer calibration and atmospheric density and wind computations. With the growing number of high-quality space-borne GPS receivers, this method could be applied to a large range of satellites. In this paper an extensive simulation study has been carried out, based on real accelerometer data from the GRACE mission, in order to determine the optimal processing strategy and the resulting accuracy of the estimated non-gravitational accelerations. It is shown that the optimal processing strategy consists of a piecewise linear parameterization of the estimated empirical accelerations, together with short 6-h orbit arcs. The GPS-based accelerometry approach makes use of triple-differenced GPS observations and the impact of considering the correlated observation noise was found to be marginal in the presence of other error sources such as GPS ephemeris errors. Using a priori non-gravitational force models improves the recovery of low temporal resolution accelerations, except during huge geomagnetic storms. With this strategy, non-gravitational accelerations can be recovered during high solar activity with an accuracy of better than 10% of the total signal in along-track direction and around 25–40% in cross-track direction, at time resolutions of around 8–20 min. During solar minimum conditions, the relative recovery error will increase to approximately 50% in along-track direction and around 60–70% in cross-track direction, due to the reduced atmospheric drag signal. Unfortunately, GPS-based accelerometry is hardly sensitive in the radial direction.     

Spiders: water-driven erosive structures in the southern hemisphere of Mars

Recent data from space missions reveal that there are ongoing climatic changes and erosive processes that continuously modify surface features of Mars. We have investigated the seasonal dynamics of a number of morphological features located at Inca City, a representative area at high southern latitude that has undergone seasonal processes. By integrating visual information from the Mars Orbiter Camera on board the Mars Global Surveyor and climatic cycles from a Mars' General Circulation Model, and considering the recently reported evidence for the presence of water-ice and aqueous precipitates on Mars, we propose that a number of the erosive features identified in Inca City, among them spiders, result from the seasonal melting of aqueous salty solutions.     

A Hitch-hiker’s Guide to Stochastic Differential Equations

In this review, an overview of the recent history of stochastic differential equations (SDEs) in application to particle transport problems in space physics and astrophysics is given. The aim is to present a helpful working guide to the literature and at the same time introduce key principles of the SDE approach via “toy models”. Using these examples, we hope to provide an easy way for newcomers to the field to use such methods in their own research. Aspects covered are the solar modulation of cosmic rays, diffusive shock acceleration, galactic cosmic ray propagation and solar energetic particle transport. We believe that the SDE method, due to its simplicity and computational efficiency on modern computer architectures, will be of significant relevance in energetic particle studies in the years to come.     

Daytime Ionosphere Retrieval Algorithm for the Ionospheric Connection Explorer (ICON)

The NASA Ionospheric Connection Explorer Extreme Ultraviolet spectrograph, ICON EUV, will measure altitude profiles of the daytime extreme-ultraviolet (EUV) OII emission near 83.4 and 61.7 nm that are used to determine density profiles and state parameters of the ionosphere. This paper describes the algorithm concept and approach to inverting these measured OII emission profiles to derive the associated \(\mathrm{O}^{+}\) density profile from 150–450 km as a proxy for the electron content in the F-region of the ionosphere. The algorithm incorporates a bias evaluation and feedback step, developed at the U.S. Naval Research Laboratory using data from the Special Sensor Ultraviolet Limb Imager (SSULI) and the Remote Atmospheric and Ionospheric Detection System (RAIDS) missions, that is able to effectively mitigate the effects of systematic instrument calibration errors and inaccuracies in the original photon source within the forward model. Results are presented from end-to-end simulations that convolved simulated airglow profiles with the expected instrument measurement response to produce profiles that were inverted with the algorithm to return data products for comparison to truth. Simulations of measurements over a representative ICON orbit show the algorithm is able to reproduce hmF2 values to better than 5 km accuracy, and NmF2 to better than 12% accuracy over a 12-second integration, and demonstrate that the ICON EUV instrument and daytime ionosphere algorithm can meet the ICON science objectives which require 20 km vertical resolution in hmF2 and 18% precision in NmF2.     

Extended Pile Driving Model to Predict the Penetration of the Insight/HP<Superscript>3</Superscript> Mole into the Martian Soil

The NASA InSight mission will provide an opportunity for soil investigations using the penetration data of the heat flow probe built by the German Aerospace Center DLR. The Heat flow and Physical Properties Probe (HP³) will penetrate 3 to 5 meter into the Martian subsurface to investigate the planetary heat flow. The measurement of the penetration rate during the insertion of the HP³ will be used to determine the physical properties of the soil at the landing site. For this purpose, numerical simulations of the penetration process were performed to get a better understanding of the soil properties influencing the penetration performance of HP³. A pile driving model has been developed considering all masses of the hammering mechanism of HP³. By cumulative application of individual stroke cycles it is now able to describe the penetration of the Mole into the Martian soil as a function of time, assuming that the soil parameters of the material through which it penetrates are known. We are using calibrated materials similar to those expected to be encountered by the InSight/HP³ Mole when it will be operated on the surface of Mars after the landing of the InSight spacecraft. We consider various possible scenarios, among them a more or less homogeneous material down to a depth of 3–5 m as well as a layered ground, consisting of layers with different soil parameters. Finally we describe some experimental tests performed with the latest prototype of the InSight Mole at DLR Bremen and compare the measured penetration performance in sand with our modeling results. Furthermore, results from a 3D DEM simulation are presented to get a better understanding of the soil response.     

Experimental Investigation of InSight HP<Superscript>3</Superscript> Mole Interaction with Martian Regolith Simulant

The InSight mission launches in 2018 to characterize several geophysical quantities on Mars, including the heat flow from the planetary interior. This quantity will be calculated by utilizing measurements of the thermal conductivity and the thermal gradient down to 5 meters below the Martian surface. One of the components of InSight is the Mole, which hammers into the Martian regolith to facilitate these thermal property measurements. In this paper, we experimentally investigated the effect of the Mole’s penetrating action on regolith compaction and mechanical properties. Quasi-static and dynamic experiments were run with a 2D model of the 3D cylindrical mole. Force resistance data was captured with load cells. Deformation information was captured in images and analyzed using Digitial Image Correlation (DIC). Additionally, we used existing approximations of Martian regolith thermal conductivity to estimate the change in the surrounding granular material’s thermal conductivity due to the Mole’s penetration. We found that the Mole has the potential to cause a high degree of densification, especially if the initial granular material is relatively loose. The effect on the thermal conductivity from this densification was found to be relatively small in first-order calculations though more complete thermal models incorporating this densification should be a subject of further investigation. The results obtained provide an initial estimate of the Mole’s impact on Martian regolith thermal properties.     

Towards a Unified View of Inhomogeneous Stellar Winds in Isolated Supergiant Stars and Supergiant High Mass X-Ray Binaries

Massive stars, at least \(\sim10\) times more massive than the Sun, have two key properties that make them the main drivers of evolution of star clusters, galaxies, and the Universe as a whole. On the one hand, the outer layers of massive stars are so hot that they produce most of the ionizing ultraviolet radiation of galaxies; in fact, the first massive stars helped to re-ionize the Universe after its Dark Ages. Another important property of massive stars are the strong stellar winds and outflows they produce. This mass loss, and finally the explosion of a massive star as a supernova or a gamma-ray burst, provide a significant input of mechanical and radiative energy into the interstellar space. These two properties together make massive stars one of the most important cosmic engines: they trigger the star formation and enrich the interstellar medium with heavy elements, that ultimately leads to formation of Earth-like rocky planets and the development of complex life. The study of massive star winds is thus a truly multidisciplinary field and has a wide impact on different areas of astronomy.In recent years observational and theoretical evidences have been growing that these winds are not smooth and homogeneous as previously assumed, but rather populated by dense “clumps”. The presence of these structures dramatically affects the mass loss rates derived from the study of stellar winds. Clump properties in isolated stars are nowadays inferred mostly through indirect methods (i.e., spectroscopic observations of line profiles in various wavelength regimes, and their analysis based on tailored, inhomogeneous wind models). The limited characterization of the clump physical properties (mass, size) obtained so far have led to large uncertainties in the mass loss rates from massive stars. Such uncertainties limit our understanding of the role of massive star winds in galactic and cosmic evolution.Supergiant high mass X-ray binaries (SgXBs) are among the brightest X-ray sources in the sky. A large number of them consist of a neutron star accreting from the wind of a massive companion and producing a powerful X-ray source. The characteristics of the stellar wind together with the complex interactions between the compact object and the donor star determine the observed X-ray output from all these systems. Consequently, the use of SgXBs for studies of massive stars is only possible when the physics of the stellar winds, the compact objects, and accretion mechanisms are combined together and confronted with observations.This detailed review summarises the current knowledge on the theory and observations of winds from massive stars, as well as on observations and accretion processes in wind-fed high mass X-ray binaries. The aim is to combine in the near future all available theoretical diagnostics and observational measurements to achieve a unified picture of massive star winds in isolated objects and in binary systems.     

A Composition Analysis Tool for the Solar Wind Around Pluto (SWAP) Instrument on New Horizons

We describe the response of the Solar Wind Around Pluto (SWAP) instrument (McComas et al. in Space Sci. Rev. 140:261, 2008) to 1–40 amu ions in order to assess whether it can be used to determine plasma composition. Our goal is to enhance the scientific return on the SWAP plasma measurements obtained during the New Horizons traversal down Jupiter’s magnetotail in 2007. We present calibration data for the SWAP flight instrument and another largely flight-like SWAP sensor, dubbed “SWAP-II”. SWAP’s mass-dependent response was characterized by analyzing the count ratios from its two channel electron multipliers (CEMs). We observe significant differences in the instrument response between light (mass ≤ He) and heavy (mass > He) ions, especially for energies below ～4 keV. We attribute these differences to the mass-dependent electron emission yield from SWAP’s ultra-thin (～1 μg/cm²) carbon foil. Using these results, we develop a plasma composition analysis technique to statistically distinguish between light and heavy plasma ions measured by the instrument.     

Evidence for the lowering of the centroid of daytime thermospheric O(D) 630.0 nm emission over the magnetic equator: First results

It is shown in this paper for the first time that the intensity of the daytime thermospheric O(¹D) 630.0 nm airglow as measured by the ground-based dayglow photometer over Trivandrum (8.5°N; 77°E; dip lat. 0.5°N), a geomagnetic dip equatorial station, exhibit a direct correlation with the electron density at 180 km. This altitude is about ∼40 km lower than the believed centroid of the O(¹D) 630.0 nm dayglow emission i.e. 220 km. This observation is contrary to the understanding of the behavior of O(¹D) 630.0 nm dayglow over equatorial/low latitudes. Over these latitudes, the variations of the measured intensity of O(¹D) 630.0 nm dayglow are known to be associated with the changes in the electron density at altitudes around 220 km, the centroid of this emission. In this context, the present results indicating the lowering of the peak altitude of O(¹D) 630.0 nm emission from ∼220 to ∼180 km over the dip equator is new. Recent results on solar XUV flux indicate that this could be an important parameter that controls the O(¹D) 630.0 nm dayglow excitation rates through modulations in the neutral and ionic composition in lower thermosphere-ionosphere region. However, the lowering of the centroid of O(¹D) 630.0 nm emission, as shown in this study, has been ascribed primarily to the fountain effect associated with the equatorial ionization anomaly.     

Auroral Substorms: Search for Processes Causing the Expansion Phase in Terms of the Electric Current Approach

Auroral substorms are mostly manifestations of dissipative processes of electromagnetic energy. Thus, we consider a sequence of processes consisting of the power supply (dynamo), transmission (currents/circuits) and dissipations (auroral substorms-the end product), namely the electric current line approach. This work confirms quantitatively that after accumulating magnetic energy during the growth phase, the magnetosphere unloads the stored magnetic energy impulsively in order to stabilize itself. This work is based on our result that substorms are caused by two current systems, the directly driven (DD) current system and the unloading system (UL). The most crucial finding in this work is the identification of the UL (unloading) current system which is responsible for the expansion phase. A very tentative sequence of the processes leading to the expansion phase (the generation of the UL current system) is suggested for future discussions.

1. (1)
   The solar wind-magnetosphere dynamo enhances significantly the plasma sheet current when its power is increased above \(10^{18}~\mbox{erg}/\mbox{s}\) (\(10^{11}\) w).
    
2. (2)
   The magnetosphere accumulates magnetic energy during the growth phase, because the ionosphere cannot dissipate the increasing power because of a low conductivity. As a result, the magnetosphere is inflated, accumulating magnetic energy.
    
3. (3)
   When the power reaches \(3\mbox{--}5\times 10^{18}~\mbox{erg}/\mbox{s}\) (\(3\mbox{--}5\times 10^{11}\) w) for about one hour and the stored magnetic energy reaches \(3\mbox{--}5\times10^{22}\) ergs (\(10^{15}\) J), the magnetosphere begins to develop perturbations caused by current instabilities (the current density \({\approx}3\times 10^{-12}~\mbox{A}/\mbox{cm}^{2}\) and the total current \({\approx}10^{6}~\mbox{A}\) at 6 Re). As a result, the plasma sheet current is reduced.
    
4. (4)
   The magnetosphere is thus deflated. The current reduction causes \(\partial B/\partial t > 0\) in the main body of the magnetosphere, producing an earthward electric field. As it is transmitted to the ionosphere, it becomes equatorward-directed electric field which drives both Pedersen and Hall currents and thus generates the UL current system.
    
5. (5)
   A significant part of the magnetic energy is accumulated in the main body of the magnetosphere (the inner plasma sheet) between 4 Re and 10 Re, because the power (Poynting flux \([ \boldsymbol{E} \times \boldsymbol{B} ])\) is mainly directed toward this region which can hold the substorm energy.
    
6. (6)
   The substorm intensity depends on the location of the energy accumulation (between 4 Re and 10 Re), the closer the location to the earth, the more intense substorms becomes, because the capacity of holding the energy is higher at closer distances. The convective flow toward the earth brings both the ring current and the plasma sheet current closer when the dynamo power becomes higher.
    

This proposed sequence is not necessarily new. Individual processes involved have been considered by many, but the electric current approach can bring them together systematically and provide some new quantitative insights.