The in-situ utilization of lunar helium-3 resource is crucial to manned lunar landings and lunar base construction. Ilmenite was selected as the representative mineral which preserves most of the helium-3 in lunar soil. The implantation of helium-3 ions into ilmenite was simulated to figure out the concentration profile of helium-3 trapped in lunar ilmenite. Based on the obtained concentration profile, the thermal release model for molecular dynamics was established to investigate the diffusion and release of helium-3 in ilmenite. The optimal heating temperature, the diffusion coefficient, and the release rate of helium-3 were analyzed. The heating time of helium-3 in lunar ilmenite under actual lunar conditions was also studied using similitude analysis. The results show that after the implantation of helium-3 into lunar ilmenite, it is mainly trapped in vacancies and interstitials of ilmenite crystal and the corresponding concentration profile follows a Gaussian distribution. As the heating temperature rises, the cumulative amounts of released helium-3 increase rapidly at first and then tend to stabilize. The optimal heating temperature of helium-3 is about 1000 K and the corresponding cumulative release amount is about 74%. The diffusion coefficient and activation energy of helium-3 increase with the temperature. When the energy of helium-3 is higher than the binding energy of the ilmenite lattice, the helium-3 is released rapidly on the microscale. Furthermore, when the heating temperature increases, the heating time for thermal release of helium-3 under actual lunar conditions decreases. For the optimal heating temperature of 1000 K, the thermal release time of helium-3 is about 1 s. The research could provide a theoretical basis for in-situ helium-3 resources utilization on the moon.

Atmospheric drag provides an indirect approach for evaluating atmospheric mass density, which can be derived from the Precise Orbit Determination (POD) of Low Earth Orbit (LEO) satellites. A method was developed to estimate nongravitational acceleration, which includes the drag acceleration of the thermospheric density model and empirical force acceleration in the velocity direction from the centimeter-level reduced-dynamic POD. The main research achievements include the study of atmospheric responses to geomagnetic storms, especially after the launch of the spherical Qiu Qiu (QQ)-Satellite (QQ-Satellite) with the global navigation system satellite (GNSS) receiver onboard tracking the Global Positioning System (GPS) and Beidou System (BDS) data. Using this derivation method, the high-accuracy POD atmospheric density was determined from these data, resulting in better agreement among the QQ-Satellite-derived densities and the NRLMSISE-00 model densities. In addition, the POD-derived density exhibited a more sensitive response to magnetic storms. Improved accuracy of short-term orbit predictions using derived density was one of the aims of this study. Preliminary experiments using densities derived from the QQ-Satellite showed promising and encouraging results in reducing orbit propagation errors within 24 h, especially during periods of geomagnetic activity.

• Post-deployment, TianMu-1's first concurrent multi-LT observation captured May 2024's superstorm using 5 satellites. • Quantifying the LT dependence of thermospheric hemispheric asymmetry driven by May 2024's superstorm based on multi-satellite synchronous observations. • Directly resolved the entire evolutionary sequence of global multi-LT TADs, from generation through propagation to dissipation, including their LT-dependent modulation.

• We present the first joint observations of thermosphere density observed by the TM-1 constellation during a geomagnetic storm. • The TM-1 mass density data are compared with ESA's Swarm-B data and NRLMSISE00 simulations. • Response of thermosphere mass density at multiple different Local Times during the geomagnetic storm.

The longitudinal distribution of upper atmospheric density has been broadly studied. However, the studies mostly focused on 24 h averaged distribution. This study presents the longitudinal distribution of thermospheric density at dawn and dusk, using observations collected by the atmospheric density detector onboard the Chinese satellite APOD (Atmospheric Density Detection and Precise Orbit Determination) during low solar activity. The APOD observations show a significant relative longitudinal variation of thermospheric density with global maxima (Δρ_rmax) near the geomagnetic pole, especially in the winter hemisphere. The annual maximum of Δρ_rmax appears in the Southern Hemisphere around the June solstices and reaches 26.3% and 39.6% at dawn and dusk, respectively. The auroral heating and meridional wind might play a significant role in the longitudinal variation of thermospheric density. We further compare the APOD observations with the semi-empirical atmospheric model MSIS (Mass Spectrometer Incoherent Scatter Radar) 2.0 predictions under low solar activity conditions. The MSIS 2.0 model reproduces similar longitudinal variations to the observations, with hemispheric asymmetry. The longitudinal variation of thermospheric density from APOD should be related to the distribution of the atmospheric average molecular weight from the model. More observational data are needed to verify the results of this study further.

In this article, high spatiotemporal resolution data obtained by the atmospheric density detector carried by China’s APOD satellite are used to study the hemispheric asymmetry of thermospheric density. A detailed analysis is first performed on the dual magnetic storm event that occurred near the autumnal equinox on 8 September 2017. The results show that the enhancement ratio of atmospheric density in the southern polar region (SPR) on the duskside was approximately 1.33–1.65 times that of the northern polar region (NPR), demonstrating a strong hemispheric asymmetry of thermospheric atmospheric density response during the magnetic storm. However, the asymmetry response was smaller on the dawnside, suggesting that the hemispheric density response asymmetry is related to local time (LT). The energy injection in high-latitude regions increases local atmospheric density and forms traveling atmospheric disturbances (TADs). TADs can propagate to low-latitude regions over several hours and affect the global distribution of thermospheric atmospheric density. Similarly, the geomagnetic index fitting slope of SPR relative density difference is greater than that of NPR. The SuperDARN convection pattern indicates that the plasma convection velocity of SPR is significantly greater than that of NPR, indicating that joule heating caused by neutral friction of ions in the Southern Hemisphere may be stronger. Subsequently, an analysis of annual solar activity and seasons was carried out on the thermospheric NPR, SPR atmospheric density, and their differences from December 2015 to December 2020. The results show that thermospheric atmospheric density decreases overall as the number of sunspots decreases. The differences between the NPR and SPR atmospheric densities in the thermosphere exhibits a noticeable annual periodicity. The NPR and SPR atmospheric densities appear to have different distribution characteristics in different seasons. The NPR density peak is mainly in March or April. In particular, the “double-peak” phenomenon occurred in 2017, with peaks in March and September, while the most obvious feature of SPR atmospheric density is that its minimum value occurs in the summer months of June and July. This paper reveals the annual, seasonal, and magnetic storm response characteristics of the hemispheric asymmetry of thermospheric atmospheric density, which has significant implications for the study of multilayer energy coupling of the magnetosphere–ionosphere–thermosphere.