2026-03-09 | | Total: 10
We investigate the magnetic signature of oceanic circulation in Ganymede's subsurface ocean using kinematic induction modeling. Our approach couples zonal jet flows from rotating thermal convection simulations with magnetic field models incorporating Ganymede's internal dynamo and external contributions from Jupiter. We solve the induction equation in spherical geometry for deep-ocean (493 km) and shallow-ocean (287 km) scenarios with varying magnetic Reynolds numbers. Ocean flows generate a predominantly toroidal magnetic field through the omega-effect, with a weaker poloidal component pervading beyond the conductive ocean layer. For some, but not all, induction configurations, analysis of the time-averaged Lowes-Mauersberger spectra reveals that ocean-induced signals dominate at spherical harmonic degrees $\ell \geq 4$. Deep ocean scenarios with magnetic Reynolds numbers above unity produce surface magnetic signals up to 9 nT. Our results demonstrate that Ganymede's intrinsic magnetic field creates favorable conditions for detecting subsurface ocean dynamics, thus emphasizing the need for low-altitude orbits for the Juice probe.
This paper investigates an autonomous navigation method for spacecraft operating in the outer solar system, up to 250 AU from the Sun, using the parallactic shifts of nearby stars. These measurements enable estimation of the spacecraft trajectory while distant stars provide attitude information through conventional star-pattern matching. Stellar observation models are developed, accounting for delta light-time, parallax, and aberration effects. Navigation performance is assessed using two approaches: (1) a least-squares estimator using simultaneous multi-star measurements, and (2) a Kalman filter processing sequential single-star observations along deep-space trajectories. Monte Carlo simulations on trajectories representative of Voyager 1, Voyager 2, Pioneer 10, Pioneer 11, and New Horizons missions show sub-AU position accuracies at 250 AU, and velocity accuracies better than 0.00004 AU/day, under realistic spacecraft and instrumentation uncertainties. These values correspond to relative errors below 0.4% in position and velocity with respect to the reference trajectories. Although less precise than radiometric tracking, this performance can support navigation in the outer solar system without reliance on Earth. When ground-based navigation remains necessary, this approach can be employed during long cruising phases, lowering the number of ground contacts. The method additionally shows potential for future missions venturing farther from the Sun.
Sulfur is a redox active element that may have helped mediate an electron flow that kickstarted life and which presently is an essential element for all life on Earth. Despite current uncertainties in global sulfur fluxes, modeling sulfur's abiotic cycling through Earth's deep history is important for understanding the impact of a planet wide biosphere on sulfur geochemical cycling and availability and vice versa. We present here an open-source, dynamical box model for estimating global sulfur fluxes and concentrations among surface and deep Earth reservoirs over Earth history, allowing tracking and estimation of the sulfur distribution in planetary reservoirs over deep time in the absence of life. While the main model presented here does not take into account the abrupt evolution of redox-shunting biosynthetic pathways such as oxygenic photosynthesis, we also modeled the abiotic sulfur cycle before and after a Great Oxidation Event-like transition on Earth-like planets. Our results suggest a considerably distinct chemical makeup of sulfur content in marine sediments in the absence of life on an Earth-like planet, leading to a marine sediment sulfate content two orders of magnitude larger than on present-day Earth and a marine sediment sulfide content 4 orders of magnitude lower than on present day Earth, attributable to the lack of microbial sulfur metabolism. This model could be useful for understanding sulfur cycling on potentially habitable exoplanets.
The K-dwarf TOI-4504 hosts two giant planets in 2:1 mean-motion resonance, with orbital periods of 41.3 days (planet d) and 82.8 days (planet c). They exhibit among the largest known absolute transit-timing variations, with respective peak-to-node amplitudes up to 5 and 3 days. Newer TESS data show that the previously non-transiting planet d has now precessed into transiting, and we derive updated system parameters with significant discrepancies with the discovery paper. The revised parameters place planets d and c deep in the resonance and close to or in the fully-relaxed limit-cycle state, with the resonant and secular modes interfering nonlinearly to induce non-zero relaxed free eccentricities which precess at the same rate as the forced eccentricities and the longitude of conjunctions, in turn enabling precise measurement of the full eccentricities and apsidal angles. We discuss the predictions of linear theory and how it can be used to understand the true state of the system revealed by N-body integrations, and more generally why it is that the posteriors of systems more compact than 2:1 tend to suffer from significant eccentricity degeneracy. We show that the extraordinary dynamical states of the giant pairs orbiting TOI-4504 and the M-dwarf GJ 876 are remarkably similar, in spite of the significant difference in their host-star masses, and discuss the implications for damping timescales during the relatively gentle formation process of Type II migration.
We investigate the dynamical evolution of particles in the $β$ Pic system to determine likely formation pathways to the present-day observed exocomet populations. We aim to relate these results to similar studies recently carried out since the discovery of the inner planet $β$ Pic c. We simulate the $β$ Pic system using the non-symplectic adaptive N-body integrator IAS15 in REBOUND. We seed the system with over 100,000 mass-less test particles that evolve for 25 Myr, and adopt initial conditions and a particle distribution that closely matches similar simulations in recent literature. Using IAS15, REBOUND resolves close-encounters between test particles and the two gas giants in the system, which is crucial for understanding aspects of the dynamical evolution. Planet-disk interactions rapidly clear most of the system within 35 AU apart from a region within the orbit of $β$ Pic c, and a region between 20 and 25 AU. After 10 Myr, exocomets can be sourced continuously from these regions, as well as from the inner edge of the region beyond ~35 AU where particles are stable on longer timescales. From the region interior to $β$ Pic c, the exocomets are formed by excitation via mean-motion resonance with $β$ Pic c, obtaining a narrow distribution of radial velocities, consistent with spectroscopic observations. Particles initialized in the outer system may enter onto stargrazing orbits due to disruption by the two gas giants, causing a wider radial velocity distribution, and we propose that this population corresponds to a second dynamical family previously observed via spectroscopy. These particles typically undergo chaotic dynamical evolution for $10^2$ to $10^3$ years after passing the water sublimation limit at ~8 AU until reaching the sublimation distance of calcium near 0.4 AU, implying that the two families of exocomets may have different volatile contents.
High-resolution scattered-light imaging has revealed complex morphologies in protoplanetary and circumstellar disks. Measuring the vertical height of the scattering surface is key to understanding disk structure, evolution, and the properties of embedded dust. We develop a methodology for fitting elliptical shapes to scattered-light images of protoplanetary disks in order to extract vertical height profiles of the dust scattering surface across a large and morphologically diverse disk sample. The dataset consists of 92 near-infrared polarimetric images obtained with VLT/SPHERE. The aim is to identify trends in vertical structure across different disk morphologies and test for correlations with stellar mass, age, and disk dust mass, as well as to investigate the implications of the derived height profiles for the masses of potential embedded planets. We implement a structure extraction and ellipse fitting (SEEF) algorithm that uses edge detection and Gaussian fitting to locate disk structures. Ellipse fitting reveals spatial offsets between the ellipse centre and the stellar position, which are interpreted as vertical height assuming circular ring geometry. Disk inclination, position angle, and the aspect ratio h/r are also derived. The method yields vertical height measurements for 92 disks, showing profiles consistent with flared disk geometries. However, the full sample cannot be described by a single power-law relation. Subdivision by morphology shows no strong correlations for most disk classes, except for extended disks with outer radii larger than about 150 au, which exhibit a clear power-law flaring trend. The lack of strong correlations with other system properties suggests that either different morphologies exhibit distinct vertical structures or that additional physical factors influence disk flaring.
Imminent impactors are natural bodies discovered in space before impacting the Earth. They provide a rare opportunity to characterize individual near-Earth objects (NEOs) in great detail as asteroids in space, meteors in Earth's atmosphere and meteorites on the ground. The Vera C. Rubin Observatory's upcoming Legacy Survey of Space and Time (LSST) is expected to transform our understanding of the NEO population. In this work, we evaluate LSST's expected discovery performance for imminent impactors using $343$ meter-size objects previously recorded in NASA's CNEOS database as fireballs impacting Earth's atmosphere. We simulate pre-impact observations of these CNEOS impactors with the Sorcha survey simulator under LSST's default three-night discovery strategy and a one-night strategy for fast-moving objects that relies on matching aligned streaks in two exposures on the same night. We estimate that LSST will discover $\sim1-2$ meter-size and larger imminent impactors per year, representing $\sim4\%$ of all Earth impactors $\gtrsim1$ m in diameter and almost doubling the current discovery rate of imminent impactors. The median time of discovery and median time of first observation for impactors discovered in our simulations are $\sim1.57$ and $\sim3.06$ days before impact, respectively. The spatial distribution of the 11 previously discovered imminent impactors is biased towards the Northern Hemisphere, where the observatories that discovered them are located. We find a similar trend towards Southern Hemisphere impacts in our simulated LSST detections of the CNEOS impactors, suggesting Rubin will provide a powerful counterpart to existing asteroid surveys primarily located in the Northern Hemisphere.
We present a novel machine learning method that is capable of rapidly and accurately producing dust-continuum model images and spectral energy distributions from training sets created using a detailed radiative transfer code. We create a training set that encompasses the parameter space for protoplanetary discs, and then couple the trained machine learning method with a Bayesian optimisation algorithm. We then simultaneously fitted 1.3 mm ALMA ODISEA survey images of protostellar discs in rho Oph, and their spectral energy distributions, in order to determine fundamental discs parameters such as dust masses and radii. We find that good simultaneous fits may be found for the Class II objects in the survey, although the spectral fits are poorer for the Class I and flat spectrum sources. We find that the dust mass distributions of discs is broader and shallower than that predicted from 1.3 mm flux dust mass estimates, substantially increasing the numbers of objects with high-mass and low-mass discs. We show that this is due to a combination of optical depth and dust temperature effects, which are strongly related to the disc size and inclination constraints provided by the imaging fits. We show that there is a significant decrease in disc scale height and disc flaring when moving from the the Class I objects, to the flat spectrum sources, and the Class II discs.
Stellar rotation is a fundamental tracer of stellar magnetic evolution, age, and activity, with broad implications for Galactic archaeology and exoplanet characterization. The Transiting Exoplanet Survey Satellite (TESS) provides high-precision time-series photometry across the sky, enabling rotation measurements for an unprecedented number of stars. We present the TESS All-Sky Rotation Survey (TARS), an all-sky catalog of stellar variability periods for 944,056 stars with T < 16 and distances within 500 pc. We estimate that 94% of these periods are rotation periods. This catalog increases the number of rotation period measurements for stars with T < 16 within 100 pc by a factor of 2.1 and within 500 pc by 3.7. We also present a method to correct half-period aliases in TESS data and show that it reliably recovers periods as long as 25 days from a single TESS sector. TARS represents the largest homogeneous catalog of stellar rotation periods to date, providing a foundation for studies of stellar ages, young associations, and Galactic structure. We make the light curves used in our analysis available as a HLSP through MAST. Beyond the default TARS catalog, we provide code that allows users to generate rotation period catalogs with adjustable completeness and reliability thresholds. This code and all rotation period measurements are available through Zenodo.
Magnetohydrodynamic (MHD) seismology uses naturally occurring MHD waves to infer plasma properties that are otherwise hard to measure, especially magnetic field strength and topology, electric currents, fine structuring, transport coefficients, and energy release. Across the solar atmosphere, heliosphere, and planetary magnetospheres, multi-wavelength remote sensing and in-situ observations of waves provide powerful diagnostics that can address major open problems including chromospheric and coronal heating, flare and eruption physics, solar wind acceleration, and space weather impacts. This White Paper sets out the case for a coordinated UK programme that couples high precision observations with advanced theory and numerical modelling, modern time-frequency methods for non-stationary signals, and machine learning approaches for detection, classification, and parameter inference from rapidly growing multi-instrument datasets. It outlines priority needs such as robust mode identification, reliable density and temperature constraints, multi line-of-sight capability, and models that include partial ionisation and non-adiabatic/collisionless effects, alongside enabling instrumentation such as next-generation spectropolarimetry, integral field units, and radio facilities including the Square Kilometre Array. The paper highlights the UK's strong track record and infrastructure, and argues that sustained investment will amplify UK scientific return through international partnerships and mission involvement, delivering transformative plasma diagnostics and downstream benefits for space weather forecasting and related applications.