-
Quasi-geostrophic convection-driven dynamos in a thick spherical shell
Authors:
Olivier Barrois,
Thomas Gastine,
Christopher C. Finlay
Abstract:
We present dynamos computed using a hybrid QG-3D numerical scheme in a thick spherical shell geometry. Our model is based on a quasi-geostrophic convection code extended with a 3D treatment of heat transport and magnetic induction. We find a collection of self-sustained, multipolar, weak field dynamos with magnetic energy one or two orders of magnitude lower than the kinetic energy. The poloidal m…
▽ More
We present dynamos computed using a hybrid QG-3D numerical scheme in a thick spherical shell geometry. Our model is based on a quasi-geostrophic convection code extended with a 3D treatment of heat transport and magnetic induction. We find a collection of self-sustained, multipolar, weak field dynamos with magnetic energy one or two orders of magnitude lower than the kinetic energy. The poloidal magnetic energy is weak and, by construction, there is a lack of equatorially anti-symmetric components in the Buoyancy and Lorentz forces. This leads to configurations where the velocity field is only weakly impacted by the magnetic field, similar to dynamos found in 3D simulations where zonal flows and the $Ω$-effect dominate. The time-dependence of these dynamos is characterised by quasi-periodic oscillations that we attribute to dynamo waves. The QG-3D dynamos found so far are not Earth-like. The inability of our setup to produce strong, dipole-dominated, magnetic fields likely points to a missing ingredient in our QG flows, and a related lack of helicity and $α$-effect. The models presented here may be more relevant for studying stellar dynamos where zonal flows are known to dominate. This study was carried out at modest control parameters, however moving to lower Ekman numbers, when smaller values of both the magnetic and hydrodynamic Prandtl numbers can be of interest, our approach will be able to gain in efficiency by using relatively coarse grids for the 3D magnetic and temperature fields and a finer grid for the QG velocity field.
△ Less
Submitted 15 December, 2023;
originally announced December 2023.
-
Trajectory Based RFI Subtraction and Calibration for Radio Interferometry
Authors:
Chris Finlay,
Bruce A. Bassett,
Martin Kunz,
Nadeem Oozeer
Abstract:
Radio interferometry calibration and Radio Frequency Interference (RFI) removal are usually done separately. Here we show that jointly modelling the antenna gains and RFI has significant benefits when the RFI follows precise trajectories, such as for satellites. One surprising benefit is improved calibration solutions, by leveraging the RFI signal itself. We present tabascal (TrAjectory BAsed RFI…
▽ More
Radio interferometry calibration and Radio Frequency Interference (RFI) removal are usually done separately. Here we show that jointly modelling the antenna gains and RFI has significant benefits when the RFI follows precise trajectories, such as for satellites. One surprising benefit is improved calibration solutions, by leveraging the RFI signal itself. We present tabascal (TrAjectory BAsed RFI Subtraction and CALibration), a new algorithm that jointly models the RFI and calibration parameters in visibilities. We test tabascal on simulated MeerKAT calibration observations contaminated by satellite-based RFI. We obtain gain estimates that are both unbiased and up to an order of magnitude better constrained compared to uncontaminated data. When combined with an ad hoc RFI subtraction scheme, tabascal solutions can be further applied to an adjacent target observation: 5 minutes of calibration data results in an image with about a third the noise achieved when using flagging alone. The recovered flux distribution of RFI subtracted data was on par with uncontaminated data. In contrast, RFI flagging alone resulted in a higher detection threshold and consistent underestimation of source fluxes. For a mean RFI amplitude of 17 Jy, using RFI subtraction leads to less than 1% loss of data compared to 75% data loss from an ideal $3σ$ flagging algorithm, a very significant increase in data available for science analysis. Although we have examined the case of satellite RFI, tabascal should work for any RFI moving on parameterizable trajectories, relative to the phase centre, such as planes and/or objects fixed to the ground.
△ Less
Submitted 28 June, 2023; v1 submitted 10 January, 2023;
originally announced January 2023.
-
Applications for CryoSat2 satellite magnetic data in studies of Earth's core field variations
Authors:
Magnus Danel Hammer,
Christopher C. Finlay,
Nils Olsen
Abstract:
We use 20 years of magnetic field measurements from the Oersted, CHAMP and Swarm satellite missions, supplemented by calibrated platform magnetometer data from the CryoSat2 satellite, to study time variations of the Earth's core field at satellite altitude and at the core-mantle boundary (CMB). From the satellite data we derive composite time series of the core field secular variation (SV) with 4m…
▽ More
We use 20 years of magnetic field measurements from the Oersted, CHAMP and Swarm satellite missions, supplemented by calibrated platform magnetometer data from the CryoSat2 satellite, to study time variations of the Earth's core field at satellite altitude and at the core-mantle boundary (CMB). From the satellite data we derive composite time series of the core field secular variation (SV) with 4month cadence, at 300 globally distributed Geomagnetic Virtual Observatories (GVO). GVO radial SV series display regional fluctuations with 5-10 years duration and amplitudes reaching 20 nT/yr, most notably at low latitudes over Indonesia (2014), over South America and the South Atlantic (2007, 2011 and 2014), and over the central Pacific (2017). Applying the Subtractive Optimally Localized Averages (SOLA) method, we map the SV at the CMB as a collection of locally averaged SV estimates. We demonstrate that using 2-year windows of CryoSat2 data, it is possible to reliably estimate the SV and its time derivative, the secular acceleration (SA), at the CMB, with a spatial resolution, corresponding to spherical harmonic degree 10. Along the CMB geographic equator, we find strong SA features under Indonesia from 2011-2014, under central America from 2015 to 2019, and sequences of SA with alternating sign under the Atlantic during 2004-2019. We find that data from CryoSat2 make a valuable contribution to the emerging picture of sub-decadal core field variations. Using 1 year windows of data from the Swarm satellites, it is possible to study SA changes at low latitudes on timescales down to 1 year, with spatial resolution corresponding to spherical harmonic degree 10. We find strong positive and negative SA features appearing side-by-side in the Pacific in 2017, and thereafter drift westward.
△ Less
Submitted 29 November, 2021;
originally announced November 2021.
-
The Hydrogen Intensity and Real-time Analysis eXperiment: 256-Element Array Status and Overview
Authors:
Devin Crichton,
Moumita Aich,
Adam Amara,
Kevin Bandura,
Bruce A. Bassett,
Carlos Bengaly,
Pascale Berner,
Shruti Bhatporia,
Martin Bucher,
Tzu-Ching Chang,
H. Cynthia Chiang,
Jean-Francois Cliche,
Carolyn Crichton,
Romeel Dave,
Dirk I. L. de Villiers,
Matt A. Dobbs,
Aaron M. Ewall-Wice,
Scott Eyono,
Christopher Finlay,
Sindhu Gaddam,
Ken Ganga,
Kevin G. Gayley,
Kit Gerodias,
Tim Gibbon,
Austin Gumba
, et al. (75 additional authors not shown)
Abstract:
The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a radio interferometer array currently in development, with an initial 256-element array to be deployed at the South African Radio Astronomy Observatory (SARAO) Square Kilometer Array (SKA) site in South Africa. Each of the 6m, $f/0.23$ dishes will be instrumented with dual-polarisation feeds operating over a frequency range of 40…
▽ More
The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a radio interferometer array currently in development, with an initial 256-element array to be deployed at the South African Radio Astronomy Observatory (SARAO) Square Kilometer Array (SKA) site in South Africa. Each of the 6m, $f/0.23$ dishes will be instrumented with dual-polarisation feeds operating over a frequency range of 400-800 MHz. Through intensity mapping of the 21 cm emission line of neutral hydrogen, HIRAX will provide a cosmological survey of the distribution of large-scale structure over the redshift range of $0.775 < z < 2.55$ over $\sim$15,000 square degrees of the southern sky. The statistical power of such a survey is sufficient to produce $\sim$7 percent constraints on the dark energy equation of state parameter when combined with measurements from the Planck satellite. Additionally, HIRAX will provide a highly competitive platform for radio transient and HI absorber science while enabling a multitude of cross-correlation studies. In this paper, we describe the science goals of the experiment, overview of the design and status of the sub-components of the telescope system, and describe the expected performance of the initial 256-element array as well as the planned future expansion to the final, 1024-element array.
△ Less
Submitted 17 January, 2022; v1 submitted 28 September, 2021;
originally announced September 2021.
-
Deep Learning improves identification of Radio Frequency Interference
Authors:
Alireza Vafaei Sadr,
Bruce A. Bassett,
Nadeem Oozeer,
Yabebal Fantaye,
Chris Finlay
Abstract:
Flagging of Radio Frequency Interference (RFI) is an increasingly important challenge in radio astronomy. We present R-Net, a deep convolutional ResNet architecture that significantly outperforms existing algorithms -- including the default MeerKAT RFI flagger, and deep U-Net architectures -- across all metrics including AUC, F1-score and MCC. We demonstrate the robustness of this improvement on b…
▽ More
Flagging of Radio Frequency Interference (RFI) is an increasingly important challenge in radio astronomy. We present R-Net, a deep convolutional ResNet architecture that significantly outperforms existing algorithms -- including the default MeerKAT RFI flagger, and deep U-Net architectures -- across all metrics including AUC, F1-score and MCC. We demonstrate the robustness of this improvement on both single dish and interferometric simulations and, using transfer learning, on real data. Our R-Net model's precision is approximately $90\%$ better than the current MeerKAT flagger at $80\%$ recall and has a 35\% higher F1-score with no additional performance cost. We further highlight the effectiveness of transfer learning from a model initially trained on simulated MeerKAT data and fine-tuned on real, human-flagged, KAT-7 data. Despite the wide differences in the nature of the two telescope arrays, the model achieves an AUC of 0.91, while the best model without transfer learning only reaches an AUC of 0.67. We consider the use of phase information in our models but find that without calibration the phase adds almost no extra information relative to amplitude data only. Our results strongly suggest that deep learning on simulations, boosted by transfer learning on real data, will likely play a key role in the future of RFI flagging of radio astronomy data.
△ Less
Submitted 12 October, 2020; v1 submitted 18 May, 2020;
originally announced May 2020.