Rhys Hawkins

Spatial retarding field energy analyzer measurements downstream of a helicon double layer plasma

Spatial ion energy measurements using a retarding field energy analyzer are performed in the exhaust of a 0.30mTorr, 250W helicon double layer plasma to investigate the divergence of the argon ion beam formed by acceleration in the double layer. Various divergence angles are computed by considering the radial distribution of beam density; the average beam ion diverging by 9°. The efficiency at which momentum is imparted parallel to the longitudinal axis of the thruster is calculated to be 98%. The results show that a few centimeters downstream of the source, the beam ions do not follow the magnetic field lines.

Direct measurement of neutral gas heating in a radio-frequency electrothermal plasma micro-thruster

Direct measurements and modelling of neutral gas heating in a radio-frequency (13.56 MHz) electrothermal collisional plasma micro-thruster have been performed using rovibrational band matching of the second positive system of molecular nitrogen (N2) for operating pressures of 4.5 Torr down to 0.5 Torr. The temperature measured with decreasing pressure for 10 W power input ranged from 395 K to 530 K in pure N2 and from 834 K to 1090 K in argon with 1% N2. A simple analytical model was developed which describes the difference in temperatures between the argon and nitrogen discharges.

REDBACK: Open‐source software for efficient noise‐reduction in plate kinematic reconstructions

Knowledge of past plate motions derived from ocean-floor finite rotations is an important asset of the Earth Sciences, because it allows linking a variety of shallow-rooted and deep-rooted geological processes. Efforts have recently been taken toward inferring finite rotations at the unprecedented temporal resolution of 1 Myr or less, and more data are anticipated in the near future. These reconstructions, like any data set, feature a degree of noise that compromises significantly our ability to make geodynamical inferences. Bayesian Inference has been recently shown to be effective in reducing the impact of noise on plate kinematics inferred from high-temporal-resolution finite-rotation data sets. We describe REDBACK, an open-source software that implements transdimensional hierarchical Bayesian Inference for efficient noise-reduction in plate kinematic reconstructions. Algorithm details are described and illustrated by means of a synthetic test.

Bayesian noise-reduction in Arabia/Somalia and Nubia/Arabia finite rotations since ∼20 Ma: Implications for Nubia/Somalia relative motion

Knowledge of Nubia/Somalia relative motion since the Early Neogene is of particular importance in the Earth Sciences, because it (i) impacts on inferences on African dynamic topography; and (ii) allows us to link plate kinematics within the Indian realm with those within the Atlantic basin. The contemporary Nubia/Somalia motion is well known from geodetic observations. Precise estimates of the past-3.2-Myr average motion are also available from paleo-magnetic observations. However, little is known of the Nubia/Somalia motion prior to ∼3.2 Ma, chiefly because the Southwest Indian Ridge spread slowly, posing a challenge to precisely identify magnetic lineations. This also makes the few observations available particularly prone to noise. Here we reconstruct Nubia/Somalia relative motions since ∼20 Ma from the alternative plate-circuit Nubia-Arabia-Somalia. We resort to trans-dimensional hierarchical Bayesian Inference, which has proved effective in reducing finite-rotation noise, to unravel the Arabia/Somalia and Arabia/Nubia motions. We combine the resulting kinematics to reconstruct the Nubia/Somalia relative motion since ∼20 Ma. We verify the validity of the approach by comparing our reconstruction with the available record for the past ∼3.2 Myr, obtained through Antarctica. Results indicate that prior to ∼11 Ma the total motion between Nubia and Somalia was faster than today. Furthermore, it featured a significant strike-slip component along the Nubia/Somalia boundary. It is only since ∼11 Ma that Nubia diverges away from Somalia at slower rates, comparable to the present-day one. Kinematic changes of some 20% might have occurred in the period leading to the present-day, but plate-motion steadiness is also warranted within the uncertainties.

Tsunami source uncertainty estimation: The 2011 Japan tsunami

This article was published in the Journal of Geophysical Research: Solid Earth [©2016 American Geophysical Union] and the definite version is available at : http://dx.doi.org/10.1002/2015JB012764 The Journals website is at: http://onlinelibrary.wiley.com/doi/10.1002/2015JB012764/abstract;jsessionid=BF5813B61924FB094F31BC8F213840A1.f02t04

Is there a link between geomagnetic reversal frequency and paleointensity? A Bayesian approach

Over the last 25 years, several studies have tested for a link between geomagnetic field intensity and reversal frequency. However, despite a large increase in the number of absolute paleointensity determinations, and improved methods for obtaining such data, two competing models have arisen. Here we employ a new tool for objectively analyzing paleomagnetic time series to investigate the possibility of a link between reversal frequency and paleointensity. Transdimensional Markov chain Monte Carlo techniques are applied to a quality-filtered version of the global paleointensity (PINT) database for the last 202 Myr to model long-term paleointensity behavior. A large ensemble of models is sampled, from which a final representative mean model is extracted. The resulting paleointensity model confirms published conclusions that the single-silicate crystal method gives significantly different results from more conventional whole rock paleointensity methods; this makes it difficult to jointly model the two data types in the same analysis. When the much larger whole rock data set is considered, a stable paleointensity of 5.46 ± 0.28 × 10 22 A/m 2 for the last 202 Myr is consistent with the 95% confidence interval of the paleointensity model. Statistical tests indicate no significant correlation between reversal frequency and field intensity at the 0.05 level. However, this result is likely due to the characteristics of the PINT database rather than being a genuine, physically representative conclusion. Given the paucity of data and general state of the global paleointensity database, concerted efforts to increase the number of high-quality, well-dated paleointensity data are required before conclusions about a link between geomagnetic field intensity and reversal frequency can be confidently drawn.

Mechanical amorphization, flash heating, and frictional melting: Dramatic changes to fault surfaces during the first millisecond of earthquake slip

The evolution of fault strength and behavior during the initial stages of slip plays an important role in driving the onset of instability and fault weakening. Using small-displacement triaxial experiments on quartz sandstone, this study highlights the rapid onset of microstructural change on fault interfaces and identifies new evidence for an evolution in physical processes with increasing slip and velocity. Pre-ground fault surfaces have been slipped over a range of velocities (0.36 μm s –1 to 18 cm s –1 ) and at normal stresses comparable to upper- to mid-crustal conditions (92–287 MPa). Microstructural analysis of the fault interfaces reveals the formation of amorphous material at displacements <170 μm and slip durations < 1 ms. Mechanical and microstructural observations have been combined with numerical modeling to present the first documented examples of a transition from mechanical amorphization to flash heating, then frictional melting, with changes in slip conditions. The sequence of processes activated during the initial stages of fault movement may provide new insights into factors that influence the onset of slip in the seismogenic crust.

Three-Dimensional Mapping of Ion Density in a Double-Layer Helicon Plasma

An investigation of the ion density produced inside the Chi Kung reactor at 250 W and 0.30 mtorr was performed through Langmuir probe measurements. Three-dimensional interpolation combined these data to produce a 3-D density volume. The resulting image was used to identify three distinct regions of differing density: power deposition, acceleration, and expansion. The ion density over these regions was found to vary from ~1.3 times 1010 to ~2.7 times 1011 cm-3.

Trans-dimensional Bayesian inversion of airborne electromagnetic data for 2D conductivity profiles

This paper presents the application of a novel trans-dimensional sampling approach to a time domain airborne electromagnetic (AEM) inverse problem to solve for plausible conductivities of the subsurface. Geophysical inverse field problems, such as time domain AEM, are well known to have a large degree of non-uniqueness. Common least-squares optimisation approaches fail to take this into account and provide a single solution with linearised estimates of uncertainty that can result in overly optimistic appraisal of the conductivity of the subsurface. In this new non-linear approach, the spatial complexity of a 2D profile is controlled directly by the data. By examining an ensemble of proposed conductivity profiles it accommodates non-uniqueness and provides more robust estimates of uncertainties. We apply a novel trans-dimensional Bayesian approach using a wavelet parameterisation to airborne electromagnetic (AEM) inversions using data from the Broken Hill region. This approach allows exploration of a range of plausible subsurface conductivity models and provides more robust uncertainty estimates while accounting for potential non-uniqueness.

Particle in cell simulation of a radiofrequency plasma jet expanding in vacuum

This research was funded by the Australian Research Council Discovery Projects DP 1096653 and DP140100571.