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.

The polar-orbiting spherical experimental satellite of China for atmospheric density detection with an altitude of ~520 km was successfully launched on 14 October 2021. Based on the dynamic inversion method for atmospheric density and the precise orbit determination data obtained by its GNSS, we inverted the orbital atmospheric density during the severe geomagnetic storm in early November 2021. In this paper, we compared the atmospheric density data obtained by the spherical satellite with the simulations of the MSISE00 and the DTM, evaluated their error distribution, and analyzed the response of the atmospheric density during the severe geomagnetic storm in the dawn–dusk orbit of 520 km altitude. The properties and the physical processes for the atmospheric density of the time evolutions in different latitudes and the global distributions during the severe geomagnetic storm were obtained. We found that the substantial disturbance enhancement and recovery of the atmospheric density of the dawn–dusk orbit have a close correlation with the geomagnetic indexes Kp and Dst. The elevation extends from the poles to the equator, and the relative variation in two hemispheres demonstrates a bimodal nearly symmetric growth structure. The maximum relative variation of the two hemispheres both occurred in the middle latitude, and, for this case, the enhancement of atmospheric density in the mid-latitude region accounted for a larger proportion. The asymmetry between the northern and southern hemispheres is demonstrated by the fact that the absolute value and absolute change in the southern hemisphere in summer are larger than those in the northern hemisphere, and the bimodal structure of the relative variation is inclined to the northern hemisphere.

Global navigation satellite system (GNSS) reflectometry (GNSS-R) developed into a promising remote sensing technique. However, few previous related studies considered the potential of its polarization. Owing to lack of sufficient in situ measurement data to support comprehensive investigation of GNSS-R polarization, this study used theoretical models and reference to our previous work to explore this topic. The commonly used microwave scattering models are employed to get the bare soil or vegetation scattering properties of GNSS-R configurations, i.e., the random surface scattering model and the first-order radiative transfer equation were improved and then employed to obtain the scattering properties of both bare soil and vegetation. Since the final output of the space-borne GNSS-R missions is a delay Doppler map (DDM), a spaceborne (DDM) simulator, oriented for the Chinese FengYun-3E (FY-3E) GNSS-R payload, was utilized to obtain the final output at different polarizations. Using the developed models (such as the bare soil and vegetation scattering models), corresponding polarization simulations were performed. That is to say, not only the commonly used LR (left hand circular polarizations (LHCP) received and the right hand circular polarizations (RHCP) received) can be presented, but also the scattering properties at RR, VR, and HR (the transmitted signals are RHCP, while the received polarizations are RHCP, vertical (V) and horizontal (H) polarizations, respectively) can be predicted by our developed models. Results reveal obvious polarization differences for the bistatic scattering and DDM. Therefore, the use of GNSS-R polarization information has potential to provide competitive and fruitful results in the future detection of land surface geophysical parameters.