More Like this

1. Electron Heating in 2D Particle-in-cell Simulations of Quasi-perpendicular Low-beta Shocks

   https://doi.org/10.3847/1538-4357/ad1f69  

   Tran, Aaron ; Sironi, Lorenzo ( April 2024 , The Astrophysical Journal)

   We measure the thermal electron energization in 1D and 2D particle-in-cell simulations of quasi-perpendicular, low-beta (β_p= 0.25) collisionless ion–electron shocks with mass ratiom_i/m_e= 200, fast Mach number${\mathcal{M}}_{\mathrm{ms}}=1$–4, and upstream magnetic field angleθ_Bn= 55°–85° from the shock normal$\hat{\mathit{n}}$. It is known that shock electron heating is described by an ambipolar,B-parallel electric potential jump, Δϕ_∥, that scales roughly linearly with the electron temperature jump. Our simulations have$\mathrm{\Delta }{\varphi }_{\parallel }/(0.5m_{\mathrm{i}}{u_{\mathrm{sh}}}^2)\sim 0.1$–0.2 in units of the pre-shock ions’ bulk kinetic energy, in agreement with prior measurements and simulations. Different ways to measureϕ_∥, including the use of de Hoffmann–Teller frame fields, agree to tens-of-percent accuracy. Neglecting off-diagonal electron pressure tensor terms can lead to a systematic underestimate ofϕ_∥in our low-β_pshocks. We further focus on twoθ_Bn= 65° shocks: a${\mathcal{M}}_{\mathrm{s}}=4$(${\mathcal{M}}_{\mathrm{A}}=1.8$) case with a long, 30d_iprecursor of whistler waves along$\hat{\mathit{n}}$, and a${\mathcal{M}}_{\mathrm{s}}=7$(${\mathcal{M}}_{\mathrm{A}}=3.2$) case with a shorter, 5d_iprecursor of whistlers oblique to both$\hat{\mathit{n}}$andB;d_iis the ion skin depth. Within the precursors,ϕ_∥has a secular rise toward the shock along multiple whistler wavelengths and also has localized spikes within magnetic troughs. In a 1D simulation of the${\mathcal{M}}_{\mathrm{s}}=4$,θ_Bn= 65° case,ϕ_∥shows a weak dependence on the electron plasma-to-cyclotron frequency ratioω_pe/Ω_ce, andϕ_∥decreases by a factor of 2 asm_i/m_eis raised to the true proton–electron value of 1836.

    

   more » « less
2. Ion‐Acoustic Waves in a Quasi‐Perpendicular Earth's Bow Shock

   https://doi.org/10.1029/2022GL098640  

   Vasko, I. Y. ; Mozer, F. S. ; Bale, S. D. ; Artemyev, A. V. ( June 2022 , Geophysical Research Letters)

   We present analysis of electrostatic waves around the ramp of a quasi‐perpendicular Earth's bow shock observed by the Magnetospheric Multiscale spacecraft. The electrostatic waves have amplitudes up to 800 mV/m, which is the largest value ever reported in the Earth's bow shock. In contrast to previous studies, the electrostatic waves have large amplitudes of the electrostatic potential, up to 20 V or 20% of local electron temperature. The wavelengths are from 150 m to 3 km, that is from 15 to 300 Debye lengths and typically from 0.4 to 1.5 thermal electron gyroradii. Importantly, these waves can propagate not only quasi‐parallel or oblique, but also almost perpendicular to local magnetic field. The electrostatic waves are interpreted in terms of ion‐acoustic waves, although the presence of electron cyclotron harmonic waves cannot be entirely ruled out. These results suggest that electrostatic waves may strongly affect the dynamics of electrons in collisionless shocks.

    

   more » « less
3. Magnetospheric Multiscale Observations of Earth's Oblique Bow Shock Reformation by Foreshock Ultralow‐Frequency Waves

   https://doi.org/10.1029/2020GL091184  

   Liu, Terry Z. ; Hao, Yufei ; Wilson, III, Lynn B. ; Turner, Drew L. ; Zhang, Hui ( January 2021 , Geophysical Research Letters)

   Collisionless shocks can be nonstationary with periodic reformation shown in many simulation results, but direct observations are still tenuous and difficult to conclusively interpret. In this study, using Magnetospheric Multiscale (MMS) observations, we report direct observational evidence of Earth's oblique bow shock reformation driven by the foreshock Ultralow‐Frequency (ULF) waves. When the four MMS spacecraft were in a string‐of‐pearls formation roughly along the bow shock normal, they observed that when each period of foreshock ULF waves encountered the bow shock, a new shock ramp formed. Meanwhile, in the magnetosheath, the old bow shock's remnants were observed periodically convecting downstream. We propose that the reformation mechanism of the oblique bow shock is the variation of the upstream conditions by the periodic ULF waves as they encounter the bow shock. We also examine the nature of reflected ions during the reformation process.

    

   more » « less
4. Energy Repartition and Entropy Generation across the Earth’s Bow Shock: MMS Observations

   https://doi.org/10.3847/1538-4357/acdb7b  

   Agapitov, O. V. ; Krasnoselskikh, V. ; Balikhin, M. ; Bonnell, J. W. ; Mozer, F. S. ; Avanov, L. ( July 2023 , The Astrophysical Journal)

   The evolution of plasma entropy and the process of plasma energy redistribution at the collisionless plasma shock front are evaluated based on the high temporal resolution data from the four Magnetospheric Multiscale spacecraft during the crossing of the terrestrial bow shock. The ion distribution function has been separated into the populations with different characteristic behaviors in the vicinity of the shock: the upstream core population, the reflected ions, the gyrating ions, the ions trapped in the vicinity of the shock, and the downstream core population. The values of ion and electron moments (density, bulk velocity, and temperature) have been determined separately for these populations. It is shown that the solar wind core population bulk velocity slows down mainly in the ramp with the electrostatic potential increase but not in the foot region as it was supposed. The reflected ion population determines the foot region properties, so the proton temperature peak in the foot region is an effect of the relative motion of the different ion populations, rather than an actual increase in the thermal speed of any of the ion population. The ion entropy evaluated showed a significant increase across the shock: the enhancement of the ion entropy occurs in the foot of the shock front and at the ramp, where the reflected ions are emerging in addition to the upstream solar wind ions, the anisotropy growing to generate the bursts of ion-scale electrostatic waves. The entropy of electrons across the shock does not show a significant change: electron heating goes almost adiabatically.

    

   more » « less
5. Temporal variations in plume flux: characterizing pulsations from tilted plume conduits in a rheologically complex mantle

   https://doi.org/10.1093/gji/ggac455  

   Neuharth, Derek ; Mittelstaedt, Eric ( November 2022 , Geophysical Journal International)

   Along age-progressive hotspot volcano chains, the emplacement rate of igneous material varies through time. Time-series analysis of changing emplacement rates at a range of hotspots finds that these rates vary regularly at periods of a few to several tens of millions of years, indicative of changing melt production within underlying mantle plumes. Many hotspots exhibit at least one period between ∼2 and 10 Myr, consistent with several proposed mechanisms for changing near-surface plume flux, and thus melting rate, such as small-scale convection, solitary waves and instability formation in tilted plume conduits. Here, we focus on quantifying instability growth within plumes tilted by overlying plate motion. Previous studies using fluids with constant or temperature-dependent viscosity suggest that such instabilities should not form under mantle conditions. To test this assertion, we use a modified version of the finite element code ASPECT to simulate 400 Myr of evolution of a whole-depth mantle plume rising through the transition zone and spreading beneath a moving plate. In a 2-D spherical shell geometry, ASPECT solves the conservation equations for a compressible mantle with a thermodynamically consistent treatment of phase changes in the mantle transition zone and subject to either a temperature- and depth-dependent linear rheology or a temperature-, depth- and strain-rate dependent non-linear rheology. Additionally, we examine plume evolution in a mantle subject to a range of Clapeyron slopes for the 410 km (1–4 MPa K–1) phase transitions. Results suggest that plume conduits tilted by >67° become unstable and develop instabilities that lead to initial pulses in the transition zone followed by repeated plume pulsing in the uppermost mantle. In these cases, pulse size and frequency depend strongly on the viscosity ratio between the plume and ambient upper mantle. Based upon our results and comparison with other studies, we find that the range of statistically significant periods of plume pulsing in our models (∼2–7 Myr), the predicted increase in melt flux due to each pulse (3.8–26 × 10−5 km3 km−1 yr−1), and the time estimated for a plume to tilt beyond 67° in the upper mantle (10–50 Myr) are consistent with observations at numerous hotspot tracks across the globe. We suggest that pulsing due to destabilization of tilted plume conduits may be one of several mechanisms responsible for modulating the melting rate of mantle plumes as they spread beneath the moving lithosphere.

    

   more » « less