Upgrading CCIR's foF2 maps using available ionosondes and genetic algorithms

We have developed a new approach towards a new database of the ionospheric parameter $f_{o}F2$. This parameter, being the frequency of the maximum of the ionospheric electronic density profile and its main modeller, is of great interest not only in atmospheric studies but also in the realm of radio propagation. The current databases, generated by CCIR (Committee Consultative for Ionospheric Radiowave propagation) and URSI (International Union of Radio Science), and used by the IRI (International Reference Ionosphere) model, are based on Fourier expansions and have been built in the 60s from the available ionosondes at that time. The main goal of this work is to upgrade the databases by using new available ionosonde data. To this end we used the IRI diurnal/spherical expansions to represent the $f_{o}F2$ variability, and computed its coefficients by means of a genetic algorithm (GA). In order to test the performance of the proposed methodology, we applied it to the South American region with data obtained by RAPEAS (Red Argentina para el Estudio de la Atmósfera Superior, i.e. Argentine Network for the Study of the Upper Atmosphere) during the years 1958–2009. The new GA coefficients provide a global better fit of the IRI model to the observed $f_{o}F2$ than the CCIR coefficients. Since the same formulae and the same number of coefficients were used, the overall integrity of IRI’s typical ionospheric feature representation was preserved. The best improvements with respect to CCIR are obtained at low solar activities, at large (in absolute value) modip latitudes, and at night-time. The new method is flexible in the sense that can be applied either globally or regionally. It is also very easy to recompute the coefficients when new data is available. The computation of a third set of coefficients corresponding to days of medium solar activity in order to avoid the interpolation between low and high activities is suggested. The same procedure as for $f_{o}F2$ can be perfomed to obtain the ionospheric parameter M(3000)F2.     

Simulating planetary wave propagation to the upper atmosphere during stratospheric warming events at different mountain wave scenarios

Parameterization schemes of atmospheric normal modes (NMs) and orographic gravity waves (OGWs) have been implemented into the mechanistic Middle and Upper Atmosphere Model (MUAM) simulating atmospheric general circulation. Based on the 12-members ensemble of runs with the MUAM, a composite of the stratospheric warming (SW) has been constructed using the UK Met Office data as the lower boundary conditions. The simulation results show that OGW amplitudes increase at altitudes above 30 km in the Northern Hemisphere after the SW event. At altitudes of about 50 km, OGWs have largest amplitudes over North American and European mountain systems before and during the composite SW, and over Himalayas after the SW. Simulations demonstrate substantial (up to 50–70%) variations of amplitudes of stationary planetary waves (PWs) during and after the SW in the mesosphere-lower thermosphere of the Northern Hemisphere. Westward travelling NMs have amplitude maxima not only in the Northern, but also in the Southern Hemisphere, where these modes have waveguides in the middle and upper atmosphere. Simulated variations of PW and NM amplitudes correspond to changes in the mean zonal wind, EP-fluxes and wave refractive index at different phases of the composite SW events. Inclusion of the parameterization of OGW effects leads to decreases in amplitudes (up to 15%) of almost all SPWs before and after the SW event and their increase (up to 40–60%) after the SW in the stratosphere and mesosphere at middle and high northern latitudes. It is suggested that observed changes in NM amplitudes in the Southern Hemisphere during SW could be caused by divergence of increased southward EP-flux. This EP-flux increases due to OGW drag before SW and extends into the Southern Hemisphere.     

Relative equilibria in quasi-homogeneous planar three body problems

In this paper we find the families of relative equilibria for the three body problem in the plane, when the interaction between the bodies is given by a quasi-homogeneous potential. The number of the relative equilibria depends on the values of the masses and on the size of the system, measured by the moment of inertia.     

Case study of inclined sporadic E layers in the Earth’s ionosphere observed by CHAMP/GPS radio occultations: Coupling between the tilted plasma layers and internal waves

We have used the radio occultation (RO) satellite data CHAMP/GPS (Challenging Minisatellite Payload/Global Positioning System) for studying the ionosphere of the Earth. A method for deriving the parameters of ionospheric structures is based upon an analysis of the RO signal variations in the phase path and intensity. This method allows one to estimate the spatial displacement of a plasma layer with respect to the ray perigee, and to determine the layer inclination and height correction values. In this paper, we focus on the case study of inclined sporadic E (E_s) layers in the high-latitude ionosphere based on available CHAMP RO data. Assuming that the internal gravity waves (IGWs) with the phase-fronts parallel to the ionization layer surfaces are responsible for the tilt angles of sporadic plasma layers, we have developed a new technique for determining the parameters of IGWs linked with the inclined E_s structures. A small-scale internal wave may be modulating initially horizontal E_s layer in height and causing a direction of the plasma density gradient to be rotated and aligned with that of the wave propagation vector k. The results of determination of the intrinsic wave frequency and period, vertical and horizontal wavelengths, intrinsic vertical and horizontal phase speeds, and other characteristics of IGWs under study are presented and discussed.     

Near-infrared study of open clusters Teutsch 10 and Teutsch 25

The astrophysical parameters have been estimated for two unstudied open star clusters Teutsch 10 and Teutsch 25 using the Two Micron All Sky Survey ($2\mathit{MASS}$) database. Radius is estimated as 4.5 arcmin for both clusters using radial density profiles. We have estimated proper motion values in both RA and DEC directions as $2.28\pm 0.3$ and $-0.38\pm 0.11$?mas?yr^?1 for Teutsch 10 and $0.48\pm 0.3$ and $3.35\pm 0.16$?mas?yr^?1 for Teutsch 25 using $\mathit{PPMXL}$¹ catalog. By estimating the stellar membership probabilities, we have identified 30 and 28 most likely members for Teutsch 10 and Teutsch 25 respectively. We have estimated the reddening as $E(B-V)=0.96\pm 0.3$?mag for Teutsch 10 and $0.58\pm 0.2$?mag for Teutsch 25, while the corresponding distances are $2.4\pm 0.2$ and $1.9\pm 0.1$?kpc. Ages of $70\pm 10$?Myr for Teutsch 10 and $900\pm 100$?Myr for Teutsch 25 are estimated using the theoretical isochrones of metallicity Z?=?0.019. The mass function slopes are derived as $1.23\pm 0.30$ and $1.09\pm 0.35$ for Teutsch 10 and Teutsch 25 respectively. Estimated mass function slope for both the clusters are close to the Salpeter value ($x=1.35$) within the errors. Estimated values of dynamical relaxation time are found to be less than cluster’s age for these objects. This concludes that both objects are dynamically relaxed. The possible reason for relaxation may be due to dynamical evolution or imprint of star formation or both.     

Midlatitude ionospheric F2-layer response to eruptive solar events-caused geomagnetic disturbances over Hungary during the maximum of the solar cycle 24: A case study

In our study we analyze and compare the response and behavior of the ionospheric F2 and of the sporadic E-layer during three strong (i.e., Dst?<??100nT) individual geomagnetic storms from years 2012, 2013 and 2015, winter time period. The data was provided by the state-of the art digital ionosonde of the Széchenyi István Geophysical Observatory located at midlatitude, Nagycenk, Hungary (IAGA code: NCK, geomagnetic latitude: 46.17° geomagnetic longitude: 98.85°). The local time of the sudden commencement (SC) was used to characterize the type of the ionospheric storm (after Mendillo and Narvaez, 2010). This way two regular positive phase (RPP) ionospheric storms and one no-positive phase (NPP) storm have been analyzed. In all three cases a significant increase in electron density of the foF2 layer can be observed at dawn/early morning (around 6:00 UT, 07:00 LT). Also we can observe the fade-out of the ionospheric layers at night during the geomagnetically disturbed time periods. Our results suggest that the fade-out effect is not connected to the occurrence of the sporadic E-layers.     

Predictive fault-tolerant control of an all-thruster satellite in 6-DOF motion via neural network model updating

The problem of controlling an all-thruster spacecraft in the coupled translational-rotational motion in presence of actuators fault and/or failure is investigated in this paper. The nonlinear model predictive control approach is used because of its ability to predict the future behavior of the system. The fault/failure of the thrusters changes the mapping between the commanded forces to the thrusters and actual force/torque generated by the thruster system. Thus, the basic six degree-of-freedom kinetic equations are separated from this mapping and a set of neural networks are trained off-line to learn the kinetic equations. Then, two neural networks are attached to these trained networks in order to learn the thruster commands to force/torque mappings on-line. Different off-nominal conditions are modeled so that neural networks can detect any failure and fault, including scale factor and misalignment of thrusters. A simple model of the spacecraft relative motion is used in MPC to decrease the computational burden. However, a precise model by the means of orbit propagation including different types of perturbation is utilized to evaluate the usefulness of the proposed approach in actual conditions. The numerical simulation shows that this method can successfully control the all-thruster spacecraft with ON-OFF thrusters in different combinations of thruster fault and/or failure.     

Study of the rheological properties of water and Martian soil simulant mixtures for engineering applications on the red planet

The rheological properties of mixtures of water and the Martian soil simulant JSC-Mars-1A are investigated by preparing and testing samples at various solids concentrations. The results indicate that the dispersion is viscoelastic and, at small timescales (～0.1?s), reacts to sudden strain as an elastic solid. At longer timescales the dispersion behaves like a Bingham fluid and exhibits a yield stress. Hysteresis loops show that rapid step-changes (2?s duration) of shear-rate result in thixotropic behaviour, but slower changes (>10?s duration) can result in rheopexy. These observations are explained with the breakdown and recovery of the packing structure under stress. The rheological information is used to generate practical tools, such as the system curve and the Moody chart that can be used for designing piping systems, and calculating pump sizes and pressure requirements.     

Evaluation of SARAL/AltiKa performance using GNSS/IMU equipped buoy in Sajafi,Imam Hassan and Kangan Ports

Performance of SARAL/AltiKa mission has been evaluated within 2016 altimeter calibration/validation framework in Persian Gulf through three campaigns conducted in the offshore waters of Sajafi, Imam Hassan and Kangan Ports, while the altimeter overflew the passes 470, 111 and 25 on 13 Feb, 7 March and 17 June 2016, respectively. As the preparation, a lightweight buoy was equipped with a GNSS receiver/choke-ring antenna and a MEMS-based IMU to measure independent datasets in the field operations. To obtain accurate sea surface height (SSH) time series, the offset of the onboard antenna from the equilibrium sea level was predetermined through surveying operations as the buoy was deploying in the onshore waters of Kangan Port. Accordingly, the double-difference carrier phase observations have been processed via the Bernese GPS Software v. 5.0 so as to provide the GNSS-derived time series at the comparison points of the calibration campaigns, once the disturbing effects due to the platform tilt and heave have been eliminated. Owing to comparing of the SSH time series and the associating altimetry 1?Hz GDR-T datasets, the calibration/validation of the SARAL/AltiKa has been performed in the both cases of radiometer and ECMWF wet troposphere corrections so as to identify potential land contamination. An agreement of the present findings in comparison with those attained in other international calibrations sites confirms the promising feasibility of Persian Gulf as a new dedicated site for calibration/validation of ongoing and future altimetry missions.     

Role of Equatorial Anomaly in Earthquake time precursive features: A few strong events over West Pacific zone

The earthquake (EQ) time coupling processes between equator-low-mid latitude ionosphere are complex due to inherent dynamical status of each latitudinal zone and qualified geomagnetic roles working in the system. In an attempt to identify such process, the paper presents temporal and latitudinal variations of ionization density (foF₂) covering 45°N to 35°S, during a number of earthquake events (M?>?5.5). The approaches adopted for extraction of features by the earthquake induced preparatory processes are discussed in the paper through identification of parameters like the ‘EQ time modification in density gradient’ defined by δ?=?(foF_{2 max}???foF_{2 min})∕τ_mm, where τ_mm – time span (in days) between EQ modified density maximum and minimum, and the Earthquake time Equatorial Anomaly, i.e. EEA, one of the most significant phenomenon which develops even during night time irrespective of epicenter position. Based on the observations, the paper presents the seismic time coupling dynamics through anomaly like manifestations between equator, low and mid latitude ionosphere bringing in the global Total Electron Content (TEC) features as supporting indices.