The Australian SKA Pathfinder (ASKAP) offers powerful new capabilities for studying the polarised and magnetised Universe at radio wavelengths. In this paper, we introduce the Polarisation Sky Survey of the Universe’s Magnetism (POSSUM), a groundbreaking survey with three primary objectives: (1) to create a comprehensive Faraday rotation measure (RM) grid of up to one million compact extragalactic sources across the southern $\sim50$% of the sky (20,630 deg$^2$); (2) to map the intrinsic polarisation and RM properties of a wide range of discrete extragalactic and Galactic objects over the same area; and (3) to contribute interferometric data with excellent surface brightness sensitivity, which can be combined with single-dish data to study the diffuse Galactic interstellar medium. Observations for the full POSSUM survey commenced in May 2023 and are expected to conclude by mid-2028. POSSUM will achieve an RM grid density of around 30–50 RMs per square degree with a median measurement uncertainty of $\sim$1 rad m$^{-2}$. The survey operates primarily over a frequency range of 800–1088 MHz, with an angular resolution of 20” and a typical RMS sensitivity in Stokes Q or U of 18 $\mu$Jy beam$^{-1}$. Additionally, the survey will be supplemented by similar observations covering 1296–1440 MHz over 38% of the sky. POSSUM will enable the discovery and detailed investigation of magnetised phenomena in a wide range of cosmic environments, including the intergalactic medium and cosmic web, galaxy clusters and groups, active galactic nuclei and radio galaxies, the Magellanic System and other nearby galaxies, galaxy halos and the circumgalactic medium, and the magnetic structure of the Milky Way across a very wide range of scales, as well as the interplay between these components. This paper reviews the current science case developed by the POSSUM Collaboration and provides an overview of POSSUM’s observations, data processing, outputs, and its complementarity with other radio and multi-wavelength surveys, including future work with the SKA.

We provide an assessment of the Infinity Two fusion pilot plant (FPP) baseline plasma physics design. Infinity Two is a four-field period, aspect ratio $A = 10$, quasi-isodynamic stellarator with improved confinement appealing to a max-$J$ approach, elevated plasma density and high magnetic fields ($ \langle B\rangle = 9$ T). Here $J$ denotes the second adiabatic invariant. At the envisioned operating point ($800$ MW deuterium-tritium (DT) fusion), the configuration has robust magnetic surfaces based on magnetohydrodynamic (MHD) equilibrium calculations and is stable to both local and global MHD instabilities. The configuration has excellent confinement properties with small neoclassical transport and low bootstrap current ($|I_{bootstrap}| \sim 2$ kA). Calculations of collisional alpha-particle confinement in a DT FPP scenario show small energy losses to the first wall (${\lt}1.5 \,\%$) and stable energetic particle/Alfvén eigenmodes at high ion density. Low turbulent transport is produced using a combination of density profile control consistent with pellet fueling and reduced stiffness to turbulent transport via three-dimensional shaping. Transport simulations with the T3D-GX-SFINCS code suite with self-consistent turbulent and neoclassical transport predict that the DT fusion power$P_{{fus}}=800$ MW operating point is attainable with high fusion gain ($Q=40$) at volume-averaged electron densities $n_e\approx 2 \times 10^{20}$ m$^{-3}$, below the Sudo density limit. Additional transport calculations show that an ignited ($Q=\infty$) solution is available at slightly higher density ($2.2 \times 10^{20}$ m$^{-3}$) with $P_{{fus}}=1.5$ GW. The magnetic configuration is defined by a magnetic coil set with sufficient room for an island divertor, shielding and blanket solutions with tritium breeding ratios (TBR) above unity. An optimistic estimate for the gas-cooled solid breeder designed helium-cooled pebble bed is TBR $\sim 1.3$. Infinity Two satisfies the physics requirements of a stellarator fusion pilot plant.

Laser-driven neutron sources are routinely produced by the interaction of laser-accelerated protons with a converter. They present complementary characteristics to those of conventional accelerator-based neutron sources (e.g. short pulse durations, enabling novel applications like radiography). We present here results from an experiment aimed at performing a global characterization of the neutrons produced using the Titan laser at the Jupiter Laser Facility (Livermore, USA), where protons were accelerated from 23 $\mathrm {\mu }$m thick plastic targets and directed onto a LiF converter to produce neutrons. For this purpose, several diagnostics were used to measure these neutron emissions, such as CR-39, activation foils, time-of-flight detectors and direct measurement of $^7{\rm Be}$ residual activity in the LiF converters. The use of these different, independently operating diagnostics enables comparison of the various measurements performed to provide a robust characterization. These measurements led to a neutron yield of $2.0\times 10^{9}$ neutrons per shot with a modest angular dependence, close to that simulated.

We present the first unbiased survey of neutral hydrogen absorption in the Small Magellanic Cloud. The survey utilises pilot neutral hydrogen observations with the Australian Square Kilometre Array Pathfinder telescope as part of the Galactic Australian Square Kilometre Array Pathfinder neutral hydrogen project whose dataset has been processed with the Galactic Australian Square Kilometre Array Pathfinder-HI absorption pipeline, also described here. This dataset provides absorption spectra towards 229 continuum sources, a 275% increase in the number of continuum sources previously published in the Small Magellanic Cloud region, as well as an improvement in the quality of absorption spectra over previous surveys of the Small Magellanic Cloud. Our unbiased view, combined with the closely matched beam size between emission and absorption, reveals a lower cold gas faction (11%) than the 2019 ATCA survey of the Small Magellanic Cloud and is more representative of the Small Magellanic Cloud as a whole. We also find that the optical depth varies greatly between the Small Magellanic Cloud’s bar and wing regions. In the bar we find that the optical depth is generally low (correction factor to the optically thin column density assumption of $\mathcal{R}_{\mathrm{HI}} \sim 1.04$) but increases linearly with column density. In the wing however, there is a wide scatter in optical depth despite a tighter range of column densities.

We present the most sensitive and detailed view of the neutral hydrogen (${\rm H\small I}$) emission associated with the Small Magellanic Cloud (SMC), through the combination of data from the Australian Square Kilometre Array Pathfinder (ASKAP) and Parkes (Murriyang), as part of the Galactic Australian Square Kilometre Array Pathfinder (GASKAP) pilot survey. These GASKAP-HI pilot observations, for the first time, reveal ${\rm H\small I}$ in the SMC on similar physical scales as other important tracers of the interstellar medium, such as molecular gas and dust. The resultant image cube possesses an rms noise level of 1.1 K ($1.6\,\mathrm{mJy\ beam}^{-1}$) $\mathrm{per}\ 0.98\,\mathrm{km\ s}^{-1}$ spectral channel with an angular resolution of $30^{\prime\prime}$ (${\sim}10\,\mathrm{pc}$). We discuss the calibration scheme and the custom imaging pipeline that utilises a joint deconvolution approach, efficiently distributed across a computing cluster, to accurately recover the emission extending across the entire ${\sim}25\,\mathrm{deg}^2$ field-of-view. We provide an overview of the data products and characterise several aspects including the noise properties as a function of angular resolution and the represented spatial scales by deriving the global transfer function over the full spectral range. A preliminary spatial power spectrum analysis on individual spectral channels reveals that the power law nature of the density distribution extends down to scales of 10 pc. We highlight the scientific potential of these data by comparing the properties of an outflowing high-velocity cloud with previous ASKAP+Parkes ${\rm H\small I}$ test observations.

We present VELOCIraptor, a massively parallel galaxy/(sub)halo finder that is also capable of robustly identifying tidally disrupted objects and separate stellar halos from galaxies. The code is written in C++11, use the Message Passing Interface (MPI) and OpenMP Application Programming Interface (API) for parallelisation, and includes python tools to read/manipulate the data products produced. We demonstrate the power of the VELOCIraptor (sub)halo finder, showing how it can identify subhalos deep within the host that have negligible density contrasts to their parent halo. We find a subhalo mass-radial distance dependence: large subhalos with mass ratios of ≳10−2 are more common in the central regions than smaller subhalos, a result of dynamical friction and low tidal mass loss rates. This dependence is completely absent in (sub)halo finders in common use, which generally search for substructure in configuration space, yet is present in codes that track particles belonging to halos as they fall into other halos, such as hbt+. VELOCIraptor largely reproduces the dependence seen without tracking, finding a similar radial dependence to hbt+ in well-resolved halos from our limited resolution fiducial simulation.

Petrographic study of the Frewens sandstone, Upper Cretaceous Frontier Formation, documents reservoir-scale diagenetic heterogeneity. Iron-bearing calcite cement occurs as large concretions that generally follow bedding and are most common near the top of the sandstone. Median thickness of the concretions is 0.6 m, length 4.5 m, and width 5.7 m; median volume is 5.2 m3. Concretions comprise 12% of the sandstone.

The minus-cement porosity of concretion samples is low, indicating that the calcite precipitated near maximum burial depth. Isotopic and burial history data suggest that the calcite precipitated at ~54°C from evolved meteoric water enriched in 18O or from a mixed meteoric±marine pore-water. Shell-bearing transgressive shales above the Frewens sandstone are interpreted to be the source of calcium carbonate. Concretions of this size and distribution would influence fluid flow in a reservoir and would reduce the amount of hydrocarbons in place.

The impact of spatial and temporal variations in the surface albedo and aerodynamic roughness length on the surface energy balance of Haut Glacier d’Arolla, Switzerland, was examined using a semi-distributed surface energy-balance model (Arnold and others, 1996). The model was updated to incorporate the glacier-wide effects of albedo and aerodynamic roughness-length variations using parameterizations following Brock (1997). After the model’s performance was validated, the glacier-wide patterns of the net shortwave, turbulent and melt energy fluxes were examined on four days, representative of surface conditions in late May, June July and August. In the model, meteorological conditions were held constant on each day in order that the impact of albedo and aerodynamic roughness-length variations could be assessed independently. A late-summer snowfall event was also simulated. Albedo and aerodynamic roughness-length variations, particularly those associated with the migration of the transient snowline and the decay of the winter snowpack, were found to exert a strong influence on the magnitude of the surface energy fluxes The importance of meteorological conditions in suppressing the surface energy fluxes and melt rate following a fresh snowfall was highlighted

Modelling the hydrology of the Greenland ice sheet, including the filling and drainage of supraglacial lakes, requires melt inputs generated at high spatial and temporal resolution. Here we apply a high spatial (100 m) and temporal (1 hour) mass-balance model to a 450 km2 subset of the Paakitsoq region, West Greenland. The model is calibrated by adjusting the values for parameters of fresh snow density, threshold temperature for solid/liquid precipitation and elevation-dependent precipitation gradient to minimize the error between modelled output and surface height and albedo measurements from three Greenland Climate Network stations for the mass-balance years 2000/01 and 2004/05. Bestfit parameter values are consistent between the two years at 400 kg m-3, 2°C and +14% (100 m)-1, respectively. Model performance is evaluated, first, by comparing modelled snow and ice distribution with that derived from Landsat-7 ETM+ satellite imagery using normalized-difference snow index classification and supervised image thresholding; and second, by comparing modelled albedo with that retrieved from the MODIS sensor M0D10A1 product. Calculation of mass-balance components indicates that 6% of surface meltwater and rainwater refreezes in the snowpack and does not become runoff, such that refreezing accounts for 31% of the net accumulation.

We use a simple energy-conservation model and a model based on Röthlisberger’s theory for steady-state water flow in a subglacial conduit to model water movement between lakes in the Adventure subglacial trench region of East Antarctica during a 1996–98 jökulhlaup. Using available field evidence to constrain the models suggests that water flow would likely be accommodated in a tunnel with a cross-sectional area of 36 m2 and a value for k (the reciprocal of Manning’s roughness parameter) larger than the 12.5 m1/3 s−1 previously calculated. We also use Nye’s theory for time-dependent conduit water flow to model the temporal evolution of conduit discharge, cross-sectional area, water pressure and lake draining and filling during the flood. We initially assume one source and one sink lake. We perform sensitivity tests on the input parameter set, matching modeled source- and sink-lake depth changes with measured surface elevation data. Using a simple function for vertical ice deformation in which surface deformation scales linearly to the lake depth change, we find the scaling factor is of the order 4 × 10−3 of the ice thickness. The most likely value of k lies in the range 55–68 m1/3 s−1, and the ratio of source to sink-lake radii is approximately 1 : 1.4. Finally, we experiment using Nye’s theory to model water movement between one source and three sink lakes. The model fails to produce the observed patterns of water movement as indicated by the surface deformation data.

During the summers of 1987 and 1988, 15 dye-tracer tests from a total of eight injection points were conducted in the ablation area of Midtdalsbreen, a northern outlet of Hardangerjokulen, southern Norway. The spatial and temporal patterns of water discharge, shapes of the dye-return curves, through-flow velocities, dye-recovery rates, dispersivities, and velocity/discharge relationships suggest the existence of distinct catchments beneath the eastern and western halves of the glacier which are characterized by different types of drainage sytem. Experiments on the eastern side were associated with high melt-water discharges and produced short-lived and highly peaked dye-return curves, fast through-flow velocities, high dye-recovery rates, low dispersivity values which decreased through the melt season, and a velocity/discharge relationship with an exponent of 1.0. Experiments on the western side were associated with low melt-water discharges and produced flat, extended dye-return curves which often displayed secondary peaks, slow through-flow velocities, low dye-recovery rates, high dispersivity values which increased during the melt season, and a velocity/discharge relationship with an exponent of 0.6. Comparison of observed through-flow velocities with values calculated theoretically using various hypothetical drainage-system structures suggests that water flows in a major sinuous conduit beneath the eastern half of the glacier and in a system of linked cavities beneath the western half. A model for the seasonal evolution of the whole drainage network is postulated which has important implications for temporal variations in subglacial water pressures and glacier-sliding velocity.

A 3029-m-long deep ice core extending nearly to bedrock has been drilled at the very top of the Greenland ice sheet (Summit) by the Greenland Ice-core Project (GRIP), an international European joint effort organized by the European Science Foundation. The ice core reaches back to 250,000 yr B.P. according to dating based partly on stratigraphic methods and partly on ice-flow modeling. A continuous and detailed stable isotope (δ18O) profile along the entire core depicts dramatic temperature changes in Greenland through the last two glacial cycles, including abrupt climatic shifts during the Eem/Sangamon Interglaciation, which is elsewhere recorded as a warm and stable period. The stratigraphic continuity of the Eemian layers has therefore been scrutinized. New ice core studies, comprising cloudy band observations, deconvolution, and frequency analyses, lead to the conclusion that the climate instability suggested during the Eem Interglaciation in Greenland is likely to be real, though no conclusive evidence is available. Whereas latitudinal displacements of the North Atlantic Ocean current are considered the immediate cause of the glacial climate instability, longitudinal displacements may be the immediate cause of the Eemian instability. If so, the Eemian climate changes will be much subdued outside the Arctic region and will probably only be recognizable in sedimentary sequences of high sensitivity and temporal resolution.

Nickel–molybdenum (Ni–Mo) materials are widely used functional oxide catalysts for the hydrogen evolution reaction. In this work, we investigate the high activity of Ni–Mo by depositing size-controlled Ni nanocrystals (NCs) onto Mo substrates. We observe a synergistic increase in catalytic activity that does not scale with the Ni–Mo interface length. This evidence points to a bulk electronic interaction of the two metals that is separate from the mechanism of enhancement seen in conventionally co-deposited Ni–Mo electrocatalysts. In addition to elucidating the catalytic behavior of the Ni–Mo system, this work offers a general NC-based paradigm for investigating fundamental interactions and synergistic effects in electrocatalytic materials.

Antarctica Bedmap2 datasets are used to calculate subglacial hydraulic potential and the area, depth and volume of hydraulic potential sinks. There are over 32 000 contiguous sinks, which can be thought of as predicted lakes. Patterns of subglacial melt are modelled with a balanced ice flux flow model, and water fluxes are cumulated along predicted flow pathways to quantify steady-state fluxes from the main basin outlets and from known subglacial lakes. The total flux from the continent is ~21 km3 a−1. Byrd Glacier has the greatest basin flux of ~2.7 km3 a−1. Fluxes from subglacial lakes range from ~1 × 10−4 to ~1.5 km3 a−1. Lake turnover times are calculated from their volumes and fluxes, and have median values of ~100 a for known ‘active’ lakes and ~500 a for other lakes. Recurrence intervals of a 0.25 km3 flood range from ~2 months to ~2000 a (median ≈130 a) for known ‘active’ lakes and from ~2 to ~2400 a (median ≈ 360 a) for other lakes. Thus, several lakes that have recently been observed to fill and drain may not do so again for many centuries; and several lakes that have not, so far, been observed to fill and drain have the potential to do so, even at annual to decadal timescales.