Search

We present a new S-wave velocity tomographic model of the upper mantle beneath the Australian Plate and its boundaries that we call Aus22. It includes azimuthal anisotropy and was constrained by waveforms from 0.9 million vertical-component seismograms, with the densest data sampling in the hemisphere centred on the Australian continent, using all available data covering this hemisphere. Waveform inversion extracted structural information from surface waves, S- and multiple S-waves and constrained S- and P-wave speeds and S-wave azimuthal anisotropy of the crust and upper mantle, down to the 660-km discontinuity. The model was validated by resolution tests and, for particular locations in Australia with notable differences from previous models, independent inter-station measurements of surface-wave phase velocities. Aus22 can be used to constrain the structure and evolution of the Australian Plate and its boundaries in fine detail at the regional scale. Thick, high-velocity (and, by inference, cold) cratonic lithosphere occupies nearly all of western and central Australia but shows substantial lateral heterogeneity. It extends up to the northern edge of the plate, where it collides with island arcs, without subducting. Diamondiferous kimberlites and lamproite deposits are underlain by cratonic lithosphere, except for the most recent diamondiferous lamproites in the King Leopold Orogen. The rugged eastern boundary of the cratonic lithosphere resolved by the model provides a lithospheric definition of the Tasman Line. Just east of the Tasman Line, an area of intermediate-thick lithosphere is observed in the southern part of the continent. The eastern part of Australia is underlain by thin, warm lithosphere, evidenced by low seismic velocities. All the sites of Cenozoic intraplate volcanism in eastern Australia are located on thin lithosphere. A low-velocity anomaly is present in the mantle transition zone (410–660 km depths) beneath the Lord Howe and Tasmanid hotspots, indicative of anomalously high temperature and consistent with a deep mantle upwelling feeding these hotspots and, possibly, also the East Australia hotspot. High seismic velocities at 200–410 km depth below New Guinea indicate the presence of slab fragments, probably linked to the subduction of the Australian Plate. High seismic velocities are observed in the transition zone below northeast Australia and indicate the presence of subducted lithospheric fragments trapped in the transition zone, possibly parts of the former northern continental margin of Australia. [ABSTRACT FROM AUTHOR]

In recent years, there has been a considerable expansion of deployments of portable seismic stations across Australia, which have been analysed by receiver function or autocorrelation methods to extract estimates of Moho depth. An ongoing program of full-crustal reflection profiles has now provided more than 25 000 km of reflection transects that have been interpreted for Moho structure. The Moho data set is further augmented by extensive marine reflection results. These new data sources have been combined with earlier refraction and receiver function results to provide full continental coverage, though some desert areas remain with limited sampling. The dense sampling of the Moho indicates the presence of rapid changes in Moho depth, and so, the Moho surface has been constructed using an approach that allows different weighting and spatial influence depending on the nature of the estimate. The inclusion of Moho results from continental-wide gravity inversion with low weighting helps to resolve the continent-ocean transition and to provide additional control in the least sampled zones. The refined distribution indicates the presence of widespread smaller-scale variations in Moho structure. Strong lateral contrasts in crustal thickness remain, but some have become more subdued with improved sampling of critical areas. The main differences from earlier results lie in previously poorly sampled regions around the Lake Eyre Basin, where additional passive seismic results indicate somewhat thicker crust though still with a strong contrast in crustal thickness to the cratonic zone to the west. [ABSTRACT FROM AUTHOR]

Most studies of sea level budget only indicate the stationary causes for the global mean sea level (GMSL) rise over the course of a time span of interest, providing limited information on temporary changes in the GMSL budget. In this contribution, we present an instantaneous budget of the GMSL for the period of 2003–2015, which offers new insights on the time evolution of the GMSL budget. We use a space-state model to compute the instantaneous rates of GMSL and their contributing sources (barystatic and steric contributions), to investigate the causes of the GMSL instantaneous budget, which accounts for and quantifies low-frequency variations in GMSL rise. Combining the barystatic GMSL instantaneous rates (estimated from satellite gravimetry) along with a particular choice of steric GMSL instantaneous rates (inferred from an Argo product) achieves closure of the GMSL instantaneous budget with a standard deviation of 0.4 mm yr−1. We find that the barystatic GMSL rate is primarily responsible for the striking fluctuations in the GMSL instantaneous rates between 2009 and 2015. For example, over only 20 months (July 2010 to February 2012), the GMSL (barystatic) rate increases from 0.32 (0.71) to 6.12 (4.55) mm yr−1. These strong fluctuations mainly result from hydrology and mountain glacier mass variations rather than mass loss in ice sheet, for instance, Australia slowed ∼1.2 mm yr−1 equivalent contribution to GMSL rise over the year 2010, which is approximately four times the linear rate of Antarctica mass loss. [ABSTRACT FROM AUTHOR]

Gravity gradient tensors (GGTs) are used to investigate the density of subsurface structures in the Earth's crust and can reduce ambiguity during data interpretation. However, costs and research area restrictions often prevent their application during surveys, thereby limiting their utility. To address this limitation, a matrix equation based on the Taylor series expansion that uses the gravity vector and its neighbours was formulated to obtain the GGTs. Higher-order derivatives of the gravity vector were utilized to constrain the calculation, which improved the accuracy of the transformation. Synthetic data were used to demonstrate that the proposed approach improved accuracy when the radial component of the gravity vector was transformed into GGTs. This approach was also applied to gravity data from the Otway Basin in Australia. Compared with the measured GGT, the results obtained using the proposed approach had a relative error of 0.46. [ABSTRACT FROM AUTHOR]

The equivalent source technique is a powerful and widely used method for processing gravity and magnetic data. Nevertheless, its major drawback is the large computational cost in terms of processing time and computer memory. We present two techniques for reducing the computational cost of equivalent source processing: block-averaging source locations and the gradient-boosted equivalent source algorithm. Through block-averaging, we reduce the number of source coefficients that must be estimated while retaining the minimum desired resolution in the final processed data. With the gradient-boosting method, we estimate the sources coefficients in small batches along overlapping windows, allowing us to reduce the computer memory requirements arbitrarily to conform to the constraints of the available hardware. We show that the combination of block-averaging and gradient-boosted equivalent sources is capable of producing accurate interpolations through tests against synthetic data. Moreover, we demonstrate the feasibility of our method by gridding a gravity data set covering Australia with over 1.7 million observations using a modest personal computer. [ABSTRACT FROM AUTHOR]

Observations from different disciplines have shown that our planet is highly heterogeneous at multiple scale lengths. Still, many seismological Earth models tend not to include any small-scale heterogeneity or lateral velocity variations, which can affect measurements and predictions based on these homogeneous models. In this study, we describe the lithospheric small-scale isotropic heterogeneity structure in terms of the intrinsic, diffusion and scattering quality factors, as well as an autocorrelation function, associated with a characteristic scale length (a) and RMS fractional velocity fluctuations (ε). To obtain this characterization, we combined a single-layer and a multilayer energy flux models with a new Bayesian inference algorithm. Our synthetic tests show that this technique can successfully retrieve the input parameter values for 1- or 2-layer models and that our Bayesian algorithm can resolve whether the data can be fitted by a single set of parameters or a range of models is required instead, even for very complex posterior probability distributions. We applied this technique to three seismic arrays in Australia: Alice Springs array (ASAR), Warramunga Array (WRA) and Pilbara Seismic Array (PSAR). Our single-layer model results suggest intrinsic and diffusion attenuation are strongest for ASAR, while scattering and total attenuation are similarly strong for ASAR and WRA. All quality factors take higher values for PSAR than for the other two arrays, implying that the structure beneath this array is less attenuating and heterogeneous than for ASAR or WRA. The multilayer model results show the crust is more heterogeneous than the lithospheric mantle for all arrays. Crustal correlation lengths and RMS velocity fluctuations for these arrays range from ∼0.2 to 1.5 km and ∼2.3 to 3.9 per cent, respectively. Parameter values for the upper mantle are not unique, with combinations of low values of the parameters (a < 2 km and ε < ∼2.5 per cent) being as likely as those with high correlation length and velocity variations (a > 5 km and ε > ∼2.5 per cent, respectively). We attribute the similarities in the attenuation and heterogeneity structure beneath ASAR and WRA to their location on the proterozoic North Australian Craton, as opposed to PSAR, which lies on the archaean West Australian Craton. Differences in the small-scale structure beneath ASAR and WRA can be ascribed to the different tectonic histories of these two regions of the same craton. Overall, our results highlight the suitability of the combination of an energy flux model and a Bayesian inference algorithm for future scattering and small-scale heterogeneity studies, since our approach allows us to obtain and compare the different quality factors, while also giving us detailed information about the trade-offs and uncertainties in the determination of the scattering parameters. [ABSTRACT FROM AUTHOR]