Correlation of IRTAM and FPMU data confirming the application of IRTAM to support ISS Program safety

A “Real-Time” plasma hazard assessment process was developed to support International Space Station (ISS) Program real-time decision-making providing solar array constraint relief information for Extravehicular Activities (EVAs) planning and operations. This process incorporates real-time ionospheric conditions, ISS solar arrays’ orientation, ISS flight attitude, and where the EVA will be performed on the ISS. This assessment requires real-time data that is presently provided by the Floating Potential Measurement Unit (FPMU) which measures the ISS floating potential (FP), along with ionospheric electron number density (Ne) and electron temperature (Te), in order to determine the present ISS environment. Once the present environment conditions are correlated with International Reference Ionosphere (IRI) values, IRI is used to forecast what the environment could become in the event of a severe geomagnetic storm. If the FPMU should fail, the Space Environments team needs another source of data which is utilized to support a short-term forecast for EVAs. The IRI Real-Time Assimilative Mapping (IRTAM) model is an ionospheric model that uses real-time measurements from a large network of digisondes to produce foF2 and hmF2 global maps in 15?min cadence. The Boeing Space Environments team has used the IRI coefficients produced in IRTAM to calculate the Ne along the ISS orbital track. The results of the IRTAM model have been compared to FPMU measurements and show excellent agreement. IRTAM has been identified as the FPMU back-up system that will be used to support the ISS Program if the FPMU should fail.     

Mapping mean lake surface from satellite altimetry and GPS kinematic surveys

Lake water height is a key variable in water cycle and climate change studies, which is achievable using satellite altimetry constellation. A method based on data processing of altimetry from several satellites has been developed to interpolate mean lake surface (MLS) over a set of 22 big lakes distributed on the Earth. It has been applied on nadir radar altimeters in Low Resolution Mode (LRM: Jason-3, Saral/AltiKa, CryoSat-2) in Synthetic Aperture Radar (SAR) mode (Sentinel-3A), and in SAR interferometric (SARin) mode (CryoSat-2), and on laser altimetry (ICESat). Validation of the method has been performed using a set of kinematic GPS height profiles from 18 field campaigns over the lake Issykkul, by comparison of altimetry’s height at crossover points for the other lakes and using the laser altimetry on ICESat-2 mission. The precision reached ranges from 3 to 7 cm RMS (Root Mean Square) depending on the lakes. Currently, lake water level inferred from satellite altimetry is provided with respect to an ellipsoid. Ellipsoidal heights are converted into orthométric heights using geoid models interpolated along the satellite tracks. These global geoid models were inferred from geodetic satellite missions coupled with absolute and regional anomaly gravity data sets spread over the Earth. However, the spatial resolution of the current geoid models does not allow capturing short wavelength undulations that may reach decimeters in mountaineering regions or for rift lakes (Baikal, Issykkul, Malawi, Tanganika). We interpolate in this work the geoid height anomalies with three recent geoid models, the EGM2008, XGM2016 and EIGEN-6C4d, and compare them with the Mean Surface of 22 lakes calculated using satellite altimetry. Assuming that MLS mimics the local undulations of the geoid, our study shows that over a large set of lakes (in East Africa, Andean mountain and Central Asia), short wavelength undulations of the geoid in poorly sampled areas can be derived using satellite altimetry. The models used in this study present very similar geographical patterns when compared to MLS. The precision of the models largely depends on the location of the lakes and is about 18 cm, in average over the Earth. MLS can serve as a validation dataset for any future geoid model. It will also be useful for validation of the future mission SWOT (Surface Water and Ocean Topography) which will measure and map water heights over the lakes with a high horizontal resolution of 250 by 250 m.     

Impact of rain on wave retrieval from Sentinel-1 synthetic aperture radar images in tropical cyclones

The main objective of our work was to investigate the impact of rain on wave observations from C-band (~5.3 GHz) synthetic aperture radar (SAR) in tropical cyclones. In this study, 10 Sentinel-1 SAR images were available from the Satellite Hurricane Observation Campaign, which were taken under cyclonic conditions during the 2016 hurricane season. The third-generation wave model, known as Simulating WAves Nearshore (SWAN) (version 41.31), was used to simulate the wave fields corresponding to these Sentinel-1 SAR images. In addition, rainfall data from the Tropical Rainfall Measuring Mission satellite passing over the spatial coverage of the Sentinel-1 SAR images were collected. The simulated results were validated against significant wave heights (SWHs) from the Jason-2 altimeter and European Centre for Medium-Range Weather Forecasts data, revealing a root mean square error (RMSE) of ~0.5 m with a 0.25 scatter index. Winds retrieved from the VH-polarized Sentinel-1 SAR images using the Sentinel-1 Extra Wide-swath Mode Wind Speed Retrieval Model after Noise Removal were taken as prior information for wave retrieval. It was discovered that rain did indeed affect the SAR wave retrieval, as evidenced by the 3.21-m RMSE of SWHs between the SAR images and the SWAN model, which was obtained for the ~1000 match-ups with raindrops. The raindrops dampened the wave retrieval when the rain rate was < ~5 mm/hr; however, they enhanced wave retrieval for higher rain rates. It was also found that the portion of the rain-induced ring wave with a wave number > 0.05 rad/m (~125 m wavelength) was clearly observed in the SAR-derived wave spectra.     

Comparison of GPS derived TEC with the TEC predicted by IRI 2012 model in the southern Equatorial Ionization Anomaly crest within the Eastern Africa region

We have compared the TEC obtained from the IRI-2012 model with the GPS derived TEC data recorded within southern crest of the EIA in the Eastern Africa region using the monthly means of the 5 international quiet days for equinoxes and solstices months for the period of 2012 – 2013. GPS-derived TEC data have been obtained from the Africa array and IGS network of ground based dual-frequency GPS receivers from four stations (Kigali (1.95°S, 30.09°E; Geom. Lat. 11.63°S), Malindi (2.99°S, 40.19°E; Geom. Lat. 12.42°S), Mbarara (0.60°S, 30.74°E; Geom. Lat. 10.22°S) and Nairobi (1.22°S, 36.89°E; Geom. Lat. 10.69°S)) located within the EIA crest in this region. All the three options for topside Ne of IRI-2012 model and ABT-2009 for bottomside thickness have been used to compute the IRI TEC. Also URSI coefficients were considered in this study. These results are compared with the TEC estimated from GPS measurements. Correlation Coefficients between the two sets of data, the Root-Mean Square Errors (RMSE) of the IRI-TEC from the GPS-TEC, and the percentage RMSE of the IRI-TEC from the GPS-TEC have been computed. Our general results show that IRI-2012 model with all three options overestimates the GPS-TEC for all seasons and at all stations, and IRI-2001 overestimates GPS-TEC more compared with other options. IRI-Neq and IRI-01-corr are closely matching in most of the time. The observation also shows that, GPS TEC are underestimated by TEC from IRI model during noon hours, especially during equinoctial months. Further, GPS-TEC values and IRI-TEC values using all the three topside Ne options show very good correlation (above 0.8). On the other hand, the TEC using IRI-Neq and IRI-01- corr had smaller deviations from the GPS-TEC compared to the IRI-2001.     

Identification of seismic activity sources on the subsatellite track by ionospheric plasma disturbances detected with the Sich-2 onboard probes

Using a cylindrical Langmuir probe and the authors’ proprietary two-channel pressure transducer, ionospheric plasma parameter distributions along the orbit of the Sich-2 satellite (Ukraine, 2011–2012) were measured. This paper is concerned with identifying the space–time location of ionospheric plasma disturbance sources, including the epicenters of actual earthquakes (before or during the satellite flyover) and incipient earthquakes on the subsatellite track, from the measured distributions of the electron density and temperature and the neutral particle temperature along the satellite orbit. To do this, the measured ionospheric plasma parameter distributions are connected to the coordinates on the subsatellite track.It is shown that local disturbances in the electron density and temperature and neutral particle temperature distributions in the satellite orbit in the ionosphere may serve as indicators of seismic activity on the subsatellite track. The epicenters of incipient earthquakes may be set off from other plasma parameter disturbance sources associated with seismic activity using information provided by special monitoring and survey centers that monitor the current seismic situation.     

Regolith X-Ray Imaging Spectrometer (REXIS) Aboard the OSIRIS-REx Asteroid Sample Return Mission

The Regolith X-ray Imaging Spectrometer (REXIS) is the student collaboration experiment proposed and built by an MIT-Harvard team, launched aboard NASA’s OSIRIS-REx asteroid sample return mission. REXIS complements the scientific investigations of other OSIRIS-REx instruments by determining the relative abundances of key elements present on the asteroid’s surface by measuring the X-ray fluorescence spectrum (stimulated by the natural solar X-ray flux) over the range of energies 0.5 to 7 keV. REXIS consists of two components: a main imaging spectrometer with a coded aperture mask and a separate solar X-ray monitor to account for the Sun’s variability. In addition to element abundance ratios (relative to Si) pinpointing the asteroid’s most likely meteorite association, REXIS also maps elemental abundance variability across the asteroid’s surface using the asteroid’s rotation as well as the spacecraft’s orbital motion. Image reconstruction at the highest resolution is facilitated by the coded aperture mask. Through this operation, REXIS will be the first application of X-ray coded aperture imaging to planetary surface mapping, making this student-built instrument a pathfinder toward future planetary exploration. To date, 60 students at the undergraduate and graduate levels have been involved with the REXIS project, with the hands-on experience translating to a dozen Master’s and Ph.D. theses and other student publications.     

ExoMars 2016 Schiaparelli Module Trajectory and Atmospheric Profiles Reconstruction

On 19^th October 2016 Schiaparelli module of the ExoMars 2016 mission flew through the Mars atmosphere. After successful entry and descent under parachute, the module failed the last part of the descent and crashed on the Mars surface. Nevertheless the data transmitted in real-time by Schiaparelli during the entry and descent, together with the entry state vector as initial condition, have been used to reconstruct both the trajectory and the profiles of atmospheric density, pressure and temperature along the traversed path.The available data-set is only a small sub-set of the whole data acquired by Schiaparelli, with a limited data rate (8 kbps) and a large gap during the entry because of the plasma blackout on the communications.This paper presents the work done by the AMELIA (Atmospheric Mars Entry and Landing Investigations and Analysis) team in the exploitation of the available inertial and radar data. First a reference trajectory is derived by direct integration of the inertial measurements and a strategy to overcome the entry data gap is proposed. First-order covariance analysis is used to estimate the uncertainties on all the derived parameters. Then a refined trajectory is computed incorporating the measurements provided by the on-board radar altimeter.The derived trajectory is consistent with the events reported in the telemetry and also with the impact point identified on the high-resolution images of the landing site.Finally, atmospheric profiles are computed tacking into account the aerodynamic properties of the module. Derived profiles result in good agreement with both atmospheric models and available remote sensing observations.     

Detection of Molecular Scattering Field from a Polarization Analysis of the Sky Background during Transitive Twilight and Temperature Measurements near the Stratopause

The simplest version of the method of detecting the single molecular scattering field based on the polarization measurements of the twilight sky background by all-sky cameras has been considered. The method can be used during transitive twilight (with solar zenith angles of 94°–98°), when effective single scattering occurs in the upper stratosphere and lower mesosphere. The long-term measurements conducted using this method in the Moscow region and Apatity make it possible to determine the temperature of these atmospheric layers and estimate the contribution and properties of multiple scattering during the transitive twilight.