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.     

A novel framework for liquid propellant engine’s cooling system design by sensitivity analysis based on RSM and multi-objective optimization using PSO

One of the challenges of combustion chamber and nozzle design in a Liquid Propellant Engine (LPE) is to predict the behavior and performance of the cooling system. Therefore, while designing, the optimization of the cooling system is always of great importance. This paper presents the multi-objective optimization of the LPE’s cooling system. To this end, a novel framework has been developed, resulting from the application of the Response Surface Method (RSM) and the correlation coefficients matrix, sensitivity analysis and the The Particle Swarm Optimization (PSO). based on this method, the input variables, constraints, objective functions, and their surfaces were identified. In terms of multi-optimization algorithms, RSM and PSO are utilized to get global optimum. In conclusion, the methodology capability is to optimize the LPE’s cooling system, 6 percentage increase in total heat transfer and 7 bar decrease cooling system pressure loss, which resulted in a 1.2-seconds increase in the specific impulse of the engine.     

Deorbiter CubeSat mission design

This paper presents the mission design for a CubeSat-based active debris removal approach intended for transferring sizable debris objects from low-Earth orbit to a deorbit altitude of 100 km. The mission consists of a mothership spacecraft that carries and deploys several debris-removing nanosatellites, called Deorbiter CubeSats. Each Deorbiter is designed based on the utilization of an eight-unit CubeSat form factor and commercially-available components with significant flight heritage. The mothership spacecraft delivers Deorbiter CubeSats to the vicinity of a predetermined target debris, through performing a long-range rendezvous maneuver. Through a formation flying maneuver, the mothership then performs in-situ measurements of debris shape and orbital state. Upon release from the mothership, each Deorbiter CubeSat proceeds to performing a rendezvous and attachment maneuver with a debris object. Once attached to the debris, the CubeSat performs a detumbling maneuver, by which the residual angular momentum of the CubeSat-debris system is dumped using Deorbiter’s onboard reaction wheels. After stabilizing the attitude motion of the combined Deorbiter-debris system, the CubeSat proceeds to performing a deorbiting maneuver, i.e., reducing system’s altitude so much so that the bodies disintegrate and burn up due to atmospheric drag, typically at around 100 km above the Earth surface. The attitude and orbital maneuvers that are planned for the mission are described, both for the mothership and Deorbiter CubeSat. The performance of each spacecraft during their operations is investigated, using the actual performance specifications of the onboard components. The viability of the proposed debris removal approach is discussed in light of the results.     

A novel satellite image encryption algorithm based on hyperchaotic systems and Josephus problem

With the rapid growth of the number of Earth observation satellite (EOS) supporting critical applications, it is required to improve the security techniques to protect the sensitive data and images during the transmission between the satellites and the ground stations. This paper introduces a new satellite image encryption algorithm based on the Linear Feedback Shift Register (LFSR) generator, SHA 512 hash function, hyperchaotic systems, and Josephus problem. LFSR generates a matrix that is used to construct the 512-bits value of the hash function. These bits are used to set the initial values and parameters of the proposed encryption algorithm. Firstly, the six dimensions (6-D) hyperchaotic system is divided into three parts, where every two equations are considered as one part. Secondly, the 1-D hyperchaotic logistic-tent system is considered as the controller to select one part. The selected part is used to generate a matrix that is XORed with the original image. Thirdly, the scrambling operation by Josephus sequences is applied to the output of the previous step by scrambling the rows and the columns according to the selected part to produce the pre-encrypted image. Finally, if the number of iterations is less than the required number which is considered as a parameter of the secret key, the previous operations will be repeated in the pre-encrypted image; otherwise, the pre-encrypted image is considered as the final cipher image. Experimental and analyses results show that the proposed algorithm has good performance in terms of high level of security, large enough key-space, tolerance to Single Event Upsets (SEU) as well as low time complexity.     

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.     

Spacecraft survivability in the natural debris environment near the stable Earth-Moon Lagrange points

The theoretical analysis of the motion of natural space debris near the stable Earth-Moon Lagrange Points, $L_4$ and $L_5$, is presented with a focus on the potential debris risks to spacecraft operating near these points. Specifically, the research formulates a debris propagation model using four-body dynamics, then applies candidate probabilistic survivability models to a notional spacecraft operating at the $L_4$ and $L_5$ Lagrange points to quantify the collision risks to the spacecraft from natural debris particles. Of the survivability models implemented, the natural debris collision risks to spacecraft survivability are found to be incredibly low, but mitigation strategies to reduce the risk further are identified in this study. Overall, research into stable Lagrange point natural debris propagation improves understanding of the collision risks posed by the naturally occurring Kordylewski clouds and enhances operational planning for Lagrange point space missions.     

Equatorial anomaly according to the Interkosmos-19 data and IRI model: A comparison

The comparison of the IRI model with the foF2 distribution in the equatorial anomaly region obtained by topside sounding onboard the Interkosmos-19 satellite has been carried out. The global distribution of foF2 in terms of LT-maps was constructed by averaging Intercosmos-19 data for summer, winter, and equinox. These maps, in fact, represent an empirical model of the equatorial anomaly for high solar activity F10.7 ~ 200. The comparison is carried out for the latitudinal foF2 profiles in the characteristic longitudinal sectors of 30, 90, 210, 270, and 330°, as well as for the longitudinal variations in foF2 over the equator. The largest difference between the models (up to 60%) for any season was found in the Pacific longitudinal sector of 210°, where there are a few ground-based sounding stations. Considerable discrepancies, however, are sometimes observed in the longitudinal sectors, where there are many ground-based stations, for example, in the European or Indian sector. The discrepancies reach their maximum at 00 LT, since a decay of the equatorial anomaly begins before midnight in the IRI model and after midnight according to the Interkosmos-19 data. The discrepancies are also large in the morning at 06 LT, since in the IRI model, the foF2 growth begins long before sunrise. In the longitudinal variations in foF2 over the equator at noon, according to the satellite data, four harmonics are distinguished in the June solstice and at the equinox, and three harmonics in the December solstice, while in the IRI model only two and one harmonics respectively are revealed. In diurnal variations in foF2 and, accordingly, in the equatorial anomaly intensity, the IRI model does not adequately reproduce even the main, evening extremum.     

Atmospheric reentry hemisphere prediction for prograde orbits using logical disjunction

A unique logic-based algorithm for atmospheric reentry hemisphere prediction is presented for spacecraft in low-eccentricity, prograde low Earth orbits at altitudes of 300 km and lower. Using two-line element (TLE) data for initial orbit conditions, coupled with coarse estimates for spacecraft aerodynamic characteristics, the algorithm relies on logical disjunction operations based on a dual analysis of histogram and two-weighted Gaussian probability density function (PDF) fits of predicted reentry latitude data. The algorithm requires the execution of a series of parametric simulations to determine the reentry hemisphere for variations in spacecraft aerodynamic coefficients and drag reference area. When implemented, the algorithm yields accurate hemisphere predictions on average 15 days from reentry as demonstrated by historical reentry cases from 1979 to 2018. All reentry cases were selected to demonstrate the algorithm’s ability to deliver accurate reentry hemisphere predictions for spacecraft with varying physical size and mass, and reentering during different periods of solar cycle activity.     

Observations of a quasi-90-day oscillation in the MLT winds and tides over an equatorial station using meteor radar winds

An evident signature of a least studied quasi-90-day oscillation is found in the winds and tides in the MLT from an equatorial station, São João do Cariri (7.4°S, 36.5°W). The oscillation is found to appear mainly in certain intervals with small but appreciable seasonal (fourth harmonic of annual oscillation) contribution. The maximum amplitude of the oscillation is found to be around 10 m/s in the zonal wind. The enhancement peak of the oscillation exhibits downward movement indicating a plausible role of upward moving waves/tides in carrying its imprint from below to the MLT. Similar oscillation feature in the tropospheric zonal wind and ozone may imply its lower atmospheric origin as a component of the intraseasonal oscillation (ISO) that moves upward by modulating the tides. Subsequently, the propagating tides (mainly semidiurnal) are enhanced by the ozone in the stratosphere through absorption of solar UV radiation and finally manifest the oscillation in the MLT. Consistency of the present findings with the past investigations are observed in some aspects of the oscillation, whereas existing mismatches in others are believed to be due to geophysical variability depending on space and time among various locations on the globe.