  This work proposes the AutoHVSR algorithm that allows for fully-automated
processing of horizontal-to-vertical spectral ratio (HVSR) measurements,
including those with zero, one, or multiple clear resonances. The AutoHVSR
algorithm integrates robust signal processing and computational methods with
state-of-the-art machine-learning models trained using a diverse dataset of
1109 HVSR measurements. The AutoHVSR algorithm demonstrates excellent
performance by correctly determining the number of HVSR resonances for 1099 of
the 1109 HVSR measurements (>99%) and predicting the mean resonant frequency of
the correctly identified resonances with a root mean square error (RMSE) of
0.05 Hz. Furthermore, the AutoHVSR algorithm was able to produce these
predictions in 13 minutes (including HVSR processing time) compared to the 30
hours required for traditional processing (a speed up of 138). The AutoHVSR
algorithm is further demonstrated on a challenging dataset from Canterbury, New
Zealand that included many HVSR curves with multiple and/or ambiguous
resonances. The AutoHVSR algorithm was capable of correctly determining the
number of HVSR resonances for 113 of the 129 HVSR measurements (>87%) and
predicting the mean resonant frequency of the correctly identified resonances
with a RMSE of 0.10 Hz. The AutoHVSR algorithm produced these predictions under
2 minutes (including HVSR processing time) compared to the 4 hours required for
traditional processing (a speed up of 120). Finally, while the AutoHVSR
algorithm was developed using microtremor measurements where the horizontal
components were combined using the geometric mean, it is shown to extend
without modification to microtremor HVSR measurements where the two horizontal
components are rotated azimuthally and to HVSR measurements from earthquake
recordings. The AutoHVSR algorithm has been made publicly available in v0.3.0
of HVSRweb.
