Regularizing the quark-level σ\sigma model

We show that the finite difference, $-i\pi^2 m^2$, between quadratic and logarithmic divergent integrals $\int d^4p[m^2(p^2-m^2)^{-2}-(p^2-m^2)^{-1}]$, as encountered in the linear $\sigma$ model, is in fact regularization independent.Comment: 9 pages, 2 figures, Latex, to appear in Mod. Phys. Lett.

Zinc bioglasses regulate mineralization in human dental pulp stem cells

China Scholarship Council (CSC

Inherited crustal deformation along the East Gondwana margin revealed by seismic anisotropy tomography

Acknowledgments We thank Mallory Young for providing phase velocity measurements in mainland Australia and Tasmania. Robert Musgrave is thanked for making available his tilt-filtered magnetic intensity map. In the short term, data may be made available by contacting the authors (S.P. or N.R.). A new database of passive seismic data recorded in Australia is planned as part of a national geophysics data facility for easy access download. Details on the status of this database may be obtained from the authors (S.P., N.R., or A.M.R.). There are no restrictions on access for noncommercial use. Commercial users should seek written permission from the authors (S.P. or N.R.). Ross Cayley publishes with the permission of the Director of the Geological Survey of Victoria.Peer reviewedPublisher PD

Evidence of micro-continent entrainment during crustal accretion

Peer reviewedPublisher PD

Insights into the structure and dynamics of the upper mantle beneath Bass Strait, southeast Australia, using shear wave splitting

© 2019 Elsevier B.V. We investigate the structure of the upper mantle using teleseismic shear wave splitting measurements obtained at 32 broadband seismic stations located in Bass Strait and the surrounding region of southeast Australia. Our dataset includes ∼366 individual splitting measurements from SKS and SKKS phases. The pattern of seismic anisotropy from shear wave splitting analysis beneath the study area is complex and does not always correlate with magnetic lineaments or current N-S absolute plate motion. In the eastern Lachlan Fold Belt, fast shear waves are polarized parallel to the structural trend (∼N25E). Further south, fast shear wave polarization directions trend on average N25–75E from the Western Tasmania Terrane through Bass Strait to southern Victoria, which is consistent with the presence of an exotic Precambrian microcontinent in this region as previously postulated. Stations located on and around the Neogene-Quaternary Newer Volcanics Province in southern Victoria display sizeable delay times (∼2.7 s). These values are among the largest in the world and hence require either an unusually large intrinsic anisotropy frozen within the lithosphere, or a contribution from both the lithospheric and asthenospheric mantle. In the Eastern Tasmania Terrane, nearly all observed fast directions are approximately NW-SE. Although part of our data set strongly favours anisotropy originating from “fabric” frozen in the lithospheric mantle, a contribution from the asthenospheric flow related to the present day plate motion is also required to explain the observed splitting parameters. We suggest that deviation of asthenospheric mantle flow around lithospheric roots could be occurring, and so variations in anisotropy related to mantle flow may be expected. Alternatively, the pattern of fast polarisation orientations observed around Bass Strait may be consistent with radial mantle flow associated with a plume linked to the recently discovered Cosgrove volcanic track. However, it is difficult to characterise the relative contributions to the observed splitting from the lithospheric vs. asthenospheric upper mantle due to poor backazimuthal coverage of the data

Automatic differentiation in geophysical inverse problems

Automatic differentiation (AD) is the technique whereby output variables of a computer code evaluating any complicated function (e.g. the solution to a differential equation) can be differentiated with respect to the input variables. Often AD tools take the form of source to source translators and produce computer code without the need for deriving and hand coding of explicit mathematical formulae by the user. The power of AD lies in the fact that it combines the generality of finite difference techniques and the accuracy and efficiency of analytical derivatives, while at the same time eliminating 'human' coding errors. It also provides the possibility of accurate, efficient derivative calculation from complex 'forward' codes where no analytical derivatives are possible and finite difference techniques are too cumbersome. AD is already having a major impact in areas such as optimization, meteorology and oceanography. Similarly it has considerable potential for use in non-linear inverse problems in geophysics where linearization is desirable, or for sensitivity analysis of large numerical simulation codes, for example, wave propagation and geodynamic modelling. At present, however, AD tools appear to be little used in the geosciences. Here we report on experiments using a state of the art AD tool to perform source to source code translation in a range of geoscience problems. These include calculating derivatives for Gibbs free energy minimization, seismic receiver function inversion, and seismic ray tracing. Issues of accuracy and efficiency are discussed

Exploiting seismic signal and noise in an intracratonic environment to constrain crustal structure and source parameters of infrequent earthquakes

In many regions of the world characterized by a relatively low rate of seismicity, the determination of local and regional seismic source parameters is often restricted to an analysis of the first onsets of P waves (or first motion analysis) due to incomplete information about Earth structure and the small size of the events. When rare large earthquakes occur in these regions, their waveforms can be used to model Earth structure. This, however, makes the nature of the earthquake source determination problem circular, as source information is mapped as structure. Presented here is one possible remedy to this situation, where through a two-step approach we first constrain Earth structure using data independent of the earthquake of interest. In this study, we focus on a region in Western Australia with low seismicity and minimal instrument coverage and use the CAPRA/LP temporary deployment to demonstrate that reliable structural models of the upper lithosphere can be obtained from an independent collection of teleseismic and ambient noise datasets. Apart from teleseismic receiver functions (RFs), we obtain group velocities from the cross-correlation of ambient noise and phase velocities from the traditional two-station method using carefully selected teleseismic earthquakes and station pairs. Crustal models are then developed through the joint inversion of dispersion data and RFs, and structural Green's functions are computed from a layered composite model. In the second step of this comprehensive approach, we apply full waveform inversion (three-component body and surface waves) to the 2007 M L= 5.3 Shark Bay, Western Australia, earthquake to estimate its source parameters (seismic moment, focal mechanism, and depth). We conclude that the full waveform inversion analysis provides constraints on the orientation of fault planes superior to a first motion interpretation

Development and application of the GIM code for the Cyber 203 computer

The GIM computer code for fluid dynamics research was developed. Enhancement of the computer code, implicit algorithm development, turbulence model implementation, chemistry model development, interactive input module coding and wing/body flowfield computation are described. The GIM quasi-parabolic code development was completed, and the code used to compute a number of example cases. Turbulence models, algebraic and differential equations, were added to the basic viscous code. An equilibrium reacting chemistry model and implicit finite difference scheme were also added. Development was completed on the interactive module for generating the input data for GIM. Solutions for inviscid hypersonic flow over a wing/body configuration are also presented