Search NASA⌕ Search

Engineering topics

Liying Qian

Publications and source records attributed to Liying Qian.

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2024

Storm-Time Variability of Terrestrial Hydrogen Exosphere: Kinetic Simulation Results

Recent studies of TWINS Lyman-α observations have reported an increase in geocoronal column brightness during geomagnetic storms, indicating enhanced exospheric hydrogen atom density (NH). This suggests a complex role of exospheric neutrals in determining storm-time magnetosphere dynamics and their energy release through charge-exchange processes. We developed a Model for Analyzing Terrestrial Exosphere (MATE) to investigate storm-time exospheric behaviors and their physical drivers. MATE traces test hydrogen atoms backward in time from locations in the exosphere to a nominal exobase altitude of 500 km, employing Newtonian mechanics with gravitational force. The model then calculates the phase-space densities (PSDs) of test hydrogen atoms at the exobase using the Maxwellian distribution with physics-based exobase conditions from the TIMEGCM upper atmosphere model. MATE maps PSDs at the exobase to the exosphere using Liouville’s Theorem under collisionless assumptions and derives NH by integrating the PSDs across velocity space. We conducted MATE simulation before, during, and after a minor geomagnetic storm from 12 to 18 June 2008, and compared the model results with NH estimates from the TWINS geocorona data. MATE reproduces storm-time density enhancements soon after the minimum Dst is reached, matching well with a general trend of TWINS NH estimates. The results suggest that upper atmospheric heating during a geomagnetic storm increases the number of ballistic and escaping hydrogen atoms entering the exosphere from the exobase, thereby boosting NH. However, the magnitude of modeled NH mismatches the TWINS NH estimates. The potential mechanisms of this density discrepancy include the physics excluded in the MATE model — such as neutralneutral collisions, neutral-plasma charge exchange, solar radiation pressure, and photoionization — as well as the higher exobase hydrogen density of TIMEGCM compared to typical empirical values, which will be addressed in future.

Hyunju K Connor↗
NASA NTRS2023

Hemispheric Asymmetry of the Annual and Semiannual Variation of Thermospheric Composition

We examine hemispheric asymmetry of the annual and semiannual variation of the ratio of O and N 2 concentrations ( O/N 2 ) using observations by the Global Ultraviolet Imager (GUVI) instrument onboard the Thermosphere Ionosphere Mesosphere Energetics and Dynamics satellite and compare them with Whole Atmosphere Community Climate Model with thermosphere-ionosphere eXtension (WACCM-X) model simulations. We found that in the equatorial region, the “equinox peaks” of the observed O/N 2 are near the end of March and October, and the two annual lows are near the beginning of July and January. Compared to the equatorial region, in the northern hemisphere (NH) low latitudes, the first “equinox peak” clearly shifts toward the December solstice, whereas in the southern hemisphere (SH) low latitudes, the “equinox peaks” shift toward the June solstice (JS), forming the hemispheric asymmetry characteristics of the annual and semiannual variation. Seasonal variation of O/N 2 shows no apparent phase variation with altitude, and the annual and semiannual pattern is consistent from year to year. WACCM-X reproduces the observed annual and semiannual pattern in NH but in SH, it simulates an annual variation instead of the observed annual and semiannual variation. The largest discrepancy occurs near JS in the lower and middle thermosphere: the simulated O density has an annual high near JS in SH; the simulated N 2 density has an annual high near JS in NH but a predominant annual low near JS in SH. These are not in the GUVI data. A weaker thermospheric meridional circulation in the winter hemisphere, or a reduced summer-to-winter latitudinal gradient of neutral temperature in WACCM-X simulations would make model-data comparisons more consistent.

thermosphere composition↗
NASA NTRS2023

Statistical Analysis of Equatorial Plasma Bubbles Climatology and Multi-Day Periodicity Using GOLD Observations

This study develops a new Bubble Index to quantify the intensity of 2-D postsunset equatorial plasma bubbles (EPBs) in the American/Atlantic sector, using Global-scale Observations of the Limb and Disk (GOLD) nighttime data. A climatology and day-to-day variability analysis of EPBs is conducted based on the newly-derived Bubble Index with the following results: (a) EPBs show considerable seasonal and solar activity dependence, with stronger (weaker) intensity around December (June) solstice and high (low) solar activity years. (b) EPBs exhibit opposite geomagnetic activity dependencies during different storm phases: EPBs are intensified concurrently with an increasing Kp, but are suppressed with high Kp occurring 3–6 hr earlier. (c) For the first time, we found that EPBs' day-to-day variation exhibited quasi-3-day and quasi-6-day periods. A coordinated analysis of Ionospheric Connection Explorer (ICON) winds and ionosonde data suggests that this multi-day periodicity was related to the planetary wave modulation through the wind-driven dynamo.

Ercha Aa↗
NASA NTRS2023

Long-Term Geospace Climate Monitoring

Climate change is characterized by global surface warming associated with the increase of greenhouse gas population since the start of the industrial era. Growing evidence shows that the upper atmosphere is experiencing appreciable cooling over the last several decades. The seminal modeling study by Roble and Dickinson (1989) suggested potential effects of increased greenhouse gases on the ionosphere and thermosphere cooling which appear consistent with some observations. However, several outstanding issues remain regarding the role of CO 2 , other important contributors, and impacts of the cooling trend in the ionosphere and thermosphere: for example, (1) what is the regional variability of the trends? (2) the very strong ionospheric cooling observed by multiple incoherent scatter radars that does not fit with the prevailing theory based on the argument of anthropogenic greenhouse gas increases, why? (3) what is the effect of secular changes in Earth’s main magnetic field? Is it visible now in the ionospheric data and can it explain some of the regional variability in the observed ionospheric trends? (4) what is the impact of long-term cooling in the thermosphere on operational systems? (5) what are the appropriate strategic plans to ensure the long-term monitoring of the critical space climate?

long-term trends↗
NASA NTRS2023

Seasonality of the Migrating Semidiurnal Tide in the Tropical Upper Mesosphere and Lower Thermosphere and its Thermodynamic and Momentum Budget

This work uses the Specified Dynamics-Whole Atmosphere Community Climate Model with Ionosphere/Thermosphere eXtension (SD-WACCM-X) to determine and explain the seasonality of the migrating semidiurnal tide (SW2) components of tropical upper mesosphere and lower thermosphere (UMLT) temperature, zonal wind, and meridional wind. This work also quantifies aliasing due to SW2 in satellite-based tidal estimates. Results show that during equinox seasons, the vertical profiles of tropical UMLT temperature SW2 and zonal-wind SW2’s amplitudes have a double-peak structure while they, along with meridional-wind SW2, have a single-peak structure in their amplitudes in June solstice. Hough mode reconstruction reveals that a linear combination of five SW2 Hough modes cannot fully reproduce these features. Tendency analysis reveals that for temperature, the adiabatic term, nonlinear advection term, and linear advection term are important. For the winds, the classical terms, nonlinear advection term, linear advection term, and gravity wave drag are important. Results of our alias analysis then indicate that SW2 can induce an ∼60% alias in zonal-mean and DW1 components calculated from sampling like that of the Thermosphere–Ionosphere–Mesosphere Energetics and Dynamics satellite and the Aura satellite. This work concludes that in situ generation by wave–wave interaction and/or by gravity waves plays significant roles in the seasonality of tropical UMLT temperature SW2, zonal-wind SW2, and meridional-wind SW2. The alias analysis further adds that one cannot simply assume that SW2 in the tropical UMLT is negligible.

Tides↗
NASA NTRS2023

Aura/MLS Observes and SD-WACCM-X Simulates the Seasonality, Quasi-Biennial Oscillation and El Niño-Southern Oscillation of the Migrating Diurnal Tide Driving Upper Mesospheric CO Primarily Through Vertical Advection

This work uses 17 years of upper mesospheric carbon monoxide (CO) and temperature observations by the Microwave Limb Sounder (MLS) on-board the Aura satellite to present and explain the seasonal and interannual variability of the migrating diurnal tide (DW1) component of upper mesospheric CO. This work then compares these observations to simulations by the Specified Dynamics – Whole Atmosphere Community Climate Model with Ionosphere/Thermosphere eXtension (SD-WACCM-X). Results show that, for all seasons, MLS CO local-time perturbations peaks above 85 km and has a latitude structure resembling the (1,1) mode in temperature. On the other hand, SD-WACCM-X DW1 also peaks above 85 km and has a latitude structure resembling the (1,1) mode but it simulates 2 local maximum of the (1,1) mode between 85 km and 92 km. Despite the differences in altitude structure, a tendency analysis and the adiabatic displacement method revealed that, on seasonal and interannual timescales, observed and modelled CO’s (1,1) component can be reproduced solely using vertical advection. It was also found that both observed and modelled CO’s (1,1) component contains interannual oscillations with periodicities close to that of the Quasi-biennial Oscillation and the El Niño–Southern Oscillation. From these results, this work concludes that on seasonal and interannual timescales, the observed and modelled (1,1) mode affects the global structure of upper mesospheric CO primarily through vertical advection.

Cornelius Csar Jude H Salinas↗
NASA NTRS2023

Aura/MLS Observes, and SD-WACCM-X Simulates the Seasonality, Quasi-Biennial Oscillation and El Nino Southern Oscillation of the Migrating Diurnal Tide Driving Upper Mesospheric Co Primarily Through Vertical Advection

This work uses 17 years of upper mesospheric carbon monoxide (CO) and temperature observations by the Microwave Limb Sounder (MLS) on-board the Aura satellite to present and explain the seasonal and interannual variability of the migrating diurnal tide (DW1) component of upper mesospheric CO. This work then compares these observations to simulations by the Specified Dynamics – Whole Atmosphere Community Climate Model with Ionosphere/Thermosphere eXtension (SD-WACCM-X). Results show that, for all seasons, MLS CO local-time perturbations peaks above 85 km and has a latitude structure resembling the (1,1) mode in temperature. On the other hand, SD-WACCM-X DW1 also peaks above 85 km and has a latitude structure resembling the (1,1) mode but it simulates 2 local maximum of the (1,1) mode between 85 km and 92 km. Despite the differences in altitude structure, a tendency analysis and the adiabatic displacement method revealed that, on seasonal and interannual timescales, observed and modelled CO’s (1,1) component can be reproduced solely using vertical advection. It was also found that both observed and modelled CO’s (1,1) component contains interannual oscillations with periodicities close to that of the Quasi-biennial Oscillation and the El Nino Southern Oscillation. From these results, this work concludes that on seasonal and interannual timescales, the observed and modelled (1,1) mode affects the global structure of upper mesospheric CO primarily through vertical advection.

Cornelius Csar Jude H. Salinas↗
NASA NTRS2022

Thermospheric Nitric Oxide Cooling Responses to the 14 December 2020 Solar Eclipse

The behaviors of the nitric oxide (NO) cooling in the lower thermosphere during the 14 December 2020 solar eclipse are studied using Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) measurements and WACCM-X simulations. We found that NO cooling rate decreases during the solar eclipse in both SABER measurements and WACCM-X simulations. The maximum decrease of the NO cooling is 40% in SABER measurements and 25% in WACCM-X simulations. The NO cooling process is initiated almost entirely through the collisions with atomic oxygen (O) which depends linearly on NO and O densities and non-linearly on the neutral temperature. During the eclipse, the NO concentration and temperature decreases are larger than that of O concentration. Consequently, the eclipse-time NO concentration and temperature decreases are the major drivers of the NO cooling rate decrease. The decreases of the temperature and the NO concentration contribute comparably to the eclipse-time NO cooling rate decrease.

solar eclipse↗
NASA NTRS2022

The Lower Thermospheric Winter-To-Summer Meridional Circulation: 1: Driving Mechanism

In this study, the mechanism driving the narrow lower-thermospheric winter-to-summer meridional circulation is thoroughly investigated for the first time using the Specified Dynamics configuration runs of the Whole Atmosphere Community Climate Model eXtended (SD-WACCMX) simulations and the TIMED Doppler Interferometer (TIDI) observations. The mean meridional circulation in the SD-WACCMX is qualitatively consistent with the TIDI measurements, though the magnitude in the SD-WACCMX is about 50% weaker. The lower-thermospheric winter-to-summer circulation is mainly driven by the resolved wave forcing, including the tides and internally generated inertia gravity waves (GWs). The momentum forcing from the parameterized sub-grid scale GWs is not as significant as the resolved wave forcing in driving the lower-thermospheric meridional circulation. The GW parameterization scheme in the SD-WACCMX only includes GWs with phase velocities in the range of ±45 m/s, which might result in most of the parameterized sub-grid GWs dissipating and breaking in the mesosphere and hardly impacting the lower thermosphere. Only including slow GWs in the SD-WACCMX parameterization could potentially lead to the underestimation of the meridional wind in the model. Analysis also indicates the lower-thermospheric meridional circulation is stronger in the summer hemisphere, which is attributed to the hemispheric asymmetry in the resolved wave momentum forcing. This study underlines the importance of the whole atmosphere coupling through wave propagation and dissipation. This understanding can guide the model development with an accurate representation of underlying physical processes in the mesosphere and lower thermosphere which drives the lower-thermospheric circulation as well as the overall dynamics of this region.

Jack C. Wang↗
NASA NTRS2022

Estimating the Migrating Diurnal Tide Component of Mesospheric Water Vapor

This work presents a method for estimating the migrating diurnal tide (DW1) component of mesospheric H2O from observations of the temperature tide and zonal-mean H2O made by the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) instrument. This work first shows that a 2D least-squares fit on SABER H2O yields an erroneous DW1 due to retrieval algorithm problems. Consequently, the adiabatic displacement method is used. Applying the method to SABER observations and comparing it with Microwave Limb Sounder (MLS) H2O observations reveals that the method reproduces the MLS-observed H2O DW1 component best in March and June over low-latitudes in the altitude range between 65 and 75 km where errors range between +5% and +20%. Applying the method for simulations using the Specified Dynamics—Whole Atmosphere Community Climate Model with Ionosphere/Thermosphere eXtension (SD-WACCM-X) indicates that in the model, it is best only in March and not in June. Model simulations further showed that in March, theoretical errors due to the vertical advection assumption and aliasing are best at less than +13% for the low-latitudes and altitudes between 65 and 75 km. These results, therefore, lead to the conclusion that despite the errors in the vertical advection assumption and aliasing, the adiabatic displacement method performs best in estimating mesospheric H2O DW1 with SABER observations during both March equinox and June solstice over the low latitudes between 65 and 75 km. Results also show that SD-WACCM-X satisfactorily simulates the H2O DW1 magnitude and mechanisms only in the March equinox period.

h20↗
NASA NTRS2022

Seasonal Variation of Thermospheric Composition Observed by NASA GOLD

We examine characteristics of the seasonal variation of thermospheric composition using column number density ratio ∑O/N(2) observed by the NASA Global Observations of Limb and Disk (GOLD) mission from low-mid to mid-high latitudes. We also use ∑O/N(2) derived from the Global Ultraviolet Imager (GUVI) limb measurements onboard the Thermosphere Ionosphere Mesosphere Energetics and Dynamics (TIMED) satellite and estimated by the NRLMSISE-00 empirical model to aid our investigation. We found that the ∑O/N(2) seasonal variation is hemispherically asymmetric: in the southern hemisphere, it exhibits the well-known annual and semiannual pattern, with highs near the equinoxes, and primary and secondary lows near the solstices. In the northern hemisphere, it is dominated by an annual variation, with a minor semiannual component with the highs shifting toward the wintertime. We also found that the durations of the December and June solstice seasons in terms of ∑O/N(2) are highly variable with longitude. Our hypothesis is that ion-neutral collisional heating in the equatorial ionization anomaly region, ion drag, and auroral Joule heating play substantial roles in this longitudinal dependency. Finally, the rate of change in ∑O/N(2) from one solstice season to the other is dependent on latitude, with more dramatic changes at higher latitudes.

thermosphere↗
NASA NTRS2022

Climatology of Mesosphere and Lower Thermosphere Residual Circulations and Mesopause Height Derived From SABER Observations

In the mesosphere and lower thermosphere (MLT) region, residual circulations driven by gravity wave breaking and dissipation significantly impact constituent distribution and the height and temperature of the mesopause. The distribution of CO2 can be used as a proxy for the residual circulations. Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) CO2 volume mixing ratio (VMR) and temperature measurements from 2003 to 2020 are used to study the monthly climatology of MLT residual circulations and the mesopause height. Our analyses show that (a) mesopause height strongly correlates with the CO2 VMR vertical gradient during solstices; (b) mesopause height has a discontinuity at midlatitude in the summer hemisphere, with a lower mesopause height at mid-to-high latitudes as a result of adiabatic cooling driven by strong adiabatic upwelling; (c) the residual circulations have strong seasonal variations at mid-to-high latitudes, but they are more uniform at low latitudes; and (d) the interannual variability of the residual circulations and mesopause height is larger in the Southern Hemisphere (SH; 4–5 km) than in the Northern Hemisphere (NH; 0.5–1 km).

SABER↗
NASA NTRS2021

A Comparison of the CIR- and CME-Induced Geomagnetic Activity Effects on Mesosphere and Lower Thermospheric Temperature

Neutral temperature responses in the mesosphere and lower thermosphere (MLT) to severe geomagnetic storms induced by coronal mass ejections (CMEs) are of growing interest to the space science research community. Recently, it was found that geomagnetic activities produced by the corotating interaction regions (CIRs) caused comparable effects on the Earth's upper atmosphere. In this work, we carried out a comparative study of the temperature responses in the MLT region to these two types of geomagnetic activities, using the temperature measured by the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) instruments onboard the Thermosphere, Ionosphere, Mesosphere Energetics and Dynamics (TIMED) satellite. Our results demonstrate that CIR-induced geomagnetic activity produced temperature variations in the MLT region and that this effect can penetrate downward to ∼100 km at high latitudes in both hemispheres. Temperature enhancements penetrated deeper during CME-induced geomagnetic activities, but the heating effects lasted longer during CIR-induced geomagnetic activities. There is a hemispherical asymmetry in the geomagnetical activity induced temperature changes in the MLT region. The temperature enhancements are stronger in the southern hemisphere than in the northern hemisphere during CME events.

geomagnetic activity↗
NASA NTRS2019

Annual and Semiannual Oscillations of Thermospheric Composition in TIMED/GUVI Limb Measurements

The Global UltraViolet Imager (GUVI) onboard the Thermosphere Ionosphere Mesosphere Energetics and Dynamics (TIMED) satellite provides a data set of vertical thermospheric composition (O, N2, and O2 densities) and temperature profiles from 2002–2007. Even though GUVI sampling is limited by orbital constraint, we demonstrated that the GUVI data set can be used to derive the altitude profiles of the amplitudes and phases of annual oscillation (AO) and semiannual oscillation (SAO), thereby providing important constraints on models seeking to explain these features. We performed a seasonal and interannual analysis of GUVI limb O, O2, and N2 densities and volume number density ratio O/N2 at constant pressure levels. These daytime observations of O and O/N2 in the lower thermosphere show a strong AO at midlatitudes and a clear SAO at lower latitudes. The global mean GUVI O/N2 number density ratio shows the AO, with slightly larger values in January than in July and a SAO with O/N2 greater during equinoxes than at the solstices. O and N2 densities on fixed pressure levels in the upper thermosphere are anticorrelated with solar extreme ultraviolet flux. On the other hand, O/N2 is smaller during solar minimum and larger during solar maximum. The thermospheric AO and SAO in composition have a constant phase with altitude throughout the thermosphere.

thermosphere composition↗