Search

Airless bodies are subjected to space-weathering effects that modify the first few microns of their surface. Therefore, understanding their impact on the optical properties of asteroids is key to the interpretation of their color variability and infrared reflectance observations. The recent Hayabusa2 sample return mission to asteroid Ryugu offers the first opportunity to study these effects, in the case of the most abundant spectral type among the main-asteroid belt, C-type objects. This study employs vibrational electron energy-loss spectroscopy in the transmission electron microscope to achieve the spatial resolution required to measure the distinct mid-infrared spectral signature of Ryugu's space-weathered surface. The comparison with the spectrum of the pristine underlying matrix reveals the loss of structural -OH and C-rich components in the space-weathered layers, providing direct experimental evidence that exposure to the space environment tends to mask the optical signatures of phyllosilicates and carbonaceous matter. Our findings should contribute to rectifying potential underestimations of water and carbon content of C-type asteroids when studied through remote sensing with new-generation telescopes.

The solar forcing prepared for Phase 6 of the Coupled Model Intercomparison Project (CMIP6) has been used extensively in climate model experiments and has been tested in various intercomparison studies. Recently, an International Space Science Institute (ISSI) working group has been established to revisit the solar forcing recommendations, based on the lessons learned from CMIP6, and to assess new datasets that have become available, in order to define a road map for building a revised and extended historical solar forcing dataset for the upcoming Phase 7 of CMIP. This paper identifies the possible improvements required and outlines a strategy to address them in the planned new solar forcing dataset. Proposed major changes include the adoption of the new Total and Spectral Solar Irradiance Sensor (TSIS-1) solar reference spectrum for solar spectral irradiance and an improved description of top-of-the-atmosphere energetic electron fluxes, as well as their reconstruction back to 1850 by means of geomagnetic proxy data. In addition, there is an urgent need to consider the proposed updates in the ozone forcing dataset in order to ensure a self-consistent solar forcing in coupled models without interactive chemistry. Regarding future solar forcing, we propose consideration of stochastic ensemble forcing scenarios, ideally in concert with other natural forcings, in order to allow for realistic projections of natural forcing uncertainties. © Author(s) 2024.

The protonated form of hydrogen cyanide, HCNH+, holds significant importance in astrochemistry, serving as an intermediate species in ion-neutral reactions occurring in the cold molecular clouds. Although it plays a crucial role in the chemistry of HCN and HNC, the excitation rate coefficients of this molecular cation by the dominant interstellar colliders have not been thoroughly investigated, leading to limitations in the radiative transfer models used to derive its abundance. In this work, we present the first hyperfine-resolved excitation rate coefficients for HCNH+ induced by collisions with both He and H2 at low temperatures, addressing a crucial requirement for precise modeling of HCNH+ abundance in typical cold dense molecular clouds. Using non-local thermodynamic equilibrium (non-LTE) radiative transfer calculations, we reproduced the 1 → 0 and 2 → 1 observational spectra of HCNH+ fairly well and derived updated molecular column densities. For the TMC-1 molecular cloud, the new HCNH+ abundance is twice as large as suggested by previous LTE modeling, whereas the column density of this molecular cation is improved only by 10% in the case of the L483 proto-star. The factor of two in the case of TMC-1 most likely arises from an error in the early analysis of observational spectra rather than an effect of the LTE assumption, given that the HCNH+ lines are predominantly thermalized at densities higher than 2 × 104 cm−3. For multiline studies of clouds of moderate densities, we strongly recommend using the collisional rate coefficients reported in this work.

Most low-mass stars form in stellar clusters that also contain massive stars, which are sources of far-ultraviolet (FUV) radiation. Theoretical models predict that this FUV radiation produces photodissociation regions (PDRs) on the surfaces of protoplanetary disks around low-mass stars, which affects planet formation within the disks. We report James Webb Space Telescope and Atacama Large Millimeter Array observations of a FUV-irradiated protoplanetary disk in the Orion Nebula. Emission lines are detected from the PDR; modeling their kinematics and excitation allowed us to constrain the physical conditions within the gas. We quantified the mass-loss rate induced by the FUV irradiation and found that it is sufficient to remove gas from the disk in less than a million years. This is rapid enough to affect giant planet formation in the disk.

Clouds of gas and dust in the Galaxy are nurseries in which stars and planetary systems are born. During their journey from the diffuse interstellar medium to the protoplanetary disks, molecular solids accumulate on cold dust grains by accretion and surface chemistry. These so-called icy grains will continuously evolve, notably by collision and aggregation processes, modifying their sizes. Our ‘Ice Age’ James Webb Space Telescope observations of the dense Chamaeleon I cloud reveal that this growth starts early, before the protostellar phase, substantially modifying the ice band profiles in the spectra. Spectral analysis confirms that the grains reach micrometre sizes, implying myriad changes in local microphysics, including mass transfer from small to large grains, reduction in the grain surface available for chemistry and modification of the penetration and propagation of radiation fields. Deformation of the observed profiles complicates the determination of chemical abundance. Observing the extensive icy grain growth in dense clouds quantitatively constrains the grain size evolution before star and planet formation.

[Aims] In this work, we determine the expected yield of detections of solar-like oscillations for the targets of the foreseen PLATO ESA mission. Our estimates are based on a study of the detection probability, which takes into account the properties of the target stars, using the information available in the PIC 1.1.0, including the current best estimate of the signal-to-noise ratio (S/N). The stellar samples, as defined for this mission, include those with the lowest noise level (P1 and P2 samples) and the P5 sample, which has a higher noise level. For the P1 and P2 samples, the S/N is high enough (by construction) that we can assume that the individual mode frequencies can be measured. For these stars, we estimate the expected uncertainties in mass, radius, and age due to statistical errors induced by uncertainties from the observations only., [Methods] We used a formulation from the literature to calculate the detection probability. We validated this formulation and the underlying assumptions with Kepler data. Once validated, we applied this approach to the PLATO samples. Using again Kepler data as a calibration set, we also derived relations to estimate the uncertainties of seismically inferred stellar mass, radius, and age. We then applied those relations to the main sequence stars with masses equal to or below 1.2 M⊙ belonging to the PLATO P1 and P2 samples and for which we predict a positive seismic detection., [Results] We found that we can expect positive detections of solar-like oscillations for more than 15 000 FGK stars in one single field after a two-year observation run. Among them, 1131 main sequence stars with masses of ≤1.2 M⊙ satisfy the PLATO requirements for the uncertainties of the seismically inferred stellar masses, radii, and ages. The baseline observation programme of PLATO consists of observing two fields of similar size (one in the southern hemisphere and one in the northern hemisphere) for two years apiece. Accordingly, the expected seismic yields of the mission amount to over 30 000 FGK dwarfs and subgiants, with positive detections of solar-like oscillations. This sample of expected solar-like oscillating stars is large enough to enable the PLATO mission’s stellar objectives to be amply satisfied., [Conclusions] The PLATO mission is expected to produce a catalog sample of extremely well seismically characterized stars of a quality that is equivalent to the Kepler Legacy sample, but containing a number that is about 80 times greater, when observing two PLATO fields for two years apiece. These stars are a gold mine that will make it possible to make significant advances in stellar modelling.

Observations of low-mass stars have frequently shown a disagreement between observed stellar radii and radii predicted by theoretical stellar structure models. This ‘radius inflation’ problem could have an impact on both stellar and exoplanetary science. We present the final results of our observation programme with the CHaracterising ExOPlanet Satellite (CHEOPS) to obtain high-precision light curves of eclipsing binaries with low-mass stellar companions (EBLMs). Combined with the spectroscopic orbits of the solar-type companions, we can derive the masses, radii, and effective temperatures of 23 M-dwarf stars. We use the PYCHEOPS data analysis software to analyse their primary and secondary occultations. For all but one target, we also perform analyses with Transiting Exoplanet Survey Satellite (TESS) light curves for comparison. We have assessed the impact of starspot-induced variation on our derived parameters and account for this in our radius and effective temperature uncertainties using simulated light curves. We observe trends in inflation with both metallicity and orbital separation. We also observe a strong trend in the difference between theoretical and observational effective temperatures with metallicity. There is no such trend with orbital separation. These results are not consistent with the idea that the observed inflation in stellar radius combines with lower effective temperature to preserve the luminosity predicted by low-mass stellar models. Our EBLM systems provide high-quality and homogeneous measurements that can be used in further studies of radius inflation.

This 1:30,000 scale geological map describes Oxia Planum, Mars, the landing site for the ExoMars Rosalind Franklin rover mission. The map represents our current understanding of bedrock units and their relationships prior to Rosalind Franklin’s exploration of this location. The map details 15 bedrock units organised into 6 groups and 7 textural and surficial units. The bedrock units were identified using visible and near-infrared remote sensing datasets. The objectives of this map are (i) to identify where the most astrobiologically relevant rocks are likely to be found, (ii) to show where hypotheses about their geological context (within Oxia Planum and in the wider geological history of Mars) can be tested, (iii) to inform both the long-term (hundreds of metres to ∼1 km) and the short-term (tens of metres) activity planning for rover exploration, and (iv) to allow the samples analysed by the rover to be interpreted within their regional geological context.

Despite the important hydrological and ecological implications of the snowpack, its real time monitoring remains challenging. This is particularly relevant in relation to the Snow Water Equivalent (SWE), as the available technologies which measure it, exhibit a number of limitations that hinder their operational implementation. In this work, we explore the potential of a new technology, Muon Scattering Radiography, to infer the SWE. We coupled snowpack simulations generated by the SNOWPACK model, with a muon scattering simulation program based on GEANT4. The SWE is modeled as a function of the muon scattering distributions. Predictions of the SWE along the year are provided showing a root-mean-square error of 12 mm for 5 hr continuous measurements. We also performed laboratory measurements using ice samples, confirming the SWE estimation capabilities and the potential of the technique to operate as a SWE monitoring tool., [Plain Language Summary] The monitoring of the seasonal snowpack is important to understand and predict the dynamics of the hydrological and ecological processes, but its continuous monitoring is still a scientific challenge. Particularly in relation to the Snow Water Equivalent (SWE). The available technologies to monitor the SWE exhibit a number of limitations that prevent its use in many real world cases. Here we explore the potential of a new technology, Muon Scattering Radiography (MSR), to quantify the SWE. MSR is a technique based in the detection of the natural and innocuous radiation of muon particles. The technique consists in the measurement of muon deviations, which are larger when the muons cross very large or dense materials. In this analysis, we simulated the snowpack evolution itself and its measurement process. Then, we determined the relation between muon deviations and SWE. Finally, we estimated the precision in the determination of SWE comparing the predictions to the ground truth in simulation. The results yielded a precision of about 12 mm. We also performed laboratory measurements with ice samples, using a 4 layer muon detector based on multiwire proportional chambers, confirming the potential of the technique to operate as a SWE monitoring tool.

This paper describes the demonstration of a regional high-resolution level-3 (L3) altimeter data unification and altimeter combination system (DUACS) developed with support from the French space agency (CNES). Deduced from full-rate (20 Hz to 40 Hz) level-2 (L2) altimeter measurements, this product provides sea level anomalies (SLA) and other essential physical variables at a spatial resolution of one sample every ~1 km over the North Atlantic Ocean. This allows us to resolve wavelengths from ~35 km to ~55 km depending on the altimeter considered. This was made possible by recent advances in radar altimeter processing for both synthetic aperture radar (SAR) and low-resolution-mode (LRM) measurements, as well as improvements made to different stages of the DUACS processing chain. Firstly, the new adaptive and low-resolution with range migration correction (LR-RMC) processing techniques were considered for Jason and Sentinel-3 (S3A), respectively. They significantly reduce errors at short wavelengths, and the adaptive processing also reduces possible land contamination near the coast. Next, up-to-date geophysical and environmental corrections were selected for this production. This includes specific corrections intended to reduce the measurement noise on LRM measurements and thus enhance the observability at short wavelengths. Compared with the 1 Hz product, the observable wavelengths reached with the demonstration high-resolution product are reduced by up to one third, or up to half in the northeast Atlantic region. The residual noises were optimally filtered from full-rate measurements, taking into consideration the different observing capabilities of the altimeters processed. A specific data recovery strategy was applied, significantly optimizing the data availability, both in the coastal and open ocean areas. This demonstration L3 product is thus better resolved than the conventional 1 Hz product, especially near the coast, where it is defined up to ~5 km against

The complexity of real-world geophysical systems is often compounded by the fact that the observed measurements depend on hidden variables. These latent variables include unresolved small scales and/or rapidly evolving processes, partially observed couplings, or forcings in coupled systems. This is the case in ocean–atmosphere dynamics, for which unknown interior dynamics can affect surface observations. The identification of computationally-relevant representations of such partially-observed and highly nonlinear systems is thus challenging and often limited to short-term forecast applications. Here, we investigate the physics-constrained learning of implicit dynamical embeddings, leveraging neural ordinary differential equation (NODE) representations. In particular, we restrict the NODE representation to linear–quadratic dynamics and enforce a global boundedness constraint, which promotes the generalization of the learned dynamics to arbitrary initial conditions. The proposed architecture is implemented within a deep learning framework, and its relevance is demonstrated with respect to state-of-the-art schemes for different case studies representative of geophysical dynamics.

Gas phase Elemental abundances in molecular CloudS (GEMS) is an IRAM 30m large program aimed at determining the elemental abundances of carbon (C), oxygen (O), nitrogen (N), and sulfur (S) in a selected set of prototypical star-forming filaments. In particular, the elemental abundance of S remains uncertain by several orders of magnitude and its determination is one of the most challenging goals of this program. We have carried out an extensive chemical modeling of the fractional abundances of CO, HCO$^+$, HCN, HNC, CS, SO, H$_2$S, OCS, and HCS$^+$ to determine the sulfur depletion toward the 244 positions in the GEMS database. These positions sample visual extinctions from A$_V$ $\sim$ 3 mag to $>$50 mag, molecular hydrogen densities ranging from a few 10$^3$~cm$^{-3}$ to 3$\times$10$^6$~cm$^{-3}$, and T$_k$ $\sim$ 10$-$35 K. Most of the positions in Taurus and Perseus are best fitted assuming early-time chemistry, t=0.1 Myr, $\zeta_{H_2}$$\sim$ (0.5$-$1)$\times$10$^{-16}$ s$^{-1}$, and [S/H]$\sim$1.5$\times$10$^{-6}$. On the contrary, most of the positions in Orion are fitted with t=1~Myr and $\zeta_{H_2}$$\sim$ 10$^{-17}$ s$^{-1}$. Moreover, $\sim$40% of the positions in Orion are best fitted assuming the undepleted sulfur abundance, [S/H]$\sim$1.5$\times$10$^{-5}$. Our results suggest that sulfur depletion depends on the environment. While the abundances of sulfur-bearing species are consistent with undepleted sulfur in Orion, a depletion factor of $\sim$20 is required to explain those observed in Taurus and Perseus. We propose that differences in the grain charge distribution in the envelopes of the studied clouds might explain these variations. The shocks associated with past and ongoing star formation could also contribute to enhance [S/H] in Orion.

We present the detection of ethanol (C2H5OH), acetone (CH3COCH3), and propanal (C2H5CHO) toward the cyanopolyyne peak of TMC-1. These three O-bearing complex organic molecules are known to be present in warm interstellar clouds, but had never been observed in a starless core. The addition of these three new pieces to the puzzle of complex organic molecules in cold interstellar clouds stresses the rich chemical diversity of cold dense cores in stages prior to the onset of star formation. The detections of ethanol, acetone, and propanal were made in the framework of QUIJOTE, a deep line survey of TMC-1 in the Q band that is being carried out with the Yebes 40m telescope. We derive column densities of (1.1 +/- 0.3)e12 cm-2 for C2H5OH, (1.4 +/- 0.6)e11 cm-2 for CH3COCH3, and (1.9 +/- 0.7)e11 cm-2 for C2H5CHO. The formation of these three O-bearing complex organic molecules is investigated with the aid of a detailed chemical model which includes gas and ice chemistry. The calculated abundances at a time around 2e5 yr are in reasonable agreement with the values derived from the observations. The formation mechanisms of these molecules in our chemical model are as follows. Ethanol is formed on grains by addition of atomic carbon on methanol followed by hydrogenation and non-thermal desorption. Acetone and propanal are produced by the gas-phase reaction between atomic oxygen and two different isomers of the C3H7 radical, where the latter follows from the hydrogenation of C3 on grains followed by non-thermal desorption. A gas-phase route involving the formation of (CH3)2COH+ through several ion-neutral reactions followed by its dissociative recombination with electrons do also contribute to the formation of acetone.

The goal of the Sea2Cloud project is to study the interplay between surface ocean biogeochemical and physical properties, fluxes to the atmosphere, and ultimately their impact on cloud formation under minimal direct anthropogenic influence. Here we present an interdisciplinary approach, combining atmospheric physics and chemistry with marine biogeochemistry, during a voyage between 41° and 47°S in March 2020. In parallel to ambient measurements of atmospheric composition and seawater biogeochemical properties, we describe semicontrolled experiments to characterize nascent sea spray properties and nucleation from gas-phase biogenic emissions. The experimental framework for studying the impact of the predicted evolution of ozone concentration in the Southern Hemisphere is also detailed. After describing the experimental strategy, we present the oceanic and meteorological context including provisional results on atmospheric thermodynamics, composition, and flux measurements. In situ measurements and flux studies were carried out on different biological communities by sampling surface seawater from subantarctic, subtropical, and frontal water masses. Air–Sea-Interface Tanks (ASIT) were used to quantify biogenic emissions of trace gases under realistic environmental conditions, with nucleation observed in association with biogenic seawater emissions. Sea spray continuously generated produced sea spray fluxes of 34% of organic matter by mass, of which 4% particles had fluorescent properties, and which size distribution resembled the one found in clean sectors of the Southern Ocean. The goal of Sea2Cloud is to generate realistic parameterizations of emission flux dependences of trace gases and nucleation precursors, sea spray, cloud condensation nuclei, and ice nuclei using seawater biogeochemistry, for implementation in regional atmospheric models.

We report on the Imaging X-ray Polarimetry Explorer (IXPE) observation of the closest and X-ray brightest Compton-thick active galactic nucleus (AGN), the Circinus galaxy. We find the source to be significantly polarized in the 2–6 keV band. From previous studies, the X-ray spectrum is known to be dominated by reflection components, both neutral (torus) and ionized (ionization cones). Our analysis indicates that the polarization degree is 28 ± 7 per cent (at 68 per cent confidence level) for the neutral reflector, with a polarization angle of 18° ± 5°, roughly perpendicular to the radio jet. The polarization of the ionized reflection is unconstrained. A comparison with Monte Carlo simulations of the polarization expected from the torus shows that the neutral reflector is consistent with being an equatorial torus with a half-opening angle of 45°–55°. This is the first X-ray polarization detection in a Seyfert galaxy, demonstrating the power of X-ray polarimetry in probing the geometry of the circumnuclear regions of AGNs, and confirming the basic predictions of standard Unification Models. © 2022 The Author(s) Published by Oxford University Press on behalf of Royal Astronomical Society.

Harmful algal blooms (HABs) have severe environmental and economic impacts worldwide. Improving HAB detection is crucial because massive blooms are likely to increase in both frequency and amplitude in the next decades due to global warming and escalating coastal eutrophication. While satellite remote sensing has proved useful to detect red tides and support HAB monitoring, the discrimination of the dominant bloom-forming species is still a challenge, all the more as the observation of highly concentrated phytoplankton patches can be hampered by a too coarse spatial resolution. Moreover, the majority of HAB studies generally had a regional focus, and a limited number of species were separately documented so far. Here, we provide a broader perspective for red tides remote sensing to better resolve HAB optical and taxonomical diversity. The main objective of the present study was to identify how many optical bloom types could be recognized with the high spatial resolution Sentinel-2 (S2) satellite mission. For that purpose, an extensive database of massive, nearly monospecific blooms, both documented in situ and using synchronous S2 observation was compiled. More than 100 S2 images of various red tides were selected worldwide. Altogether, the S2 database covered the typical reflectance spectra of 27 red tide forming species. The remote-sensing reflectance of each red tide was analysed to evaluate S2 ability to distinguish the dominant species of the bloom. A hierarchical clustering analysis suggested that six optical bloom types could be identified: (1) surface accumulation of cyanobacteria or of green Noctiluca scintillans, (2) surface accumulation of red N. scintillans (a purely heterotrophic plankton devoid of chlorophyll a), (3) red tides of Mesodinium rubrum (a phycoerythrin-bearing ciliate), (4) green seawater discolorations of Lepidodinium chlorophorum (a dinoflagellate with unusual carotenoids), (5) blooms dominated by a dinoflagellate such as Prorocentrum, Gy

The transient X-ray source MAXI J1810−222 was discovered in 2018 and has been active ever since. A long combined radio and X-ray monitoring campaign was performed with Australia Telescope Compact Array and Swift, respectively. It has been proposed that MAXI J1810−222 is a relatively distant black hole X-ray binary, albeit showing a very peculiar outburst behaviour. Here, we report on the spectral study of this source making use of a large sample of NICER observations performed between 2019 February and 2020 September. We detected a strong spectral absorption feature at ∼1 keV, which we have characterized with a physical photoionization model. Via a deep scan of the parameters space, we obtained evidence for a spectral-state dependent outflow, with mildly relativistic speeds. In particular, the soft and intermediate states point to a hot plasma outflowing at 0.05–0.15c. This speeds rule-out thermal winds and hence, they suggest that such outflows could be radiation pressure or (most likely) magnetically driven winds. Our results are crucial to test current theoretical models of wind formation in X-ray binaries.

[Context] The so-called action variables are specific functions of the positions and velocities that remain constant along the stellar orbit. The astrometry provided by Gaia Early Data Release 3 (EDR3), combined with the velocities inferred from the Radial Velocity Spectrograph (RVS) spectra of Gaia DR3, allows for the estimation of these actions for the largest volume of stars to date., [Aims] We explore such actions with the aim of locating structures in the Galactic disc., [Methods] We computed the actions and the orbital parameters of the Gaia DR3 stars, assuming an axisymmetric model for the Milky Way. Using Gaia DR3 photometric data, we also selected a subset of giant stars with better astrometry as a control sample., [Results] We find that the maps of the percentiles of the radial action JR reveal arc-like segments. We found a high JR region centered at R ≈ 10.5 kpc of 1 kpc width, as well as three arc-shape regions dominated by circular orbits at inner radii. We also identified the spiral arms in the overdensities of the giant population., [Conclusions] For Galactic coordinates (X, Y, Z), we find good agreement with the literature in the innermost region for the Scutum-Sagittarius spiral arms. At larger radii, the low JR structure tracks the Local arm at negative X, while for the Perseus arm, the agreement is restricted to the X < 2 kpc region, with a displacement with respect to the literature at more negative longitudes. We detected a high JR area at a Galactocentric radii of ∼10.5 kpc, consistent with some estimations of the Outer Lindblad Resonance location. We conclude that the pattern in the dynamics of the old stars is consistent in several places with the spatial distribution of the spiral arms traced by young populations, with small potential contributions from the moving groups.

This work provides evidence of the influence of large volcanic eruptions on Sahel rainfall relying on PMIP4/past1000 multi-model simulations, covering the last millennium. A classification of volcanic eruptions in the last millennium according to the meridional symmetry of the associated radiative forcing reveals different mechanisms of the West African Monsoon response at inter-annual timescale. In all cases, these simulated changes result in Sahel drying up to 2 years after an eruption. Besides, we add evidence of a role of varying volcanic activity across the past millennium in the Sahel precipitation variability at multi-decadal to secular timescales.

We report on an X-ray polarimetric observation of the high-mass X-ray binary LMC X-1 in the high/soft state, obtained by the Imaging X-ray Polarimetry Explorer (IXPE) in 2022 October. The measured polarization is below the minimum detectable polarization of 1.1 per cent (at the 99 per cent confidence level). Simultaneously, the source was observed with the Neutron Star Interior Composition Explorer (NICER), Nuclear Spectroscopic Telescope Array (NuSTAR), and Spectrum-Rontgen-Gamma (SRG)/Astronomical Roentgen Telescope – X-ray Concentrator (ART-XC) instruments, which enabled spectral decomposition into a dominant thermal component and a Comptonized one. The low 2–8 keV polarization of the source did not allow for strong constraints on the black hole spin and inclination of the accretion disc. However, if the orbital inclination of about 36 is assumed, then the upper limit is consistent with predictions for pure thermal emission from geometrically thin and optically thick discs. Assuming the polarization degree of the Comptonization component to be 0, 4, or 10 per cent, and oriented perpendicular to the polarization of the disc emission (in turn assumed to be perpendicular to the large-scale ionization cone orientation detected in the optical band), an upper limit to the polarization of the disc emission of 1.0, 0.9, or 0.9 per cent, respectively, is found (at the 99 per cent confidence level). © The Author(s) 2023.

We report on a comprehensive analysis of simultaneous X-ray polarimetric and spectral data of the bright atoll source GX 9+9 with the Imaging X-ray Polarimetry Explorer (IXPE) and NuSTAR. The source is significantly polarized in the 4- 8 keV band, with a degree of 2.2% ± 0.5% (uncertainty at the 68% confidence level). The NuSTAR broad-band spectrum clearly shows an iron line, and is well described by a model including thermal disc emission, a Comptonized component, and reflection. From a spectro-polarimetric fit, we obtain an upper limit to the polarization degree of the disc of 4% (at the 99% confidence level), while the contribution of Comptonized and reflected radiation cannot be conclusively separated. However, the polarization is consistent with resulting from a combination of Comptonization in a boundary or spreading layer, plus reflection off the disc, which significantly contributes in any realistic scenario. The Authors 2023.

Context. Extreme ultraviolet (EUV) observations of the quiet solar atmosphere reveal extended regions of weak emission compared to the ambient quiescent corona. The magnetic nature of these coronal features is not well understood. Aims. We study the magnetic properties of the weakly emitting extended regions, which we name coronal voids. In particular, we aim to understand whether these voids result from a reduced heat input into the corona or if they are associated with mainly unipolar and possibly open magnetic fields, similar to coronal holes. Methods. We defined the coronal voids via an intensity threshold of 75% of the mean quiet-Sun (QS) EUV intensity observed by the high-resolution EUV channel (HRIEUV) of the Extreme Ultraviolet Imager on Solar Orbiter. The line-of-sight magnetograms of the same solar region recorded by the High Resolution Telescope of the Polarimetric and Helioseismic Imager allowed us to compare the photospheric magnetic field beneath the coronal voids with that in other parts of the QS. Results. The coronal voids studied here range in size from a few granules to a few supergranules and on average exhibit a reduced intensity of 67% of the mean value of the entire field of view. The magnetic flux density in the photosphere below the voids is 76% (or more) lower than in the surrounding QS. Specifically, the coronal voids show much weaker or no network structures. The detected flux imbalances fall in the range of imbalances found in QS areas of the same size. Conclusions. We conclude that coronal voids form because of locally reduced heating of the corona due to reduced magnetic flux density in the photosphere. This makes them a distinct class of (dark) structure, different from coronal holes. © The Authors 2023

Context. The brightness of faculae and network depends on the angle at which they are observed and the magnetic flux density. Close to the limb, the assessment of this relationship has until now been hindered by the increasingly lower signal in magnetograms. Aims. This preliminary study aims to highlight the potential of using simultaneous observations from different vantage points to better determine the properties of faculae close to the limb. Methods. We used data from the Solar Orbiter/Polarimetric and Helioseismic Imager (SO/PHI), and the Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI), recorded at an -60 angular separation of their lines of sight at the Sun. We used the continuum intensity observed close to the limb by SO/PHI and complemented it with the co-observed BLOS from SDO/HMI, originating closer to disc centre (as seen by SDO/HMI), thus avoiding the degradation of the magnetic field signal near the limb. Results. We derived the dependence of facular brightness in the continuum on disc position and magnetic flux density from the combined observations of SO/PHI and SDO/HMI. Compared with a single point of view, we were able to obtain contrast values reaching closer to the limb and to lower field strengths. We find the general dependence of the limb distance at which the contrast is maximum on the flux density to be at large in line with single viewpoint observations, in that the higher the flux density is, the closer the turning point lies to the limb. There is a tendency, however, for the maximum to be reached closer to the limb when determined from two vantage points. We note that due to the preliminary nature of this study, these results must be taken with caution. Conclusions. Our analysis shows that studies involving two viewpoints can significantly improve the detection of faculae near the solar limb and the determination of their brightness contrast relative to the quiet Sun. © The Authors 2023.

Contact. Spectropolarimetric reconstructions of the photospheric vector magnetic field are intrinsically limited by the 180 ambiguity in the orientation of the transverse component. So far, the removal of such an ambiguity has required assumptions about the properties of the photospheric field, which makes disambiguation methods model-dependent. Aims. The successful launch and operation of Solar Orbiter have made the removal of the 180 ambiguity possible solely using observations of the same location on the Sun obtained from two different vantage points. Methods. The basic idea is that the unambiguous line-of-sight component of the field measured from one vantage point will generally have a nonzero projection on the ambiguous transverse component measured by the second telescope, thereby determining the "true"orientation of the transverse field. Such an idea was developed and implemented as part of the stereoscopic disambiguation method (SDM), which was recently tested using numerical simulations. Results. In this work we present a first application of the SDM to data obtained by the High Resolution Telescope (HRT) on board Solar Orbiter during the March 2022 campaign, when the angle with Earth was 27 degrees. The method was successfully applied to remove the ambiguity in the transverse component of the vector magnetogram solely using observations (from HRT and from the Helioseismic and Magnetic Imager) for the first time. Conclusions. The SDM is proven to provide observation-only disambiguated vector magnetograms that are spatially homogeneous and consistent. A discussion on the sources of error that may limit the accuracy of the method, and strategies to remove them in future applications, is also presented. © 2023 The Authors 2023.

Arch-like loop structures filled with million Kelvin hot plasma form the building blocks of the quiet-Sun corona. Both high-resolution observations and magnetoconvection simulations show the ubiquitous presence of magnetic fields on the solar surface on small spatial scales of ∼100 km. However, the question of how exactly these quiet-Sun coronal loops originate from the photosphere and how the magnetic energy from the surface is channeled to heat the overlying atmosphere is a long-standing puzzle. Here we report high-resolution photospheric magnetic field and coronal data acquired during the second science perihelion of Solar Orbiter that reveal a highly dynamic magnetic landscape underlying the observed quiet-Sun corona. We found that coronal loops often connect to surface regions that harbor fleeting weaker, mixed-polarity magnetic field patches structured on small spatial scales, and that coronal disturbances could emerge from these areas. We suggest that weaker magnetic fields with fluxes as low as 10 Mx and/or those that evolve on timescales less than 5 minutes are crucial to understanding the coronal structuring and dynamics. © 2023. The Author(s). Published by the American Astronomical Society., Full list of authors: Chitta, L. P.; Solanki, S. K.; del Toro Iniesta, J. C.; Woch, J.; Calchetti, D.; Gandorfer, A.; Hirzberger, J.; Kahil, F.; Valori, G.; Orozco Suárez, D.; Strecker, H.; Appourchaux, T.; Volkmer, R.; Peter, H.; Mandal, S.; Aznar Cuadrado, R.; Teriaca, L.; Schühle, U.; Berghmans, D.; Verbeeck, C.; Zhukov, A. N.; Priest, E. R.

Visual navigation has become a standard in robotic applications with the emergence of robust and versatile algorithms. In particular, Visual Odometry (VO) has proven to be the most reliable navigation solution for space missions to estimate an unmanned vehicle's motion and state. Lava Tubes exploration is one of the recent challenges in this field of applied robotics. VO in this scenario requires more robustness to poor lighting conditions while keeping a low computational cost. We propose investigating an indirect bi-monocular VO based on sliding-window optimization in such a context. It focuses on maintaining the sparsity of the problem while keeping the information of the marginalized frames to reduce the computational burden. Different sparse graph topologies are studied to encode information from the past and are evaluated on accuracy and computation load. The best method retained is then compared to state-of-the-art systems on real data under extreme illumination conditions and reaches similar accuracy results at a lower computational cost.

The extreme wave height distribution in the Caribbean Sea is studied using a new method based on the maximum basin-wide aggregate of significant wave height, Hs, values per month. Besides, by means of the Self-Organizing Maps (SOM) technique, we identify coherent geographical regions with similar extreme wave height variability in the Caribbean Sea. Our findings revealed three primary regions: the eastern side with comparatively lower values, the central region with intermediate values, and the western side with the highest extreme wave heights. The study also examines the wind forcing conditions driving the spatial and temporal variability of the extreme waves, highlighting the influence of the low-pressure belt dynamics as well as the role played by the Caribbean Low-Level Jet (CLLJ) index, and the impact of cold fronts and hurricanes on extreme wave heights. Additionally, we explore the relationship between the extreme wave height distribution and climatic indices, such as the Atlantic Multidecadal Oscillation (AMO), the Atlantic Meridional Mode (AMM), the North Atlantic Oscillation (NAO) and the Oceanic Niño (ONI). The results reveal that the spatial distribution of extreme wave heights in the Caribbean Sea is mostly ruled by the influence of the CLLJ, with correlations close to 80%. In addition, significant correlations were observed between the extreme wave heights and the ENSO in the central Caribbean, as well as positive correlations between the extreme wave heights and NAO in the eastern part of the basin, and significant values of correlation with the negative phases of AMO and AMM in the whole basin. We show that, unlike conventional (or broadly used) methods deployed to identify extreme wave height, such as percentile 99th, (Formula presented.), our methodology allows a further assessment of the wind and climate forcing conditions associated with the extreme wave events. Although, we acknowledge that the method here presented has limitations to capture e

We present a study of the tilt-to-length coupling noise during the LISA Pathfinder mission and how it depended on the system’s alignment. Tilt-to-length coupling noise is the unwanted coupling of angular and lateral spacecraft or test mass motion into the primary interferometric displacement readout. It was one of the major noise sources in the LISA Pathfinder mission and is likewise expected to be a primary noise source in LISA. We demonstrate here that a recently derived and published analytical model describes the dependency of the LISA Pathfinder tilt-to-length coupling noise on the alignment of the two freely falling test masses. This was verified with the data taken before and after the realignments performed in March (engineering days) and June 2016, and during a two-day experiment in February 2017 (long cross-talk experiment). The latter was performed with the explicit goal of testing the tilt-to-length coupling noise dependency on the test mass alignment. Using the analytical model, we show that all realignments performed during the mission were only partially successful and explain the reasons why. In addition to the analytical model, we computed another physical tilt-to-length coupling model via a minimizing routine making use of the long cross-talk experiment data. A similar approach could prove useful for the LISA mission.

(Abridged) Massive stars can determine the evolution of molecular clouds with their strong ultraviolet (UV) radiation fields. Moreover, UV radiation is relevant in setting the thermal gas pressure in star-forming clouds, whose influence can extend from the rims of molecular clouds to entire star-forming galaxies. Probing the fundamental structure of nearby molecular clouds is therefore crucial to understand how massive stars shape their surrounding medium and how fast molecular clouds are destroyed, specifically at their UV-illuminated edges, where models predict an intermediate zone of neutral atomic gas between the molecular cloud and the surrounding ionized gas whose size is directly related to the exposed physical conditions. We present the highest angular resolution (~$0.5$", corresponding to $207$ au) and velocity-resolved images of the molecular gas emission in the Horsehead nebula, using CO J=3-2 and HCO$^+$ J=4-3 observations with ALMA. We find that CO and HCO$^+$ are present at the edge of the cloud, very close to the ionization (H$^+$/H) and dissociation fronts (H/H$_2$), suggesting a very thin layer of neutral atomic gas (<$650$ au) and a small amount of CO-dark gas ($A_V=0.006-0.26$ mag) for stellar UV illumination conditions typical of molecular clouds in the Milky Way. The new ALMA observations reveal a web of molecular gas filaments with an estimated thermal gas pressure of $P_{\mathrm{th}} = (2.3 - 4.0) \times 10^6$ K cm$^{-3}$, and the presence of a steep density gradient at the cloud edge that can be well explained by stationary isobaric PDR models with pressures consistent with our estimations. However, in the HII region and PDR interface, we find $P_{\mathrm{th,PDR}} > P_{\mathrm{th,HII}}$, suggesting the gas is slightly compressed. Therefore, dynamical effects cannot be completely ruled out and even higher angular observations will be needed to unveil their role.

The horizontal-to-vertical (H/V) spectral ratio inversion is a traditional technique for deriving the local subsurface structure on Earth. We calculated the H/V from the ambient vibrations at different wind levels at the InSight landing site, on Mars, and also computed the H/V from the S-wave coda of the martian seismic events (marsquakes). Different H/V curves were obtained for different wind periods and from the marsquakes. From the ambient vibrations, the recordings during low-wind periods are close to the instrument self-noise level. During high-wind periods, the seismic recordings are highly contaminated by the interaction of the lander with the wind and the martian ground. Therefore, these recordings are less favourable for traditional H/V analysis. Instead, the recordings of the S-wave coda of marsquakes were preferred to derive the characteristic H/V curve of this site between 0.4 and 10 Hz. The final H/V curve presents a characteristic trough at 2.4 Hz and a strong peak at 8 Hz. Using a full diffuse wavefield approach as the forward computation and the Neighbourhood Algorithm as the sampling technique, we invert for the 1-D shear wave velocity structure at the InSight landing site. Based on our inversion results, we propose a strong site effect at the InSight site to be due to the presence of a shallow high-velocity layer (SHVL) over low-velocity units. The SHVL is likely placed below a layer of coarse blocky ejecta and can be associated with Early Amazonian basaltic lava flows. The units below the SHVL have lower velocities, possibly related to a Late Hesperian or Early Amazonian epoch with a different magmatic regime and/or a greater impact rate and more extensive weathering. An extremely weak buried low velocity layer (bLVL) between these lava flows explains the data around the 2.4 Hz trough, whereas a more competent bLVL would not generate this latter feature. These subsurface models are in good agreement with results from hammering experiment and compl

We present timing solutions for 21 pulsars discovered in 350 MHz surveys using the Green Bank Telescope (GBT). All were discovered in the Green Bank North Celestial Cap pulsar survey, with the exception of PSR J0957−0619, which was found in the GBT 350 MHz Drift-scan pulsar survey. The majority of our timing observations were made with the GBT at 820 MHz. With a spin period of 37 ms and a 528 days orbit, PSR J0032+6946 joins a small group of five other mildly recycled wide binary pulsars, for which the duration of recycling through accretion is limited by the length of the companion's giant phase. PSRs J0141+6303 and J1327+3423 are new disrupted recycled pulsars. We incorporate Arecibo observations from the NANOGrav pulsar timing array into our analysis of the latter. We also observed PSR J1327+3423 with the Long Wavelength Array, and our data suggest a frequency-dependent dispersion measure. PSR J0957−0619 was discovered as a rotating radio transient, but is a nulling pulsar at 820 MHz. PSR J1239+3239 is a new millisecond pulsar (MSP) in a 4 days orbit with a low-mass companion. Four of our pulsars already have published timing solutions, which we update in this work: the recycled wide binary PSR J0214+5222, the noneclipsing black widow PSR J0636+5128, the disrupted recycled pulsar J1434+7257, and the eclipsing binary MSP J1816+4510, which is in an 8.7 hr orbit with a redback-mass companion.

The shallowest intracrustal layer (extending to 8 ± 2 km depth) beneath the Mars InSight Lander site exhibits low seismic wave velocity, which is likely related to a combination of high porosity and other lithological factors. The SsPp phase, an SV- to P-wave reflection on the receiver side, is naturally suited for constraining the seismic structure of this top crustal layer since its prominent signal makes it observable with a single station without the need for stacking. We have analyzed six broadband and low-frequency seismic events recorded on Mars and made the first coherent detection of the SsPp phase on the red planet. The timing and amplitude of SsPp confirm the existence of the ∼8 km interface in the crust and the large wave speed (or impedance) contrast across it. With our new constraints from the SsPp phase, we determined that the average P-wave speed in the top crustal layer is between 2.5 and 3.2 km/s, which is a more precise and robust estimate than the previous range of 2.0–3.5 km/s obtained by receiver function analysis. The low velocity of Layer 1 likely results from the presence of relatively low-density lithified sedimentary rocks and/or aqueously altered igneous rocks that also have a significant amount of porosity, possibly as much as 22%–30% by volume (assuming an aspect ratio of 0.1 for the pore space). These porosities and average P-wave speeds are compatible with our current understanding of the upper crustal stratigraphy beneath the InSight Lander site.

Forty years ago, it was proposed that gas-phase organic chemistry in the interstellar medium can be initiated by the methyl cation CH3+ (refs. 1–3), but so far it has not been observed outside the Solar System 4,5. Alternative routes involving processes on grain surfaces have been invoked 6,7. Here we report James Webb Space Telescope observations of CH3+ in a protoplanetary disk in the Orion star-forming region. We find that gas-phase organic chemistry is activated by ultraviolet irradiation.

Context. The availability of large bandwidth receivers for millimeter radio telescopes allows the acquisition of position-position-frequency data cubes over a wide field of view and a broad frequency coverage. These cubes contain much information on the physical, chemical, and kinematical properties of the emitting gas. However, their large size coupled with inhomogenous signal-to-noise ratio (SNR) are major challenges for consistent analysis and interpretation.Aims. We search for a denoising method of the low SNR regions of the studied data cubes that would allow to recover the low SNR emission without distorting the signals with high SNR.Methods. We perform an in-depth data analysis of the 13 CO and C 17 O (1 -- 0) data cubes obtained as part of the ORION-B large program performed at the IRAM 30m telescope. We analyse the statistical properties of the noise and the evolution of the correlation of the signal in a given frequency channel with that of the adjacent channels. This allows us to propose significant improvements of typical autoassociative neural networks, often used to denoise hyperspectral Earth remote sensing data. Applying this method to the 13 CO (1 -- 0) cube, we compare the denoised data with those derived with the multiple Gaussian fitting algorithm ROHSA, considered as the state of the art procedure for data line cubes.Results. The nature of astronomical spectral data cubes is distinct from that of the hyperspectral data usually studied in the Earth remote sensing literature because the observed intensities become statistically independent beyond a short channel separation. This lack of redundancy in data has led us to adapt the method, notably by taking into account the sparsity of the signal along the spectral axis. The application of the proposed algorithm leads to an increase of the SNR in voxels with weak signal, while preserving the spectral shape of the data in high SNR voxels.Conclusions. The proposed algorithm that combines a detailed ana

[Aims] DE Boo is a unique system, with an edge-on view through the debris disk around the star. The disk, which is analogous to the Kuiper belt in the Solar System, was reported to extend from 74 to 84 AU from the central star. The high photometric precision of the Characterising Exoplanet Satellite (CHEOPS) provided an exceptional opportunity to observe small variations in the light curve due to transiting material in the disk. This is a unique chance to investigate processes in the debris disk., [Methods] Photometric observations of DE Boo of a total of four days were carried out with CHEOPS. Photometric variations due to spots on the stellar surface were subtracted from the light curves by applying a two-spot model and a fourth-order polynomial. The photometric observations were accompanied by spectroscopic measurements with the 1m RCC telescope at Piszkésteto and with the SOPHIE spectrograph in order to refine the astrophysical parameters of DE Boo., [Results] We present a detailed analysis of the photometric observation of DE Boo. We report the presence of nonperiodic transient features in the residual light curves with a transit duration of 0.3–0.8 days. We calculated the maximum distance of the material responsible for these variations to be 2.47 AU from the central star, much closer than most of the mass of the debris disk. Furthermore, we report the first observation of flaring events in this system., [Conclusions] We interpreted the transient features as the result of scattering in an inner debris disk around DE Boo. The processes responsible for these variations were investigated in the context of interactions between planetesimals in the system.

Knowledge of Martian crust and uppermost mantle aid us studying the planet's evolution. NASA's InSight mission provides seismic data being used to reveal the interior structure. Most studies have focused on the crustal structure beneath InSight lander, but the seismic structure of other regions has remained poorly known. We use surface-wave data to investigate the crustal structure of a large region along the Medusa Fossae Formation and the dichotomy. We adopt the largest-magnitude marsquake (S1222a) that has been recorded, which provides both Rayleigh- and Love-wave signals. We measure and jointly invert these surface-wave fundamental-mode group velocities from ∼15 to 40 s to estimate the average 1D isotropic velocity models. These models includes a high-velocity layer at ∼7-km depth, which could be due to a regional basaltic activity or regional stress. Our models also indicate that a common intra-crustal structure (∼12–40 km depth) may exist in this region along the dichotomy.

The European Space Agency's Euclid mission will provide high-quality imaging for about 1.5 billion galaxies. A software pipeline to automatically process and analyse such a huge amount of data in real time is being developed by the Science Ground Segment of the Euclid Consortium; this pipeline will include a model-fitting algorithm, which will provide photometric and morphological estimates of paramount importance for the core science goals of the mission and for legacy science. The Euclid Morphology Challenge is a comparative investigation of the performance of five model-fitting software packages on simulated Euclid data, aimed at providing the baseline to identify the best-suited algorithm to be implemented in the pipeline. In this paper we describe the simulated dataset, and we discuss the photometry results. A companion paper is focussed on the structural and morphological estimates. We created mock Euclid images simulating five fields of view of 0.48 deg2 each in the IE band of the VIS instrument, containing a total of about one and a half million galaxies (of which 350 000 have a nominal signal-to-noise ratio above 5), each with three realisations of galaxy profiles (single and double Sérsic, and 'realistic' profiles obtained with a neural network); for one of the fields in the double Sérsic realisation, we also simulated images for the three near-infrared YE, JE, and HE bands of the NISP-P instrument, and five Rubin/LSST optical complementary bands (u, g, r, i, and z), which together form a typical dataset for an Euclid observation. The images were simulated at the expected Euclid Wide Survey depths. To analyse the results, we created diagnostic plots and defined metrics to take into account the completeness of the provided catalogues, as well as the median biases, dispersions, and outlier fractions of their measured flux distributions. Five model-fitting software packages (DeepLeGATo, Galapagos-2, Morfometryka, ProFit, and SourceXtractor++) were compared, a

The Euclid mission will conduct an extragalactic survey over 15 000 deg2 of the extragalactic sky. The spectroscopic channel of the Near-Infrared Spectrometer and Photometer (NISP) has a resolution of R ~ 450 for its blue and red grisms that collectively cover the 0.93–1.89 µm range. NISP will obtain spectroscopic redshifts for 3 × 107 galaxies for the experiments on galaxy clustering, baryonic acoustic oscillations, and redshift space distortion. The wavelength calibration must be accurate within 5 Å to avoid systematics in the redshifts and downstream cosmological parameters. The NISP pre-flight dispersion laws for the grisms were obtained on the ground using a Fabry-Perot etalon. Launch vibrations, zero gravity conditions, and thermal stabilisation may alter these dispersion laws, requiring an in-flight recalibration. To this end, we use the emission lines in the spectra of compact planetary nebulae (PNe), which were selected from a PN database. To ensure completeness of the PN sample, we developed a novel technique to identify compact and strong line emitters in Gaia spectroscopic data using the Gaia spectra shape coefficients. We obtained VLT/X-shooter spectra from 0.3 to 2.5 µm for 19 PNe in excellent seeing conditions and a wide slit, mimicking Euclid’s slitless spectroscopy mode but with a ten times higher spectral resolution. Additional observations of one northern PN were obtained in the 0.80–1.90 µm range with the GMOS and GNIRS instruments at the Gemini North Observatory. The collected spectra were combined into an atlas of heliocentric vacuum wavelengths with a joint statistical and systematic accuracy of 0.1 Å in the optical and 0.3 Å in the near-infrared. The wavelength atlas and the related 1D and 2D spectra are made publicly available.

Material outgassing in a vacuum leads to molecular contamination, a well-known problem in spaceflight. Water is the most common contaminant in cryogenic spacecraft, altering numerous properties of optical systems. Too much ice means that Euclid’s calibration requirements cannot be met anymore. Euclid must then be thermally decontaminated, which is a month-long risky operation. We need to understand how ice affects our data to build adequate calibration and survey plans. A comprehensive analysis in the context of an astrophysical space survey has not been done before. In this paper we look at other spacecraft with well-documented outgassing records. We then review the formation of thin ice films, and find that for Euclid a mix of amorphous and crystalline ices is expected. Their surface topography – and thus optical properties – depend on the competing energetic needs of the substrate-water and the water-water interfaces, and they are hard to predict with current theories. We illustrate that with scanning-tunnelling and atomic-force microscope images of thin ice films. Sophisticated tools exist to compute contamination rates, and we must understand their underlying physical principles and uncertainties. We find considerable knowledge errors on the diffusion and sublimation coefficients, limiting the accuracy of outgassing estimates. We developed a water transport model to compute contamination rates in Euclid, and find agreement with industry estimates within the uncertainties. Tests of the Euclid flight hardware in space simulators did not pick up significant contamination signals, but they were also not geared towards this purpose; our in-flight calibration observations will be much more sensitive. To derive a calibration and decontamination strategy, we need to understand the link between the amount of ice in the optics and its effect on the data. There is little research about this, possibly because other spacecraft can decontaminate more easily, quenching the ne

The Aneto glacier, although it may be considered a very small glacier (<0.5 km2), is the largest glacier in the Pyrenees. Its surface and thickness loss have been continuous in recent decades, and there have been signs of accelerated melting in recent years. In this study, thickness and surface losses of the Aneto glacier from 1981 to 2022 are investigated using historical aerial imagery, airborne lidar point clouds and unoccupied aerial vehicle (UAV) imagery. A ground-penetrating radar (GPR) survey conducted in 2020, combined with data from photogrammetric analyses, allowed us to reconstruct the current ice thickness and also the existing ice distribution in 1981 and 2011. Over the last 41 years, the total glacierised area has decreased by 64.7 %, and the ice thickness has decreased, on average, by 30.5 m. The mean remaining ice thickness in autumn 2022 was 11.9 m, as against the mean thickness of 32.9, 19.2 and 15.0 m reconstructed for 1981 and 2011 and observed in 2020, respectively. The results demonstrate the critical situation of the glacier, with an imminent segmentation into two smaller ice bodies and no evidence of an accumulation zone. We also found that the occurrence of an extremely hot and dry year, as observed in the 2021–2022 season, leads to a drastic degradation of the glacier, posing a high risk to the persistence of the Aneto glacier, a situation that could extend to the rest of the Pyrenean glaciers in a relatively short time.

The Laser Interferometer Space Antenna (LISA) has two scientific objectives of cosmological focus: to probe the expansion rate of the universe, and to understand stochastic gravitational-wave backgrounds and their implications for early universe and particle physics, from the MeV to the Planck scale. However, the range of potential cosmological applications of gravitational-wave observations extends well beyond these two objectives. This publication presents a summary of the state of the art in LISA cosmology, theory and methods, and identifies new opportunities to use gravitational-wave observations by LISA to probe the universe.

The unprecedented precision of satellite laser altimetry data from the NASA ICESat-2 mission and the increasing availability of high-resolution elevation datasets open new opportunities to measure snow depth in mountains, a critical variable for ecosystem and water resource monitoring. We retrieved snow depth over the upper Tuolumne basin (California, USA) for 3 years by differencing ICESat-2 ATL06 snow-on elevations and various snow-off digital elevation models. Snow depth derived from ATL06 data only (snow-on and snow-off) offers a poor temporal and spatial coverage, limiting its potential utility. However, using a digital terrain model from airborne lidar surveys as the snow-off elevation source yielded a snow depth accuracy of ∼ 0.2 m (bias) and precision of ∼ 1 m (random error) across the basin, with an improved precision of 0.5 m for low slopes (< 10∘), compared to eight reference airborne lidar snow depth maps. Snow depths derived from ICESat-2 ATL06 and a satellite photogrammetry digital elevation model have a larger bias and reduced precision, partly induced by increased errors in forested areas. These various combinations of repeated ICESat-2 snow surface elevation measurements with satellite or airborne products will enable tailored approaches to map snow depth and estimate water resource availability in mountainous areas with limited snow depth observations.

In order to make accurate inferences about the solar interior using helioseismology, it is essential to understand all the relevant physical effects on the observations. One effect to understand is the (complex-valued) ratio of the horizontal to vertical displacement of the p- and f-modes at the height at which they are observed. Unfortunately, it is impossible to measure this ratio directly from a single vantage point, and it has been difficult to disentangle observationally from other effects. In this paper we attempt to measure the ratio directly using 7.5 h of simultaneous observations from the Polarimetric and Helioseismic Imager on board Solar Orbiter and the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory. While image geometry problems make it difficult to determine the exact ratio, it appears to agree well with that expected from adiabatic oscillations in a standard solar model. On the other hand it does not agree with a commonly used approximation, indicating that this approximation should not be used in helioseismic analyses. In addition, the ratio appears to be real-valued. © The Authors 2023

The impact of the Double Asteroid Redirection Test spacecraft into Dimorphos, moon of the asteroid Didymos, changed Dimorphos’s orbit substantially, largely from the ejection of material. We present results from 12 Earth-based facilities involved in a world-wide campaign to monitor the brightness and morphology of the ejecta in the first 35 days after impact. After an initial brightening of ∼1.4 mag, we find consistent dimming rates of 0.11-0.12 mag day−1 in the first week, and 0.08-0.09 mag day−1 over the entire study period. The system returned to its pre-impact brightness 24.3-25.3 days after impact though the primary ejecta tail remained. The dimming paused briefly eight days after impact, near in time to the appearance of the second tail. This was likely due to a secondary release of material after re-impact of a boulder released in the initial impact, though movement of the primary ejecta through the aperture likely played a role. © 2023. The Author(s). Published by the American Astronomical Society.

The phase- and energy-resolved polarization measurements of accreting X-ray pulsars (XRPs) allow us to test different theoretical models of their emission, and they also provide an avenue to determine the emission region geometry. We present the results of the observations of the XRP GX 301-2 performed with the Imaging X-ray Polarimetry Explorer (IXPE). A persistent XRP, GX 301-2 has one of the longest spin periods known: ∼680 s. A massive hyper-giant companion star Wray 977 supplies mass to the neutron star via powerful stellar winds. We did not detect significant polarization in the phase-averaged data when using spectro-polarimetric analysis, with the upper limit on the polarization degree (PD) of 2.3% (99% confidence level). Using the phase-resolved spectro-polarimetric analysis, we obtained a significant detection of polarization (above 99% confidence level) in two out of nine phase bins and a marginal detection in three bins, with a PD ranging between ∼3% and ∼10% and a polarization angle varying in a very wide range from ∼0 to ∼160. Using the rotating vector model, we obtained constraints on the pulsar geometry using both phase-binned and unbinned analyses, finding excellent agreement. Finally, we discuss possible reasons for a low observed polarization in GX 301-2. © 2023 EDP Sciences. All rights reserved.

Data assimilation techniques that integrate available observations with snow models have been proposed as a viable option to simultaneously help constrain model uncertainty and add value to observations by improving estimates of the snowpack state. However, the propagation of information from spatially sparse observations in high-resolution simulations remains an under-explored topic. To remedy this, the development of data assimilation techniques that can spread information in space is a crucial step. Herein, we examine the potential of spatio-temporal data assimilation for integrating sparse snow depth observations with hyper-resolution (5 m) snow simulations in the Izas central Pyrenean experimental catchment (Spain). Our experiments were developed using the Multiple Snow Data Assimilation System (MuSA) with new improvements to tackle the spatio-temporal data assimilation. Therein, we used a deterministic ensemble smoother with multiple data assimilation (DES-MDA) with domain localization. Three different experiments were performed to showcase the capabilities of spatio-temporal information transfer in hyper-resolution snow simulations. Experiment I employed the conventional geographical Euclidean distance to map the similarity between cells. Experiment II utilized the Mahalanobis distance in a multi-dimensional topographic space using terrain parameters extracted from a digital elevation model. Experiment III utilized a more direct mapping of snowpack similarity from a single complete snow depth map together with the easting and northing coordinates. Although all experiments showed a noticeable improvement in the snow patterns in the catchment compared with the deterministic open loop in terms of correlation (r=0.13) and root mean square error (RMSE = 1.11 m), the use of topographical dimensions (Experiment II, r=0.63 and RMSE = 0.89 m) and observations (Experiments III, r=0.92 and RMSE = 0.44 m) largely outperform the simulated patterns in Experiment I (r=0.3

Photometric galaxy surveys probe the late-time Universe where the density field is highly non-Gaussian. A consequence is the emergence of the super-sample covariance (SSC), a non-Gaussian covariance term that is sensitive to fluctuations on scales larger than the survey window. In this work, we study the impact of the survey geometry on the SSC and, subsequently, on cosmological parameter inference. We devise a fast SSC approximation that accounts for the survey geometry and compare its performance to the common approximation of rescaling the results by the fraction of the sky covered by the survey, fSKY, dubbed ‘full-sky approximation’. To gauge the impact of our new SSC recipe, that we call ‘partial-sky’, we perform Fisher forecasts on the parameters of the (w0, wa)-CDM model in a 3 × 2 point analysis, varying the survey area, the geometry of the mask, and the galaxy distribution inside our redshift bins. The differences in the marginalised forecast errors –with the full-sky approximation performing poorly for small survey areas but excellently for stage-IV-like areas– are found to be absorbed by the marginalisation on galaxy bias nuisance parameters. For large survey areas, the unmarginalised errors are underestimated by about 10% for all probes considered. This is a hint that, even for stage-IV-like surveys, the partial-sky method introduced in this work will be necessary if tight priors are applied on these nuisance parameters. We make the partial-sky method public with a new release of the public code PySSC.

We present a machine learning framework to simulate realistic galaxies for the Euclid Survey, producing more complex and realistic galaxies than the analytical simulations currently used in Euclid. The proposed method combines a control on galaxy shape parameters offered by analytic models with realistic surface brightness distributions learned from real Hubble Space Telescope observations by deep generative models. We simulate a galaxy field of 0.4 deg2 as it will be seen by the Euclid visible imager VIS, and we show that galaxy structural parameters are recovered to an accuracy similar to that for pure analytic Sérsic profiles. Based on these simulations, we estimate that the Euclid Wide Survey (EWS) will be able to resolve the internal morphological structure of galaxies down to a surface brightness of 22.5 mag arcsec−2, and the Euclid Deep Survey (EDS) down to 24.9 mag arcsec−2. This corresponds to approximately 250 million galaxies at the end of the mission and a 50% complete sample for stellar masses above 1010.6 M⊙ (resp. 109.6 M⊙) at a redshift z ∼ 0.5 for the EWS (resp. EDS). The approach presented in this work can contribute to improving the preparation of future high-precision cosmological imaging surveys by allowing simulations to incorporate more realistic galaxies.