Saturn’s magnetosphere is not the clean, empty vacuum that early space-age scientists might have imagined. It is a crowded arena of electrons, ions, and microscopic dust particles, many of them shed from the planet’s vast E ring and from the geysers of Enceladus. What makes this environment truly remarkable is that the dust grains themselves can carry electric charges of either sign. Depending on their size, composition, and exposure to sunlight and plasma currents, some grains accumulate negative charge while others can become positively charged. This bipolar dust population, coexisting with ordinary plasma species, creates a medium whose wave behavior is far richer than that of a conventional electron-ion plasma.

The new model focuses specifically on the region between Dione and Rhea, a stretch of the middle magnetosphere where Cassini-era measurements have documented complex plasma populations. The authors construct a five-component fluid description: two distinct dust fluids representing positively and negatively charged grains, cold Maxwellian electrons and ions, and a population of hot suprathermal electrons. The suprathermal electrons are particularly important because spacecraft observations have repeatedly shown that Saturn’s magnetospheric electron distributions deviate from thermal equilibrium, possessing energetic tails that a simple Maxwellian cannot capture. By treating these hot electrons with a kappa distribution, the model stays faithful to what instruments such as Cassini’s Radio and Plasma Wave Science experiment actually measured.

The mathematical heart of the study is the reductive perturbation method, a classical technique in nonlinear plasma physics that has been used since the 1970s to tame the full complexity of fluid equations. The idea is elegant: rather than solving the entire nonlinear system at once, one assumes that the wave amplitude is small and stretches the space and time coordinates in a carefully chosen way. The original equations are then expanded order by order in a small parameter, and at the lowest significant order one recovers a linear wave mode, in this case the dust-acoustic wave, a compressive oscillation in which the massive, charged dust grains move collectively while the much lighter electrons and ions simply respond to restore charge balance.

Carrying the expansion to the next order yields a Korteweg-de Vries equation, one of the most celebrated equations in all of nonlinear science. The KdV equation famously describes systems in which nonlinearity, which steepens waves, and dispersion, which spreads them out, strike an exact balance. When these two effects cancel each other, the result is a solitary wave: a localized pulse that propagates without changing shape, behaving almost like a particle. First observed as a shallow water canal wave in the nineteenth century, solitons have since been identified in optical fibers, ocean currents, and, as this work shows, potentially in the dusty plasma surrounding Saturn.

One of the study’s most striking findings is the character of the solitary solutions. Unlike many laboratory dusty plasma experiments, which typically produce compressive solitons, the Saturnian configuration analyzed here supports only rarefactive solitary waves, localized depletions in the electrostatic potential rather than humps. The authors show that this behavior is governed decisively by the equilibrium negative dust charge. Because the negatively charged grains dominate the charge balance in this region, their abundance and charge state control the sign and magnitude of the nonlinearity coefficient in the KdV equation, and therefore determine both the polarity and the amplitude of the resulting solitary structures. Small changes in the dust charge parameter translate into significant changes in the wave profile, making the equilibrium dust charge a powerful diagnostic knob for the entire system.

The numerical analysis of the stationary solutions yields concrete, testable predictions. The electrostatic potential amplitudes of these solitary structures reach the order of a few millivolts, a small but physically meaningful value in the tenuous magnetospheric plasma. More intriguingly, the associated electric field exhibits a bipolar signature, swinging positive and negative in sequence, with a magnitude in the millivolt-per-meter range and a characteristic period of a few hundred seconds. This places the structures firmly in the ultra-low-frequency band, a regime where planetary magnetospheres are known to host a variety of wave phenomena but where dusty plasma effects have been difficult to isolate.

The bipolar electric field signature is particularly significant for observational purposes. Spacecraft such as Cassini, which spent thirteen years orbiting Saturn, routinely detected electrostatic solitary waves and broadband wave bursts in the planet’s inner magnetosphere. Previous studies documented such structures inside ten Saturn radii and near Enceladus, but the middle magnetosphere between Dione and Rhea remained comparatively underexplored from a theoretical standpoint. The new predictions give mission scientists and future data analysts a concrete template: if ultra-low-frequency, bipolar electric field pulses with periods of hundreds of seconds and amplitudes in the millivolt-per-meter range appear in archived or future data from this region, they could be the fingerprint of dust-acoustic solitons shaped by bipolar dust populations.

Beyond its immediate application to Saturn, the work speaks to a broader question in space plasma physics: how do charged dust grains reshape wave dynamics in planetary environments generally? Dusty plasmas are now recognized as ubiquitous, appearing in cometary comas, planetary rings, the lunar plasma environment, and interstellar clouds. The presence of dust introduces new low-frequency modes, modifies existing dispersion relations, and, as this study demonstrates, can even determine the polarity of nonlinear structures. The finding that a single parameter, the equilibrium dust charge, exerts decisive control over the nonlinear behavior suggests that similar diagnostics could be applied to other dusty magnetospheres, including those of Jupiter and the dusty plasmas near Mars and comets.

The authors are careful to ground their model in realistic Saturnian conditions, drawing on decades of spacecraft measurements of plasma densities, electron temperatures, and dust properties compiled from the Voyager era through the Cassini mission. Their order-by-order perturbation analysis, presented in full mathematical detail, shows that the secondary electron emission current, while physically present, remains roughly an order of magnitude smaller than the photoemission current and three orders of magnitude smaller than the dominant electron and ion collection currents under the considered conditions, justifying its neglect without altering the main conclusions. This kind of careful bookkeeping strengthens the case that the predicted rarefactive solitons are not mathematical artifacts but genuine features of the Saturnian middle magnetosphere. As planetary scientists continue to mine the Cassini archive and contemplate future missions to the outer solar system, this theoretical framework provides a ready-made lens for recognizing the subtle, slow, and silent waves that charged dust may be sending through the realm of the ringed planet.

Subject of Research: Nonlinear dust-acoustic solitary waves in the bipolar dusty plasma of Saturn's middle magnetosphere

Article Title: A perturbative approach to investigate nonlinear coherent structures in complex bipolar Saturnian plasmas

Article References: Bezbaruah, T., & Karmakar, P. K. (2026). A perturbative approach to investigate nonlinear coherent structures in complex bipolar Saturnian plasmas. Astrophysics and Space Science, 371(9), Article 106. https://doi.org/10.1007/s10509-026-04637-6

Image Credits: AI Generated

DOI: 10.1007/s10509-026-04637-6

Keywords: Saturn, dusty plasma, dust-acoustic waves, solitary waves, Korteweg-de Vries equation, magnetosphere, suprathermal electrons, nonlinear waves, Dione, Rhea, planetary rings, plasma physics

The Carrington storm takes its name from the English astronomer Richard Carrington, who on 1 September 1859 sketched a singular, brilliant white-light flare on the Sun, an observation independently confirmed by Richard Hodgson. Within roughly seventeen hours, the coronal mass ejection associated with that flare slammed into Earth’s magnetosphere, igniting auroras visible as far south as the Caribbean and Hawaii and disrupting telegraph communications across the world. For decades, estimates of the storm’s intensity have rested heavily on a single near-equatorial record from the Colaba Observatory in Bombay, and scientists have repeatedly cautioned that the intensity estimates derived from that record carry substantial uncertainty. The new Trivandrum data provide an independent equatorial vantage point, located close to the magnetic equator, and thereby offer a rare second witness to one of the most extreme geomagnetic disturbances in recorded history.

The instrument at the heart of the discovery is the bifilar magnetometer, a device in which a magnet is suspended by two parallel wires so that changes in the horizontal component of Earth’s magnetic field cause measurable rotations of the suspended system. John Allan Broun, the Scottish physicist who directed the Trivandrum Observatory under the patronage of the Maharajah of Travancore, was among the pioneers of this technique and published foundational work on the instrument’s errors and corrections in the Proceedings of the Royal Society of Edinburgh in 1862. His meticulous handwritten logs, correspondence, and papers, catalogued under accession 10064/1-172 at the National Library of Scotland, record hourly readings of the horizontal magnetic intensity at Trivandrum through the critical weeks of 1859. By digitizing and calibrating these archival entries, the research team converted the historical readings into modern nanotesla units, producing a quantitative record of geomagnetic disturbance at a latitude where the equatorial electrojet, an intense east-west current flowing in the ionosphere above the magnetic equator, strongly shapes the local magnetic field.

The analysis yielded a striking and unexpected result: the disturbance recorded at Trivandrum on 29 August 1859 was nearly as severe as the celebrated Carrington storm itself. The team measured a depression in the horizontal magnetic field intensity of approximately 1,851 nanotesla for the 29 August event, compared with approximately 1,926 nanotesla for the 2 September superstorm. In other words, the pre-Carrington storm of late August was of comparable magnitude to the event that has long defined the upper limit of space weather severity. This finding positions the 29 August disturbance as a newly calibrated benchmark for coronal mass ejection activity preceding the main storm, and it suggests that the Sun was already in an extraordinarily eruptive state days before Carrington pointed his telescope at the sunspot region that produced the historic flare.

The chronology that emerges from the Trivandrum records is one of repeated blows rather than a single punch. The data reveal geomagnetic disturbances in response to likely coronal mass ejection impacts on 16 August and 29 August, preceding the main event of 2 September. After the superstorm, the magnetosphere did not immediately settle. The records show further magnetospheric disturbances on 5 September, 13 September, and 15 September, which the researchers attribute to weaker coronal mass ejection events and/or high-speed solar wind streams arriving in the aftermath of the main activity. This pattern of successive impacts is consistent with the modern understanding that active solar regions can launch multiple eruptions over days or weeks, each compressing Earth’s magnetic shield and driving currents through the magnetosphere and ionosphere.

A particularly telling signature in the archival data is the disruption of the normal diurnal periodicity of the horizontal magnetic field intensity on 16 August, 5 September, 13 September, and 15 September. Under quiet conditions, the horizontal field at a near-equatorial station follows a regular daily cycle driven largely by the equatorial electrojet, which strengthens during daylight hours as solar heating and tidal winds drive ionospheric currents. When interplanetary coronal mass ejections strike the magnetosphere, magnetic reconnection at the magnetopause, the process by which interplanetary magnetic field lines merge with terrestrial field lines and transfer solar wind energy into the magnetosphere, injects particles and energy into the ring current, a torus of charged particles encircling Earth. The intensification of this ring current weakens the horizontal magnetic field at low latitudes, producing the characteristic negative excursions in the H component that magnetometers record as storms. The Trivandrum logs show exactly this kind of breakdown of the regular daily rhythm, underscoring the magnetospheric fluctuations induced by successive interplanetary coronal mass ejection impacts.

The physical framework for interpreting these observations draws on decades of magnetospheric research. The concept of geomagnetic storms as ring current phenomena was formalized in the widely cited definition by Gonzalez and colleagues in the Journal of Geophysical Research in 1994, and subsequent studies of ring current energy input and decay by Kozyra and Liemohn, and of ring current ion energy density evolution based on Van Allen Probes measurements by Zhao and colleagues, have clarified how storm-time currents build and dissipate. Reconnection physics, from Dungey’s seminal 1961 model linking the interplanetary magnetic field to the auroral zones through modern analyses of magnetopause reconnection location, explains how the orientation of the arriving solar wind magnetic field governs the efficiency of energy transfer. The Trivandrum data, though recorded with nineteenth-century instrumentation, capture the integrated low-latitude response of this entire coupled system, which is precisely why independent equatorial records are so valuable to space weather historians.

The significance of the new work extends beyond historical curiosity. Modern technological civilization is far more vulnerable to extreme space weather than the telegraph age was: power grids, satellite constellations, GPS navigation, aviation routes over the poles, and radio communications all depend on conditions in near-Earth space. Estimates of the upper limits of geomagnetic storm intensity feed directly into engineering standards and risk assessments for critical infrastructure. Recent scholarship, including Love and colleagues’ 2024 reassessment in the Journal of Space Weather and Space Climate, has emphasized that the intensity of the 1859 Carrington storm remains uncertain, and that the field needs more independent records to constrain the true worst case. Every newly recovered historical dataset narrows that uncertainty. The Trivandrum measurements, by confirming a near-Carrington-class storm on 29 August 1859, demonstrate that storms approaching the most extreme known intensity may occur more frequently than single-event statistics suggest, a conclusion with direct implications for how often infrastructure planners should expect once-in-a-century events.

The recovery of the records is itself a story of scientific detective work and international collaboration. The researchers acknowledge the support of the National Library of Scotland in retrieving the handwritten notes of John Allan Broun, and the visit was facilitated by Dr. Anna Szolucha of Jagiellonian University under a project supported by the Aryabhatta Research Institute of Observational Sciences. The work was funded by the Indian National Science Academy through project HS/RC/233. The study builds on the team’s earlier research, published in Advances in Space Research in 2025, which first unraveled the detection of the Carrington storm from Trivandrum’s historical magnetic declination observations, and on their documentation of the observatory’s 162-year engagement with the magnetic equator. Together, these efforts are restoring the Travancore Observatory, one of the colonial-era scientific institutions of India, to its rightful place in the history of geomagnetism.

As space agencies deploy fleets of monitors at the Lagrange points and improve models of coronal mass ejection propagation, the 1859 event remains the touchstone against which all space weather preparedness is measured. The Trivandrum chronology now shows that the Carrington storm was a multi-phase affair, a two-week siege of the terrestrial magnetosphere bracketed by disturbances nearly as fierce as the main event. In the delicate ink strokes of a Victorian magnetometer log, preserved by chance in a Scottish archive, scientists have found a warning written across a century and a half: the Sun is capable of sustained, storm-after-storm violence, and Earth’s defenses are tested not once but repeatedly when the star enters its most eruptive moods.

Subject of Research: A multi-phase reconstruction of the 1859 Carrington geomagnetic storm using archival bifilar magnetometer records from Trivandrum Observatory

Article Title: A multi-phase chronology of the Carrington storm of 1859 based on archival records of Trivandrum observatory

Article References: Jayakrishnan, R., C K, F., & A, A. (2026). A multi-phase chronology of the Carrington storm of 1859 based on archival records of Trivandrum observatory. Astrophysics and Space Science, 371(9), Article 105. https://doi.org/10.1007/s10509-026-04634-9

Image Credits: AI Generated

DOI: 10.1007/s10509-026-04634-9

Keywords: Carrington storm, space weather, geomagnetic storm, Trivandrum observatory, bifilar magnetometer, coronal mass ejection, magnetosphere, ring current, solar flare, magnetic equator, history of science, John Allan Broun

Gravitational wave detection has transformed our understanding of the universe since the historic first detection by LIGO in 2015. Ground-based instruments like LIGO, Virgo and KAGRA listen for waves in the 10 Hz to 1 kHz band, catching the final violent moments of neutron star and stellar-mass black hole mergers. But an entire cosmos of quieter, slower phenomena—supermassive black hole binaries spiraling together over millennia, extreme mass ratio inspirals of compact objects plunging into giant black holes—emits signals in the low-frequency band from 0.1 millihertz to 1 hertz. These can only be heard from space, which is why the European Space Agency’s LISA, China’s Taiji and TianQin programs, and Japan’s DECIGO concept are all racing to build giant laser interferometers among the stars.

TianQin, proposed in 2014 by Chinese scientists, is one of the most ambitious of these efforts. Its three satellites, deployed in a 100,000-kilometer geocentric orbit with an orbital period of roughly 3.9 days, will use laser interferometry to measure pico-scale changes in the distances between freely floating test masses housed within the spacecraft. For the mission to succeed, the acceleration noise on those test masses must be suppressed to a staggering level—on the order of 1×10⁻¹⁵ m s⁻² Hz⁻¹/² in the mission’s sensitive frequency band. Anything above that threshold risks drowning out the faint chirps of gravitational waves reaching Earth from across the universe.

Herein lies the problem. The TianQin orbit, unlike that of LISA, spends most of its time inside Earth’s magnetosphere, a region where the magnetic environment is far more complex and turbulent than interplanetary space. Meanwhile, energetic particles from galactic cosmic rays and the solar wind continually bombard the spacecraft, penetrating its outer walls and striking the test masses. This bombardment leaves the test masses electrically charged. A charged object moving through a magnetic field experiences the Lorentz force, and the residual magnetic moment that inevitably survives even the most careful manufacturing couples with the local field and its gradients to produce additional forces. All of these effects translate directly into acceleration noise that threatens the mission’s ability to sense gravitational waves.

In a study published in Results in Physics, Cheng-Long Yu of Sun Yat-sen University and colleagues have performed a comprehensive numerical simulation of this magnetic noise environment, modeling the full 100,000-kilometer TianQin orbit using the TA16 empirical magnetospheric model developed by Nikolai Tsyganenko and Valeria Andreeva. The TA16 model is a sophisticated piece of geophysical engineering: it divides the magnetospheric magnetic field into toroidal and polar components, each expanded as a weighted sum of radial basis functions calibrated against decades of satellite data from missions including Geotail, Polar, Cluster, THEMIS/ARTEMIS and the Van Allen Probes. The model responds to solar wind dynamic pressure, the interplanetary magnetic field, and geomagnetic perturbation indices, allowing researchers to reconstruct the magnetic field along the TianQin orbit at any point in history.

The team’s approach was unusually thorough. Rather than modeling a snapshot of space conditions, they computed magnetic field data continuously from 1997 to 2021—a span covering roughly two complete 11-year solar cycles. Because the TA16 model’s native time resolution of five minutes is insufficient to capture the spectrum up to 0.1 hertz, the researchers used cubic spline interpolation to refine the data to one-second resolution, then calculated the magnetic field gradient using a fourth-order central difference method enhanced with Richardson extrapolation, pushing the numerical accuracy to sixth order. The result was a high-fidelity map of both the magnetic field and its spatial gradients along the orbit, expressed as power spectral densities across frequencies from 1×10⁻⁵ hertz to 0.1 hertz.

From these spectra, the researchers calculated the acceleration noise from three distinct mechanisms: the force of the magnetic field acting on the residual magnetic moment of the test mass, the force of magnetic field gradients on both residual and induced magnetic moments, and the Lorentz force acting on the charged test mass. Using realistic parameters drawn largely from the LISA Pathfinder mission—including a test mass of 1.928 kilograms, a residual magnetic moment of about 0.14 nanoampere square meters along the sensitive axis, a magnetic susceptibility of −3.3723×10⁻⁵, and a test mass charge of 10⁻¹² coulombs—they found that the highest total acceleration noise remained below 1×10⁻¹⁴ m s⁻² Hz⁻¹/² at low frequencies and below 1×10⁻¹⁶ m s⁻² Hz⁻¹/² in the mid- and high-frequency bands. Crucially, when they defined a ratio β comparing the calculated noise to TianQin’s strict requirement curve, the maximum value of β across all analyzed cases was 0.7591—meaning the noise never exceeded the mission’s noise budget.

The findings also carry practical implications for instrument design. The calculations revealed that the noise from magnetic field gradients coupling to the test mass’s magnetization is negligible compared to the contributions from the magnetic field strength itself and the Lorentz force. Moreover, when the team examined the effect of solar activity—comparing the stormy solar maximum from June 1998 to June 2003 with the quiet minimum from January 2006 to December 2010—they found that noise levels rose noticeably during the active period, reflecting how enhanced solar wind pressure compresses the sun-facing side of Earth’s magnetosphere and strengthens the local magnetic field. In fact, during periods of extreme solar activity, the noise exceeds the current design budget, which means TianQin’s magnetic shielding must be improved beyond current assumptions.

The parametric sweeps performed by the team quantified these safety margins precisely. The noise ratio remains below one, indicating compliance with the mission requirement, as long as the test mass charge stays below about 7.47×10⁻¹² coulombs—a condition that can be maintained through established charge management techniques such as ultraviolet light discharge systems, which have already been demonstrated on LISA Pathfinder. Similarly, the magnetic shielding leakage coefficient, conservatively assumed to be 0.1 in the baseline calculations, must stay below 0.579; if shielding degrades beyond that threshold, additional magnetic shielding materials would be required to protect the test mass. Encouragingly, the team also showed that even if the magnetic moment within the test mass is not perfectly uniform—a possibility they probed by dividing the mass into a billion randomly magnetized elements—the resulting noise contribution is on the order of 10⁻²² m s⁻² Hz⁻¹/², far too small to matter.

For a mission that demands the most silent environment humanity has ever engineered, these results are deeply reassuring. TianQin’s constellation will fly through one of the more challenging magnetic environments imaginable for a gravitational wave observatory, buffeted by solar storms, cosmic rays and the roiling currents of Earth’s magnetosphere. Yet the new analysis shows that with realistic parameters and adequate shielding, the magnetic noise can be kept safely within budget under normal conditions—and that the residual risks are identifiable, quantifiable and correctable. As space-based gravitational wave astronomy moves from blueprint to hardware, studies like this one provide the confidence that when the first low-frequency gravitational wave whispers arrive at TianQin’s test masses, the universe’s signal will not be lost in the hum of the magnetosphere.

Subject of Research: Simulation of magnetic-field-induced acceleration noise on the test masses of China's TianQin space gravitational wave detector

Article Title: Integrated effect of the cosmic space magnetic field on the acceleration noise of the TQ gravitational wave detection program

Article References: Yu, C.-L., Yan, J., Ji, L., Shi, W.-K., Zhang, Y., Lau, Y. K., & Huo, H.-Q. (2026). Integrated effect of the cosmic space magnetic field on the acceleration noise of the TQ gravitational wave detection program. Results in Physics, Article 108755. https://doi.org/10.1016/j.rinp.2026.108755

Image Credits: AI Generated

DOI: 10.1016/j.rinp.2026.108755

Keywords: gravitational waves, TianQin, magnetic field noise, space-based interferometry, Lorentz force, test mass, magnetosphere, solar activity, Tsyganenko model, acceleration noise, charge management, magnetic shielding

Subpulse drifting is one of the most striking phenomena in radio pulsar astronomy. Rather than emitting a perfectly repeating pulse, many pulsars produce subpulses, discrete blobs of radio emission within the broader pulse profile, that shift systematically in spin longitude from one rotation to the next. When these subpulses are tracked across successive pulses, they appear to march, or drift, through the pulse window in a quasi-regular pattern. The prevailing interpretation, rooted in the polar-cap spark model developed by Michael Ruderman and Peter Sutherland in 1975 and refined by Jaroslaw Gil and colleagues, is that the drifting reflects a carousel of spark-associated plasma columns rotating around the star’s magnetic pole. As these sparks circulate, the radio beams they generate sweep across our line of sight at slightly different longitudes each rotation, producing the characteristic drifting pattern. The rate at which subpulses drift, and the spacing between drift bands, encode information about the geometry of the emission region and the conditions in the pulsar’s inner acceleration region.

For PSR J1059-5742, the new analysis shows that this carousel behavior dominates the star’s output to a remarkable degree. The team found that 96.2 percent of the pulsar’s radio radiation occurs in the drifting state, with subpulses drifting toward earlier longitudes across successive rotations. Quantitatively, the drift pattern is characterized by a vertical drift period of P3 equal to 5.7 plus or minus 0.5 times the pulsar’s spin period, meaning a full drift cycle repeats roughly every six rotations, and a horizontal spacing of P2 equal to 5.6 plus or minus 0.7 degrees of spin longitude. These parameters place J1059-5742 among the growing catalog of pulsars whose drifting behavior can be precisely characterized, joining objects studied in surveys such as the Meterwavelength Single-pulse Polarimetric Emission Survey and the Thousand-Pulsar-Array programme on MeerKAT, which have catalogued subpulse modulation across more than a thousand pulsars.

But the real surprise came from the small fraction of time when the drifting stopped. The Parkes single-pulse observation revealed two distinct non-drifting state events, each lasting 107 and 108 consecutive rotation periods respectively. Together these two episodes accounted for just 3.8 percent of the total observing time, making them genuinely rare interruptions to an otherwise relentlessly periodic pattern. During these events, the subpulses did not march across the pulse window at all; instead, the emission remained anchored at fixed longitudes, rotation after rotation, for stretches lasting on the order of a hundred spins. Given that pulsars rotate with clockwork regularity, often hundreds of times per second or, in slower cases, roughly once per second, a hundred-rotation interval represents a brief but unmistakable episode of magnetospheric reorganization.

What makes these events even more intriguing is that the emission did not merely stop drifting; it changed character. Single pulses recorded during the non-drifting state were brighter than those in the drifting state, and the peak of the average pulse profile for the non-drifting state led that of the drifting state by about 2 degrees in spin longitude. That longitudinal shift is small in absolute terms, a tiny sliver of the pulse phase, but it is systematic and measurable, and it signals that something physical changed in the emission geometry. The authors interpret this as evidence that changes may have occurred in the polar-cap discharge process or in the magnetospheric configuration itself. In other words, the carousel of sparks that normally produces the drifting pattern may have temporarily halted its rotation, or the entire discharge mechanism underpinning the radio emission may have switched to a different operating mode.

This interpretation connects J1059-5742 to a broader family of pulsar state-switching phenomena. Pulsars are known to alternate between emission modes, a behavior called mode changing, in which the average pulse profile abruptly switches between two or more stable shapes. Some pulsars also exhibit nulling, in which the radio emission ceases entirely for stretches of time, a phenomenon first reported by Donald Backer in 1970. More dramatically, the landmark 2010 study by Andrew Lyne and colleagues published in Science showed that switched magnetospheric states in several pulsars correlate with changes in the rate at which the stars spin down, implying that the magnetosphere and the star’s rotational evolution are intimately linked. The occasional cessation of drifting in J1059-5742 may represent a milder cousin of these phenomena, a mode change confined to the polar-cap discharge rather than a global magnetospheric reconfiguration.

The energy analysis of individual pulses added another layer to the picture. By examining the energies of all detected single pulses, the team found that the pulsar may have a nulling fraction of 13 plus or minus 1 percent, meaning that roughly one pulse in eight simply fails to appear. Notably, the study found that pulse nulling occurs only within the drifting state. The non-drifting state, despite its rarity, never nulls; when the pulsar enters this mode, it shines steadily and brightly. This asymmetry suggests that the two states represent genuinely different plasma conditions above the magnetic pole. In the drifting state, the spark discharge is active but intermittent, occasionally failing to produce detectable emission. In the non-drifting state, the discharge appears to operate in a more stable, luminous configuration, one in which the emission is both stronger and more reliable.

The theoretical stakes here are considerable. In the Ruderman-Sutherland framework, the inner acceleration region above the pulsar’s polar cap sustains a large electric field that tears charges from the stellar surface in discrete spark discharges. The pattern of these sparks, their number, their spacing, and their circulation speed around the magnetic axis, determines the observed drifting parameters. A halt in the carousel motion implies that the conditions governing the spark plasma’s E-cross-B drift, the electric and magnetic field interplay that carries the plasma around the pole, changed abruptly and then reverted. Alternatively, the magnetospheric configuration, the large-scale geometry of field lines and currents threading the light cylinder, may have shifted into a different stable state. Distinguishing between these possibilities requires comparing J1059-5742 with other pulsars that show similar behavior, such as PSR B1918+19, which displays drifting, moding, and nulling in combination, or PSR J1326-6700, a well-studied mode switcher.

For now, J1059-5742 joins a short list of pulsars in which drifting has been observed to stop entirely while emission continues. Each such object constrains the models in a different way, and the fact that the non-drifting state here is brighter and longitudinally offset rather than dimmer makes it a particularly clean diagnostic. The authors suggest that the occasional cessation of drifting in PSR J1059-5742 may thus reveal a mode change of the polar-cap discharge, a conclusion that, if confirmed by longer observations and multi-frequency follow-up, would add an important piece to the puzzle of how pulsars generate their beams. As single-pulse studies with instruments like Parkes, FAST, and MeerKAT accumulate ever larger samples, rare events like these two hundred-rotation episodes are shifting from curiosities into quantitative tests of the physics of neutron-star magnetospheres, reminding astronomers that even the most regular clocks in the universe occasionally change their tune.

Subject of Research: A single-pulse radio study of the subpulse-drifting pulsar PSR J1059-5742 revealing a rare brighter non-drifting emission state and pulse nulling.

Article Title: A rare non-drifting state in the subpulse drifting pulsar J1059-5742

Article References: Shang, L. H., & Ma, H. X. (2026). A rare non-drifting state in the subpulse drifting pulsar J1059-5742. Astrophysics and Space Science, 371(9), Article 108. https://doi.org/10.1007/s10509-026-04640-x

Image Credits: AI Generated

DOI: 10.1007/s10509-026-04640-x

Keywords: pulsar, PSR J1059-5742, subpulse drifting, nulling, Parkes telescope, magnetosphere, polar cap, mode changing, neutron star, radio emission, spark model, drift period

Solar energetic particles are among the most hazardous phenomena in the solar system. Unlike the steady wind of charged particles that flows constantly from the Sun, these episodic bursts can accelerate electrons, protons and heavier ions to a substantial fraction of the speed of light, driven by solar flares and the shock waves ahead of coronal mass ejections. At Earth, our planet’s strong dipolar magnetic field and thick atmosphere protect life and infrastructure from the worst of these events, although intense storms can still disrupt satellites, radio communications and power grids. At Mercury, the picture is radically different. The planet has only a weak magnetic field, roughly one percent of Earth’s surface field strength, and essentially no atmosphere, so any particle that penetrates its magnetosphere strikes bare rock directly.

BepiColombo is ideally placed to study this environment. The mission, a collaboration between the European Space Agency and the Japan Aerospace Exploration Agency, is en route to orbit Mercury with two scientific spacecraft, the ESA-led Mercury Planetary Orbiter and the JAXA-led Mercury Magnetospheric Orbiter. To reach its destination, the combined stack has been performing a series of planetary flybys, using gravity assists at Earth, Venus and Mercury to shed energy and adjust its trajectory. During one of its close Mercury flybys, the spacecraft passed through the planet’s magnetospheric environment at a time when the Sun was unusually active, turning a routine trajectory-correction maneuver into an unplanned natural experiment.

As the spacecraft approached the planet, its particle and field instruments registered the signature of a solar energetic particle event in progress. Outside the magnetosphere, the sensors recorded intense fluxes of energetic electrons arriving from the Sun. But as BepiColombo crossed into the region dominated by Mercury’s magnetic field, the picture changed dramatically. The fluxes of solar energetic particles dropped sharply, revealing that the planet’s magnetosphere, despite its modest strength, was acting as a genuine shield, carving out a protective cavity in the otherwise hostile radiation environment. This planetary shielding effect had been predicted by models, but direct in situ evidence of it during a major solar particle event had remained scarce.

The shielding, however, is not complete, and the second half of the story is arguably more important. The measurements showed that while most of the incoming particle population was deflected around the planet, a significant fraction of the energetic particles found their way to low altitudes and precipitated onto the surface. Charged particles can enter a magnetosphere along open magnetic field lines that connect the polar and near-polar regions to interplanetary space, and they can also be scattered into loss orbits by wave-particle interactions and by the complex, asymmetric structure of Mercury’s field, which is offset from the planet’s center and strongly compressed on the sunward side.

This surface precipitation matters because it is a primary engine of space weathering on Mercury. When energetic protons and electrons strike the regolith, they sputter atoms out of surface minerals, break chemical bonds and alter the optical and chemical properties of the soil. Over geological time, this bombardment contributes to the darkening and modification of Mercury’s surface that scientists observe remotely. It also plays a role in generating the planet’s tenuous exosphere, the wispy envelope of atoms knocked loose from the ground, including species such as sodium, calcium and magnesium that telescopes have long observed glowing around the planet. Quantifying how many energetic particles actually reach the surface during a storm therefore places a direct constraint on how fast these processes operate.

The event also carries implications for understanding Mercury’s magnetosphere as a system. Unlike Earth’s relatively stable dipole, Mercury’s magnetosphere is small, comparable in scale to the planet itself, and highly dynamic. Its boundary, the magnetopause, sits close to the surface, and during periods of high solar wind pressure it can be pushed down until the planet’s dayside field lines are compressed against the crust. Reconnection between the interplanetary magnetic field and Mercury’s own field opens the magnetosphere to the solar wind, funneling energy and particles into the magnetotail and down to the surface. Observing the system during a solar energetic particle event, rather than during quiet conditions, reveals how this delicate balance responds to extreme forcing.

The BepiColombo observations provide a benchmark for testing magnetospheric models under exactly these extreme conditions. Simulations of Mercury’s space environment must reproduce both the shielding of the bulk of the solar energetic particle population and the leakage that allows a subset of particles to reach the surface. Getting this right is essential not only for Mercury science but also for the broader problem of planetary habitability and radiation environments. Mercury represents an end-member case: a planet with a magnetic field too weak to fully protect it and no atmosphere to absorb what gets through. Studying how much shielding such a field provides informs our understanding of exoplanets around active stars, where close-orbiting rocky worlds may face relentless particle bombardment.

There is also a practical dimension for the mission itself. BepiColombo’s cruise phase, with its series of flybys, exposes the spacecraft to the full spectrum of solar activity near the innermost planet, and the solar energetic particle event during the flyby subjected the spacecraft’s electronics to elevated radiation doses. Understanding the particle environment around Mercury is therefore a matter of operational survival as well as science. When the two orbiters separate and begin their science phase in orbit around the planet, they will encounter these storms regularly, and the flyby measurements offer an early, detailed preview of what they will face.

For now, the flyby observations stand as a vivid demonstration that even a small planetary magnetosphere leaves a measurable imprint on a solar storm, deflecting most of the incoming radiation while simultaneously channeling a damaging fraction of it down to the airless ground. As BepiColombo continues its journey and ultimately settles into orbit, scientists expect many more opportunities to watch Mercury’s magnetic shield at work, transforming a fleeting flyby measurement into a long-term record of how the innermost planet weathers the Sun’s worst moods.

Subject of Research: In situ observations of planetary shielding and surface precipitation of solar energetic particles in Mercury's magnetosphere during a BepiColombo close flyby.

Article Title: Planetary shielding and surface precipitation of solar energetic particles during BepiColombo’s close Mercury flyby

Article References: Kilpua, E. K. J., Vainio, R., Grande, M., Edwards, L., Laurenza, M., Esko, E., Lehtolainen, A., Palmroos, C., Gieseler, J., Oleynik, P., Ho, G., Lawrence, D. J., Liu, S. J., Massetti, S., Heyner, D., Pump, K., Sanchez-Cano, B., & Huovelin, J. (2026). Planetary shielding and surface precipitation of solar energetic particles during BepiColombo’s close Mercury flyby. Nature Astronomy. https://doi.org/10.1038/s41550-026-02914-6

Image Credits: AI Generated

DOI: 10.1038/s41550-026-02914-6

Keywords: BepiColombo, Mercury, solar energetic particles, magnetosphere, planetary shielding, space weathering, surface precipitation, solar storms, ESA, JAXA, exosphere, space weather

The review, led by Minna Palmroth of the University of Helsinki and her colleagues, traces the physics and numerics of the Vlasov equation, the mathematical backbone of collisionless plasma theory. Plasma, the collectively behaving soup of charged particles that makes up most of the visible matter in the universe, is described by a distribution function that encodes how many particles occupy each point in a six-dimensional phase space combining three dimensions of ordinary space and three dimensions of velocity. In collisionless space plasmas, where particles interact primarily through long-range electromagnetic forces rather than frequent collisions, this distribution function evolves according to the Vlasov equation, coupled self-consistently to Maxwell’s equations for the electric and magnetic fields. Because nearly every measurable plasma quantity can be derived from the distribution function, it is, in many respects, the very core of plasma physics.

Alternative simulation strategies make different compromises. Magnetohydrodynamics, the workhorse fluid approach, treats plasma as a single thermalised fluid and is computationally cheap, but it assumes a Maxwellian, single-temperature plasma that simply does not exist in space, where the absence of collisions leaves particle velocity distributions multi-temperature and non-Maxwellian. The popular particle-in-cell (PIC) method propagates vast numbers of simulated particles and reconstructs the distribution from their statistics, but the resulting distributions are noisy, which can obscure the delicate physical processes at stake. The hybrid-Vlasov approach, by contrast, solves the Vlasov equation directly for ions on a full six-dimensional grid while treating electrons as a massless, charge-neutralising fluid. The decisive advantage is that the distribution function is evolved as an entity, without statistical noise, allowing sharp gradients and subtle kinetic signatures to be resolved with confidence.

The price of that fidelity is staggering. A straightforward Eulerian discretisation of Earth’s entire magnetosphere out to the lunar orbit, resolving the solar wind ion inertial length in space and the solar wind thermal speed in velocity, would require roughly 10 to the power of 18 phase-space cells, corresponding to a minimum of four exbibytes of memory. Even the most powerful supercomputers on Earth cannot hold such a dataset. The trick that makes global Vlasov simulations feasible at all is sparsity: large portions of velocity space contain essentially no plasma, so the code stores and propagates the distribution only where its density exceeds a threshold, retaining buffer regions for accurate transport. Combined with adaptive mesh refinement in ordinary space, this strategy cuts the computational burden by many orders of magnitude, bringing realistic global simulations within reach of petascale machines.

That machinery culminates in Vlasiator, developed in Finland and first proposed in 2007 to the newly established European Research Council as a high-risk, high-gain venture. Vlasiator advances the ion distribution using Strang splitting, alternating a spatial translation step with an acceleration step driven by the Lorentz force, both handled by a semi-Lagrangian solver called SLICE-3D that remaps the distribution with high-order reconstruction. Magnetic fields are propagated with a divergence-free, upwind constrained transport scheme that preserves the crucial solenoidality of the magnetic field by construction. The code couples the magnetospheric domain to a height-integrated ionosphere model that maps field-aligned currents down to 100 kilometres altitude, solves for the ionospheric electric potential, and feeds the resulting convection back into the simulation. The source code is openly available on GitHub, and the review carefully documents verification against analytical wave-dispersion solutions and hybrid-PIC benchmarks.

The physics harvest has been remarkable. In the terrestrial foreshock, the region upstream of the bow shock where reflected ions stream back toward the Sun, Vlasiator has revealed how foreshock waves and transient structures such as cavitons and spontaneous hot flow anomalies erode and reform the bow shock, and how these waves can transmit through the shock itself into the magnetosheath and even into the magnetosphere, where they are observed as Pc3 pulsations. Simulations showed that magnetosheath high-speed jets, fast plasma bursts that can hammer the magnetopause, can be launched when steepened foreshock waves strike the shock like bullets. In the magnetotail, a 3D breakthrough published in Nature Geoscience in 2023 demonstrated for the first time that magnetic reconnection and ion-kinetic instabilities operate simultaneously during the explosive eruptions that release plasmoids, a paradigm-shifting result that fluid models could never deliver because it requires resolving small and large scales in the same simulated volume.

The consequences for space weather are tangible. Energetic particles from the Sun disturb radio communications at high latitudes, and sudden magnetic changes induce currents in pipelines, railways and power grids; in 2022, a moderately stormy day cost the Starlink company 38 satellites, and worst-case estimates of an extreme event run to enormous economic damage. By reproducing observed ion distribution functions, auroral proton precipitation fluxes and magnetopause reconnection signatures in striking agreement with in situ spacecraft data from missions such as MMS and DMSP, hybrid-Vlasov simulation provides the physically grounded foundation on which reliable geospace prediction must ultimately be built. The review notes that even coarse spatial resolutions, far from the ion gyroradius, still yield genuine kinetic physics, because the high-energy ions that dominate global dynamics have large gyroradii and are faithfully represented on such grids.

Technological evolution has been inseparable from the science. Vlasiator runs on three levels of parallelisation: domain decomposition across thousands of supercomputer tasks with dynamic load balancing through the Zoltan library, shared-memory threading within each node, and vectorised processing of velocity-space blocks on each core. The latest frontier is graphics processing units. Early CUDA-based experiments date back more than a decade, but the current semi-Lagrangian code has been ported through performance-portable frameworks that translate to CUDA or HIP on demand, and the authors report that careful kernel fusion and redesigned memory handling should soon permit exascale 3D-3V simulations of the full terrestrial magnetosphere. New computational ideas are also entering the field, from low-rank tensor-train representations of the distribution function, cousins of the matrix-product techniques now ubiquitous in machine learning, to quantum algorithms for the Vlasov equation and physics-informed neural networks that already reproduce Vlasov-Poisson solutions with a few percent error.

Perhaps most striking is how far the Vlasov frontier now extends beyond Earth. The review highlights growing astrophysical applications, including relativistic Vlasov solvers for black-hole accretion coronae where radiation fields and quantum-electrodynamic processes matter, Vlasov-Poisson descriptions of galactic stellar dynamics, and emerging targets ranging from Mercury’s small magnetosphere to Mars, comets and the plasma wakes of airless bodies. The unifying lesson, the authors emphasise, is that everything affects everything: scale coupling between microscopic kinetic processes and global dynamics is the engine of space plasma behaviour, and only a method that treats both honestly in one simulation can capture it. For scientists seeking to understand everything from tomorrow’s geomagnetic storm to the eruptions on the Sun and the jets of distant galaxies, the noiseless, distribution-resolving Vlasov approach, once dismissed as computationally impossible, has become an indispensable window onto the kinetic universe.

Subject of Research: Vlasov-based numerical methods for kinetic plasma simulations in space physics and astrophysics

Article Title: Vlasov methods in space physics and astrophysics

Article References: Palmroth, M., Ganse, U., Pfau-Kempf, Y., Battarbee, M., Alho, M., Nättilä, J., Zaitsev, I., Cozzani, G., Papadakis, K., Kotipalo, L., Zhou, H., Turc, L., Hoilijoki, S., Grandin, M., Pänkäläinen, L., Sandroos, A., & von Alfthan, S. (2025). Vlasov methods in space physics and astrophysics. Living Reviews in Computational Astrophysics, 11(1), Article 3. https://doi.org/10.1007/s41115-025-00024-0

Image Credits: AI Generated

DOI: 10.1007/s41115-025-00024-0

Keywords: Vlasov equation, space plasma, hybrid-Vlasov simulation, Vlasiator, magnetosphere, space weather, magnetic reconnection, collisionless shocks, high-performance computing, kinetic physics, astrophysical plasmas, GPU computing

Pulsars are the rapidly rotating, magnetized remnants of massive stars, sweeping radio beams across the Galaxy with clockwork regularity. Ever since their discovery in 1967, astronomers have known that the radio flux of a pulsar generally follows a power law as a function of frequency: the star is brighter at low radio frequencies and fades steadily toward higher ones. The steepness of that decline, encoded in the spectral index, is not a mere detail of bookkeeping. It carries information about the energies of the charged particles spiraling in the star’s magnetic field and about the radiation mechanism that converts their motion into the coherent radio waves we detect on Earth.

What makes the new work powerful is its treatment of time within the pulse. A pulsar’s average profile, built up from thousands of individual rotations, is not a single blob of emission. It often shows a central core component flanked by one or two cone-shaped components, a geometry long described by empirical classification schemes and by theoretical models such as inverse Compton scattering. Each component is thought to trace a different emission region at a different height above the magnetic pole. If the spectral index changes across the profile, it means the particle energy distribution changes with location in the magnetosphere. Phase-resolved spectroscopy therefore acts as a tomographic tool, slicing the emission beam into radial and longitudinal layers.

Ultra-wideband receivers make this kind of analysis far more efficient than it used to be. Instead of observing a pulsar separately at several discrete frequencies, Parkes’ wideband system captures a broad contiguous stretch of the radio spectrum simultaneously, allowing multi-band pulse profiles to be derived from a single observation with consistent calibration. Li and Huang exploited this capability to construct spectra at many pulse phases for each of their nine targets, tracing how the spectral index evolves from the leading edge of the profile, through its peaks, and out along the trailing edge.

The result is a taxonomy. When the authors inspected the shapes of the phase-resolved spectra, the curves fell naturally into four categories: W-type, normal V-type, inverted V-type, and irregular. W-type spectra show the spectral index dipping and rising again in a double-valley pattern across the profile, while V-type spectra display a single pronounced valley, either in the conventional sense or inverted. Irregular spectra, as the name suggests, resist simple description. The key insight is that the category is not random: it correlates with the morphology of the pulse profile itself.

Specifically, pulsars whose profiles show poorly resolved single or multiple emission peaks tend to produce the orderly W-type and V-type phase-resolved spectra. Sources with clearly distinct, well-separated multiple peaks, by contrast, generally yield irregular spectra. This pattern suggests that the smooth, overlapping components of a simple beam geometry produce systematic spectral behavior, whereas a profile composed of sharply separated components reflects a more complicated emission region where different radiation zones contribute in ways that do not blend into a tidy curve. In other words, the shape of the spectrum across the pulse is a fingerprint of the beam’s internal architecture.

To quantify these shapes, the authors fitted the phase-resolved spectra with single-valley and multi-valley functions, and found that all nine pulsars were well described by such models. That success matters because a reliable analytic fit is what allows the next step in the chain of inference: inverting the observed spectral behavior to recover the distribution of particle energies in the emission regions. Under radiation models such as inverse Compton scattering, the spectral index at a given pulse phase maps onto the energy of particles radiating from a particular location along the field lines. A well-fitted spectrum is therefore a proxy for a map of the magnetosphere’s particle population.

The study builds on a long observational tradition. Pulsar spectra were first characterized decades ago, and multi-frequency studies of individual objects, such as detailed analyses of PSR B1133+16 and PSR B0329+54, previously demonstrated how profile components evolve with frequency and how those changes constrain emission height and particle energy. Scatter broadening by the interstellar medium, which smears pulses in a frequency-dependent way, has been measured for large samples of pulsars at meter wavelengths and must be accounted for when comparing profiles across a wide band. The Parkes Observatory Pulsar Data Archive, maintained by the Australia Telescope National Facility, provides the archival backbone for much of this work, and the authors acknowledge the facility for providing the data used in the study.

Why should anyone beyond radio astronomy care? Pulsars are foundational tools across modern astrophysics. Their timing stability underpins tests of general relativity, the search for low-frequency gravitational waves through pulsar timing arrays, and even proposals for spacecraft navigation based on pulsar clocks. Yet the fundamental question of how a pulsar makes its radio beam, the coherent emission mechanism, remains only partially solved after more than half a century. Every constraint on where in the magnetosphere different frequencies are emitted, and how particle energies vary across the emission beam, narrows the space of viable theories. Phase-resolved spectra are among the sharpest such constraints available, because they connect spectral behavior directly to pulse geometry.

The Parkes results also carry practical implications for future surveys. Instruments such as the Square Kilometre Array and its precursors will detect vast numbers of pulsars across enormous bandwidths. Understanding how spectral index varies with pulse phase, and how that variation relates to profile class, will help astronomers interpret wideband measurements correctly, avoid biases in flux calibration and population studies, and select the best targets for precision timing. A catalog of spectral behavior by morphological type, even one beginning with nine objects, provides a template for scaling to thousands.

There is also a conceptual payoff. The finding that profile morphology predicts spectral morphology supports the picture in which a pulsar’s beam is organized into a core-plus-cone structure, with each component occupying a distinct region of the magnetosphere and hosting its own particle energy distribution. When components overlap in our line of sight, their spectral signatures blend into smooth W or V shapes; when they are distinct, the spectrum becomes irregular because the observer samples genuinely different emission zones in quick succession. The spectra, in effect, let researchers read the beam the way a geologist reads strata.

Li and Huang’s work, funded by the Guizhou Provincial Basic Research Program and science programs of the Xinjiang Uygur Autonomous Region, demonstrates how a single well-instrumented telescope can continue to extract new physics from familiar objects. The nine pulsars studied are ordinary radio pulsars by most standards, yet their phase-resolved spectra reveal structured, classifiable behavior that single-frequency observations would entirely miss. As wideband receivers spread to telescopes worldwide, the technique demonstrated here is likely to become standard practice, turning every calibrated pulsar observation into a spectrum-rich dataset and bringing the community closer to answering the oldest question in pulsar science: what, exactly, is shining inside the beam.

Subject of Research: Phase-resolved radio spectra and emission beam geometry of nine pulsars observed with the Parkes 64-m Radio Telescope

Article Title: Phase-resolved spectra of nine pulsars from ultra-wideband observations with the Parkes 64-m Radio Telescope

Article References: Li, Z., & Huang, A. (2026). Phase-resolved spectra of nine pulsars from ultra-wideband observations with the Parkes 64-m Radio Telescope. Astrophysics and Space Science, 371(9), Article 98. https://doi.org/10.1007/s10509-026-04630-z

Image Credits: AI Generated

DOI: 10.1007/s10509-026-04630-z

Keywords: pulsars, radio astronomy, Parkes Observatory, phase-resolved spectra, spectral index, pulse profiles, emission beam geometry, ultra-wideband receiver, magnetosphere, radiation mechanisms, Astrophysics and Space Science, Phase-resolved

The particle and X-ray detector known as SIXS, designed and built in Finland, recorded the entire encounter. Now the measurements have been published in the journal Nature Astronomy, offering scientists their first close-range view of how solar energetic particles penetrate a planetary magnetosphere as small and feeble as Mercury’s and precipitate onto the planet itself. The results matter far beyond one small, scorched world. They provide a natural laboratory for understanding what happens when the most violent space weather strikes a magnetosphere, knowledge that directly informs how researchers assess the threat such storms pose to Earth’s own technological infrastructure and atmosphere.

Emilia Kilpua, principal investigator of SIXS and professor of space physics at the University of Helsinki, describes the flyby as an extraordinary stroke of fortune. The spacecraft came much closer to the surface than it will on its final science orbit, and a significant particle eruption occurred on the Sun at precisely that moment. Observing a powerful solar energetic particle event from a few hundred kilometers above Mercury, with instruments built to measure exactly such particles, is the kind of convergence mission scientists rarely dare to hope for. Instead of a quiet pass used mainly for steering and instrument testing, the team obtained a dataset of genuine scientific value that would have been impossible to schedule.

What SIXS revealed is striking. During the eruption, a large number of highly energetic particles reached Mercury’s surface across a wide area, sweeping over the planet rather than being funneled into narrow polar cusps alone. Mercury lacks a substantial atmosphere, so there is no cushion of air to absorb or slow the incoming radiation the way Earth’s atmosphere shields the surface. The planet’s magnetic field, moreover, is far weaker than Earth’s, and its magnetosphere is correspondingly much smaller, leaving it far more exposed to the full force of solar outbursts. The measurement therefore captures, in real time and at close range, how effectively a planetary magnetic field can or cannot shield the ground beneath it when the Sun unleashes its energy.

The physics of what happens next on the surface is central to why these observations matter. When energetic electrons and protons slam into the regolith, the loose layer of dust and rock that coats Mercury, they knock atoms and molecules loose from the surface material in a process called sputtering. The same impacts also generate X-ray radiation, as particles decelerate and excite atoms in the crust. Both the ejected particles and the induced X-rays carry chemical fingerprints of the surface composition. By measuring them, scientists can deduce what Mercury’s crust is made of and how the planet has been processed by billions of years of exposure to solar radiation and micrometeorite impacts. In orbit, several BepiColombo instruments will study these processes together, cross-referencing particle measurements from SIXS with X-ray observations of the surface from the MIXS instrument, in which Finland is also heavily involved, to build a coherent picture of the planet’s chemistry and evolution.

Rami Vainio, co-principal investigator of SIXS and professor of space physics at the University of Turku, emphasizes the comparative dimension of the findings. Mercury’s magnetic field is weaker than Earth’s and its magnetosphere much smaller, but the conditions at Mercury during this event resemble what Earth would experience if a powerful solar storm compressed its own magnetosphere. Under extreme conditions, Earth’s protective magnetic bubble can be squeezed inward dramatically, allowing energetic particles to penetrate far deeper into near-Earth space than usual, threatening satellites, disrupting radio communications, and in the most severe cases depositing energy into the upper atmosphere. Mercury during this flyby effectively demonstrated the end-member scenario, a magnetosphere overwhelmed by particle flux, showing how such radiation propagates and where it ultimately lands.

This is why the SIXS observations are expected to feed directly into space weather research focused on our own planet. Understanding how destructive particle radiation penetrates a planetary near-space environment and atmosphere during the most powerful space storms helps researchers refine models of radiation dose at aviation altitudes, satellite orbital predictions, and the behavior of Earth’s magnetosphere under stress. The data will also be used in the Center of Excellence in Space Resilience, funded by the Research Council of Finland for the period 2026 to 2033, in which both Kilpua’s and Vainio’s research groups participate. That program brings together space physicists, engineers, and resilience researchers to anticipate and mitigate the effects of extreme space weather on society, a concern that has grown as solar cycle activity has intensified and as modern infrastructure has become ever more dependent on satellites and sensitive electronics.

The broader context of the measurement is the remarkable journey that made it possible. BepiColombo was launched in October 2018 and has spent eight years traveling to Mercury, using a series of flybys of Earth, Venus, and Mercury itself to shed the orbital energy that a direct trajectory to the Sun’s neighborhood would make prohibitively expensive in fuel. In early September, the mission reached a milestone: the two orbiters, the Mercury Planetary Orbiter built by ESA and the Mercury Magnetospheric Orbiter built by JAXA, separated from the transfer module that had carried them across the solar system. The orbiters will enter Mercury’s orbit in November and separate from each other in December, beginning the dedicated science phase in complementary orbits. In total, the two spacecraft carry sixteen scientific instruments, forming the most comprehensive payload ever sent to the smallest and least explored of the solar system’s eight planets.

Mercury holds a special place in that exploration history, and a humble one. It is the planet closest to the Sun, which makes it exceptionally strongly affected by solar activity, yet only three spacecraft have ever studied it at close range. Mariner 10 flew past three times in the mid-1970s, revealing a heavily cratered world with a surprising magnetic field. MESSENGER orbited the planet from 2011 to 2015, mapping its surface, confirming water ice in permanently shadowed polar craters, and revealing an offset magnetic field whose origin remains debated. BepiColombo is only the third mission to visit and Europe’s first, carrying the ambition not merely to revisit these discoveries but to understand the planet as a coupled system: its interior, its surface, its thin exosphere, and its magnetosphere, all responding continuously to the ferocious output of the nearby Sun.

The fourth flyby measurements now published add an unexpected early dividend to that mission. Long before the orbiters begin their systematic observations, SIXS has demonstrated how solar energetic particles behave at Mercury, how widely they spread across the surface, and what that means for both planetary science and space weather forecasting. The lucky timing of the September 2024 encounter has effectively given the BepiColombo team a preview of the storm-driven physics they will study in detail from orbit, and given space weather scientists on Earth a benchmark dataset drawn from the most sun-exposed planet in the solar system. As the Sun approaches and passes through the active phase of its cycle, the lessons learned from Mercury’s bombardment will resonate directly in models that protect satellites, astronauts, and power grids far from the scorched world where they were first measured.

Subject of Research: Solar energetic particle precipitation onto Mercury's surface measured by the BepiColombo spacecraft during its close flyby

Article Title: BepiColombo measured Mercury's particle bombardment up close

Article References: BepiColombo measured Mercury's particle bombardment up close. (n.d.). Original publication

Image Credits: AI Generated

DOI: Not provided

Keywords: BepiColombo, Mercury, solar energetic particles, space weather, magnetosphere, SIXS instrument, ESA, JAXA, solar storms, X-ray radiation, planetary shielding, Nature Astronomy