Search for: All records

Abstract We implement a nudging module into the Thermosphere Ionosphere Electrodynamics General Circulation Model (TIEGCM) to identify effective techniques for incorporating global‐scale tides and medium‐scale gravity waves (GWs) that induce ionospheric variability. Nudging the full fields of basic state variables minimizes contamination from spectral aliasing and mode coupling, ensuring the most accurate reproduction of each tidal component. In contrast, nudging solely diurnal tides has substantial spectral leakage into semidiurnal tides, leading to underestimations of their own amplitudes and day‐to‐day variabilities (DTDVs). Nudging both diurnal and semidiurnal tides mitigates such underestimations, establishing a minimal requirement for reproducing tidal dynamics and ionospheric DTDVs. Lower boundary forcing (LBF) causes significant deviations of tidal amplitudes and DTDVs near the boundary, but only a ∼10% underestimation above it. The DTDV of vertical ion drift gradually increases with more wave components incorporated and shows a ∼10% underestimation with LBF. Constraining geopotential height (Z*) is critical in TIEGCM to properly add GWs at lower levels. Model runs withZ* constrained exhibit reduced sensitivity to nudging levels: one‐level nudging and LBF runs show 20%–30% underestimations of TID magnitudes compared to a four‐scale‐height nudging run. Conversely, whenZ* is unavailable and onlyU,V,Tare constrained, one‐level nudging and LBF lead to 80%–90% underestimations of TIDs, with LBF entirely missing wave features. Therefore, multi‐level nudging, especially withZ* unconstrained, is recommended to incorporate GWs. Overall, nudging provides a powerful tool to realistically incorporate observed or simulated waves across medium to global scales into ionosphere‐thermosphere models, offering a data‐driven perspective of variability for lower boundary conditions. 

Abstract The 12‐year continuous observation of gravity wave momentum fluxes (GWMFs) estimated by the Mohe meteor radar (53.5°N, 122.3°E) revealed prominent intraseasonal variability around the extratropical mesopause (82–94 km) during boreal winters. Composite analysis of the December‒January‒February (DJF) season according to the Madden‒Julian Oscillation (MJO) phases revealed that the zonal GWMFs notably increased in MJO Phase 4 (P4) by ∼2–4 m²/s², and a Monte Carlo test was designed to examine the statistical significance. The response in zonal winds lags behind the GWMF response by two MJO phases (i.e., 1/2π), indicating a “force‒response” interaction between them. Additionally, time‐lagged composites revealed that strengthened westward GWMFs occurred ∼25–35 days after MJO P4, coincident with the MJO impact on the zonal winds in the stratosphere. The analysis results also suggested that the mechanism of MJO by which the MJO influences the stratospheric circulation might involve poleward propagating effects of stationary planetary waves with zonal wavenumber one. This work emphasizes the importance of GW intraseasonal variability, which impacts tropical sources from the troposphere to the extratropical mesopause. 

Abstract We use the TIEGCM‐NG nudged by MAGIC gravity waves to study the impacts of a severe thunderstorm system, with a hundred tornado touchdowns, on the ionospheric and thermospheric disturbances. The generated waves induce a distinct concentric ring pattern on GNSS TIDs with horizontal scales of 150–400 km and phase speeds of 150–300 m/s, which is well simulated by the model. The waves show substantial vertical evolution in period, initially dominated by 0.5 hr at 200 km, shifting to 0.25 hr and with more higher‐frequency waves appearing at higher altitudes (∼400 km). The TADs reach amplitudes of 100 m/s, 60 m/s, 80 K, and 10% in horizontal winds, vertical wind, temperature, and relative neutral density, respectively. Significantly perturbations in electron density cause dramatic changes in its nighttime structure around 200 km and near the EIA crest. The concentric TIDs are also simulated in ion drifts and mapped from the Tornado region to the conjugate hemisphere likely due to neutral wind‐induced electric field perturbations. The waves manage to impact the ionosphere at altitudes of ICON and COSMIC‐2, which pass through the region of interest on a total of 8 separate orbits. In situ ion density observations from these spacecrafts reveal periodic fluctuations that frequently show good agreement with the TIEGCM‐NG simulation. The O⁺fraction observations from ICON indicate that the density fluctuations are the result of vertical transport of the ions in this region, which could result from either direct forcing by neutral winds or electrodynamic coupling. 

Abstract The statistics of day‐to‐day tidal variability within 35‐day running mean windows is obtained from Michelson Interferometer for Global High‐Resolution Thermospheric Imaging (MIGHTI)/Ionospheric Connection Explorer (ICON) observations in the 90–107 km height region for the year 2020. Temperature standard deviations for 18 diurnal and semidiurnal tidal components, and for four quasi‐stationary planetary waves are presented, as function of latitude, altitude, and day‐of‐year. Our results show that the day‐to‐day variability (DTDV) can be as large as 70% of the monthly mean amplitudes, thus providing a significant source of variability for the ionospheric E‐region dynamo and hence for the F‐region plasma. We further validate our results with COSMIC‐2 ionospheric observations and present an approach to extend the MIGHTI/ICON results to all latitudes using Hough Mode Extension fitting, to produce global tidal fields and their statistical DTDV that are suitable as lower boundary conditions for nudging and ensemble modeling of TIE‐GCM. In the future, this will likely help to establish a data‐driven perspective of space weather variability caused by the tidal weather of the lower atmosphere. 

Abstract We provide observational evidence that the stability of the stratospheric Polar vortex (PV) is a significant driver of sub‐seasonal variability in the thermosphere during geomagnetically quiet times when the PV is anomalously strong or weak. We find strong positive correlations between the Northern Annular Mode (NAM) index and subseasonal (10–90 days) Global Observations of the Limb and Disk (GOLD) O/N₂perturbations at low to mid‐northern latitudes, with a largest value of +0.55 at ∼30.0°N when anomalously strong or weak (NAM >2.5 or < −2.1) vortex times are considered. Strong agreement for O/N₂variability and O/N₂‐NAM correlations is found between GOLD observations and the Whole Atmosphere Community Climate Model with thermosphere‐ionosphere eXtension (WACCM‐X) simulations, which is then used to delineate the global distribution of O/N₂‐NAM correlations. We find negative correlations between subseasonal variability in WACCM‐X O/N₂and NAM at high northern and southern latitudes (as large as −0.54 at ∼60.0°S during anomalous vortex times). These correlations suggest that PV driven upwelling at low latitudes is accompanied by corresponding downwelling at high latitudes in the lower thermosphere (∼80–120 km), which is confirmed using calculations of residual mean meridional circulation from WACCM‐X. 

Abstract Using 17 years of Modern‐Era Retrospective analysis for Research and Applications, Version 2 (MERRA‐2) data, significant responses of gravity wave (GW) variances, zonal winds and parameterized GW drag to the Madden‐Julian Oscillation (MJO) are identified globally during boreal winter, and their relations are examined. The relative anomalies of GW variances range from −4% (phase 7) to 8% (phase 4) in tropics, and −20% (phase 1) to 20% (phase 5) in the northern polar region (NPR). The anomalies of zonal winds are from −3–3 m/s and −4–8 m/s in tropics and NPR, respectively. The vertical and latitudinal structures of MJO signals in GW, wind and GW drag show coherent patterns. Further analysis implies that in the NPR, the eastward wind leads to westward momentum flux carried by the GWs. This flux leads to westward drag, which drives that of zonal winds and imprint the MJO signal in GWs to the wind. 

Abstract We develop a new methodology for the multi‐resolution assimilation of electric fields by extending a Gaussian process model (Lattice Kriging) used for scalar field originally to vector field. This method takes the background empirical model as “a priori” knowledge and fuses real observations under the Gaussian process framework. The comparison of assimilated results under two different background models and three different resolutions suggests that (a) the new method significantly reduces fitting errors compared with the global spherical harmonic fitting (SHF) because it uses range‐limited basis functions ideal for the local fitting and (b) the fitting resolution, determined by the number of basis functions, is adjustable and higher resolution leads to smaller errors, indicating that more structures in the data are captured. We also test the sensitivity of the fitting results to the total amount of input data: (a) as the data amount increases, the fitting results deviate from the background model and become more determined by data and (b) the impacts of data can reach remote regions with no data available. The assimilation also better captures short‐period variations in local PFISR measurements than the SHF and maintains a coherent pattern with the surrounding. The multi‐resolution Lattice Kriging is examined via attributing basis functions into multiple levels with different resolutions (fine level is located in the region with observations). Such multi‐resolution fitting has the smallest error and shortest computation time, making the regional high‐resolution modeling efficient. Our method can be modified to achieve the multi‐resolution assimilation for other vector fields from unevenly distributed observations. 

Abstract The 17‐year SABER‐observed gravity wave (GW) temperature variances reveal significant responses of GWs to the Madden‐Julian Oscillation (MJO) over the middle atmosphere (30–100 km) in tropics and extratropics (45°S to 45°N) for boreal winter. The responses vary significantly with latitude but barely with altitude. From 20°S to 45°N, strong positive anomalies are found for MJO Phases 3–5, while negative anomalies for Phases 7–8. From 45–20°S, these patterns are reversed. The peak‐to‐peak differences (positive‐to‐negative anomalies) are ~6–16% relative to the seasonal mean. Comparison with MJO modulations on tropical convection and polar vortex suggests that GW responses in tropics may result from the modulation of GW source, while responses in northern extratropics may result from the modulation of polar vortex, which in turn modulates GW activities. These results highlight the importance of GWs to imprint the tropical MJO signals vertically to the middle atmosphere and horizontally to extratropical regions. 

Abstract A statistical study of 18 years of diurnal temperature tides observed by the SABER instrument on board the TIMED satellite reveals a substantial response of the tides in the upper atmosphere (>60 km) to the Madden‐Julian Oscillation (MJO) in the tropical troposphere. Nonmigrating tidal amplitudes are modulated at the intraseasonal MJO periods up to ~25% relative to the seasonal mean, twice as much as for the migrating tides (~10%). We fully characterize the tidal response for active MJO days as a function of season and MJO location as prescribed by the MJO index. The MJO modulation of the tides was predicted by models but could not be unequivocally observed before. Our results further point to an important role of background winds that partly cause a different response for equatorial and nonequatorial tidal modes in different seasons, which has implications for the MJO imprint on the ionospheric dynamo region. 

Abstract The Madden‐Julian Oscillation (MJO), an eastward‐moving disturbance near the equator (±30°) that typically recurs every ∼30–90 days in tropical winds and clouds, is the dominant mode of intraseasonal variability in tropical convection and circulation and has been extensively studied due to its importance for medium‐range weather forecasting. A previous statistical diagnostic of SABER/TIMED observations and the MJO index showed that the migrating diurnal (DW1) and the important nonmigrating diurnal (DE3) tide modulates on MJO‐timescale in the mesosphere/lower thermosphere (MLT) by about 20%–30%, depending on the MJO phase. In this study, we address the physics of the underlying coupling mechanisms using SABER, MERRA‐2 reanalysis, and SD‐WACCMX. Our emphasis was on the 2008–2010 time period when several strong MJO events occurred. SD‐WACCMX and SABER tides show characteristically similar MJO‐signal in the MLT region. The tides largely respond to the MJO in the tropospheric tidal forcing and less so to the MJO in tropospheric/stratospheric background winds. We further quantify the MJO response in the MLT region in the SD‐WACCMX zonal and meridional momentum forcing by separating the relative contributions of classical (Coriolis force and pressure gradient) and nonclassical forcing (advection and gravity wave drag [GWD]) which transport the MJO‐signal into the upper atmosphere. Interestingly, the tidal MJO‐response is larger in summer due to larger momentum forcing in the MLT region despite the MJO being most active in winter. We find that tidal advection and GWD forcing in MLT can work together or against each other depending on their phase relationship to the MJO‐phases.