Solar Wind Protons Research Articles

The Alpha-proton Differential Flow in the Alfvénic Young Solar Wind: from Sub-Alfvénic to Super-Alfvénic

Data obtained from Parker Solar Probe (PSP) since 2021 April have shown the first in situ observation of the solar corona, where the solar wind is formed and accelerated. Here, we investigate the alpha-proton differential flow and its characteristics across the critical Alfvén surface (CAS) using data from PSP during encounters 8–10 and 12–13. We first show the positive correlation between the alpha-proton differential velocity and the bulk solar wind speed at PSP encounter distances. Then we explore how the characteristics of the differential flow vary across the CAS and how they are affected by Alfvénic fluctuations including switchbacks. We find that the differential velocity below the CAS is generally smaller than that above the CAS, and the local Alfvén speed well limits the differential speed both above and below the CAS. The deviations from the alignment between the differential velocity and the local magnetic field vector are accompanied by large-amplitude Alfvénic fluctuations and decreases in the differential speed. Moreover, we observe that V α p increases from M A < 1 to M A ≃ 2 and then starts to decrease, which suggests that alphas may remain preferentially accelerated well above the CAS. Our results also reveal that in the sub-Alfvénic solar wind both protons and alphas show a strong correlation between their velocity fluctuations and magnetic field fluctuations, with a weaker correlation for alphas. By contrast, in the super-Alfvénic regime the correlation remains high for protons, but is reduced for alphas.

Intermittent Electron-only Reconnection at Lunar Mini-magnetospheres

Lunar crustal magnetic anomalies (LCMA) are sub-ion-gyroradius structures that have been shown to stand off the solar wind (SW) plasma from the Moon’s surface, forming shock-like discontinuities and reflecting incident SW protons. In this Letter, the results of high-resolution, 2D fully kinetic simulations show a bursty electron-only magnetic reconnection in the SW-LCMA interaction region, characterized by the quasiperiodic formation and ejection of magnetic islands and strong parallel electron flows along the X-point separator lines. The islands are observed to modify the magnetic pressure pileup and Hall electric field above the LCMA, leading to sharp increases in reflected protons that drive electromagnetic fluctuations downstream and short distances upstream in the SW.

A Silicon-Photo-Multiplier-Based Camera for the Terzina Telescope on Board the Neutrinos and Seismic Electromagnetic Signals Space Mission

NUSES is a pathfinder satellite project hosting two detectors: Ziré and Terzina. Ziré focuses on the study of protons and electrons below 250 MeV and MeV gamma rays. Terzina is dedicated to the detection of Cherenkov light produced by ultra-high-energy cosmic rays above 100 PeV and ultra-high-energy Earth-skimming neutrinos in the atmosphere, ensuring a large exposure. This work mainly concerns the description of the Cherenkov camera, composed of SiPMs, for the Terzina telescope. To increase the data-taking period, the NUSES orbit will be Sun-synchronous (with a height of about 550 km), thus allowing Terzina to always point toward the dark side of the Earth’s limb. The Sun-synchronous orbit requires small distances to the poles, and as a consequence, we expect an elevated dose to be received by the SiPMs. Background rates due to the dose accumulated by the SiPM would become a dominant contribution during the last two years of the NUSES mission. In this paper, we illustrate the measured effect of irradiance on SiPM photosensors with a variable-intensity beam of 50 MeV protons up to a 30 Gy total integrated dose. We also show the results of an initial study conducted without considering the contribution of solar wind protons and with an initial geometry with Geant4. The considered geometry included an entrance lens as one of the options in the initial design of the telescope. We characterize the SiPM output signal shape with different μ-cell sizes. We describe the developed parametric SiPM simulation, which is a part of the full Terzina simulation chain.

Skewness and kurtosis of solar wind proton distribution functions: The normal inverse-Gaussian model and its implications

Context.In the solar wind (SW), the particle distribution functions are generally not Gaussian. They present nonthermal features that are related to underlying acceleration and heating processes. These processes are critical in the overall dynamics of this expanding astrophysical fluid.Aims.The Proton Alpha Sensor (PAS) on board Solar Orbiter commonly observes skewed proton distributions, with a more populated high-energy side in the magnetic field direction than the Gaussian distribution. Our objectives are: (1) to identify a theoretical statistical function that adequately models the observed distributions and (2) to use its statistical interpretation to constrain the acceleration and heating processes.Methods.We analyzed the 3D velocity distribution functions (VDFs) measured by PAS and compared them to model statistical functions.Results.We show that the normal inverse Gaussian (NIG), a type of hyperbolic statistical distribution, provides excellent fits of skewed and leptokurtic proton distributions. NIG can model both the core distribution and the beam, if present. We propose an interpretation that is inspired by the mathematical formulation of the NIG. It assumes that the acceleration or heating mechanism can be modeled as a drifting diffusion process in velocity space, controlled (or subordinated) by the time of interaction of the particles with “accelerating structures”. The probability function of the interaction time is an inverse Gaussian (IG), obtained by considering a random drift across structures of a given size. The control of the diffusion by interaction times that follow an IG probability function formally defines the NIG distribution. Following this model, we show that skewness and kurtosis can be used to estimate the kinetic and thermal energy gains provided by the interaction with structures. For example, in the case studies presented here, the analyzed populations would have gained kinetic energy representing approximately two to four times their thermal energy, with an increase in velocity – due to acceleration – of from one-tenth to one-third of the observed flow velocity. We also show that the model constrains the initial temperature of the populations.Conclusions.Overall, the NIG model offers excellent fits of the observed proton distributions. Combining the skewness and the kurtosis, it also leads to constraints in the part of acceleration and heating due to the interactions with structures in the formation of the proton populations. We suggest that these effects add to the classical thermal evolution of the bulk velocity and temperature resulting from SW expansion.

Dependences of Energetic Neutral Atoms Energy on the Solar Wind Energy and Solar Zenith Angle Observed by the Chang’E-4 Rover

Solar wind can directly interact with the lunar surface and bring a space weathering effect. Some solar wind protons can be scattered as energetic neutral atoms (ENAs), which include rich information of the solar wind–surface interaction. However, people still know little about the ENA truth on the lunar ground due to the lack of in situ measurements. Different from the previous in-orbit measurements, here we present the first ground-based ENA measurements by the Chang’E-4 rover and find a good correlation between the mean ENA energy and the solar wind energy. Moreover, the loss rate of ENA energy can strongly depend on both the solar wind energy and the solar zenith angle (SZA), in which the energy loss rate can be enhanced by 73% when the solar wind energy increases from 400 to 1400 eV and can be reduced by 32% when the SZA increases from 57° to 71°. Combined with numerical simulations by SDTrimSP code, we propose that the solar wind protons can penetrate deeper into the lunar surface with a longer path length when the solar wind energy is higher or the SZA is lower, which results in a larger energy loss rate for the scattered ENAs. Our results provide an important constraint for the solar wind–surface research and have general implications in studying the surficial space weathering of the Moon and other airless bodies.

The nature of the mesoscale field-aligned currents in the auroral oval for positive IMF BZ: More frequent occurrence in the dawnside sector than in the duskside sector

To understand the similarities and differences between the duskside mesoscale quasistatic field-aligned currents (FACs) and the dawnside mesoscale FACs, we examined the magnetic field data obtained through the constellation measurements of Swarm satellites for a four-year period. The automatic event identification method developed in the previous study (Yokoyama Y et al., 2021) identified events of quasistatic mesoscale FACs, which are embedded in the diminished dawnside and duskside Region 1/2 current systems, in 774 passes out of 4001 passes on the dawnside, and 443 out of 3755 passes on the duskside, respectively. The dawnside and duskside mesoscale FACs have similar occurrence tendencies such that both have a relatively high occurrence ratio for positive IMF BZ, occur regardless of sunlight conditions, and have the current density increasing with the solar wind proton density. On the other hand, one notable difference was found; the occurrence ratio of the dawnside mesoscale FACs is approximately 1.7 times higher than that of the duskside mesoscale FACs. This difference is discussed in terms of the dawn–dusk asymmetry of the occurrence of Kelvin Helmholtz instabilities on the magnetopause.

Periodic Mesoscale Density Structures Comprise a Significant Fraction of the Solar Wind and Are Formed at the Sun

AbstractMesoscale density structures in the solar wind are often periodic, with f ∼ 0.1–5 mHz. They are trains of advected density structures with radial length scales of ∼100–10,000 Mm. While studies have shown that these periodic density structures (PDSs) are often formed at the Sun and released into the solar wind, it is unknown what percent of the solar wind at 1 AU is comprised of PDSs from the Sun, as opposed to periodicities formed through dynamics en route. We expand on Kepko et al. (2020, https://doi.org/10.1029/2020ja028037) which analyzed 25 years of in situ solar wind proton data, and include here alphas to examine the compositional characteristics of PDSs. Compositional changes, such as the alpha‐to‐proton ratio (α/p), are frozen into the solar wind plasma low in the corona, and so do not evolve as the solar wind advects and fills the Heliosphere. We find a broad occurrence enhancement in both the proton and α/p distributions between 1 and 3 mHz, centered near ∼2.1 mHz, and demonstrate that this distribution can be modeled assuming ∼30% of the solar wind segments contain a PDS. We find a distinct distribution below 1 mHz, with markedly different α/p characteristics. The α/p indicates that both populations are from the Sun, with likely different generation mechanisms. We conclude by summarizing mechanisms at the Sun that could produce periodic mass release, namely, periodic magnetic reconnection. The lower frequency PDSs likely involve reconnection at S‐web arcs including the heliospheric current sheet, while the higher frequency PDSs may be driven by p‐mode‐related transverse coronal oscillations.

About the magnetic barrier of Venus

The magnetized supersonic solar wind, when flowing around planets, forms a magnetic barrier near the streamlined surface. The main feature of the magnetic barrier is that the magnetic pressure prevails over the plasma pressure. The Hall-MHD model is used to simulate the magnetic barrier in the case of solar wind flow around the atmosphere of Venus. The obtained numerical results are compared with an analytical approximation of the magnetic barrier thickness, which expresses the dependence of the magnetic barrier on solar wind parameters. Particular attention is paid to the physical reasons for the asymmetry of the magnetic barrier caused by the Hall effects, which are mainly concentrated in the boundary layer near the ionopause, where the electric current has a maximum strength. An additional source of asymmetry is also considered, which acts in the same direction and is associated with the influence of the normal component of the electric field on the specific behavior of new atmospheric ions. It is shown that solar wind protons are loaded by new atmospheric ions mainly in the E+ hemisphere. In the case of more intense loading, the boundary of the magnetic barrier and the shock wave are located farther from the ionopause.

Analysis of the first coronagraphic multi-band observations of a sungrazing comet

Context. Between 24 and 25 December 2021 a sungrazing comet (SOHO-4341) approached the Sun, being observed by “classical” visible light (VL) coronagraphs on board the SOHO and STEREO missions, and also by the innovative Metis coronagraph on board the ESA-NASA Solar Orbiter mission in the VL and ultraviolet (UV H I Lyman-α) band. Aims. We show how VL data acquired by the Metis coronagraph can be combined with those provided by other space-based coronagraphs to reconstruct the comet orbit, but also to provide information on the dust composition from the polarized VL emission. Moreover, we show how the UV emission can be employed to measure local plasma parameters of the ambient solar wind. Methods. By using the comet positions tracked with VL Metis images (with spatial resolution that is four times better than UV), the UV images (with a time cadence that is five times faster than VL) have been coaligned to maximize the signal-to-noise ratio in the UV band. The local electron density ne was measured from the observed exponential decay of the UV Lyman-α intensity along the tail, while the solar wind speed vwind was measured from the UV Lyman-α tail inclination with respect to the cometary orbital path deprojected in 3D. Moreover, the proton kinetic temperature Tk was also obtained by the aperture angle of the UV Lyman-α tail. Results. When the comet was at an average heliocentric distance of 14.3 R⊙, the comet had a radial speed of 155 km s−1 and a tangential speed of 59 km s−1. The comet had a UV Lyman-α tail extending in the anti-solar direction over more than 1.5 R⊙. From the analysis of the tail shape in UV we obtained the local solar wind speed (vwind = 190 km s−1), electron density (ne = 1.5 × 104 cm−3), and proton temperature (Tk = 1.2 × 106 K). Moreover, theoretical analysis of the measured UV Lyman-α intensity allowed us to estimate the radius of the cometary nucleus (Rcom = 65 m) and the water outgassing rate (QH2O = 4.8 × 1028 molec s−1). Conclusions. These results show that sungrazing comets are unique “local probes” for the ambient coronal plasma, providing measurements that are not as affected by the line-of-sight integration effects as those provided by remote sensing instruments, in regions of the Heliosphere that are not explored in situ by the ongoing space missions.

Properties of Quiet Magnetotail Plasma Sheet at Lunar Distances

AbstractThe Earth's magnetotail at lunar distances (R ≈ 60 RE) serves as a unique laboratory to study plasma dynamics in a weak, highly fluctuating magnetic field with a strong magnetic field gradient. In particular, studies of the of quiet‐time plasma in the lunar‐distant magnetotail can inform us about plasma entry to the magnetosphere and sources of magnetospheric plasma populations. We use the data collected by the two lunar‐orbiting Acceleration, Reconnection, Turbulence and Electrodynamics of Moon's Interaction with the Sun (ARTEMIS) spacecraft during its 2013–2019 magnetotail traversals to infer the average thermodynamic and spectral properties of plasma populations in the quiet‐time (i.e., low geomagnetic activity and absence of fast plasma flows) magnetotail at lunar distances. We found that plasma temperature and density in the quiet‐time magnetotail at R ≈ 60 RE are organized by the magnetic field. Three distinct regions with plasma β ≫ 1, β ∼ 1, and β ≪ 1, the central plasma sheet (CPS), outer plasma sheet (OPS), and lobes are sampled. We found that temperatures and energy spectra of ion populations in CPS, OPS, and lobes regimes are different: the hotter CPS temperatures scale with the kinetic energy of solar wind protons; cold OPS/lobe ions are, likely, of ionospheric origin. The ion and electron particle spectra in CPS, OPS, and lobes are nonthermal and reasonably well fitted by the Kappa function, with κ exponent varying between 2.5 and 3.5.

Extreme Rarefaction of Solar Wind: A Study on Origin and Characteristics Using Ulysses Observations

From the Ulysses observation of the solar wind between the heliocentric distance r h of ∼1.0 and ∼5.4 au during 1990–2009, we identified 53 intervals when the solar wind exhibited extreme rarefaction, ∼2 orders of magnitude decreases in the solar wind proton density N p from their ambient values. These extremely low-density solar wind (ELDSW) events, characterized by an average (median) N p of ∼0.28 ± 0.09 (∼0.30) cm−3, ram pressure of ∼0.07 ± 0.04 (∼0.07 nPa) and mass flux of ∼166 ± 84 (∼159) 10−22 kg cm−2 s−1 all normalized to 1 au, have an average (median) duration of ∼6.0 ± 3.5 days (∼5.5 days), and radial extent of ∼1.9 ± 1.1 au (∼1.9 au). A clear hemispheric asymmetry is noted in their solar/interplanetary origin, with 70% being identified in the south hemisphere, and 30% in the north hemisphere of the heliosphere. About 23% of the events were encountered between r h of 2.25 and 4 au, and 77% at r h > 4 au, indicating that these are not intrinsic properties of the Sun/solar corona but are created by the evolution of the solar wind with increasing radial distance from the Sun. The majority (49%) of the events occurred during magnetic clouds, 34% in solar wind high-speed stream (HSS) tails, 11% during the proper HSSs, and 6% during interplanetary sheaths. The identification of ELDSWs will have important consequences for their interaction with the magnetospheres of Jupiter and Saturn.

The Effects of Nonequilibrium Velocity Distributions on Alfvén Ion-cyclotron Waves in the Solar Wind

In this work, we investigate how the complex structure found in solar wind proton velocity distribution functions (VDFs), rather than the commonly assumed two-component bi-Maxwellian structure, affects the onset and evolution of parallel-propagating microinstabilities. We use the Arbitrary Linear Plasma Solver, a numerical dispersion solver, to find the real frequencies and growth/damping rates of the Alfvén modes calculated for proton VDFs extracted from Wind spacecraft observations of the solar wind. We compare this wave behavior to that obtained by applying the same procedure to core-and-beam bi-Maxwellian fits of the Wind proton VDFs. We find several significant differences in the plasma waves obtained for the extracted data and bi-Maxwellian fits, including a strong dependence of the growth/damping rate on the shape of the VDF. By applying the quasilinear diffusion operator to these VDFs, we pinpoint resonantly interacting regions in velocity space where differences in VDF structure significantly affect the wave growth and damping rates. This demonstration of the sensitive dependence of Alfvén mode behavior on VDF structure may explain why the Alfvén ion-cyclotron instability thresholds predicted by linear theory for bi-Maxwellian models of solar wind proton background VDFs do not entirely constrain spacecraft observations of solar wind proton VDFs, such as those made by the Wind spacecraft.

Annual Variations in the Near-Earth Solar Wind

Earth’s orbit and rotation produces systematic variations in geomagnetic activity, most notably via the changing orientation of the dayside magnetospheric magnetic field with respect to the heliospheric magnetic field (HMF). Aside from these geometric effects, it is generally assumed that the solar wind in near-Earth is uniformly sampled. But systematic changes in the intrinsic solar wind conditions in near-Earth space could arise due to the annual variations in Earth heliocentric distance and heliographic latitude. In this study, we use 24 years of Advanced Composition Explorer data to investigate the annual variations in the scalar properties of the solar wind, namely the solar wind proton density, the radial solar wind speed and the HMF intensity. All parameters do show some degree of systematic annual variation, with amplitudes of around 10 to 20%. For HMF intensity, the variation is in phase with the Earth’s heliocentric distance variation, and scaling observations for distance largely explains the observed variation. For proton density and solar wind speed, however, the phase of the annual variation is inconsistent with Earth’s heliocentric distance. Instead, we attribute the variations in speed and density to Earth’s heliographic latitude variation and systematic sampling of higher speed solar wind at higher latitudes. Indeed, these annual variations are most strongly ordered at solar minimum. Conversely, combining scalar solar wind parameters to produce estimates of dynamic pressure and potential power input to the magnetosphere results in solar maximum exhibiting a greater annual variation, with an amplitude of around 40%. This suggests Earth’s position in the heliosphere makes a significant contribution to annual variations in space weather, in addition to the already well-studied geometric effects.

Energetic Neutral Atom (ENA) Emission Characteristics at the Moon and Mercury From 3D Regolith Simulations of Solar Wind Reflection

AbstractThe reflection of solar wind protons as energetic neutral atoms (ENAs) from the lunar surface has regularly been used to study the plasma‐surface interaction at the Moon. However, there still exists a fundamental lack of knowledge of the scattering process. ENA emission from the surface is expected to similarly occur at Mercury and will be studied by BepiColombo. Understanding this solar wind backscattering will allow studies of both Mercury's plasma environment as well as properties of the hermean surface itself. Here, we expand on previous simulation studies of the solar‐wind‐regolith interaction with 3D grains in SDTrimSP‐3D to compare the predicted scattering energies and angles to ENA measurements from the Moon by the Chandrayaan‐1 and IBEX missions. The simulations reproduce a backward emission toward the Sun, which can be connected to the geometry of the regolith grain stacking. In contrast, the ENA energy distribution and its Maxwellian shape is mostly connected to the solar wind velocity. Our simulations also correctly describe a lunar ENA albedo between 10% and 20% and support its decrease with solar wind velocity. We further expand our studies to illustrate how BepiColombo will be able to observe ENAs at Mercury using hybrid simulations of Mercury's magnetosphere as an input for the complex surface precipitation patterns. We demonstrate that the variable ion precipitation will directly influence ENA emission from the surface. The orbits of BepiColombo's Mercury Planetary Orbiter and Mercury Magnetospheric Orbiter/Mio spacecraft are shown to be suitable to observe ENA emission patterns both on a local and a global scale.

New insights on the behaviour of solar wind protons and alphas in the stream interaction region in solar cycle 23 and 24

New insights on the behaviour of solar wind protons and alphas in the stream interaction region in solar cycle 23 and 24

Kinetic study of ion acoustic waves in Venusian ionosphere

Kinetic dispersion of the ion acoustic waves has been explored for an unmagnetized five component plasma system comprising of Venusian protons, Venusian oxygen ions, Venusian electrons, solar wind protons, and kappa electrons. The solar wind protons and electrons are assumed to be streaming along the ambient magnetic field. The plasma parameters for this study have been obtained from Lundin et al. [Icarus 215(2), 751–758 (2011)] for the dawn dusk meridian of Venus Express with the data from the ASPERA-4 ion mass analyzer. Our analysis revealed that two modes, viz., ion acoustic mode and beam driven mode, are excited for the considered plasma parameters. The ion acoustic mode exists due to the Venusian ions, and its growth rate is influenced by the solar wind beam electrons. The beam driven mode's existence and its growth rate depend on the solar wind beam protons. We conjecture that the ion acoustic mode and the beam driven mode could be useful in explaining the electrostatic noise in the Venusian ionosphere in the range of several hundreds Hz to 1 kHz and several tens kHz, respectively.

Quantifying the Energy Budget in the Solar Wind from 13.3 to 100 Solar Radii

A variety of energy sources, ranging from dynamic processes, such as magnetic reconnection and waves, to quasi-steady terms, such as plasma pressure, may contribute to the acceleration of the solar wind. We utilize a combination of charged particle and magnetic field observations from the Parker Solar Probe (PSP) to attempt to quantify the steady-state contribution of the proton pressure, the electric potential, and the wave energy to the solar wind proton acceleration observed by PSP between 13.3 and ∼100 solar radii (R ☉). The proton pressure provides a natural kinematic driver of the outflow. The ambipolar electric potential acts to couple the electron pressure to the protons, providing another definite proton acceleration term. Fluctuations and waves, while inherently dynamic, can act as an additional effective steady-state pressure term. To analyze the contributions of these terms, we utilize radial binning of single-point PSP measurements, as well as repeated crossings of the same stream at different distances on individual PSP orbits (i.e., fast radial scans). In agreement with previous work, we find that the electric potential contains sufficient energy to fully explain the acceleration of the slower wind streams. On the other hand, we find that the wave pressure plays an increasingly important role in the faster wind streams. The combination of these terms can explain the continuing acceleration of both slow and fast wind streams beyond 13.3 R ☉.

Slow solar wind modeling of the Metis/Solar Orbiter – Parker Solar Probe quadrature

In January 2021, Metis/SolO and PSP formed a quadrature from which the slow solar wind was able to be measured from the extended solar corona (3.5 – 6.3 R⊙) to the very inner heliosphere (23.2 R⊙). Metis/SolO remotely measured the coronal solar wind, finding a speed of 96 – 201 kms−1, and PSP measured the solar wind in situ, finding a speed of 219.34 kms−1. Similarly, the normalized cross-helicity and the normalized residual energy measured by PSP are 0.96 and -0.07. In this manuscript, we study the evolution of the proton entropy and the turbulence cascade rate of the outward Elsässer energy during this quadrature. We also study the relationship between solar wind speed, density and temperature, and their relationship with the turbulence energy, the turbulence cascade rate, and the solar wind proton entropy. We compare the theoretical results with the observed results measured by Metis/SolO and PSP.

Electrostatic Solitary Waves in the Venusian Ionosphere Pervaded by the Solar Wind: A Theoretical Perspective

Electrostatic solitary waves (ESWs) in the Venusian ionosphere that are impinged by the solar wind are investigated using a homogeneous, collisionless, and magnetized multicomponent plasma consisting of Venusian H+ and O+ ions, Maxwellian Venusian electrons and streaming solar wind protons, and suprathermal electrons following κ − distribution. The model supports the propagation of positive potential slow O+ and H+ ion-acoustic solitons. The evolution and properties of the solitons occurring in two sectors, viz., dawn-dusk and noon-midnight sector of the Venus ionosphere at an altitude of (200–2000) km, are studied. The theoretical model predicts positive potential solitons with amplitude ∼(0.067–56) mV, width ∼(1.7–53.21) m, and velocity ∼(1.48–8.33) km s−1. The bipolar soliton electric field has amplitude ∼(0.03–27.67) mV m−1 with time duration ∼(0.34–22) ms. These bipolar electric field pulses when Fourier transformed to the frequency domain occur as a broadband electrostatic noise, with frequency varying in the range of ∼9.78 Hz–8.77 kHz. Our results can explain the observed electrostatic waves in the frequency range of 100 Hz–5.4 kHz in the Venus ionosphere by the Pioneer Venus Orbiter mission. The model can also be relevant in explaining the recent observation of ESWs in the Venus magnetosheath by the Solar Orbiter during its first gravity assist maneuver of Venus.

Interaction of Interplanetary Shocks With the Moon: Hybrid Simulations and ARTEMIS Observations

AbstractWhen in the solar wind, the Moon is exposed to its embedded discontinuities such as interplanetary (IP) shocks. In this study we utilize 3‐D electromagnetic hybrid (kinetic ions, fluid electrons) simulations and observations of two events, by THEMIS/ARTEMIS spacecraft to understand the interaction of IP shocks with the Moon. Simulation parameters are based on one of the ARTEMIS events and include the presence of supra‐thermal ions in the solar wind. The results show that the absorption of the cold solar wind protons on the Moon's dayside leads to a density hole in the shock front behind the Moon. This density hole refills and recovers as the IP shock moves down the lunar tail. Penetration of the interplanetary magnetic field through Moon's body leads to the survival of a steepened magnetic structure associated with the shock surface. However, due to the lower pressures in the lunar tail the structure broadens in space and is expected to steepen when the IP shock is beyond the tail. The IP shock is also found to interact with the energetic ions in the lunar tail resulting in their acceleration to higher energies. Comparing the densities measured by ARTEMIS in the solar wind and lunar tail shows the absence of a shock front in density in the tail consistent with the formation of a density hole at the shock front. Similarly, comparing the magnetic field profiles in the solar wind and lunar tail shows the expected broadening of the magnetic field shock front.