The Student Dust Counter on the New Horizons Mission

The Student Dust Counter (SDC) experiment of the New Horizons Mission is an impact dust detector to map the spatial and size distribution of dust along the trajectory of the spacecraft across the solar system. The sensors are thin, permanently polarized polyvinylidene fluoride (PVDF) plastic films that generate an electrical signal when dust particles penetrate their surface. SDC is capable of detecting particles with masses m>10^?12 g, and it has a total sensitive surface area of about 0.1 m², pointing most of the time close to the ram direction of the spacecraft. SDC is part of the Education and Public Outreach (EPO) effort of this mission. The instrument was designed, built, tested, integrated, and now is operated by students.     

In-flight testing of the injection of the LISA Pathfinder test mass into a geodesic

LISA Pathfinder is a technology demonstrator space mission, aimed at testing key technologies for detecting gravitational waves in space. The mission is the precursor of LISA, the first space gravitational waves observatory, whose launch is scheduled for 2034. The LISA Pathfinder scientific payload includes two gravitational reference sensors (GRSs), each one containing a test mass (TM), which is the sensing body of the experiment. A mission critical task is to set each TM into a pure geodesic motion, i.e. guaranteeing an extremely low acceleration noise in the sub-Hertz frequency bandwidth. The grabbing positioning and release mechanism (GPRM), responsible for the injection of the TM into a geodesic trajectory, was widely tested on ground, with the limitations imposed by the 1-g environment. The experiments showed that the mechanism, working in its nominal conditions, is capable of releasing the TM into free-fall fulfilling the very strict constraint imposed on the TM residual velocity, in order to allow its capture on behalf of the electrostatic actuation.However, the first in-flight releases produced unexpected residual velocity components, for both the TMs. Moreover, all the residual velocity components were greater than maximum value set by the requirements. The main suspect is that unexpected contacts took place between the TM and the surroundings bodies. As a consequence, ad hoc manual release procedures had to be adopted for the few following injections performed during the nominal mission. These procedures still resulted in non compliant TM states which were captured only after impacts. However, such procedures seem not practicable for LISA, both for the limited repeatability of the system and for the unmanageable time lag of the telemetry/telecommand signals (about 4400 s). For this reason, at the end of the mission, the GPRM was deeply tested in-flight, performing a large number of releases, according to different strategies. The tests were carried out in order to understand the unexpected dynamics and limit its effects on the final injection. Some risk mitigation maneuvers have been tested aimed at minimizing the vibration of the system at the release and improving the alignment between the mechanism and the TM. However, no overall optimal release strategy to be implemented in LISA could be found, because the two GPRMs behaved differently.     

How the astronomical aspects of climate science were settled? On the Milankovitch and Bacsák anniversaries,with lessons for today

It was 100 years ago (on August 7, 1920), that the comprehensive mathematical foundations of climate change research, written by a Serbian researcher, Milutin Milankovitch, were published. A later interpreter and developer of his results, Georg (in Hungarian: György) Bacsák (Pozsony/Pressburg/Bratislava, June 5, 1870 - Fonyód, March 4, 1970) was born 150 years ago and died at the age of one hundred, half a century ago. In this commemorative paper we look back to special circumstances in revealing the secrets of ice ages that had puzzled scientists for at least several centuries. Recently, after 100 years, the Milankovitch theory, including related short-term forcings (ranging from interannual, multidecadal to millennial timescales) has not only been confirmed, but its climate forcing mechanism has also been identified and proposed. Owing to the uniqueness of the problem, the science of the orbital forcing of climate change can be proclaimed to be essentially settled.     

Encounter trajectories for deep space mission ASTER to the triple near Earth asteroid 2001-SN263. The laser altimeter (ALR) point of view

In cooperation with Russia, the Brazilian deep space mission ASTER plans to send a small spacecraft to investigate the triple asteroid 2001-SN263. The nearest launch opportunities for this project include June 2022 and June 2025. One main exploration campaign is being planned with focus on the largest asteroid (Alpha). Among the instruments under development, a laser altimeter (named ALR) was preliminarily designed and presented in 2010–2011. Many studies to define mission and instruments requirements were performed aiming at the characterization of important issues for the successful realization of the mission. Among them, the identification of a suitable trajectory that could be followed by the ASTER spacecraft in the encounter phase, when the main campaign will take place. This paper describes the effort undertaken with focus on the laser altimeter operation. Possible encounter trajectories were modelled and simulated to identify suitable approach parameters and conditions allowing the accomplishment of the intended investigation. The simulation also involves the instrument operation, considering approach geometry, attitude, relative motion, time/date, and the dynamics of the main asteroid. From the laser altimeter point of view, keeping in mind the desired coverage results (50% minimum surface coverage of asteroid Alpha, complying with horizontal and vertical resolution requirements), results point out crucial features for the encounter trajectory, like the need for a small inclination (10^-6 degrees; with respect to the asteroid's orbit), the most favourable spacecraft positioning (between the Sun and the asteroid) and pointing condition (back to the Sun), the minimum amount of achievable surface coverage (58%, focused on central areas), and the most proper time to conduct the main campaign (January 2025). Concerning the instrument, results offer refined values for divergence angle (500 to 650 μrad, half-cone), pulse repetition frequencies (from 1/20 to 1 Hz), and consequent data generation rates. A simulation tool that can use any 3D generated trajectories as input data was created for the analyses presented here. Although created for the ALR in this mission, this simple analysis tool can be adapted to other instruments in this or other missions.     

Dispersive shock waves in partially ionised plasmas

Compressional waves propagating in the partially ionised solar lower atmospheric plasmas can easily steepen into nonlinear waves, including shocks. Here we investigate the effect of weak dispersion generated by Hall currents perpendicular to the ambient magnetic field on the characteristics of shock waves. Our study will also focus on the interplay between weak dispersion and partial ionisation of the plasma. Using a multiple scale technique we derive the governing equation in the form of a Korteweg-de Vries-Burgers equation. The effect of weak dispersion on shock waves is obtained using a perturbation technique. The secular behaviour of second order terms is addressed with the help of a renormalization technique. Our results show that dispersion modifies the characteristics of shock waves and this change is dependent also on the ionisation degree of the plasma. Dispersion can create short lived oscillations in the shocked plasma. The shock fronts become wider with the increase in the number of neutrals in the plasma.     

Horizontal and vertical thermospheric cross-wind from GOCE linear and angular accelerations

Thermospheric wind measurements obtained from linear non-gravitational accelerations of the Gravity field and steady-state Ocean Circulation Explorer (GOCE) satellite show discrepancies when compared to ground-based measurements. In this paper the cross-wind is derived from both the linear and the angular accelerations using a newly developed iterative algorithm. The two resulting data sets are compared to test the validity of wind derived from angular accelerations and quantify the uncertainty in accelerometer-derived wind data. In general the difference is found to be less than 50?m/s vertically after high-pass filtering, and 100?m/s horizontally. A sensitivity analysis reveals that continuous thrusting is a major source of uncertainty in the torque-derived wind, as are the magnetic properties of the satellite. The energy accommodation coefficient is identified as a particularly promising parameter for improving the consistency of thermospheric cross-wind data sets in the future. The algorithm may be applied to obtain density and cross-wind from other satellite missions that lack accelerometer data, provided the attitude and orbit are known with sufficient accuracy.     

Spatio-temporal variability of the wet component of the troposphere – Application to satellite altimetry

Due to the presence of water vapour and cloud liquid water in the atmosphere, the wet component of the troposphere is responsible for a delay in the propagation of the altimeter signals, the Wet Path Delay (WPD). The high space–time variability of the water vapour distribution makes the modelling of WPD difficult, its effect still being one of the main error sources in satellite altimetry applications, e.g. in the estimation of Mean Sea Level (MSL). The understanding and the quantification of the WPD variability on various spatial and temporal scales are the main purposes of this study, in view to improve the MSL error budget. The dominant timescales of WPD variability and its correlation with Sea Level Anomaly (SLA) are examined. In these analyses, the atmospheric reanalysis ERA-Interim model from the European Centre for Medium-Range Weather Forecasts (ECMWF) is used to derive a global dataset of daily grids of WPD, spanning a 28-year period from January 1988 to December 2015. The Seasonal-Trend decomposition procedure based on Loess (STL) is used to extract precise WPD annual and interannual signals. Linear trends have been derived from the interannual time series and the contribution of each STL component was mapped globally, allowing the understanding of the WPD variability in spatial terms. The correlation between SLA and WPD is mapped and decomposed into seasons using monthly mean grids, for a period of 21-years, from January 1993 to December 2013.Aiming at inspecting the sensitivity of the results to the used data set, the WPD temporal analysis is extended to the data set provided by the Special Sensor Microwave Imager (SSM/I) and SSM/I Sounder (SSM/IS) Sensors. The WPD from SSM/I(S) is compared against those from the ERA-Interim and from the National Centers for Environmental Prediction (NCEP).Results show that climate phenomena, especially the El Niño Southern Oscillation (ENSO) are the cause for this high variability, since they affect the water vapour and temperature. The observed trends from ERA-Interim, computed globally and over ocean regions only, allow concluding that WPD is increasing with time by approximately 0.1?mm per year, and the maximum trends are observed for the Pacific North and Indian Oceans. High correlation between WPD and SLA is found over the western tropical Pacific.The comparison between WPD from SSM/I(S) and from ERA-Interim and NCEP, allows concluding that the trends computed using only the SSM/I(S) measurement points are substantially larger.     

Pre-earthquake ionospheric anomalies before three major earthquakes by GPS-TEC and GIM-TEC data during 2015–2017

The dual-frequency satellite-based measurements from Global Positioning System (GPS) may provide feasible ways of studying and potentially detecting of earthquake (EQ) related anomalies in the ionosphere. In this paper, GPS based Total Electron Content (TEC) data are studied for three major M?>?7.0 EQs in Nepal and Iran-Iraq border during 2015–2017 by implementing statistical procedures on temporal and spatial scale. Previous studies presented different time interval of pre-seismic ionospheric anomalies, however, this study showed that EQs ionospheric precursors may occur within 10?days. Furthermore, the ionospheric anomalies on the suspected day occurred during UT?=?08:00–12:00?h before the main shock. The Global Ionospheric Map TEC (GIM-TEC) data retrieved over the epicenter of M7.8 (Nepal EQ) showed a significant increase of 6 TECU on April 24, 2015 (one day before the main shock), which is recorded by the ground GPS station data of Islamabad (station lies within the EQ preparation zone). Furthermore, the spatial GIM-TEC result imply significant anomalies over the epicenter during the time interval between UT?=?08:00–12:00?h (LT?=?13:00–17:00). For M7.3 (Nepal EQ), the TEC anomalies were detected on May 10, 2015 (2?days before the event) in the temporal data. The spatial TEC data imply the huge clouds over the epicenter at about UT?=?08:00–12:00?h on May 10, 2015, that may be associated with this EQ in the quiet geomagnetic storm conditions. Similarly, temporal and spatial TEC showed anomaly on November 3, 2017, during UT?=?08:00–12:00 (9?days before the Iran-Iraq border EQ) after implementing the statistical method on it. Conversely, there exists a short-term but low magnitude TEC anomaly synchronized with a geomagnetic storm on November 7–8, 2017 (4 to 5?days prior to M7.3 Iran-Iraq border EQ). The diurnal and hourly GIM-TEC and VTEC data also imply the execution of ionospheric anomalies within 10?days prior to all events. All these positive anomalies in TEC may be due to the existence of a huge energy from the epicenter during the EQ preparation period.