AutoHVSR: A machine-learning-supported algorithm for the fully-automated processing of horizontal-to-vertical spectral ratio measurements

The horizontal-to-vertical spectral ratio (HVSR) has seen widespread use as a tool for measuring a site's resonant behavior. However, the processing of HVSR with traditional methods can be tedious and time-consuming when the HVSR is complex and exhibits multiple resonances. This work proposes the AutoHVSR algorithm that allows for fully-automated processing of HVSR measurements, including those with zero, one, or multiple clear reso-nances. The AutoHVSR algorithm accepts microtremor or earthquake recordings and user-defined HVSR pro-cessing parameters as input and returns the computed HVSR from each time window or earthquake record, the mean/median HVSR curve, and statistics on each automatically identified HVSR resonance in terms of frequency and amplitude. The AutoHVSR algorithm integrates robust signal processing and computational methods with state-of-the-art machine-learning models trained using a diverse dataset of 1108 HVSR measurements. The AutoHVSR algorithm demonstrates excellent performance by correctly determining the number of HVSR reso-nances for 1098 of the 1108 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 min (including HVSR processing time) compared to the 30 h 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. Despite the challenging nature of the Canterbury dataset, 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 min (including HVSR processing time) compared to the 4 h required for traditional processing (a speed up of 120). Finally, while the AutoHVSR al-gorithm 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, it can be accessed at https://hvsrweb.designsafe-ci.org/.