Variation in Rayleigh wave ellipticity as a possible indicator of earthflow mobility: a case study of Sobradinho landslide compared with pile load testing

Rainfall-induced landslides pose a significant risk to communities and infrastructures. To improve the prediction of such events, it is imperative to adequately investigate the rainfall-dependent dynamics (leading to fluidization) and any associated internal sliding along shear planes within clayey slopes. Therefore, the present study adopted ambient noise analysis based on the Horizontal to Vertical Spectral Ratio (HVSR) method, to measure the seasonal variation of Rayleigh wave ellipticity as an indicator for the internal deformation and transition in the material state. The methodology was applied to an existing landslide, where variations in soil stiffness and internal sliding were expected to occur in response to rainfall. To improve the interpretation of the HVSR results (and hence the prediction of landslide’ reactivation by rainfall), HVSR measurements were also conducted on a field-scale pile load test. The pile test allowed a comparison of the seismic data generated by the soil movement along shear planes. The HVSR curves of this field test showed two frequency peaks with no changes in the resonance. In comparison with the data obtained from the landslide, the resultant HVSR curves showed three frequency patterns: ubiquitous (2Hz), landslide (4-8Hz), and flat (no peak). However, the HVSR curves did not show any response to the expected seasonally induced variations in the landslide mass because of the relatively short data acquisition. Nevertheless, time-lapse HVSR is a promising technique that can complement other geophysical methods for improving landslide monitoring.


Introduction
Landslides pose a threat to the global terrestrial hazards affecting many countries, for example: the Sidoarjo mudflow in Indonesia in 2006 (Wibisana, 2013), the Beichuan landslide in Chania in 2008 (Tang et al., 2011), and the Montecito mudflows in California in 2018 (Donnellan et al., 2018). Other examples include the 2008 Mw 7.9 Wenchuan earthquake causing about 200,000 coseismic landslides in China (Yunus et al., 2020), the 2004 Mw 6.8 Chuetsu and the 2018 Mw 6.6 Hokkaido earthquake-induced over 10000 landslides in Japan (Dou et al., 2020), respectively. Also, dense rainfall in mountainous terrain can seriously trigger numerous landslides in hill slopes in different countries and regions, such as China, Taiwan, Japan, Vietnam, etc. (Chang et al., 2019;Li et al., 2020;Pham et al., 2019, Fan et al., 2018. Due to the substantial human and economic losses accounted for the disastrous landslides, scientists worldwide focus on the fact of monitoring locations of future landslide events for mitigating such disasters. In the case of Brazil, a significant number of shallow rainfall-triggered landslides have occurred (Hussain et al., 2017;Rosi et al., 2019). From 1988 to 2015, there were 733 landslides events in 243 Brazilian municipalities, with 4,000 fatalities. Economic losses and homelessness are the other social impacts of these hazardous events in the states of Rio de Janeiro, São Paulo, Minas Gerais, and Santa Catarina. The tragedies of 1966The tragedies of , 1967The tragedies of , 1985The tragedies of , 1988The tragedies of , 1995The tragedies of , 2008The tragedies of , 2009The tragedies of , 2010The tragedies of , 2011The tragedies of and 2014 are historical reminders of the atrocities caused by landslides in Rio de Janeiro (Martins et al., 2018, Hussain et al., 2019a. The Sobradinho landslide (considered in this study) is an example of such events, which has superficially developed in lacustrine sediments, and whose movements were mostly controlled by water infiltration and fluvial dynamics (Hussain et al., 2017;Rosi et al., 2019).
Landslides involving sudden mobility of earthflows have been associated with two mutually inclusive mechanisms: (i) excess pore-water pressures generated along shear plane (Baum et al., 2003;van Asch and Malet, 2009;Hussain et al., 2019b) and (ii) loss of shear strength because of the increased degree of saturation (Pastor et al., 2010;Jongmans et al., 2015;Hamza and Bellis, 2008;Zhou et al., 2013Zhou et al., , 2016. Accordingly, the mechanical behavior of the earthflows is considered a grouping of Coulomb plastic solids where the mobility is associated with distributed internal shearing that causes flow-like appearance instead of mass liquefaction (Hungr et al., 2001).
The pore water pressure assumption explains the mobility of landslide in terms of the reduction of shear strength along a slip surface (shear plane). However, in the case of clayey material, the undrained shear strength (cohesion) is not expected to reach zero. In such cases, the mobility can be better explained as a transition of a material state from solid to plastic and from plastic to fluid (i.e. solid-plastic-fluid) combined with internal sliding along shear planes. This process is called fluidization, which differs from liquefaction (Lowe, 1976;Bertello et al., 2018). However, the identification of such material transition and internal sliding is a challenging task, which requires continuous monitoring in order to improve landslides' prediction. For the monitoring, a continuous source of energy at high spatiotemporal resolution is required which can be provided by ambient noise (vibrations in the absence of any single source such as an earthquake).
The time-lapse spectral analysis and ambient noise interferometry (ANI), which measure Rayleigh wave ellipticity and Rayleigh wave velocity, respectively, are the two commonly adopted ambient noise-based techniques for the dynamic monitoring of landslides. In theory, both techniques allow time-lapse monitoring to identify temporal changes of these parameters (ellipticity and wave velocity) (Mainsant et al., 2012;Köhler and Weidle, 2019). However, very few attempts have been conducted in the literature, where seismic waves (Rayleigh wave) are utilized for the demarcation of such transition. In these studies (Mainsant et al., 2012;Bertello et al., 2018), a substantial drop of Rayleigh wave velocity has been observed in response to the earthflow acceleration, which is followed by a slow return to pre-disturbance of Rayleigh velocities as the landslide decelerates. However, these studies applied expensive, time-consuming and labor extensive techniques for the estimation of Rayleigh wave velocities such as time-lapse interferometry (Mainsant et al., 2012;Hussain et al., 2019c) and time-lapse active Multi-Channel Analysis of Surface Waves (MASW) (Bertello et al., 2018).
Rayleigh wave ellipticity calculated from an HVSR curve (Horizontal to Vertical Spectral Ratio) provides an alternative choice. Recently, few studies have interpreted seasonal changes in HVSR ellipticity curves as a result of the thaw-freeze cycles of an active permafrost layer (Abbott et al., 2016;Kula et al., 2018;Köhler and Weidle, 2019) and shallow rainfall-triggered landslide (Hussain et al., 2019c). The results of these studies have suggested that surface wave ellipticity (HVSR peak) might constitute a good indicator for the internal deformation and transition in material state (fluidization process) and can be used as a precursor for predicting earthflow. Therefore, the HVSR method could bear the potential to become a low-cost, passive, and non-invasive method for long-term monitoring and site characterization of landslides with high temporal resolution.
The aim of the present study is to improve the prediction of rainfallinduced landslides by understanding the rainfall-dependent dynamics (leading to fluidization) and any associated internal sliding along shear planes within clayey slopes. The study applied the HVSR monitoring on Sobradinho landslide, where variations in soil stiffness and internal sliding (fluidization process) would be expected to occur in response to rainfall events. To improve the interpretation of the HVSR results (and hence the prediction of landslide reactivation by rainfall), HVSR was also conducted to monitor a field-scale pile load test. This prototype experiment allows a similar mechanism (i.e., the development of shear planes leading to soil failure) to be observed under an increasing compression load. The pile load test would allow a comparison of the seismic data generated by the soil movement along shear planes. Moreover, the study included site characterization (landslide geometry and its stratigraphy) as well as directional analysis in order to identify any possible directional effects. This study provides a good basis for future investigation on landslides' early warning system.

Study area
The geology of the Federal District can be summarized in four lithological boundaries, according to regional stratigraphy, delineated as (base to top): (i) Paranoá (metasedimentary rocks -Meso-Neoproterozoic), Canastra (phyllites -Meso-Neoproterozoic); (ii) Araxá (schists -Neoproterozoic); (iii) Bambuí (clayed metasiltites rolled, clay and metasiltites banks -Meso-Neoproterozoic); (iv) Groups of soil or shallow colluvial deposits (pedimentary type -Quaternary) (Freitas- Silva and Campos, 1998). According to the Brazilian soil classification system, the region of the study area, Brasilia (Federal District), is comprised of a mantle of Tertiary-Quaternary age detritus-lateritic soil. The soil is composed mainly of red-yellow latosols and has high porosity and a cementitious bond type. Therefore, it has a highly unstable structure when subjected to an increased moisture or changes in the state of stresses. This soil almost always presents a volumetric change as high as the variation of these factors (referred to as a collapsible structure) (Hussain et al., 2019b).
On the slope of the study area, the channel depositions in the basin are distinguished as (i) landslide debrites due to recent rotational failure mechanism; (ii) fluvial deposits on recent riverbanks as alluvial fan deposits, colluvial, and alluvial. The detailed geology and geomorphology of the area can be found elsewhere (Hussain et al., 2019c and references therein).

Sobradinho landslide
The area of landslide is part of the Maranhão River Basin. The hydrographic network in this region is composed of a high-grade segment of the Contagem stream, which is a direct tributary of the Rio Maranhão. The main drainage system is controlled by the structural lineaments (faults and fractures) and low sinuosity, while the secondary network is typically dendritic with high sinuosity. The main watercourse is the Ribeirão da Contagem, but there are also small intermittent or perennial tributaries with flow rates of less than 0.005 m 3 / second in the dry months of the year. The selected landslide is E-W trending roto-translational slide (Hungr et al. 2014), approximately 150 m long and about 70 m wide, located in a small vicinity called 'Rua do Mato' (Figure 1). At the top of the landslide mass (along the main scarp), there is a small scarp in the middle created by the erosion caused by the release effects of the Contagem River (according to in-situ inspections), which is quite high during the rainy season (Hussain et al., 2019b). The Sobradinho landslide has been extensively studied with geophysical, geotechnical, and geomorphological methods. Figure 1 shows the location of the study area, the distribution of the HVSR stations and some aspects of the landslide.
The mechanical behavior of early-stage earthflows in clayey formation might be considered as Coulomb plastic solids whose mobility is developed by sliding (as described above). In this manner, the material flows by internal deformation (internal shears) (Berti et al., 2019). So, the sliding can be resolved into several occurrences of propagation of sliding plane or soil failure (Vouillamoz et al., 2018;Hussain et al., 2019b). A similar mechanism (i.e., where shear planes are developed leading to soil failure) can be observed when a concrete pile is subjected to an increased compression force leading to soil bearing failure. However, this mechanism differs in the type and level of applied stresses, which are very small in the landslide under the influence of triggering factors (e.g. rainfall).

Material and Methods
Pile load test (prototype experiment) was conducted in this study to artificially generate a failure mechanism with sliding of soil along shear planes at field-scale. Figure 2-A and 2-B show the location of the testing site and its soil profile. By conducting this test, our aim was to aid the interpretation of the results obtained from Sobradinho landslide. This hypothesis was based on the potential similarity of some aspects of the failure mechanisms in both cases (i.e. the landslide and the pile load test) as discussed below.

E F
The fieldwork of the pile load test was conducted in several stages. In the first stage of this prototype experiment, ground characterization of the testing site was carried out using Standard Penetration Test (SPT-N). In the second stage, a concrete pile of 0.8 m diameter and 12 m depth was cast-in-situ and allowed to harden for 4 weeks. The pile was then subjected to progressive axial load while measuring the vertical displacement of the pile head. According to Brazilian standards (Hussain et al., 2019b), soil failure occurs at a displacement greater than 8 mm, where the base resistance is likely to be mobilized i.e. when shaft resistance has been or is close to being fully mobilized.
As the pile load was gradually increased at 30-minute intervals, the HVSR data was obtained from each 30-minute long ambient noise record. The pile loading mechanism was supported by two 8m deep reaction piles that balanced the applied load on the main pile. A load cell was placed below the steel rods, and a hydraulic jack was attached to the load cell. This formed a good base for applying a compression force on the pile by a manually operated oil pump. The hydraulic jack was attached to two smaller metallic rods on which four extensometers were installed for monitoring the pile vertical displacement ( Figure 2C and D).

Ambient noise data acquisition at landslide area
The experimental data was acquired using Sercel L-4A-3D short period seismometers having a natural frequency response of 2 Hz. The threecomponent seismometer records three-dimensional ground motion in the northsouth (NS), east-west (EW), and vertical (Z) directions as required by the single station HVSR analysis. For continuous data acquisition, DAS-130 REFTEK dataloggers were used that can record up to 1000 samples per second (SPS). The time and position of each measurement point were provided by the GPS-130 locks. The seismometer was aligned with the geographic north determined by a magnetic compass. The sensors were also leveled and, in order to avoid environmental interference, buried in the 50cm deep hole. The ambient noise measurements were automatically stored on two memory cards installed in the datalogger.
In the landslide experiment, the data acquisitions were conducted in two stages. In the first stage, single station HVSR measurements were taken both at the landslide mass and in the surrounding areas ( Figure 1C) in order to delineate the landslide's boundary. Ambient noise at each point was recorded for 30-minutes length in accordance with the guidelines provided by SESAME (2004). In the second stage, for continuous ambient noise recording, three sensors were installed on the landslide in a triangular geometry, as adopted by Hussain et al. (2019c). The continuous data used for the seasonal impact evaluation was divided into three acquisition campaigns: (i) dry period, from 02-11-2016 to 08-11-2016; (ii) intermediate, from 12-04-2017 to 15-04-2017; and (iii) saturated period e.g. from 10-12-2017 to 17-12-2017.

HVSR technique
In the single-station HVSR technique, ambient noise recording is conducted at a single three-component station (N-S, E-W and Up-Down). The HVSR curve is interpreted in terms of its peak frequency and amplitude. The HVSR curve peak frequency is traditionally interpreted as a representative of the Rayleigh wave ellipticity or shear wave amplification over horizontally stratified layers. The physical significance of the HVSR peak amplitude is still dubious and open to debate (Fäh et al., 2001). Many authors (e.g., Bonnefoy-Claudet et al., 2008) agree on a relationship between the HVSR amplitude and the impedance contrast. Few of the documented interpretations of the peak amplitude are the reflection of ambient noise wavefield typology (Rayleigh, Love, and body waves). However, the environmental conditions such as sensor ground coupling, water saturation Gosar and Lenart, 2010;Pazzi et al., 2018a;2018b), weather nearby structures, underground structures, and transient noise can possibly affect the results obtained from HVSR (Chatelain et al., 2008). In a porous material, an increase of water content can influence (increase) the P-wave velocity more than Vs, thus causing an increase of Poisson ratio, which, in turn, induces an increase of Rayleigh wave ellipticity (i.e., of the H/V ratio of the elliptical particle motion). This makes the practical application of HVSR for the dynamic (fluidization) monitoring of landslide possible. In this study, the HVSR survey was carried out in accordance with the approach provided by SESAME (2004). In this approach, the duration of ambient noise acquisition is dependent on the ambient noise levels in the area; usually, the record length should be higher than 15 minutes.
The HVSR measurements yield important insights into the resonance behavior of particular lithomorphologies and geomorphological features, while the polarization analysis of the ambient noise records can provide information about the particle motion and directivity properties, which are associated with the mechanical vibrational behavior of structures (Burjánek et al., 2012;Zare et al., 2017). The HVSR has been applied in several application including the determination of the resonance of soil deposit over the basement by Nogoshi and Igarashi (1971), earthquake hazards assessments through microzonation (Di Giulio et al., 2014), bedrock depth estimation (García-Jerez et al., 2006;Gosar and Lenart, 2010), archeology studies (Kalil et al., 2016), estimation of the ice sheet thickness and of the basal properties of a glacier (Picotti et al., 2017) as well as for the dynamic monitoring of glacier (Köhler et al., 2015) and landslide assessment (Imposa et al., 2017;Bertello et al., 2018, Hussain et al., 2019c. Despite these enormous applications in the last decade or so, the theoretical basis of the technique is still unclear and open to debates.

Data processing adopted in this study
The same data processing workflow was applied to the records of all experiments. The data was acquired in REFTEK format, which is first converted to a Seismic Analysis Codes (SAC) format. Then using SAC software, the data was detrended and demeaned. In the next stage, the three-component records of each point were uploaded to a freely available software Geopsy® (www.geopsy.org), where the HVSR analysis was conducted. The recorded transient events, due to any earthquakes and intermittent anthropogenic noises, were removed by choosing automatic window selection in the software. The entire length of the data was divided into short windows while performing the calculation. In the last step, the spectra were smoothened by Konno-Ohmachi function with a constant value of 40 and a cosine taper of 10% width. The entire HVSR calculation processing workflow is presented in Figure 3.

Directional analysis
While it is generally accepted that the peak frequency is reliably diagnostic of site subsoil conditions, several researchers have observed that amplitude can be considerably influenced by factors independent of site conditions e.g. the noise wavefield composition including different portions of the body, Love and Rayleigh waves (Bonnefoy-Claudet et al., 2006). The objective of the directional analysis is to understand how the landslide body receives a contribution of ambient noise energy along the different axis. If the amplification at two horizontal components is unequal with larger motion at a specific azimuthal, this is known as directional amplification (Bonamassa and Vidale, 1991). The presence of fractures and cracks as well as fault zones may show directionality. The directional amplification is suggested to be due to the presence of largescale open cracks (Moore et al., 2011) or microcracks (Pischiutta et al., 2017). In this way, links between the distribution of noise sources and the subsurface conditions are evaluated. The analysis of directional effects in the site response is conducted by rotating the components of the spectral ratios obtained at each measurement site by steps of 10•, starting from 0•N to 180•N, and plotting the contours of spectral amplitude as a function of frequency and direction of motion (Panzera et al., 2012). The 0° stands for the North direction marked on seismometer and 180° stands for South (from 0° to 360° results are symmetric). The directional effects of landslide-affected areas have been shown in previous studies (Burjánek et al., 2010;Imposa et al., 2017).

Pile load test
The results of the pile load test showed three peaks on HVSR curves: first at 5 Hz (deep interface), second at 32 Hz (not considered), and third at 55Hz (shallow interface). Both peaks (shallow and deep) represent the depth of soil horizons, i.e., shallow and deep interfaces. Following an 1D simple modeling, a site resonance frequency fo was related to the thickness h of surface material having velocity Vs through the following equation: The Vs value is calculated from SPT-N using Imai et al. (1977) relation for all soils as 550 m/s, which for fo=5 Hz would imply a thickness of less than 30 m. Assuming a value of 55Hz for the relevant peak, the thickness of the material involved in resonance would be 2.5 m, i.e. a thin surface soil layer (Hussain et al., 2019d). Both peaks show variations in their amplitudes with the increase of an applied load with constant resonance frequency. According to the results shown in Figure 4 (A), the variations of the mean amplitude measured on successive time intervals of 30 minutes were very small and did not follow a trend correlated to the load increase. These variations are most likely irrelevant to the changes in soil mechanical properties; it seems that they depend on the variations in ambient noise wavefield composition recorded over time. There was another HVSR peak at a frequency range of 32Hz, having peak amplitude rather low, not reaching an HVSR factor of 2, which is generally considered a minimum threshold to infer the occurrence of significant resonance phenomena (see SESAME-2004 guidelines); thus, this peak was not considered for the analysis.
The geometrical similarities of the curves revealed that the peak frequency of HVSR curves at all time steps were similar. This might be attributed to the presence of identical subsurface stratigraphy (the same impedance contrast), which did not change with the applied load. The shear sliding (internal sliding) in the soil was too small to mark any difference in the HVSR curve peak, although other studies have also documented such changes in peak ambient under different phenomena (e.g., Köhler et al. 2006;Bonnefoy-Claudet et al., 2008).
The recorded low-frequency resonance peak (5Hz) is almost azimuth independent (as shown in A-H in Figure 5), denoting the absence of 2-D effects. The stratigraphic peak did not show any variation with azimuth, while the higher frequency peak showed azimuthal dependent variations ( Figure 5). The prominent azimuth range was 20-60 and 100-180, implying the nonuniform distribution of ambient noise sources surrounding the sensor.

Ubiquitous frequency
The ubiquitous frequency was the frequency of saprolite (highly weathered rock), the soft layer over the seismic bedrock, and it was generally recorded throughout the studied area having high amplitude relative to the landslide peak ( Figures 5A, 6 and 7). This peak was not affected by the seasonal variations (stiffness) produced by meteorological changes (degree of saturation). However, there were slight variations in peak amplitude recorded at different points in the study area (Figures 6 and 7) that could be linked with the composition of ambient noise wavefield. Based on the active MASW results of a study carried out by Hussain et al. (submitted), the value of shear wave velocity was considered as 160 m/s. This value is low because of the high porosity of material and disturbance brought by the onset of the landslide. Using the equation (1) with Vs=160 m/s and h = the thickness of the sediment package over the seismic bedrock, it was found that the thickness of the saprolite plus clay layer is 20m deep. These results are consistent with the results obtained by joint inversion of HVSR and f-k curves, where the depth of Saprolite was found equal to 23 m (Hussain et al., 2019d).
The ubiquitous frequency did not show the amplitude variations of the HVSR along the azimuths, mainly because it represented the stratigraphy and not the fractured landslide mass that led to trap seismic energy at a particular direction (azimuth) (Zare et al., 2017;Rezaei et al., 2018). We observed that the rotated HVSR obtained in the landslide fractured zones showed clear directional effects, which were quite complex and could not be interpreted easily. There were two groups on the HVSR curves: 1) where the amplitude is similar in almost all directions (2Hz peak), and 2) where the higher amplitudes of the landslide peak (higher peak) are between 40° and 150° ( Figure 6). It is possible that can provide some more conformation about the ambient source distributions or the presence of fractures linked with the landslide boundary.

Landslide frequency
In Figures 6 and 7 along the high amplitude stratigraphic peak (2 Hz), there are peaks at higher frequencies having low amplitudes, which presumably indicate the presence of shallower velocity boundaries within the sediments. The amplitudes of these side peaks were, in most cases, considerably lower than the amplitude of the main peak; only at P5 and P11, the amplitudes were higher. This relatively low amplitude over the HVSR curve secondary peak is mainly related to the presence of thin and soft soil covers (Stanko et al., 2017). The presence of such peaks has also been reported in previous studies (Gosar and Lenort, 2010). The variation in amplitude is possibly related to the ambient noise sources in the region and the contribution of Love waves that affect the amplitude of the H/V peak, as described in section 4.1. These higher frequency peaks 6Hz and 8Hz represent landslide depth (slip surface) at 6.6 m and 5 m, respectively. The depth of the slip surface is relatively consistent with the usual depth surface of the landslides in the region.
The higher resonance frequencies represent the detached soil blocks created by the seasonal erosion of the Contagem River flowing at the bottom. The river crosses the landslide at the bottom, which is a possible source of disturbance (dynamism) in the landslide mass. As a result, fractures are produced; some of these fractures are very deep, supported by visual field inspections -conducted by the authors. Each block has its own degree of freedom and has a different resonance frequency, which may explain the emergence of secondary peaks of different frequencies depending on the location along the landslide (Figure 6, 7). The emergency of stratigraphic and landslide peaks has been reported in a previous study by Rezaei et al. (2018) at a landslide in Tehran, Iran.

Flat HVSR
This class of HVSR curves lacked any frequency peak according to the reliability criteria of clear peak (amplitude greater than 2). There are two possible reasons for the absence of a peak on HVSR curves (i) bedrock are exposed or (ii) low-quality data. At some locations, the secondary peak also appeared with very low amplitude because of the presence of coarser grain material towards the floodplain of the river that lies behind the landslide. These curves were not considered for the analysis. The results are shown in Figure 8. Figure 6. HVSR curves of group-A Peak at 2 Hz and 4-6 Hz. The two dashed lines represent the HVSR standard deviation, while the gray areas represent the peak frequency standard deviation, which quantifies the experimental error associated with the average peak frequency value. Below we show the contour plots of HVSR curves as a function of frequency (x-axis) and rotation angle (y-axis) with respect to the geographic north. The color scale represents the HVSR amplitude. Figure 7. HVSR curves of group-B peak at 2 Hz and 8-10Hz (P refers to HVSR station, Figure 1). Legend is the same as of Figure 5. The color bar presents the HVSR amplitude.

Figure 8.
No peak/flat response of the strata on the HVSR curves. HVSR standard deviation, while the gray areas represent the peak frequency standard deviation, which quantifies the experimental error associated with the average peak frequency value.

Degree of saturation
The effects of rainfall infiltration are possibly related to changes in the Poisson ratio that affects Rayleigh wave ellipticity (HVSR). The changes in Rayleigh wave ellipticity may also be related to the rise in pore pressures in the landslide mass, internal sliding, and changes in the state of landslide mass (as described in section1). In addition to the reduction of the soil stiffness due to the pore water pressure, the filtered water can reduce the shear wave velocity through increasing the degree of saturation. These phenomena have extensively been studied in the landslides triggered by rainfall and earthquakes (Mainsant et al. 2012;Harba and Pilecki, 2017;Hussain et al., 2019b;2019c).
The results of the present study did not show changes in the ubiquitous and higher frequency peaks and their amplitudes in response to rainfall dependent degree of saturation (Figure 9, 10). One possible reason as documented by Hussain et al. (2019c), is the insufficient amounts of rainfall recorded during the shorter data acquisition campaigns. These conditions did not let the highly porous landslide mass to reach high saturation levels.
At station DF01, no secondary peak was observed, and only stratigraphic peak appeared. This station was located in a stable place. The peak showed no variations in response to meteorological changes ( Figure 9). The records of station DF02 showed a well-defined HVSR peak at higher frequency (landslide peak). This higher peak was present throughout the entire length of the experiment irrespective of the day-night ambient noise variabilities ( Figure 10). However, there was no evidence of the changes in the resonance of landslide with rain, which showed the absence of dynamism (e.g., fluidization or internal sliding) in the landslide mass.

Conclusions
The HVSR results obtained from the landslide showed several peaks. In particular, the ubiquitous peak at low frequency (2Hz) was not related to landslide dynamics. Instead, it represented the site response of seismically lowvelocity soil layer overlying higher velocity Saprolite (20m deep). Whereas the landslide peaks at high frequency were related to the detached blocks at the landslide mass and their amplitudes reflected the degree of impedance contrasts created by either the presence of thin, soft sediments or the composition of ambient noise wavefield. Both peaks (ubiquitous and landslide) did not show any response to the effects of rainfall, mainly because of the relatively small degree of soil saturation. In comparison, the pile load test showed two stratigraphic peaks, for every increase in soil stress there was a small variation in the amplitudes of the HVSR peaks, which might be related to the composition of ambient noise wavefield, Overall, these observations suggest that HVSR amplitude variabilities were complex and created by the structural change in the shallow sub-surface and the ambient noise wavefield composition. The variability reflected local noise conditions at each individual station affected by topography, vicinity to streams (Contagem River), exposure to wind, and extend and timing of degrading instrument coupling. Therefore, our data did not allow precise identification of the relationship between rainfall-induced dynamic (fluidization), displacement rate, and Rayleigh ellipticity. However, it was still possible to conduct the landslide site characterization using limited ambient noise data.
The present study provides a sound basis for a more comprehensive application of seismic noise measurements, a potential method in the disaster diagnosis and prediction, with an emphasis in densely populated mountainous areas. Further investigation is recommended where greater deformations and larger scale fractures (centimeters) can be created in the soil in order to identify the possible effects of soil deformations on HVSR. The study can also benefit from the development of an improved physical simulation of the emergence and propagation of fractures as well as the fall of soil blocks in order to determine the possible effects of the landslide -before, during and after the event. Moreover, these results can be validated by conducting remote sensing of the landslide surface to confirm the actual movement of the soil mass.