**Research article**
11 Nov 2021

**Research article** | 11 Nov 2021

# An update on techniques to assess normal-mode behavior of rock arches by ambient vibrations

Mauro Häusler Paul Richmond Geimer Riley Finnegan Donat Fäh and Jeffrey Ralston Moore

^{1},

^{2},

^{2},

^{1},

^{2}

**Mauro Häusler et al.**Mauro Häusler Paul Richmond Geimer Riley Finnegan Donat Fäh and Jeffrey Ralston Moore

^{1},

^{2},

^{2},

^{1},

^{2}

^{1}Swiss Seismological Service, ETH Zurich, Zurich, 8092, Switzerland^{2}Department of Geology and Geophysics, University of Utah, Salt Lake City, UT 84112, USA

^{1}Swiss Seismological Service, ETH Zurich, Zurich, 8092, Switzerland^{2}Department of Geology and Geophysics, University of Utah, Salt Lake City, UT 84112, USA

**Correspondence**: Mauro Häusler (mauro.haeusler@sed.ethz.ch)

**Correspondence**: Mauro Häusler (mauro.haeusler@sed.ethz.ch)

Received: 16 Apr 2021 – Discussion started: 06 May 2021 – Revised: 07 Oct 2021 – Accepted: 12 Oct 2021 – Published: 11 Nov 2021

Natural rock arches are rare and beautiful geologic landforms with important cultural value. As such, their management requires periodic assessment of structural integrity to understand environmental and anthropogenic influences on arch stability. Measurements of passive seismic vibrations represent a rapid and non-invasive technique to describe the dynamic properties of natural arches, including resonant frequencies, modal damping ratios, and mode shapes, which can be monitored over time for structural health assessment. However, commonly applied spectral analysis tools are often limited in their ability to resolve characteristics of closely spaced or complex higher-order modes. Therefore, we investigate two techniques well-established in the field of civil engineering through application to a set of natural arches previously characterized using polarization analysis and spectral peak-picking techniques. Results from enhanced frequency domain decomposition and parametric covariance-driven stochastic subspace identification modal analyses showed generally good agreement with spectral peak-picking and frequency-dependent polarization analyses. However, we show that these advanced techniques offer the capability to resolve closely spaced modes including their corresponding modal damping ratios. In addition, due to preservation of phase information, enhanced frequency domain decomposition allows for direct and convenient three-dimensional visualization of mode shapes. These techniques provide detailed characterization of dynamic parameters, which can be monitored to detect structural changes indicating damage and failure, and in addition have the potential to improve numerical models used for arch stability assessment. Results of our study encourage broad adoption and application of these advanced modal analysis techniques for dynamic analysis of a wide range of geological features.

Natural rock arches form by erosion (Bruthans et al., 2014; Ostanin et al., 2017) and are major tourist attractions worldwide. However, ongoing weathering can lead to partial or complete collapse, posing a hazard to visitors; prominent examples include collapse of London Bridge (Australia) in 1990 (Woodroffe, 2002), rockfall from Landscape Arch (USA) in 1991 and 1995 above a hiking trail (Deseret News, 1991), and collapse of the Azure Window in Malta in 2017 (Satariano and Gauci, 2019). As arches occur in a variety of forms and settings, simple tools for stability assessment do not exist, and current practices often include site-specific geomechanical characterization and numerical modeling (Budetta et al., 2019). In recent decades, the stability of engineered structures, such as buildings and bridges, has been increasingly analyzed using measurements of their vibrational properties associated with resonance. Understanding this dynamic response to ambient loading forms the basis for the field of structural health monitoring (SHM; Doebling et al., 1996). More recently, SHM concepts have been applied at natural rock arches and other geological formations to improve site characterization and hazard assessment associated with failure of these features (e.g., Bottelin et al., 2013; Burjánek et al., 2018; Iannucci et al., 2020; Kleinbrod et al., 2019; Mercerat et al., 2021; Moore et al., 2018). Passive seismic measurements then provide a non-invasive means to monitor dynamic behavior and evaluate stability in the presence of natural or anthropogenic stimuli, which is especially valuable at culturally important sites where more invasive or destructive monitoring techniques (e.g., taking rock samples) may not be permitted.

Passive stability assessment often involves repeated or continuous measurements of a structure to monitor deviations in baseline structural dynamic
behavior. This dynamic behavior is characterized by natural frequencies, corresponding mode shapes (i.e., structural deflection at those frequencies),
and damping ratios (e.g., Chopra, 2015). While damping describes internal energy dissipation and radiation out of the system, resonant frequencies
are primarily a function of stiffness and mass. As the mass of a rock landform is approximately constant over time (in the absence of mass-wasting
events), variations in resonant frequencies arise primarily due to changes in rock mass stiffness, which can be associated with rock damage and
environmental influences, such as temperature and moisture (Colombero et al., 2021; Bottelin et al., 2013). As internal crack growth accumulates
during progressive failure, stiffness and thus frequencies are anticipated to decrease; for example, Lévy et al. (2010) reported a drop in
resonant frequency of about 20 % in less than 2 weeks prior to collapse of a 21 000 m^{3} rock column, which they attributed to
progressive failure. More quantitative assessments of stress conditions prior to failure require individual features to be numerically modeled with
realistic values for rock density and Young's modulus. With density constrained by rock samples, Young's modulus can be derived from dynamic
properties by minimizing the error between observed and modeled resonance attributes (Moore et al., 2018, 2020; Geimer et al., 2020a). Such model
validation facilitates the estimation of the three-dimensional stress field, which was used by Moore et al. (2020) to identify arches that may be more prone to
tensile crack growth and structural failure.

Modal analysis techniques used in structural health monitoring of geological features rely primarily on identification of spectral attributes from in situ ambient vibration data. Power spectra visualizations provide a means for first interpretation, often leading to identification of resonant frequencies that can be confirmed through numerical modeling (Moore et al., 2018), while site-to-reference spectral ratios may be used to eliminate source and path effects in order to identify and track resonant frequencies (e.g., Weber et al., 2018). Selecting the maximum peak directly on the power spectra to determine the resonant frequency is usually referred to as “peak picking”. By additionally selecting the frequencies left and right of the resonant peak, for which the power drops by 3 dB, a simple estimate of the modal damping ratio can be obtained (see Sect. 3.1 for details).

Frequency-dependent polarization analysis (PA) provides a tool to estimate the modal deflection at resonance (Burjánek et al., 2012; Geimer et al., 2020a). However, these spectral analysis techniques fall short when applied to more complex systems, such as cases with closely spaced and overlapping modes, which have identical or similar frequencies but different mode shapes. In addition, phase information is not preserved across separate recording stations, impeding precise determination of mode shapes for higher modes. Thus, new techniques are necessary for refined modal analysis supporting structural health monitoring of rock landforms and providing accurate input parameters for stability assessment using numerical models. Among these, enhanced frequency domain decomposition (EFDD; Brincker et al., 2001a, b) is a promising approach to identify resonant frequencies, damping, and polarization attributes, and it is well-suited to distinguishing closely spaced modes. Covariance-driven stochastic subspace identification (SSI-COV) is an alternative time domain technique that is especially beneficial for accurate estimates of modal damping ratios (van Overschee, 1996). Since their introduction, both techniques have become standard methods for analysis of engineered structures (e.g., Brincker and Ventura, 2015) and have been compared, yielding similar results (Cheynet et al., 2017; Brincker et al., 2000). Using these complementary techniques, Bayraktar et al. (2015) found good agreement between EFDD and SSI-COV in their study on historical masonry arch bridges with resonant frequencies and damping ratios comparable to the natural rock arches and bridges studied here. Furthermore, frequency domain decomposition has been applied to natural features, such as sedimentary valleys, glaciers, and rock slope instabilities (Poggi et al., 2015; Preiswerk et al., 2019; Häusler et al., 2021, 2019; Ermert et al., 2014), while application of SSI-COV has remained restricted to engineered structures.

In this study, we analyze the modal characteristics of four natural rock arches in Utah (USA) previously investigated by Geimer et al. (2020a). As these arches exhibit various spectral complexities which complicate dynamic analyses, we apply two operational modal analysis techniques – EFDD and SSI-COV – to improve identification and characterization of normal modes. Our results highlight the value and versatility of EFDD and SSI-COV for structural characterization and monitoring in geologic hazard applications, which we propose is useful across a broad range of geomorphologic features beyond our studied landforms, such as rock slope instabilities and rock towers (Bottelin et al., 2013; Häusler et al., 2021; Moore et al., 2019).

Ambient vibration data processed in this study were collected at four natural rock arches in Utah by Geimer et al. (2020a). These consist of three
single-station measurements conducted using a Nanometrics Trillium Compact 20 s seismometer (TC 20 s, sites: Rainbow Bridge, Corona
Arch, Squint Arch) and two array measurements using three-component Fairfield Zland 5 Hz nodal geophones with synchronous recording (sites:
Squint Arch and Musselman Arch). Table 1 summarizes the arch measurements, including data acquisition length, site coordinates, and number of sensors
deployed. Prior to processing, all data were corrected using the respective instrument response (to velocity units of m s^{−1}), and the mean
and linear trend were removed.

In the study by Geimer et al. (2020a), Rainbow Bridge showed clear normal modes, although the higher-order modes are partly overlapping (Fig. 1a and b). We include this arch in our study as an example of having well-defined modes. For Corona Arch (Fig. 1c), the numerical models by Geimer et al. (2020a) predicted two modes between 5 and 6 Hz, but only one single peak could clearly be observed in the experimental data (Fig. 1d). Therefore, we include this arch as an example of having a possibly hidden mode. At Squint Arch, the opposite phenomenon was observed: while two peaks were observed in the power spectrum between 10 and 15 Hz, only one mode was predicted by the numerical model in this frequency range (Fig. 1e and f). Finally, the large array dataset acquired at Musselman Arch provides the possibility to test the techniques presented here using a dense sensor array, highlighting the value of retained phase information.

## 3.1 Peak-picking and polarization analyses

In previous studies of the dynamic response of natural rock arches, the resonant frequencies of the landform were determined by selecting the local
maxima of the power spectra of the recordings, which is so-called peak picking (Starr et al., 2015; Moore et al., 2018, e.g., *f*_{1} in Fig. 2a). The
corresponding modal damping ratio can be estimated using the half-power bandwidth technique, whereby the frequencies left and right of the resonant
frequency *f*_{n} are selected (*f*_{A} and *f*_{B}, respectively) as those for which power has decreased to $\mathrm{1}/\sqrt{\mathrm{2}}$ (or
approximately −3 dB; see Fig. 2a). The damping ratio *ζ* is then obtained by

Mode shape information can be retrieved by polarization analysis (PA), for example, using the approach by Koper and Hawley (2010) as applied to rock arches by Moore et al. (2016) and Geimer et al. (2020a). These single-station techniques are easy to use and provide reliable modal parameters in the case of well-separated modes. However, they fall short in the case of closely spaced or overlapping modes, as the mode bells are not visible or do not correspond to the underlying resonant mode (e.g., Papagiannopoulos and Hatzigeorgiou, 2011; Wang et al., 2012). This is illustrated for the example of three single-degree-of-freedom (SDOF; see Appendix A) systems in Fig. 2a: one well-defined mode at 9 Hz is damped at 2 %, whereas two closely spaced modes at 10.5 and 10.7 Hz are damped with 3 % and 4 %, respectively. The superposition of the three SDOF systems provides the resulting response of the multi-degree-of-freedom (MDOF) system (black line in Fig. 2a), which is observed in the power spectrum. Analysis of the well-separated fundamental mode is straightforward, as the peak corresponds to the resonant frequency and applying the half-power bandwidth technique provides the correct damping ratio of 2 %. In contrast, the peaks of both higher modes merge to one single mode bell at 10.5 Hz with an apparent damping estimate of 4.5 %. Therefore, the superposition of the two modes results in broadening of the mode bell and consequently overestimation of damping. Furthermore, the third mode cannot be detected in the power spectra. In addition to damping overestimation by close and hidden modes, the half-power bandwidth technique tends to overestimate damping due to spectral leakage (Seybert, 1981).

The most direct estimate of modal damping ratios is obtained by active source experiments in which the structure studied is excited artificially and
energy dissipation is measured, for example, in the time domain by the logarithmic decrement *δ* (Fig. 2b):

Here, *y* represents the amplitude of the measured quantity (e.g., acceleration, velocity) at time *t* and *T*_{D} refers to the damped
natural period (i.e., the inverse of the resonant frequency). For small damping (< 20 %; Chopra, 2015), this can be approximated and solved
for *ζ*:

Active source experiments can be considered to provide good estimates of damping ratios, but their application is restricted to structures that can be excited artificially (without inducing damage, Magalhães et al., 2010). In contrast, passive (i.e., ambient vibration) experiments can be applied to a broad range of structures but are subject to larger uncertainties (up to 20 % is possible, e.g., Au et al., 2021; Döhler et al., 2013; Gersch, 1974; Griffith and Carne, 2007).

## 3.2 Enhanced frequency domain decomposition

We processed three-component ambient vibration data using enhanced frequency domain decomposition (EFDD) (Brincker et al., 2001a, b; Brincker and Ventura, 2015; Michel et al., 2010). The method first computes the cross-power spectral density between all input traces and for every discrete frequency. Next, singular value decomposition for each frequency provides the singular values and singular vectors. The singular values can be understood as the collection of virtual SDOF systems of the structure, which enables detection of close and hidden modes that are not visible in the power spectrum. The first singular value shows peaks at the dominant natural frequencies of the system. If present, overlapping secondary (i.e., non-dominant) modes result in elevated higher singular values. Resonant frequencies are then determined from analysis of the singular value plot and the singular vector at the identified frequencies gives the three-dimensional modal vector (i.e., mode shape) of the chosen mode, with higher singular vectors representing the mode shape of non-dominating modes. These processing steps represent the frequency domain decomposition method described by Brincker et al. (2001b). The half-power bandwidth technique could now be applied to the singular values to estimate damping, as the bias by modal superposition is now addressed. However, spectral leakage may still broaden the mode bell.

A more accurate technique to estimate modal damping is the enhanced FDD (EFDD) technique introduced by Brincker et al. (2001a). Here the mode bell
around each resonant frequency is picked manually and transformed to the time domain, providing the impulse response function (see Fig. 2b). Energy
decay in the linear part of the impulse response function is expressed by the damping ratio *ζ*, which can be determined using the logarithmic
decrement technique. Linear regression of the zero-crossing times within the linear part of the decay curve additionally provides an updated estimate
of the resonant frequency. The advantage of EFDD over the half-power bandwidth method is that the damping estimate is not based on only three picks
but on a curve fitting approach, which reduces errors introduced by noise. However, EFDD still tends to overestimate damping due to spectral leakage
(e.g., Bajric et al., 2015). A detailed description of the EFDD processing workflow applied in this study can be found in Häusler et al. (2019) and
Häusler et al. (2021), who applied the technique to unstable rock slopes.

## 3.3 Covariance-driven stochastic subspace identification

The second technique used in this study is the covariance-driven stochastic subspace identification (SSI-COV) method (Peeters and De Roeck, 1999; Van Overschee and De Moor, 1993; van Overschee, 1996). Like EFDD, SSI-COV is a modal analysis technique frequently used in civil and mechanical engineering. Contrary to EFDD, SSI-COV is a time domain parametric technique, which searches for the best set of modal parameters (resonant frequencies and modal damping) representing the observed structural response in a mathematical manner, i.e., minimizing misfit between modeled and observed data. Because it is a time domain approach, overestimation of damping from spectral leakage is avoided. The most important processing parameter is the maximum lag time between two time samples used for computing the covariance matrices, which should be 2 to 6 times larger than the longest eigenperiod of the structure. Other user-controlled parameters include the number of possible modes, the accuracy threshold for modal frequency and damping, the maximum spectral distance inside a cluster, and the variation of the minimum modal assurance criterion (e.g., Allemang and Brown, 1982), which is a measure of the similarity of the mode shape at neighboring frequencies. We applied the SSI-COV algorithm software by Cheynet (2020), which is based on the implementation by Magalhães et al. (2009) and was used for comparison to EFDD on long suspension bridges (Cheynet et al., 2017). We followed the parameter suggestions by Cheynet (2020) and chose the parameters in a trial-and-error approach such that the first three resonant modes were reproduced (see Table A1). As SSI-COV establishes a mathematical model of the structure studied, the dynamic response can be defined by poles and zeros (in the sense of mathematics of complex numbers). Therefore, the term “pole” can be used as representative for “resonant mode” and is used hereafter to be in line with the terminology in the field. Since SSI-COV is a parametric method, its resulting resonant frequencies should be verified by a frequency domain technique to prevent misinterpretation by model overfitting.

Results from SSI-COV (and other SSI variants) are illustrated using stability diagrams (e.g., Fig. 3c). Initially, the structure's response is modeled with a low number of modes (poles), which is continuously increased to the maximum number of poles defined by the user. The maximum number of poles should be chosen to be significantly larger than the expected number of modes in order to establish an overdetermined mathematical model. The resulting resonant frequencies for each mode at every model run are plotted in the stability diagram (blue crosses in Fig. 3c). Repeated poles, i.e., identical or very similar values for resonant frequencies, damping, and mode shape, represent stable poles and can be identified as vertical stacks of poles in the stability diagram (red circles in Fig. 3c). Poles not fulfilling the user-defined accuracy criteria are not interpreted as stable poles and are scattered at arbitrary values as a result of noise fitting. Stable poles are clustered using hierarchical clustering, thereby grouping poles with similar characteristics to the final resonant modes of the structure.

We observe the first three resonant modes of Rainbow Bridge determined by the single-station measurement at 1.1, 2.2, and 2.5 Hz (Fig. 3a–c).
While the fundamental mode (*f*_{1} at 1.1 Hz) is distinctly separated from other spectral peaks, the second and third modes (*f*_{2},
*f*_{3}) occur at closely spaced frequencies but are clearly identified by the elevated second singular value. Damping is estimated at between
0.6 % and 1.3 % for all three modes (Table 1). For the fundamental mode, we estimate damping at 0.9 % and 0.6 % using EFDD and
SSI-COV, respectively, which is significantly lower than estimated by Geimer et al. (2020a) using the half-power bandwidth method (2.4 %). Modal
vectors (i.e., azimuth and incidence angle with a lower hemisphere projection) derived by EFDD are similar to the PA results of Geimer et al. (2020a)
with some minor differences for *f*_{3} and a 180^{∘} ambiguity in the azimuth of the nearly horizontally polarized mode *f*_{2}. Note that
Geimer et al. (2020a) allowed a polarity flip for mode shapes with sub-horizontal incidence angles equal to or larger than 85^{∘} in order to
compare to numerical models. As SSI-COV and EFDD provide similar results, we only compare values from EFDD to PA in Table 2 and provide SSI-COV
results in Table A2.

The singular value plot of Corona Arch reveals two distinct spectral peaks at 2.7 and 5.3 Hz (Fig. 3e). However, the second singular value
also peaks at ∼ 5.3 Hz, indicating the presence of a closely spaced mode at that frequency. Therefore, we confirm the two close modes
proposed by Geimer et al. (2020a). However, EFDD and SSI-COV suggest nearly identical frequencies for *f*_{2} and *f*_{3} (5.3 Hz), while Geimer
et al. (2020a) selected more separated frequencies (5.0 and 5.4 Hz) based on numerical modeling and PA. For EFDD, the mode shape vector of the
non-dominant mode *f*_{3} is determined by the second singular vector. Modal vectors for *f*_{1} and *f*_{2} resolved by PA and EFDD are in good
agreement; however, azimuth and incidence differ for *f*_{3}. While EFDD and SSI-COV gave similar values for incidence of *f*_{3} (44 and 54^{∘},
respectively), PA estimated incidence at 73^{∘}. Damping is estimated between 0.9 % and 2.0 % for all three modes, with 0.9 % and
1.4 % for the fundamental mode (via SSI-COV and EFDD, respectively). These values are again slightly lower than the half-power bandwidth estimates
of Geimer et al. (2020a, 1.9 %). Damping ratios for *f*_{2} and *f*_{3} are between 1.5 % and 2.0 %, with EFDD and SSI-COV providing similar
damping ratios within the expected uncertainty range.

For Squint Arch, we observe two closely spaced modes at 11.5 and 12.5 Hz, as well as a third mode at 19.9 Hz (Fig. 3h and i). Geimer
et al. (2020a) interpreted the first two peaks as one mode as it could not be confirmed as a separate mode by numerical models. Our analysis of the
second mode suggests the modal vector has a steeper incidence angle compared to the first mode (49^{∘}) and is therefore oriented 60^{∘}
from *f*_{1}. If the two spectral peaks are analyzed separately by PA, the match between PA and EFDD is very good (see values in brackets in
Table 2). The damping ratio of 1.6 % determined by Geimer et al. (2020a) is in perfect agreement with the estimation by SSI-COV but differs
slightly from the EFDD result (2.4 %).

We demonstrate the ability of EFDD to retrieve the full-length normal-mode shapes at Squint Arch, where data acquired by a nodal geophone array during
a separate experiment are available (Fig. 4b; raw power spectra are shown in Fig. A1 a–c). We note that modal frequencies for *f*_{1} and *f*_{2}
increased by about 1 Hz compared to the single-station measurement. No other higher modes are visible on the singular value plot during this
measurement. Modal vectors for the first two modes at all stations resulting from EFDD analysis are shown in Fig. 4c and d. We were not able to define
a set of SSI-COV parameters that could successfully reproduce the observed modes.

We performed EFDD and SSI-COV modal analyses on geophone array data acquired at Musselman Arch, revealing the first four resonant modes at 3.4, 4.2, 5.6, and 6.6 Hz (Fig. 5b and c; raw power spectra are shown in Fig. A1 d and f). The resonant frequencies and mode shapes are in good agreement with results of the single-component cross-correlation analysis by Geimer et al. (2020a). Visualization of the three-dimensional modal vectors for each station determined by the first singular vector derived by EFDD is shown in Fig. 5d–g. The first two modes are full-span, first-order bending modes in the horizontal and vertical directions. The third mode is a nearly symmetric second-order vertical bending mode with the node point at the center of the arch. Mode four is a slightly asymmetric second-order horizontal bending mode with the node point shifted towards the eastern abutment. Modal damping ratios for the first three modes are estimated at 1.3 %, 1.0 %, and 1.9 % with EFDD and 1.3 %, 1.1 %, and 1.6 % with SSI-COV. Note that the second singular value is elevated at each of the resonant frequencies. However, in the case of Musselman Arch, this phenomenon is caused by an anomalous sensor component and not by close or hidden modes, as discussed in Sect. 5.1.

## 5.1 Modal identification

At Rainbow Bridge, the first three spectral peaks are well-separated. Therefore, simple peak picking on the power spectra provides the same resonant
frequencies as obtained by EFDD and SSI-COV. All techniques provide comparable values for azimuth and incidence of the modal vector. The largest
discrepancies between EFDD/SSI-COV and PA are 6 and 4^{∘} for the azimuth and incidence, respectively. These discrepancies are in the range of the
misfit between field observations and the numerical models presented by Geimer et al. (2020a). Therefore, we conclude that, for well-separated modes,
all techniques provide similar results within the range of uncertainties.

At Corona Arch, only two resonant peaks are observed in the power spectra between 1 and 7 Hz. The first includes both horizontal components, while the second also includes the vertical component. The numerical model by Geimer et al. (2020a) predicted two close modes at the location of the second peak, which supported their interpretation of a close hidden mode. However, the same model also showed significant misfit between the model prediction and observed data. With EFDD and SSI-COV, the presence of two close modes can be verified (note the elevated second singular value in Fig. 4e). For the example of Corona Arch, we demonstrate that EFDD and SSI-COV are valuable techniques to detect close and hidden modes.

At Squint Arch, the power spectra show two resonant peaks between 10 and 14 Hz, both including all three components. However, the numerical models by Geimer et al. (2020a) only predicted one mode. Consequently, they interpreted this doublet peak as the signature of one mode alone. In contrast, EFDD and SSI-COV independently indicate two closely spaced modes. Therefore, these two techniques help to identify the resonant modes at Squint Arch.

At Musselman Arch, all four resonant modes between 1 and 10 Hz are observed in the power spectra, as the resonant frequencies are well-separated. However, the large array dataset results in 3 × 32 power spectra to be analyzed, which requires extensive effort (Fig. A1). In addition, the power spectra do not provide direct evidence of whether the peaks correspond to one single mode or if there are additional close modes. Here, EFDD combines all input traces in one single plot, providing a direct illustration of the resonant modes and indicating that no close modes are present. In addition, the mode shapes can be directly plotted by evaluating the singular vectors. Therefore, the EFDD analysis of the large Musselman Arch dataset provides a demonstration of the user-friendliness and simplicity of the EFDD technique.

We observed that the second singular value at Musselman Arch is elevated to about 85 dB at each of the resonant peaks, which corresponds to
the value of the first singular value at frequencies at which no resonance occurs (Fig. 6a). The difference in the overall noise floor, i.e., to the
level at which the singular values remain flat (∼ 70 dB), is about 15 dB. Such nearly constant gaps between singular values across
a broad frequency range are indicative of an erroneous sensor component or sensor coupling issue, which increases the noise floor to the noise level of
the faulty component (e.g., Brincker, 2014). In our survey, the anomalous component was identified to be the transverse horizontal component of
sensor 000, which is evident from the elevated power spectra in Fig. A1d and the anomalous mode shape vector in Fig. A2g. By omitting sensor 000 in
our EFDD analysis, the overall noise floor drops to 70 dB and removes the elevated higher singular values at the third and fourth resonant
frequency (Fig. 6b). The second singular value is still elevated between the first two resonant modes, indicating that these modes have an overlapping
bandwidth. In contrast to Musselman Arch, the higher singular values at Corona Arch were only elevated at the second resonant frequency and not at the
first. In addition, the peak at the second singular value (*f*_{3}) exceeds the noise floor of the system, i.e., power of the first singular value
outside the mode bells. This shows that the hidden mode *f*_{3} at Corona Arch is indeed a structural mode and not an artifact induced by an erroneous
sensor component.

## 5.2 Damping estimates

For most resonant modes, EFDD and SSI-COV provide comparable damping results within the anticipated range of uncertainty. However, we observe that EFDD results in 30 % to 35 % higher damping ratios for the fundamental modes of Rainbow Bridge, Corona Arch, and Squint Arch, as well as for the first higher mode of Squint Arch. We interpret this observation as an effect of damping overestimation through broadening of the resonant peak caused by spectral leakage (e.g., Bajric et al., 2015).

Damping ratios obtained by the half-power bandwidth technique at Rainbow Bridge and Corona Arch are 75 % and 53 % higher than those estimated by SSI-COV, and they are 63 % and 26 % higher than estimated by EFDD. This is illustrated in Fig. 7, where we compare damping ratios obtained by the various techniques for each arch. The resonant frequency and damping ratio derived are used to model an SDOF system, which is superimposed on the singular value plots. The amplitude of the modeled SDOF is normalized to the maximum amplitude of the first singular value. The mode bell of the fundamental mode of Rainbow Bridge is well-reproduced by EFDD and SSI-COV, but damping is overestimated by the half-power bandwidth technique (Fig. 7a).

Rainbow Bridge had the lowest fundamental frequency in the study by Geimer et al. (2020a), who used the same settings to compute the power spectra for all arches, including arches with higher resonant frequencies such as Squint Arch. Therefore, it is likely that these parameter settings were not ideal to resolve the low resonant frequencies of Rainbow Bridge with sufficient resolution. Therefore, we interpret the discrepancy between the half-power bandwidth and EFDD and SSI-COV as a result of strong spectral leakage for the half-power bandwidth technique. However, neither SSI-COV nor EFDD is able to perfectly reproduce the mode bell due to its slightly asymmetric shape, likely reflecting the oversimplified assumption of an SDOF system. At Corona Arch, SSI-COV is capable of reproducing the mode bell of the fundamental mode, and the half-power bandwidth again overestimates damping (Fig. 7b). SSI-COV and the half-power bandwidth technique yielded identical damping ratios for the fundamental mode of Squint Arch (Fig. 7c). For higher-order modes, the discrepancy between EFDD and SSI-COV is smaller on all arches.

## 5.3 Differences between surveys at Squint Arch

The normal-mode analysis of Squint Arch resulted in different resonant frequencies for the single-station broadband measurement in February 2018 and
the geophone array measurement in April 2018. We attribute this shift in frequency to seasonal variations and mainly temperature differences
(11.5 ^{∘}C for the single-station and 16 ^{∘}C for the nodal measurement; see also Starr et al.,
2015). Seasonal effects are also expected to influence modal damping ratios and mode shapes (Häusler et al., 2021). Another difference between the
two surveys is that only two modes can be detected by the geophone array. This is likely an effect of the higher self-noise of the Zland geophones
compared to the broadband TC 20 s seismometer, which may be greater than the excitation level of the higher modes (e.g., Brincker and Larsen,
2007). We were also not able to find a set of SSI-COV parameters that could reliably reproduce the resonant frequencies, which is again attributed to
the lower signal-to-noise ratio of the geophone data (e.g., Brincker, 2014; Liu et al., 2019).

We applied enhanced frequency domain decomposition (EFDD) and covariance-driven stochastic subspace identification (SSI-COV) modal analyses to a set of four natural rock arches previously analyzed by Geimer et al. (2020a) using frequency-dependent polarization analysis (PA). Our results show that EFDD and SSI-COV are able to determine the natural frequencies, damping ratios, and mode shapes of these landforms, including close, hidden, and higher resonant modes. For well-separated resonant modes, these techniques reproduce the results by Geimer et al. (2020a). In the case of more complex spectra, EFDD and SSI-COV are able to extract additional modal details not resolved with PA. EFDD facilitated identification and interpretation of closely spaced (i.e., spectrally overlapping) and hidden modes at Corona and Squint arches. EFDD additionally combines information for all input traces in a single plot, allowing rapid analysis of the dynamic response, especially when compared to picking the resonant peaks and determining polarization information on every station spectrum individually. The singular vectors resulting from EFDD can be directly interpreted as the three-dimensional modal deflection vector at each station, providing rapid and convenient visualization of mode shapes.

While modal analysis via peak picking and subsequent PA has been shown to be satisfactory for adequately spaced spectral peaks and strongly amplified resonant frequencies, here we demonstrate that more sophisticated modal analysis techniques can provide refined modal characterization for more complex dynamic systems. Improving the accuracy and our understanding of resonance properties could in turn help generate more refined numerical models, facilitating more accurate arch stability assessment. Future efforts in modeling the dynamic response of rock arches (and other geological features) should additionally involve calibration of the modal damping ratio, as we have shown that this parameter can be measured on complex structures.

Our results encourage widespread application of EFDD and SSI-COV modal analysis techniques, which are commonly used in civil engineering, to complement existing seismological techniques for dynamic analysis of geological features. Both techniques are well-suited for future near-real-time monitoring of the structural integrity of geological features beyond rock arches: for example, rock slope instabilities, unstable glaciers, and freestanding rock towers. EFDD and SSI-COV are only two methods out of many other available algorithms for modal analysis, including other SSI variants and curve fit FDD (Peeters and De Roeck, 2001; Jacobsen, 2008). Therefore, future research could explore the potential of these techniques for applications involving modal analyses and monitoring of Earth surface landforms.

The frequency response *H*(*ω*) of a single-degree-of-freedom (SDOF) system is given by

with *ω* being the angular frequency and *ω*_{n} being the angular resonant frequency. *ζ* refers to the modal damping ratio
and *i* is the imaginary unit (see, for example, Chopra, 2015).

Data used in this study and originating from the study by Geimer et al. (2020b) are available at https://doi.org/10.7278/S50D-G31E-NFW2. Data for the nodal deployment at Squint Arch are available at https://doi.org/10.7914/SN/5P_2013 (Moore, 2013).

The paper was written by MH with significant contributions from all co-authors. PG, RF, and JM acquired seismic data. MH carried out data processing and software development of EFDD modal analysis. PG performed data curation and validation of results. All authors reviewed and approved the paper.

The contact author has declared that neither they nor their co-authors have any competing interests.

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

This research has been supported by the Eidgenössische Technische Hochschule Zürich (grant no. 0-20361-17) and the US National Science Foundation under grant no. EAR-1424896.

This paper was edited by Jens Turowski and reviewed by two anonymous referees.

Allemang, R. J. and Brown, D. L.: A Correlation Coefficient for Modal Vector Analysis, in: Proceedings of the 1st International Modal Analysis Conference, Orlando, USA, 1982.

Au, S.-K., Brownjohn, J. M. W., Li, B., and Raby, A.: Understanding and managing identification uncertainty of close modes in operational modal analysis, Mech. Syst. Signal Pr., 147, 107018, https://doi.org/10.1016/j.ymssp.2020.107018, 2021.

Bajric, A., Brincker, R., and Thöns, S.: Evaluation of damping estimates in the presence of closely spaced modes using operational modal analysis techniques, in: Proceedings of the 6th International Operational Modal Analysis Conference, Gijón, Spain, 2015.

Bayraktar, A., Türker, T., and Altunişik, A. C.: Experimental frequencies and damping ratios for historical masonry arch bridges, Const. Build. Mater., 75, 234–241, https://doi.org/10.1016/j.conbuildmat.2014.10.044, 2015.

Bottelin, P., Lévy, C., Baillet, L., Jongmans, D., and Guéguen, P.: Modal and thermal analysis of Les Arches unstable rock column (Vercors massif, French Alps), Geophys. J. Int., 194, 849–858, https://doi.org/10.1093/gji/ggt046, 2013.

Brincker, R.: Some Elements of Operational Modal Analysis, Shock Vib., 2014, 325839, https://doi.org/10.1155/2014/325839, 2014.

Brincker, R. and Larsen, J. A.: Obtaining and Estimating Low Noise Floors in Vibration Sensors, in: Proceedings of the 24th International Modal Analysis Conference, Orlando, USA, 2007.

Brincker, R. and Ventura, C.: Introduction to operational modal analysis, John Wiley & Sons Inc, Chichester, UK, 360 pp., 2015.

Brincker, R., Frandsen, J. B., and Andersen, P.: Ambient Response Analysis of the Great Belt Bridge, in: Proceedings of the 18th International Modal Analysis Conference, San Antonio, USA, 2000.

Brincker, R., Ventura, C., and Andersen, P.: Damping estimation by frequency domain decomposition, in: Proceedings of the 19th International Modal Analysis Conference, Orlando, USA, 2001a.

Brincker, R., Zhang, L., and Andersen, P.: Modal identification of output-only systems using frequency domain decomposition, Smart Mater. Struct., 10, 441–445, https://doi.org/10.1088/0964-1726/10/3/303, 2001b.

Bruthans, J., Soukup, J., Vaculikova, J., Filippi, M., Schweigstillova, J., Mayo, A. L., Masin, D., Kletetschka, G., and Rihosek, J.: Sandstone landforms shaped by negative feedback between stress and erosion, Nat. Geosci., 7, 597–601, https://doi.org/10.1038/ngeo2209, 2014.

Budetta, P., De Luca, C., Simonelli, M. G., and Guarracino, F.: Geological analysis and stability assessment of a sea arch in Palinuro, southern Italy, Eng. Geol., 250, 142–154, https://doi.org/10.1016/j.enggeo.2019.01.009, 2019.

Burjánek, J., Moore, J. R., Yugsi Molina, F. X., and Fäh, D.: Instrumental evidence of normal mode rock slope vibration, Geophys. J. Int., 188, 559–569, https://doi.org/10.1111/j.1365-246X.2011.05272.x, 2012.

Burjánek, J., Gischig, V., Moore, J. R., and Fäh, D.: Ambient vibration characterization and monitoring of a rock slope close to collapse, Geophys. J. Int., 212, 297–310, https://doi.org/10.1093/gji/ggx424, 2018.

Cheynet, E.: Operational modal analysis with automated SSI-COV algorithm, Zenodo [code], https://doi.org/10.5281/ZENODO.3774061, 2020.

Cheynet, E., Jakobsen, J. B., and Snæbjörnsson, J.: Damping estimation of large wind-sensitive structures, Procedia Engineer., 199, 2047–2053, https://doi.org/10.1016/j.proeng.2017.09.471, 2017.

Chopra, A. K.: Dynamics of structures: theory and applications to earthquake engineering, fourth edn., Prentice-Hall international series in civil engineering and engineering mechanics, Pearson Prentice Hall, Boston, 2015.

Colombero, C., Jongmans, D., Fiolleau, S., Valentin, J., Baillet, L., and Bièvre, G.: Seismic Noise Parameters as Indicators of Reversible Modifications in Slope Stability: A Review, Surv. Geophys., 42, 339–375, https://doi.org/10.1007/s10712-021-09632-w, 2021.

DeseretNews: Slabs fall from landscape arch, available at: https://www.deseret.com/1991/9/7/18939827/slabs-fall-from-landscape-arch (last access: 10 December 2020), 1991.

Doebling, S. W., Farrar, C. R., Prime, M. B., and Shevitz, D. W.: Damage identification and health monitoring of structural and mechanical systems from changes in their vibration characteristics: A literature review, Technical Report, LANL – Los Alamos National Laboratory, Los Alamos, USA, https://doi.org/10.2172/249299, 1996.

Döhler, M., Hille, F., Mevel, L., and Rücker, W.: Estimation of modal parameters and their uncertainty bounds from subspace-based system identification, in: IRIS Industrial Safety and Life Cycle Engineering – Technologies/Standards/Applications, edited by: Margit, K., VCE, Vienna, Austria, 91–106, 2013.

Ermert, L., Poggi, V., Burjánek, J., and Fäh, D.: Fundamental and higher two-dimensional resonance modes of an Alpine valley, Geophys. J. Int., 198, 795–811, https://doi.org/10.1093/gji/ggu072, 2014.

Geimer, P. R., Finnegan, R., and Moore, J. R.: Sparse Ambient Resonance Measurements Reveal Dynamic Properties of Freestanding Rock Arches, Geophys. Res. Lett., 47, e2020GL087239, https://doi.org/10.1029/2020GL087239, 2020a.

Geimer, P. R., Finnegan, R., and Moore, J. R.: Data for: Sparse Ambient Resonance Measurements Reveal Dynamic Properties of Freestanding Rock Arches, University of Utah Research Data Repository (Hive) [data set], https://doi.org/10.7278/S50D-G31E-NFW2, 2020b.

Gersch, W.: On the achievable accuracy of structural system parameter estimates, J. Sound Vib., 34, 63–79, https://doi.org/10.1016/S0022-460X(74)80355-X, 1974.

Griffith, D. T. and Carne, T. G.: Experimental Uncertainty Quantification of Modal Test Data, in: Proceedings of the 25th International Modal Analysis Conference, Orlando, USA, 2007.

Häusler, M., Michel, C., Burjánek, J., and Fäh, D.: Fracture Network Imaging on Rock Slope Instabilities Using Resonance Mode Analysis, Geophys. Res. Lett., 46, 6497–6506, https://doi.org/10.1029/2019GL083201, 2019.

Häusler, M., Michel, C., Burjánek, J., and Fäh, D.: Monitoring the Preonzo rock slope instability using resonance mode analysis, J. Geophys. Res.-Earth, 126, e2020JF005709, https://doi.org/10.1029/2020JF005709, 2021.

Iannucci, R., Martino, S., Paciello, A., D'Amico, S., and Galea, P.: Investigation of cliff instability at Ghajn Hadid Tower (Selmun Promontory, Malta) by integrated passive seismic techniques, J. Seismol., 24, 897–916, https://doi.org/10.1007/s10950-019-09898-z, 2020.

Jacobsen, N.-J. and Andersen, P.: Operational Modal Analysis on Structures with Rotating Parts, in: Proceedings of the International Conference on Noise and Vibration Engineering, Leuven, Belgium, 2008.

Kleinbrod, U., Burjánek, J., and Fäh, D.: Ambient vibration classification of unstable rock slopes: A systematic approach, Eng. Geol., 249, 198–217, https://doi.org/10.1016/j.enggeo.2018.12.012, 2019.

Koper, K. D. and Hawley, V. L.: Frequency dependent polarization analysis of ambient seismic noise recorded at a broadband seismometer in the central United States, Earthquake Science, 23, 439–447, https://doi.org/10.1007/s11589-010-0743-5, 2010.

Lévy, C., Baillet, L., Jongmans, D., Mourot, P., and Hantz, D.: Dynamic response of the Chamousset rock column (Western Alps, France), J. Geophys. Res.-Earth, 115, F04043, https://doi.org/10.1029/2009JF001606, 2010.

Liu, F., Wu, J., Gu, F., and Ball, A. D.: An Introduction of a Robust OMA Method: CoS-SSI and Its Performance Evaluation through the Simulation and a Case Study, Shock Vib., 2019, 6581516, https://doi.org/10.1155/2019/6581516, 2019.

Magalhães, F., Cunha, Á., and Caetano, E.: Online automatic identification of the modal parameters of a long span arch bridge, Mech. Syst. Signal Pr., 23, 316–329, https://doi.org/10.1016/j.ymssp.2008.05.003, 2009.

Magalhães, F., Cunha, Á., Caetano, E., and Brincker, R.: Damping estimation using free decays and ambient vibration tests, Mech. Syst. Signal Pr., 24, 1274–1290, https://doi.org/10.1016/j.ymssp.2009.02.011, 2010.

Mercerat, E. D., Payeur, J. B., Bertrand, E., Malascrabes, M., Pernoud, M., and Chamberland, Y.: Deciphering the dynamics of a heterogeneous sea cliff using ambient vibrations: case study of the Sutta–Rocca overhang (southern Corsica, France), Geophys. J. Int., 224, 813–824, https://doi.org/10.1093/gji/ggaa465, 2021.

Michel, C., Guéguen, P., Lestuzzi, P., and Bard, P. Y.: Comparison between seismic vulnerability models and experimental dynamic properties of existing buildings in France, B. Earthq. Eng., 8, 1295–1307, https://doi.org/10.1007/s10518-010-9185-7, 2010.

Moore, J. R.: Structural health monitoring of rock arches and towers, International Federation of Digital Seismograph Networks (FDSN) [data set], https://doi.org/10.7914/SN/5P_2013, 2013.

Moore, J. R., Thorne, M. S., Koper, K. D., Wood, J. R., Goddard, K., Burlacu, R., Doyle, S., Stanfield, E., and White, B.: Anthropogenic sources stimulate resonance of a natural rock bridge, Geophys. Res. Lett., 43, 9669–9676, https://doi.org/10.1002/2016GL070088, 2016.

Moore, J. R., Geimer, P. R., Finnegan, R., and Thorne, M. S.: Use of Seismic Resonance Measurements to Determine the Elastic Modulus of Freestanding Rock Masses, Rock Mech. Rock Eng., 51, 3937–3944, https://doi.org/10.1007/s00603-018-1554-6, 2018.

Moore, J. R., Geimer, P. R., Finnegan, R., and Michel, C.: Dynamic Analysis of a Large Freestanding Rock Tower (Castleton Tower, Utah), Short Note, B. Seismol. Soc. Am., 109, 2125–2131, https://doi.org/10.1785/0120190118, 2019.

Moore, J. R., Geimer, P. R., Finnegan, R., and Bodtker, J.: Between a beam and catenary: Influence of geometry on gravitational stresses and stability of natural rock arches, Geomorphology, 364, 107244, https://doi.org/10.1016/j.geomorph.2020.107244, 2020.

Ostanin, I., Safonov, A., and Oseledets, I.: Natural Erosion of Sandstone as Shape Optimisation, Sci. Rep.-UK, 7, 17301, https://doi.org/10.1038/s41598-017-17777-1, 2017.

Papagiannopoulos, G. A. and Hatzigeorgiou, G. D.: On the use of the half-power bandwidth method to estimate damping in building structures, Soil Dyn. Earthq. Eng., 31, 1075–1079, https://doi.org/10.1016/j.soildyn.2011.02.007, 2011.

Peeters, B. and De Roeck, G.: Reference-based stochastic subspace identification for output-only modal analysis, Mech. Syst. Signal Pr., 13, 855–878, https://doi.org/10.1006/mssp.1999.1249, 1999.

Peeters, B. and De Roeck, G.: Stochastic System Identification for Operational Modal Analysis: A Review, J. Dyn. Syst.-T ASME, 123, 659–667, https://doi.org/10.1115/1.1410370, 2001.

Poggi, V., Ermert, L., Burjanek, J., Michel, C., and Fäh, D.: Modal analysis of 2-D sedimentary basin from frequency domain decomposition of ambient vibration array recordings, Geophys. J. Int., 200, 615–626, https://doi.org/10.1093/gji/ggu420, 2015.

Preiswerk, L. E., Michel, C., Walter, F., and Fäh, D.: Effects of geometry on the seismic wavefield of Alpine glaciers, Ann. Glaciol., 60, 112–124, https://doi.org/10.1017/aog.2018.27, 2019.

Satariano, B. and Gauci, R.: Landform Loss and Its Effect on Health and Well-being: The Collapse of the Azure Window (Gozo) and the Resultant Reactions of the Media and the Maltese Community, in: Landscapes and Landforms of the Maltese Islands, edited by: Gauci, R. and Schembri, J. A., Springer International Publishing, Cham, 289–303, 2019.

Seybert, A. F.: Estimation of damping from response spectra, J. Sound Vib., 75, 199–206, https://doi.org/10.1016/0022-460X(81)90339-4, 1981.

Starr, A. M., Moore, J. R., and Thorne, M. S.: Ambient resonance of Mesa Arch, Canyonlands National Park, Utah, Geophys. Res. Lett., 42, 6696–6702, https://doi.org/10.1002/2015GL064917, 2015.

van Overschee, P.: Subspace Identification for Linear Systems: Theory – Implementation – Applications, 1st edn., edited by: de Moor, B. L., Springer US, New York, NY, 1996.

Van Overschee, P. and De Moor, B.: Subspace algorithms for the stochastic identification problem, Automatica, 29, 649–660, https://doi.org/10.1016/0005-1098(93)90061-W, 1993.

Wang, J.-T., Jin, F., and Zhang, C.-H.: Estimation error of the half-power bandwidth method in identifying damping for multi-DOF systems, Soil Dyn. Earthq. Eng., 39, 138–142, https://doi.org/10.1016/j.soildyn.2012.02.008, 2012.

Weber, S., Fäh, D., Beutel, J., Faillettaz, J., Gruber, S., and Vieli, A.: Ambient seismic vibrations in steep bedrock permafrost used to infer variations of ice-fill in fractures, Earth Planet. Sc. Lett., 501, 119–127, https://doi.org/10.1016/j.epsl.2018.08.042, 2018.

Woodroffe, C. D.: Coasts: form, process and evolution, Cambridge University Press, Cambridge, 2002.