Geophysical Journal International

SummaryWe present and apply a pseudo trans-dimensional inversion method for 3D geometrical gravity inversion, in which the number of rock units, their geometry, and their density can vary during sampling. The method is designed for efficient exploration of the model space and to infer the presence and properties of units not directly observable but detectable with geophysical data. Sampling relies on a non-reversible Metropolis-Hastings algorithm, during which rock units can be added or removed from the model, interface geometries are perturbed using random fields, and densities are sampled from distributions informed by prior information. To visualise the space of sampled models and to aid interpretation, a workflow is proposed that combines dimensionality reduction with the clustering of models in families. The capabilities of the inversion method are evaluated using two synthetic cases. The first is a motivating example aimed at recovering an intrusion missing from the prior model. It features a horizontal layer-cake where fixed-dimensional inversion fails to adequately fit the data and sample models close to the true model, while the proposed pseudo trans-dimensional approach is much more successful. The second case investigates the recovery of two missing units and the capability to overcome prior model biases. Results show the potential of our method to infer the presence of unseen geological features such as intrusions. However, they suggest that with biased prior geological modelling, it may be challenging to infer with certainty the presence of more than two previously unknown rock units at depth.

SummaryWhen highly viscous fluids are present in a rock medium, the viscous effect of such fluids cannot be neglected in the propagation of elastic waves. In this paper, a fractional thermoporoelastic theory is newly proposed, which is a further improvement of the two-temperature generalized thermoporoelastic theory. Firstly, by introducing the Kelvin-Voigt model into the stress-strain constitutive equation, the viscous effect of highly viscous fluids is considered. Then, fractional derivatives are introduced into the heat conduction equations of the solid and fluid phase to consider the anomalous heat conduction caused by the viscous effect in rock media. Plane wave analysis method is adopted to obtain the phase velocity and attenuation factor of four longitudinal waves (P1, P2, T1, T2). Numerical results show that the introduction of fluid viscosity leads to the appearance of new relaxation peaks in the P wave at high frequencies, and the introduction of fractional derivatives causes a decrease in the phase velocity and attenuation factor of T waves. The results provide a reference for further research on the wave propagation in rock media containing highly viscous fluids.

SummaryInduced polarization (IP) effects in airborne electromagnetic (AEM) surveys have commonly been investigated in helicopter-borne systems, leaving both a bibliographic and application gap for fixed-wing configurations. This gap partly reflects the large relative number of helicopter compared to fixed wing AEM systems, but also the geometric complexity of fixed wing platforms. In these platforms, nine geometric parameters come into play: the pitch, roll, and yaw of both transmitter and receiver, plus the three-axis offsets between the coils. Shifts in these factors can distort the measured data in ways that aren’t uniquely attributable, making it hard to pinpoint whether negative recordings truly arise from IP or from geometry-related effects. The non-fixed geometry also complicates removal of the primary field, often requiring iterative processing steps that may suppress or alter spectral content linked to IP. With advances in airborne IP understanding from helicopter-borne systems, revisiting fixed-wing platforms is both timely and necessary. Part A of this two-part study addresses this issue using the TEMPEST™ fixed-wing system connecting numerical modelling with field evidence. A suite of synthetic two-layer models with variable resistivity and chargeability parameters was developed to evaluate the system’s sensitivity to polarizable structures. The experiments demonstrate that IP effects, including negative secondary field responses, can be reliably detected in fixed-wing AEM data, both in X and Z magnetic field components. The capacity of these systems to detect IP phenomena is, however, strongly dependent on the electrical conductance of the environment. For instance, both fixed-wing and helicopter-borne systems, elevated near-surface conductance enhances the amplitude of purely electromagnetic induction currents, which in turn can dominate the recorded response and obscure the comparatively weaker polarization currents. More in general, IP detectability depends on the strength of the EM response generated by induction currents flowing elsewhere, which can dominate the small reverse current flow from a polarizable target. This highlights the critical role of near-surface conductivity in controlling the expression of IP responses and underscores the need to carefully account for these factors when interpreting survey data. The synthetic results are then connected with field-scale observations from a subset of the AusAEM dataset, over 470 000 line-km of TEMPESTTM data, where negative responses align with areas of low shallow conductance, confirming the simulation results. These finding open the way to the Part B of this study, where TEMPESTTM data are inverted taking into account IP and compared with helicopter-borne results and geological information.

SummaryAcoustic signals can couple to the ground, providing an opportunity to use seismic stations to investigate airborne sources. The study of Bishop et al. (2022) used wavefield simulations in a fluid-solid medium to quantify the role of topography on the seismic (ground) recordings of a monopole source in the air. We build upon this study by linking wavefield forward modeling with the source estimation code MTUQ, which can accommodate point forces or moment tensors in a solid medium, as well as sources in the air (or water) if they are enabled by the forward-modeling solver. We perform a series of synthetic numerical experiments to demonstrate that a dipole airborne source can be estimated using ground-based receivers, even within the presence of realistic topography. We investigate the influence of receiver coverage, topography, and assumed source location on the estimated results. The established capabilities raise the prospects for future efforts to estimate dipole sources in 3D models that include heterogeneity in the air and the earth in addition to topography.