© 2019 The Author(s) . We apply virtual deep seismic sounding (VDSS) to data collected from both permanent and temporary seismic stations in Australia with the goal of examining (i) the resilience of the method to the presence of complex lithospheric structure and (ii) the effectiveness of different approaches for estimating bulk crustal properties (namely thickness and Vp). Data from the permanent station WRAB in the Northern Territory is ideal for benchmarking VDSS (large number and favourable distribution of recorded earthquakes), with the results from several approaches agreeing on a thickness of 40-42 km. Application of VDSS to data from the temporary BILBY array, a linear distribution of broadband stations that traverses central Australia, shows that strong Moho reflections can be retrieved with as few as two earthquakes even at the transition between crustal blocks of different character and in the presence of thick sedimentary basins. Crustal thickness varies between 36 and 54 km and compares well with the reflectivity character of nearby deep seismic reflection lines. Furthermore, we find that off-line estimates of crustal thickness, calculated by binning the source regions according to back-Azimuth, produce values of crustal thickness that are consistent with the regional geology. Overall, we find that VDSS is a powerful technique for estimating crustal thickness and velocity due to its insensitivity to complex short-wavelength structure and requirement of a small number earthquakes to produce a stable result. However, not all schemes tested for extracting bulk crustal properties appear to be robust and stringent data quality checking is still required during implementation.