Abstract Using the global Lagrangian version of the piecewise parabolic method‐magnetohydrodynamic (PPMLR‐MHD) model, we simulate two consecutive storms in December 2015, a moderate storm on 14–15 December and a strong storm on 19–22 December, and calculate the radial diffusion coefficients (DLL) from the simulated ultralow frequency waves. We find that even though the strong storm leads to more enhancedBzandEφpower than the moderate storm, the two storms share in common a lot of features on the azimuthal mode structure and power spectrum of ultralow frequency waves. For both storms, the totalBzandEφpower is better correlated with the solar wind dynamic pressure in the storm initial phase and more correlated withAEindex in the recovery phase.Bzwave power is shown to be mostly distributed in low mode numbers, whileEφpower spreads over a wider range of modes. Furthermore, theBzandEφpower spectral densities are found to be higher at higherLregions, with a strongerLdependence in theBzspectra. The estimatedDLLbased on MHD fields shows that inside the magnetopause, the contribution from electric fields is larger than or comparable to that from magnetic fields, and our event‐specific MHD‐basedDLLcan be smaller than some previous empiricalDLLestimations by more than an order of magnitude. At last, by validating against in situ observations from Magnetospheric Multiscale spacecraft, our MHD results are found to generally well reproduce the totalBzfields and wave power for both storms, while theEφpower is underestimated in the MHD simulations.
more »
« less
Reconstruction of Extreme Geomagnetic Storms: Breaking the Data Paucity Curse
Abstract Reconstruction of the magnetic field, electric current, and plasma pressure is provided using a new data mining (DM) method with weighted nearest neighbors (NN) for strong storms with the storm activity indexSym‐H < −300 nT, the Bastille Day event (July 2000), and the 20 November 2003 superstorm. It is shown that the new method significantly reduces the statistical bias of the original NN algorithm toward weaker storms. In the DM approach the magnetic field is reconstructed using a small NN subset of the large historical database, with the subset numberKNN ≫ 1being still much larger than any simultaneous multiprobe observation number. This allows one to fit with observations a very flexible magnetic field model using basis function expansions for equatorial and field‐aligned currents, and at the same time, to keep the model sensitive to storm variability. This also allows one to calculate the plasma pressure by integrating the quasi‐static force balance equation with the isotropic plasma approximation. For strong storms of particular importance becomes the resolution of the eastward current, which prevents the divergence of the pressure integral. It is shown that in spite of the strong reduction of the dominant NN number in the new weighted NN algorithm to capture strong storm features, it is still possible to resolve the eastward current and to retrieve plasma pressure distributions. It is found that the pressure peak for strong storms may be as close as≈2.1REto Earth and its value may exceed 300 nPa.
more »
« less
- Award ID(s):
- 1702147
- PAR ID:
- 10375242
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Space Weather
- Volume:
- 18
- Issue:
- 10
- ISSN:
- 1542-7390
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We present examples of high‐latitude field‐aligned current (FAC) and toroidal magnetic potential patterns in both hemispheres reconstructed at a 2‐min cadence using an updated optimal interpolation (OI) method that ingests magnetic perturbation data provided by the Active Magnetosphere and Planetary Electrodynamics Response Experiment (AMPERE) program. A solstice and an equinoctial event are studied to demonstrate the reconstructed patterns and to provide scientific insights into FAC response to different solar wind drivers. For the 14 June 2011 high‐speed stream event with mostly northwardBzdriving, we found persistently stronger FACs in the Northern Hemisphere. Extreme interhemispheric asymmetry is associated with the interplanetary magnetic field (IMF) direction and large dipole tilt, consistent with earlier studies. FAC asymmetries seen during an isolated substorm can be attributed to dipole tilt. During relatively low geomagnetic activity, the FAC response to IMFBxchanges is identified. For the 17–18 March 2013 period, we provide global snapshots of rapid FAC changes related to an interplanetary shock passage. We further present comparisons between instantaneous and mean behaviors of FAC for the solar wind sheath passage and interplanetary coronal mass ejection southwardBzinterval and northwardBzintervals. We show that (1) sheath passage results in strong FAC and high variation in the dayside polar cap region and pre‐midnight region, different from the typical R1/R2 currents during prolonged southwardBz; (2) four‐cell reverse patterns appear during northwardBzbut are not stable; and (3) persistent dawn‐dusk asymmetry is seen throughout the storm, especially during an extreme substorm, likely associated with a dawnside current wedge.more » « less
-
Abstract Using a two‐dimensional particle‐in‐cell simulation, we investigate the effects and roles of upper‐hybrid waves (UHW) near the electron diffusion region (EDR). The energy dissipation via the wave‐particle interaction in our simulation agrees withJ · E′measured by magnetospheric multiscale (MMS) spacecraft. It means that UHW contributes to the local energy dissipation. As a result of wave‐particle interactions, plasma parameters which determine the larger‐scale energy dissipation in the EDR are changed. They‐directional current decreases while the pressure tensorPyzincreases/decreases when the agyrotropic beam density is low/high, where(x, y, z)‐coordinates correspond the(L, M, N)‐boundary coordinates. Because the reconnection electric field comes from−∂Pyz/∂z, our result implies that UHW plays an additional role in affecting larger‐scale energy dissipation in the EDR by changing plasma parameters. We provide a simple diagram that shows how the UHW activities change the profiles of plasma parameters near the EDR comparing cases with and without UHW.more » « less
-
Abstract We present a new version of the high‐resolution Kühlungsborn Mechanistic general Circulation Model (KMCM) extended toz ∼ 450 km. This model is called HIAMCM (HI Altitude Mechanistic general Circulation Model) and explicitly simulates gravity waves (GWs) down to horizontal wavelengths ofλh ∼ 165 km. We find predominant tertiary GWs in the winter thermosphere at middle/high latitudes. These GWs typically have horizontal wavelengthsλh ∼ 300–1,100 km, ground‐based periods∼25–90 min, and intrinsic horizontal phase speedscIh ∼ 250–350 m s−1. Abovez∼ 200 km, the predominant GW horizontal propagation directions are roughly against the background winds from the diurnal tide; the GWs propagate mainly poleward at midnight, eastward at 6 local time (LT), equatorward at noon, and westward at 18 LT. Wintertime GWs atz∼ 300 km having 165 km ≤λh≤ 330 km create a large hot spot over the Southern Andes/Antarctic Peninsula that agrees well with quiet time satellite measurements. Due to cancelation effects, the time‐averaged zonal mean Eliassen‐Palm flux divergence from the resolved GWs in the thermosphere is negligible compared to that of the tides and compared to the zonal component of the time‐averaged zonal mean ion drag. We also find that the thermospheric GWs dissipate mainly from macroturbulent diffusion and, abovez∼ 200 km, from molecular diffusion, whereas the tides dissipate mainly from ion drag. The averaged dissipative heating in the thermosphere due to tides is much stronger than that due to GWs.more » « less
-
Abstract In the present study we investigate the response of the dayside ground magnetic field to the sequence of interplanetary magnetic field (IMF)BYchanges during the May 2024 geomagnetic storm. We pay particular attention to its extraordinarily large (>120 nT) and abrupt flip, and use GOES‐18 (G18) magnetic field measurements in the dayside magnetosheath as a time reference. In the dayside auroral zone, the northward magnetic component changed by as much as 4,300 nT from negative to positive indicating that the direction of the auroral electrojet changed from westward to eastward. The overall sequence was consistent with the conventional understanding of the IMFBYdriving of zonal ionospheric flows and Hall currents, which is also confirmed by a global simulation conducted for this storm. Surprisingly, however, the time delay from G18 to the ground increased significantly in time. The delay was 2–3 min for a sharpBYreduction ∼30 min prior to theBYflip, but it became as long as 10 min for the zero‐crossing of theBYflip. It is suggested that the prolonged time delay reflected the travel time from G18 to the reconnection site, which sensitively depends on the final velocity at the magnetopause, that is, the inflow velocity of the magnetic reconnection. Around theBYflip, the solar wind number density transiently exceeded 100 cm−3, and should have increased further through the bow shock crossing. It is suggested that this unusually dense plasma reduced the reconnection rate, and therefore, the solar wind‐magnetosphere energy coupling due to the extraordinary IMF.more » « less
An official website of the United States government
