Configuration space and stability analysis of solar sail near-vertical earth-trailing orbits

Regions outside the reach of traditional propulsion systems or the ones that require significant propellant, may be reached by harnessing the solar radiation pressure and leveraging coupled dynamics to maneuver a sail-based spacecraft. Earth-trailing orbits have recently been investigated for getting a unique perspective of the Sun while maintaining the spacecraft in close proximity to Earth. Vertical orbits trailing the Earth exhibit the additional capability to view the Sun from above and below the ecliptic plane. In this work, families of sail-based orbits are explored for varying Earth-trailing angles and Z amplitudes in the Sun-Earth circular restricted three-body problem. Optimization is carried out to ensure that the non-traditional vertical orbits exhibit a constant pitch angle control history, as well as symmetry across the X-Y plane. The stability of the resulting orbit families is assessed using an extension of Flouquet theory to Differential Algebraic Equations. Results indicate that sail-based Earth-trailing vertical orbits can be more stable than traditional sub-L₁ sail-based vertical orbits.     

One dimensional blood flow in a planetocentric orbit

All life on earth is accustomed to the presence of gravity. When gravity is altered, biological processes can go awry. It is of great importance to ensure safety during a spaceflight. Long term exposure to microgravity can trigger detrimental physiological responses in the human body. Fluid redistribution coupled with fluid loss is one of the effects. In particular, in microgravity blood volume is shifted towards the thorax and head. Sympathetic nervous system-induced vasoconstriction is needed to maintain arterial pressure, while venoconstriction limits venous pooling of blood prevents further reductions in venous return of blood to the heart. In this paper, we modify an existing one dimensional blood flow model with the inclusion of the hydrostatic pressure gradient that further depends on the gravitational field modified by the oblateness and rotation of the Earth. We find that the velocity of the blood flow V_B is inversely proportional to the blood specific volume d, also proportional to the oblateness harmonic coefficient J₂, the angular velocity of the Earth ω_E, and finally proportional to an arbitrary constant c. For c = −0.39073 and ξ_H = −0.5 mmHg, all orbits result to less blood flow velocities than that calculated on the surface of the Earth. From all considered orbits, elliptical polar orbit of eccentricity e = 0.2 exhibit the largest flow velocity V_B = 1.031 m/s, followed by the orbits of inclination i = 45°and 0°. The Earth’s oblateness and its rotation contribute a 0.7% difference to the blood flow velocity.     

Magnetosphere response to the 2005 and 2006 extreme solar events as observed by the Cluster and Double Star spacecraft

The four identical Cluster spacecraft, launched in 2000, orbit the Earth in a tetrahedral configuration and on a highly eccentric polar orbit (4–19.6 R_E). This allows the crossing of critical layers that develop as a result of the interaction between the solar wind and the Earth’s magnetosphere. Since 2004 the Chinese Double Star TC-1 and TC-2 spacecraft, whose payload comprise also backup models of instruments developed by European scientists for Cluster, provided two additional points of measurement, on a larger scale: the Cluster and Double Star orbits are such that the spacecraft are almost in the same meridian, allowing conjugate studies. The Cluster and Double Star observations during the 2005 and 2006 extreme solar events are presented, showing uncommon plasma parameters values in the near-Earth solar wind and in the magnetosheath. These include solar wind velocities up to ∼900 km s⁻¹ during an ICME shock arrival, accompanied by a sudden increase in the density by a factor of ∼5 and followed by an enrichment in He⁺⁺ in the secondary front of the ICME. In the magnetosheath ion density values as high as 130 cm⁻³ were observed, and the plasma flow velocity there reached values even higher than the typical solar wind velocity. These resulted in unusual dayside magnetosphere compression, detection of penetrating high-energy particles in the magnetotail, and ring current development following several successive injections of energetic particles in the inner magnetosphere, which “washed out” the previously formed nose-like ion structures.     

Intracellular molecular distributions in spacecraft experiments in orbit around Earth

It is possible that the nucleolous inside the cell plays the role of a “gravity receptor”. Furthermore, cells up to 10 μm in diameter can demonstrate some effect due to the redistribution of mitochondria or nucleolous. Effects of gravity should be present in various cell systems where larger objects such as the ribosomes move from cell to cell. In this paper we study the effects of gravity on cells. In particular, we examine the resulting intracellular molecular distribution due to Brownian motion and the ordered distribution of molecules under the action of gravity, where n₀ is the number per unit volume at certain level, and n is the number per unit volume above that level. This is an experiment that takes place at a certain orbital altitude in a spacecraft in orbit around Earth, where the acceleration due to the central field is corrected for the oblateness and also the rotation of the Earth. We found that equatorial circular and elliptical orbits have the highest n/n₀ ratios. This experiment takes place in circular and elliptical orbits, with eccentricities e = 0, 0.1 and involves a bacterial cell at an orbital altitude of 300 km. We found that n/n₀ = 1.00299 and 1.0037 respectively, which is still a 0.6–0.7 % higher than n/n₀ = 0.0996685 calculated on the surface of the Earth. Examining mitochondria in similar orbital experiments we found that equatorial orbits result to higher n/n₀ ratios. In particular, we found that n/n₀ = 8.38119, where an elliptical orbit of eccentricity e = 0.1 results to n/n₀ = 13.8525. Both are high above 100%, signifying the importance of Brownian motion over gravity. Our results are of interest to biomedical applications. Molecular concentrations are important for various processes such as the embryogenesis, positional homeostasis and its relation to cell energy expenditure, cell torque, cell deformation, and more. These results indicate that statistical molecular distributions play an important role for the recognition of a particular environment by the cell, in biological space experiment to come.     

Dynamical evolution of high area-to-mass ratio debris released into GPS orbits

A large set of simulations, including all the relevant perturbations, was carried out to investigate the long-term dynamical evolution of fictitious high area-to-mass ratio (A/M) objects released, with a negligible velocity variation, in each of the six orbital planes used by Global Positioning System (GPS) satellites. As with similar objects discovered in near synchronous trajectories, long lifetime orbits, with mean motions of about 2 revolutions per day, were found possible for debris characterized by extremely high area-to-mass ratios. Often the lifetime exceeds 100 years up to A/M ∼ 45 m²/kg, decreasing rapidly to a few months above such a threshold. However, the details of the evolution, which are conditioned by the complex interplay of solar radiation pressure and geopotential plus luni-solar resonances, depend on the initial conditions. Different behaviors are thus possible. In any case, objects like those discovered in synchronous orbits, with A/M as high as 20–40 m²/kg, could also survive in this orbital regime, with semi-major axes close to the semi-synchronous values, with maximum eccentricities between 0.3 and 0.7, and with significant orbit pole precessions (faster and wider for increasing values of A/M), leading to inclinations between 30° and more than 90°.     

The feasibility of using an L1 positioned dust cloud as a method of space-based geoengineering

In this paper a method of geoengineering is proposed involving clouds of dust placed in the vicinity of the L₁ point as an alternative to the use of thin film reflectors. The aim of this scheme is to reduce the manufacturing requirement for space-based geoengineering. It has been concluded that the mass requirement for a cloud placed at the classical L₁ point, to create an average solar insolation reduction of 1.7%, is 7.60 × 10¹⁰ kg yr⁻¹ whilst a cloud placed at a displaced equilibrium point created by the inclusion of the effect of solar radiation pressure is 1.87 × 10¹⁰ kg yr⁻¹. These mass ejection rates are considerably less than the mass required in other unprocessed dust cloud methods proposed and are comparable to thin film reflector geoengineering requirements. Importantly, unprocessed dust sourced in-situ is seen as an attractive scheme compared to highly engineered thin film reflectors. It is envisaged that the required mass of dust can be extracted from captured near Earth asteroids, whilst stabilised in the required position using the impulse provided by solar collectors or mass drivers used to eject material from the asteroid surface.     

Plasma electron temperature variability in lunar surface potential and in electric field under average solar wind conditions

The Moon is immersed in plasma environment. The most interesting challenge of the lunar plasma– field environment is that it is alternatively dominated by the extended but variable outer atmosphere of the Earth – the magnetosphere – and by the extended but highly variable solar atmosphere – the solar wind. Understanding the plasma environment and its interaction with the lunar surface will be beneficial to both manned and robotic surface exploration activities and to scientific investigations. Presented is a preliminary map of variations of lunar surface electric potential over the day side and night side using probe equations and a discussion on dust dynamics in this E-field structure using the data from Electron Reflectometer in Lunar Prospector spacecraft during 1998–1999. On the day side, potential is around 5 V and on the night side it reaches up to −82 V. On the night side region, only highly energetic electrons can overcome this large negative potential. The variation at electron temperature (T_e) strongly reflects in the surface potential. The potential reaches to a value of −82 V for T_e = 58 eV. Surface charging causes the electrostatic transport of charged dust grains. Dust grain size of 0.1 μm shows a levitation height of 4.92 m on lunar day side, 748 m on terminator region and 3.7 km on the night side. The radius of maximum sized grain to be lofted, R_max, peaks at the terminator region (R_max = 0.83 μm). At the transition region dust levitation is almost absent. This region is most suited for exploration activities as the region is free from hazards caused by lunar dust.     

Two-dimensional phase plane structure and the stability of the orbital motion for space debris in the geosynchronous ring

Based on the orbital resonance model, we study the two-dimensional phase plane structure of the motion of space debris orbiting the geosynchronous ring under the combined effects of the tesseral harmonics J₂₂, J₃₁ and J₃₃ of the Earth’s gravitational field. We present the main characteristic parameters of the two-dimensional phase plane structure. We also analyze the stability of the two-dimensional phase plane structure with numerical method. Our main findings indicate that the combined effects of the tesseral harmonics J₂₂, J₃₁ and J₃₃ fully determine the two-dimensional phase plane structure of the space debris, and it remains robust under the effect of the Earth’s actual gravitational field, the luni-solar perturbations and the solar radiation pressure with the normal area-to-mass ratios.     

High energy charged particle experiment – A payload proposal to KuaFu-B mission

High Energy Charged Particle Experiment (HECPE) is to measure the fluxes of MeV electrons and tens of MeV protons. The two satellites of KuaFu-B are in the same polar orbit with apogee 7.0R_E, perigee 1.8R_E. They can sweep large L values and pass through the inner and outer radiation belts. The high energy electrons and protons in the radiation belts are principal sources for failures of satellites and spacecrafts in the Earth orbits. The enhancements of the high energy electrons and protons, so-called energetic particle events, are important phenomena of the Space Weather. The energy ranges monitored by HECPE are 0.3–0.5 MeV, 0.5–1.0 MeV, 1.0–2.0 MeV, and E > 2.0 MeV for electrons, 5–10 MeV, 10–20 MeV, 20–40 MeV, and 40–80 MeV for protons.     

Dynamics of a spacecraft and normalization around Lagrangian points in the Neptune–Triton system

The problem of a spacecraft orbiting the Neptune–Triton system is presented. The new ingredients in this restricted three body problem are the Neptune oblateness and the high inclined and retrograde motion of Triton. First we present some interesting simulations showing the role played by the oblateness on a Neptune’s satellite, disturbed by Triton. We also give an extensive numerical exploration in the case when the spacecraft orbits Triton, considering Sun, Neptune and its planetary oblateness as disturbers. In the plane a × I (a = semi-major axis, I = inclination), we give a plot of the stable regions where the massless body can survive for thousand of years. Retrograde and direct orbits were considered and as usual, the region of stability is much more significant for the case of direct orbit of the spacecraft (Triton’s orbit is retrograde). Next we explore the dynamics in a vicinity of the Lagrangian points. The Birkhoff normalization is constructed around L₂, followed by its reduction to the center manifold. In this reduced dynamics, a convenient Poincaré section shows the interplay of the Lyapunov and halo periodic orbits, Lissajous and quasi-halo tori as well as the stable and unstable manifolds of the planar Lyapunov orbit. To show the effect of the oblateness, the planar Lyapunov family emanating from the Lagrangian points and three-dimensional halo orbits are obtained by the numerical continuation method.     

Quantifying the role of solar radiative forcing over the 20th century

The 20th century temperature anomaly record is reproduced using an energy balance model, with a diffusive deep ocean. The model takes into account all the standard radiative forcings, and in addition the possibility of a non-thermal solar component. The model is parameterized and then optimized to produce the most likely values for the climate parameters and radiative forcings which reproduce the 20th century global warming. We find that the best fit is obtained with a negligible net feedback. We also show that a non-thermal solar component is necessarily present, indicating that the total solar contribution to the 20th century global warming, of ΔT_solar = 0.27 ± 0.07 °C, is much larger than can be expected from variation in the total solar irradiance alone. However, we also find that the largest contribution to the 20th century warming comes from anthropogenic sources, with ΔT_man = 0.42 ± 0.11 °C.     

Evolution of the orbital elements for objects with high area-to-mass ratios in geostationary transfer orbits

A new population of uncatalogued objects in geosynchronous Earth orbits (GEO), with a mean motion of about 1 rev/day and eccentricities up to 0.6, has been identified recently. The first observations of this new type of objects were acquired in the framework of the European Space Agency’s (ESA) search for space debris in GEO and the geostationary transfer orbit (GTO) using the ESA 1-m telescope on Tenerife. Earlier studies have postulated that the perturbations due to the solar radiation pressure can lead to such large eccentricities for GEO objects with a high area-to-mass ratio (A/M). The simulations showed that the eccentricities of GEO objects with large A/M exhibit periodic variations with periods of about one year and amplitudes depending on the value of A/M. The findings of these studies could be confirmed by observations from the ESA 1-m telescope on Tenerife.     

Analytic and numerical treatment of motion of dust grain particle around triangular equilibrium points with post-AGB binary star and disc

This paper investigates the motion around the triangular equilibrium points, of a passively gravitating dust particle in the gravitational field of a low-mass post-AGB binary system, surrounded by circumbinary disc. The two bodies of the binary are modeled as a triaxial star and a radiating-oblate star. Due to small deviation of disc stars on circular orbits, we have assumed that the Coriolis and centrifugal forces of the stars are slightly perturbed. The triangular equilibrium points of the particle are found. These points are defined by, triaxiality of the primary star, oblateness and radiation of the secondary one and the gravitational potential from the disc mass. Further, when the disc mass increases, the particle moves nearer to the stars and farther away from the disc. In general, these equilibrium points are linearly stable when μ < μ_C; where μ is the mass ratio and μ_C is the critical mass function, defined by the parameters of the system. The effects of each of these parameters on the size of the stability region are stated, and the periodic motion around the stable points is examined. It is seen that the orbits are ellipses, and the orientation, eccentricities, lengths of the semi-major and semi-minor axes are influenced by the parameters of the problem. In particular, for our numerical linear stability analysis, we have taken an extremely depleted pulsating star, IRAS 11472-0800 as the post-AGB triaxial star, with a weakly-radiating young white dwarf star; G29-38 as the secondary. For this system, the stability result of the triangular points comes out different. Here, μ_C < μ throughout the entire range of the mass ratio and the critical mass function. Hence, the triangular equilibrium points are unstable. The stability of the orbits is tested using the Poincaré surfaces of section (Pss). The region of stability is controlled by the introduced parameters and the Jacobi constant.     

On the practical exploitation of perturbative effects in low Earth orbit for space debris mitigation

This paper presents the results of a numerical evaluation of the natural lifetime reduction in low Earth orbit, due to dynamical perturbations. The study considers two values for the area-to-mass ratio, a nominal ratio which resembles a typical value of spacecraft in orbit today, and an enhanced ratio which covers the surface augmentation. The results were obtained with two orbit propagators, one of a semi-analytical nature and the second one using non-averaged equations of motion. The simulations for both propagators were set up similarly to allow comparison. They both use the solar radiation pressure and the secular terms of the geopotential ($J_2,J_4$ and $J_6$). The atmospheric drag was turned on and off in both propagators to alternatively study the eccentricity build up and the residual lifetime. The non-averaging case also covers a validation with the full 6?×?6 geopotential. The results confirm the findings in previous publications, that is, the possibility for de-orbiting from altitudes above the residual atmosphere if a solar sail is deployed at the end-of-life, due to the combined effect of solar radiation pressure and the oblateness of the Earth. At near polar inclinations, shadowing effects can be exploited to the same end. The results obtained with the full, non-averaging propagator revealed additional de-orbiting corridors associated with solar radiation pressure which were not found by previous work on space debris mitigation. The results of both tools are compared for specific initial conditions. For nominal values of area-to-mass ratio, instead, it is confirmed that this resonance effect is negligible.The paper then puts the findings in the perspective of the current satellite catalogue. It identifies space missions which are currently close to a resonance corridor and shows the orbit evolution within the resonances with a significantly shorter residual orbital lifetime. The paper finishes with a discussion on the exploitation of these effects with regards to the long-term simulation of the space debris environment and a flux and collision probability comparison.     

Semi-analytical investigations of high area-to-mass ratio geosynchronous space debris including Earth’s shadowing effects

This paper investigates the long-term perturbations of the orbits of geosynchronous space debris influenced by direct radiation pressure including the Earth’s shadowing effects. For this purpose, we propose an extension of our homemade semi-analytical theory [Valk, S., Lemaître, A., Deleflie, F. Semi-analytical theory of mean orbital motion for geosynchronous space debris under gravitational influence. Adv. Space Res., submitted for publication], based on the method developed by Aksnes [Aksnes, K. Short-period and long-period perturbations of a spherical satellite due to direct solar radiation. Celest. Mech. Dyn. Astron. 13, 89–104, 1976] and generalized into a more convenient non-singular formalism. The perturbations accounting for the direct radiation pressure with the Earth’s shadow are computed on a revolution-by-revolution basis, retaining the original osculating Hamiltonian disturbing function. In this framework, we compute the non-singular mean longitude at shadow entry and shadow exit at every orbital revolution in opposition to classical approaches where the singular eccentric anomalies at shadow entry and shadow exit are computed. This new algorithm is developed using non-singular variables. Consequently, it is particularly suitable for both near-circular and near-equatorial orbits as well as orbits which transit periodically around null eccentricities and null inclinations.The algorithm is tested by means of numerical integrations of the equations, averaged over the short periods, including radiation pressure, J₂, the combined Moon and Sun third body attraction as well as the long-term effects of the 1:1 resonance occurring for geosynchronous objects. As an extension of [Valk, S., Lemaître, A., Anselmo, L. Analytical and semi-analytical investigations of geosynchronous space debris with high area-to-mass ratios influenced by solar radiation pressure. Adv. Space Res., doi:10.1016/j.asr.2007.10.025, 2007b], we especially apply our analysis to space debris with area-to-mass as high as 20 m²/kg. This paper provides numerical and semi-analytical investigations leading to a deep understanding of the long-term evolution of the semi-major axis. Finally, these semi-analytical investigations are compared with accurate numerical integrations of the osculating equations of motion over time scales as high as 25 years.     

An empirical model of electron density in low latitude at 600 km obtained by Hinotori satellite

An empirical model of electron density (Ne) was constructed by using the data obtained with an impedance probe on board Japanese Hinotori satellite. The satellite was in circular orbit of the height of 600 km with the inclination of 31 degrees from February 1981 to June 1982. The constructed model gives Ne at any local time with the time resolution of 90 min and between −25 and 25 degrees in magnetic latitude with its resolution of 5 degrees in the range of F_10.7 from 150 to 250 under the condition of K_p < 4. Spline interpolations are applied to the functions of day of year, geomagnetic latitude and solar local time, and linear interpolation is applied to the function of F_10.7. Longitude dependence of Ne is not taken into account. Our density model can reproduce solar local time variation of electron density at 600 km altitude better than current International Reference Ionosphere (IRI2001) model which overestimates Ne in night time and underestimates Ne in day time. Our density model together with electron temperature model which has been constructed before will enable more understanding of upper ionospheric phenomenon in the equatorial region.     

Highly Elliptical Orbits for Arctic observations: Assessment of ionizing radiation

The ionizing radiation environment was analyzed for a variety of potential Highly Elliptical Orbits (HEOs) with orbital periods ranging from 6 h to 24 h suitable to continuously monitor the Arctic region. Several models available from the ESA Space Environment Information System (SPENVIS) online tool were employed, including the new-generation AE9/AP9 model for trapped radiation. Results showed that the Total Ionizing Dose (TID) has a well-pronounced local minimum for the 14-h orbit, which is nearly identical to the overall minimum observed for the longest orbital period (24 h). The thickness of slab aluminum shielding required to keep the annual TID below 10, 5 and 3.33 krad (i.e. 150, 75 and 50 krad for 15 years of mission duration) for a 14-h orbit is 2.1, 2.7 and 3.1 mm respectively. The 16-h orbit requires an additional 0.5 mm of aluminum to achieve the same results, while the 24-h orbit requires less shielding in the order of 0.2–0.3 mm. Comparison between the AE8/AP8 and AE9/AP9 models was conducted for all selected orbits. Results demonstrated that differences ranged from −70% to +170% depending on orbit geometry.     

The effect of geocenter motion on Jason-2 orbits and the mean sea level

We compute a series of Jason-2 GPS and SLR/DORIS-based orbits using ITRF2005 and the std0905 standards ( Lemoine et al., 2010). Our GPS and SLR/DORIS orbit data sets span a period of 2 years from cycle 3 (July 2008) to cycle 74 (July 2010). We extract the Jason-2 orbit frame translational parameters per cycle by the means of a Helmert transformation between a set of reference orbits and a set of test orbits. We compare the annual terms of these time-series to the annual terms of two different geocenter motion models where biases and trends have been removed. Subsequently, we include the annual terms of the modeled geocenter motion as a degree-1 loading displacement correction to the GPS and SLR/DORIS tracking network of the POD process. Although the annual geocenter motion correction would reflect a stationary signal in time, under ideal conditions, the whole geocenter motion is a non-stationary process that includes secular trends. Our results suggest that our GSFC Jason-2 GPS-based orbits are closely tied to the center of mass (CM) of the Earth consistent with our current force modeling, whereas GSFC’s SLR/DORIS-based orbits are tied to the origin of ITRF2005, which is the center of figure (CF) for sub-secular scales. We quantify the GPS and SLR/DORIS orbit centering and how this impacts the orbit radial error over the globe, which is assimilated into mean sea level (MSL) error, from the omission of the annual term of the geocenter correction. We find that for the SLR/DORIS std0905 orbits, currently used by the oceanographic community, only the negligence of the annual term of the geocenter motion correction results in a – 4.67 ± 3.40 mm error in the Z-component of the orbit frame which creates 1.06 ± 2.66 mm of systematic error in the MSL estimates, mainly due to the uneven distribution of the oceans between the North and South hemisphere.     

GSFC DORIS contribution to ITRF2008

The NASA GSFC DORIS analysis center has provided weekly DORIS solutions from November 1992 to January 2009 (839 SINEX files) of station positions and Earth Orientation Parameters for inclusion in the DORIS contribution to ITRF2008. The NASA GSFC GEODYN orbit determination software was used to process the orbits and produce the normal equations. The weekly SINEX gscwd10 submissions included DORIS data from Envisat, TOPEX/Poseidon, SPOT-2, SPOT-3, SPOT-4, SPOT-5. The orbits were mostly seven days in length (except for weeks with data gaps or maneuvers). The processing used the GRACE-derived EIGEN-GL04S1 gravity model, updated modeling for time-variable gravity, the GOT4.7 ocean tide model and tuned satellite-specific macromodels for SPOT-2, SPOT-3, SPOT-4, SPOT-5 and TOPEX/Poseidon. The University College London (UCL) radiation pressure model for Envisat improves nonconservative force modeling for this satellite, reducing the median residual empirical daily along-track accelerations from 3.75 × 10⁻⁹ m/s² with the a priori macromodel to 0.99 × 10⁻⁹ m/s² with the UCL model. For the SPOT and Envisat DORIS satellite orbits from 2003 to 2008, we obtain average RMS overlaps of 0.8–0.9 cm in the radial direction, 2.1–3.4 cm cross-track, and 1.7–2.3 cm along-track. The RMS orbit differences between Envisat DORIS-only and SLR & DORIS orbits are 1.1 cm radially, 6.4 cm along-track and 3.7 cm cross-track and are characterized by systematic along-track mean offsets due to the Envisat DORIS system time bias of ±5–10 μs. We obtain a good agreement between the geometrically-determined geocenter parameters and geocenter parameters determined dynamically from analysis of the degree one terms of the geopotential. The intrinsic RMS weekly position repeatability with respect to the IDS-3 combination ranges from 2.5 to 3.0 cm in 1993–1994 to 1.5 cm in 2007–2008.     

Charging of dust grains by anisotropic solar wind multi-component plasma

In this paper we study the charging process of small grain particles by anisotropic multi-component solar wind plasmas (electrons, protons and heavy ions), versus two-component (electron/proton) plasmas. We are focusing attention on the important characteristics of the charging process, namely the charging time, floating potential and current content as functions of plasma parameters such as He⁺⁺/H⁺ (α/p) number density and T_α/T_p temperature ratios of alpha particles to protons, as well as plasma streaming velocity v₀. Measured statistical properties of solar wind plasma parameters at 1 AU show considerable variations in α/p-temperature ratios from 1 to 10, in α/p-number density ratio from 0.01 to 0.35, as well as in values of streaming velocity v₀ from 200 km/s to 1000 km/s and more. Periods of these variations could last for several days each, leading to significant variability in the charging process, according to newly derived general analytical expressions. Numerical calculations performed for protons/alphas plasmas showed large disparity in the charging characteristics. For example, in anisotropic plasma, grain charging time varies up to 90% depending on α/p-particles temperature and number density ratios, whereas changes in floating potential are up to 40%. In contrast, in isotropic plasma, charging characteristic for grains do not change very much for the same plasma parameters variations, with charging time varying about 12% and floating potential only varying about 4%. It is also shown that in highly anisotropic plasma, with all ballistic electrons and ions, dust grains could not hold their charges, and characteristic discharged time is calculated. We note that the analysis is equally applicable to any sized body immersed in solar wind plasma.