Analytic and numerical solutions to the seismic wave equation in continuous media.

This paper presents two approaches to mathematical modelling of a synthetic seismic pulse, and a comparison between them. First, a new analytical model is developed in two-dimensional Cartesian coordinates. Combined with an initial condition of sufficient symmetry, this provides a valuable check for the validity of the numerical method that follows. A particular initial condition is found which allows for a new closed-form solution. A numerical scheme is then presented which combines a spectral (Fourier) representation for displacement components and wave-speed parameters, a fourth-order Runge-Kutta integration method, and an absorbing boundary layer. The resulting large system of differential equations is solved in parallel on suitable enhanced performance desktop hardware in a new software implementation. This provides an alternative approach to forward modelling of waves within isotropic media which is efficient, and tailored to rapid and flexible developments in modelling seismic structure, for example, shallow depth environmental applications. Visual comparisons of the analytic solution and the numerical scheme are presented.

Evidence of micro-continent entrainment during crustal accretion.

Simple models involving the gradual outboard accretion of material along curvilinear subduction zones are often inconsistent with field-based evidence. A recent study using 3-D geodynamic modelling has shown that the entrainment of an exotic continental fragment within a simple subduction system can result in a complex phase of growth. Although kinematic models based on structural mapping and high-resolution gravity and magnetic maps indicate that the pre-Carboniferous Tasmanides in southeastern Australia may have been subjected to this process, to date there has been little corroboration from crustal scale geophysical imaging. Here, we apply Bayesian transdimensional tomography to ambient noise data recorded by the WOMBAT transportable seismic array to constrain a detailed (20 km resolution in some areas) 3-D shear velocity model of the crust beneath southeast Australia. We find that many of the velocity variations that emerge from our inversion support the recently developed geodynamic and kinematic models. In particular, the full thickness of the exotic continental block, responsible for orocline formation and the tectonic escape of the back arc region, is imaged here for the first time. Our seismic results provide the first direct evidence that exotic continental fragments may profoundly affect the development of an accretionary orogen.