Horizontal-to-vertical Spectral Ratio Method Research Articles

Identification of Subsurface Water Pipe Based on Horizontal to Vertical Spectral Ratio (HVSR)

Disruption of water distribution through underground pipe is a frequent case. As a handling, a method that can detect disturbances from subsurface pipe is needed. Microtremor is one of the geophysical methods that can provide an overview of the subsurface. Processing is done with the HVSR method with input in the form of seismic waves in the time domain which is then processed using Geopsy software to get output in the form of dominant frequency values, amplification, and H/V spectrum. The data used in this study are primary data obtained from data acquisition with a total of 60 measurement points. The distance between points is given 5 dm with a measurement time of 10 minutes for each point. The study was conducted on pipe with empty, partially filled with water, and fully filled with water with a pipe diameter of 0,9 dm. The pattern in each pipe condition from the modeling results appears at a frequency of 0-10 Hz. The difference in the range of amplification values in pipe anomalies is due to different levels of water content which can be concluded that pipe filled with water have higher amplification than empty pipe and be showed in 3D modelling.

Identification Of Shear Strain On The Surface Ground Of Wangi-Wangi Island, Southeast Sulawesi, Indonesia, Using Nakamura’s Technique and The Possibility Of Its Impacts

This research was conducted to determine the possible impact of an earthquake on the mainland of Wangi-Wangi Island based on the presence of shear strain on the surface ground (γ). The size of γ is obtained by multiplying the ground susceptibility index and the acceleration of basement ground or PGA using Nakamura’s technique. The data used are microtremor data and earthquake data from 1920 to 2020 sourced from the USGS. Microtremor data are obtained from the results of filtering ground vibration signals using a Band Pass Filter in the frequency range between 0.5 to 25 Hz. Ground vibration signals were recorded at 47 measurement points spread over the surface of Wangi-Wangi Island within 29.25 to 48.16 minutes. Furthermore, the microtremor data were processed using the HVSR (Horizontal to Vertical Spectral Ratio) method. The use of earthquake data must meet the requirements for a surface magnitude (MS) ≥ 5.0 SR and an earthquake epicenter depth (h) ≤ 45 km. The results obtained are the γ sizes of Wangi-Wangi Island in the order of 10-06 to 10-03. Based on the size distribution, it is known that the majority of the Wangi-Wangi Island area has the potential to experience cracks and land subsidence due to settlements if an earthquake occurs, and only a portion of the area is vibrating. In addition, it is also known that the mainland of Wangi-Wangi Island is not prone to landslides and liquefaction because γ<10-2.

Application of HVSR Method on Microtremor Data for Analysis of Earthquake Potential in Candipuro District, Lumajang, Indonesia

The Earthquake of 6.1 magnitude occurred on April 10, 2021, in Lumajang district. One affected area suffered building losses, road damage, and casualties in Candipuro District, Lumajang Regency, East Java. Then the research was conducted to determine the strength of the soil against the potential earthquake by micro zoning the area in Candipuro District. This research uses the microtremor method with HVSR (Horizontal to Vertical Spectral Ratio) analysis. Analysis of earthquake potential includes parameters of Ground Amplification (A0), Dominant Frequency (f0), Ground Vulnerability Index (Kg), and Ground Shear Strain (GSS). The acquisition was completed with 16 measurement points in Supiturang Village with 500 meters between points using the Portable Seismograph TDS - 303. The results obtained in Supiturang Village have the results of the amplification value factor which is classified as medium-high (2.7 - 8.7), low dominant frequency value (0.5 - 1.4 Hz), medium-high soil susceptibility index (10 - 90), VS30 value which includes medium - hard soil (187.81 - 548.38) m/s, and GSS value which shows the response of cracks and soil subsidence to shaking (2 x 10-5 - 2 x 10-4).

Earthquake Hazards Vulnerability Analysis of The Nankou-Sunhe Fault in Changping (Beijing) Using HVSR Method

During the Holocene period, three paleoseismic events (Mw=7.7-7.9) occurred on the Nankou-Sunhe Fault in Beijing. For disaster prevention and mitigation, earthquake zoning is needed to be carried out to map the areas that will experience severe damage when an earthquake occurs. This study uses the Horizontal to Vertical Spectral Ratio (HVSR) method to obtain not only the sedimentary layer thickness and fault structure according to empirical formulas, but also the spatial variations in the peak frequency and amplification values that can explain the characteristics of the geological layer beneath the surface. Based on the calculated results, the peak frequency(f) ranges from 0.2 to 4 Hz with the lowest frequency distribution in the southwest parts, which indicates a thicker sediment layer. Amplification(A) ranges from 2 to 14 with the distribution of higher values in the southeastern and northeastern parts. The K-value (A2/f) is in the range of 0-180 and the highest K-value in the southwest of the fracture indicates that this region is at high risk for damage. Our results provide a better understanding of the earthquake hazards vulnerability of the Nankou-Sunhe Fault in Changping, Beijing, and can be used as a basis for urban disaster prevention and reduction.

Mapping of Soil Dynamic Response of Palu City Based on Horizontal to Vertical Spectral Ratio (HVSR) Method

Palu is one of the earthquake prone coastal cities in Indonesia. The 28 September 2018 earthquake (Mw 7.5) destroyed buildings in Palu due to ground shaking amplification, liquefaction, and tsunami. This phenomenon implies the importance of soil dynamic properties knowledge as a basis for proper building design to reduce such an earthquake. This paper presents the results of microtremor measurements conducted at 110 locations in Palu to estimate the dynamic properties of alluvial and marine deposits for earthquake mitigation. The microtremor data were processed using the HVSR method to estimate dominant frequency and amplification factors. Based on these two parameters, the seismic vulnerability index value can be used to estimate the dynamic characteristic variations of the area. The results show that alluvial and marine sediments in Palu have dominant frequency values from 0.54 to 11.72 Hz, amplification from 1.64 to 6.15, and seismic vulnerability index of 1.49 to 56. Sediments with low dominant frequency and high amplification factors are mostly found in the western and eastern Palu and Mantikulore districts, which are associated with thick sediment layers. The areas with the highest amplification factor and seismic vulnerability index are Petobo and Balaroa, which experienced severe damage during the 2018 earthquake.

Preliminary result of dominant frequency and seismic amplification in Penajam Paser Utara and its surrounding regions using the HVSR method

East Kalimantan has been designated as the new capital city of Indonesia, precisely in Penajam Paser Utara and Kutai Kartanegara Regency. According to the historical records of seismic activity, there were several destructive earthquakes in East Kalimantan. Therefore, it is necessary to study the characteristics of the local site to assess the seismic hazard potential in the new capital city area and its surroundings. In this research, the Horizontal to Vertical Spectral Ratio (HVSR) method was applied by utilizing microtremor data from 21 seismometers installed in the research area. The results of the HVSR analysis obtained the dominant frequency (f0) and amplification factor (A0) from the H/V curve. Based on the H/V calculation, the highest f0 is 13.04 Hz and the lowest value is 0.58 Hz. Meanwhile, the highest value for A0 is 6.41, and the lowest value is 1.41. The results of the study show that the research area has a low to medium level of seismic hazard.

Amplification Analysis and Seismic Vulnerability Index of Air Tawar Barat Village of Padang City Using HVSR Method

Padang City is an area that has a very high level of seismic activity in Indonesia. The city of Padang is also an area severely affected by the earthquake in the Subduction Zone. Considerable earthquakes impact building damage and infrastructure, one of the factors that can see the vulnerability of an area with a seismic vulnerability index. Therefore, a study was conducted to determine and analyze the amplification and seismic vulnerability index in the Air Tawar Barat Village of Padang City and create a map of the natural frequency distribution, amplification, and seismic vulnerability index. Microtremor data retrieval was carried out in February-May 2022 using the Seismograph Sysmatrack MAE and S3S Sensors at Air Tawar Barat Village, Padang City. The method used in this study is the HVSR (Horizontal to Vertical Spectral Ratio) method using geopsy software to remove noise in the data and see the natural frequency value and the value of the amplification factor. Furthermore, from the data, contour maps were made using Surfer 13, where the type of research carried out was in the form of descriptive research. Based on the results of the study found that the natural frequency value (f0) in the study area ranged from 0.59141 Hz – 12.9172 Hz, the amplification (A0) ranged from 2.39825 to 13.7759, and the seismic susceptibility index (kg) was 0.4802 – 320.885. Based on the data processing results, it can be concluded that fresh western water is an area with a high seismic vulnerability index, and the location is very vulnerable when shaken by an earthquake.

Seismic microzonation based on HVSR inversion results for shear wave (Vs30) mapping and soil vulnerability in West Sulawesi and South Sulawesi Regions

The Horizontal to Vertical Spectral Ratio (HVSR) is a seismic analysis method used to obtain information about subsurface geological characteristics. This method is based on the spectral ratio analysis between the horizontal and vertical components of earthquakes or other vibration sources. HVSR has proven to be an effective tool in evaluating rock layers, layer thickness, rock hardness, and geological hazard potential. In the HVSR method, seismic data are measured using three-component seismometers, the N-S (North-South), E-W (East-West), and vertical components. In this study, the HVSR method is used for Vs30 inversion to estimate shear wave velocity at a depth of 30m. The southern and western regions of Sulawesi have various geological conditions and high seismicity levels. The results of this inversion can be used for mapping soil vulnerability indices. The obtained soil vulnerability index values indicate high vulnerability near the Palu-Koro fault, with A0 9.27 at F0 3.2 Hz. This method can also estimate shear wave velocity values at a depth of 30 m. Vs30 values range from 200 m/s to 500 m/s at various measurement points, indicating variability. According to the SNI 1726:2012 and NEHRP classifications, the research area has soil conditions ranging from medium to soft.

Identification of soil characteristics in the Central Java Region using the HVSR (Horizontal to Vertical Spectral Ratio) inversion method

This research aims as one of the efforts to mitigate earthquake disasters by using the HVSR method based on microtremor data at 28 BMKG seismic stations in the central Java Island region to estimate the dominant frequency (f0) and amplification (A0) that describe the dynamic characteristics of the soil and analyse the shear wave velocity at a depth of 30 meters (Vs30). The dominant frequency values in the study area range from 0.648 to 18.257 Hz. The amplification values scattered in the study area ranged from 0.788 - 7.496. This research also produces dominant period parameters (T0) and seismic vulnerability (Kg). The highest seismic vulnerability values are at BUJI in Buayan, Kebumen; GBJI in Gringsing, Batang; and SBJM in Sanden, Bantul. In addition, the Vs30 value is helpful to determine how strong the building foundation resists earthquakes. Based on the Vs value in the ground profile with ellipticity curve modelling, the Vs30 value ranges from 204.756 - 903.046 m/s. The distribution of Vs30 values shows that the research points in central Java are spread by a variety of hard and very dense soil types to soft rock (SC).

Peculiarities of the HVSR Method Application to Seismic Records Obtained by Ocean-Bottom Seismographs in the Arctic

The application of the horizontal-to-vertical spectral ratio (HVSR) modeling and inversion techniques is becoming more and more widespread for assessing the seismic response and velocity model of soil deposits due to their effectiveness, environmental friendliness, relative simplicity and low cost. Nevertheless, a number of issues related to the use of these techniques in difficult natural conditions, such as in the shelf areas of the Arctic seas, where the critical structures are also designed, remain poorly understood. In this paper, we describe the features of applying the HVSR modeling and inversion techniques to seismic records obtained by ocean-bottom seismographs (OBS) on the outer shelf of the Laptev Sea. This region is characterized by high seismotectonic activity, as well as sparse submarine permafrost distribution and the massive release of bubble methane from bottom sediments. The seismic stations were installed for one year and their period of operation included periods of time when the sea was covered with ice and when the sea was ice-free. The results of processing of the recorded ambient seismic noise, as well as the wave recorder data and ERA5 and EUMETSAT reanalysis data, showed a strong dependence of seafloor seismic noise on the presence of sea ice cover, as well as weather conditions, wind speed in particular. Wind-generated gravity waves, as well as infragravity waves, are responsible for the increase in the level of ambient seismic noise. The high-frequency range of 5 Hz and above is strongly affected by the coupling effect, which in turn also depends on wind-generated gravity waves and infragravity waves. The described seafloor seismic noise features must be taken into account during HVSR modeling and interpretation. The obtained HVSR curves plotted from the records of one of the OBSs revealed a resonant peak corresponding to 3 Hz, while the curves plotted from the records of another OBS did not show clear resonance peaks in the representative frequency range. Since both OBSs were located in the area of sparse distribution of submarine permafrost, the presence of a resonance peak may be an indicator of the presence of a contrasting boundary of the upper permafrost surface under the location of the OBS. The absence of a clear resonant peak in the HVSR curve may indicate that the permafrost boundary is either absent at this site or its depth is beyond the values corresponding to representative seismic sensor frequency band. Thus, HVSR modeling and inversion techniques can be effective for studying the position of submarine permafrost.

Seismic Hazard Potential in Yogyakarta Based on HVSR Curve Estimation

In 2006, the Special Region of Yogyakarta was shaken by a destructive earthquake. The United States Geological Survey (USGS) recorded that the earthquake had a magnitude of 6.3 Mw at a depth of 12.5 km below the surface that was triggered by strike slip fault activity. The severe damage occurred in Bantul district and Klaten district, where the Klaten area location were far from the earthquake epicenter. This proves that the magnitude and the earthquake source distance are not the only parameters for seismic hazard potential, but rather the presence of local site effects and building conditions. The present study aims to determine the characteristics of the soil and estimate the seismic hazard potential by using HVSR method. Horizontal to Vertical Spectral Ratio (HVSR) is one method that can be used to obtain the subsurface information from single station measurements. Furthermore, the inversion of Rayleigh wave ellipticity curve is used to obtain 1-D of shear wave velocity (Vs). The parameters from HVSR calculation and Rayleigh wave inversion were mapped to understand the subsurface structure beneath the Yogyakarta area. Based on the soil classification to Vs, Yogyakarta area is categorized to SD (Stiff soil/soft soil) and SC (Very dense soil and soft rock). The results of the mapping analysis indicated that Bantul district is an area with the highest potential of seismic hazard.

Novel Evaluation Method for Site Effect on Earthquake Ground Motion Based on Modified Horizontal to Vertical Spectral Ratio

The effect of local site conditions on the earthquake ground motion is a very important factor to be considered in engineering seismic fortification. Many methods, such as numerical simulation methods based on site analysis models and statistical empirical relation methods based on the earthquake ground motion observations and numerical simulation data, have been used to consider the site effects in actual engineering seismic fortification and earthquake disaster assessment. The statistical analysis to obtain characteristic parameters of site condition effect based on strong motion and microtremor records become an economical and practical method of determining the designed ground motion of engineering sites, especially for large survey areas and engineering sites where it is difficult to carry out a site survey. In this paper, a novel evaluation method for site effect on earthquake ground motion is proposed. The new method is based on the horizontal to vertical spectral ratio (HVSR) method, but the original HVSR is replaced by a modified HVSR considering the effect of the soil layer on the vertical ground motion. In order to build the model and determine the corresponding parameters of the modified HVSR, first, the ground motions in the bedrock below the soil layer are calculated using the one-dimensional equivalent linear method. These calculated records are independent of the influence of the downgoing wavefield, and the differences between the ground surface to bedrock spectral ratio (SBSR). The HVSR for the local sites of ground motion observation stations are analyzed using the strong ground motion records from the Kiban-Kyoshin network (KiK-net) in Japan. The statistical characteristics of the relationship between SBSR and HVSR are revealed, and then, a quantitative relationship between SBSR/HVSR and HVSR is established. The proposed evaluation method for the site effect has the advantage that the original HVSR method only requires ground motion records on the ground surface of the site, and it further considers the influence of the vertical seismic effect on the accuracy of the HVSR method. The proposed method can characterize the influence of the site conditions on ground motion more reasonably than the conventional method.

Analysis of Ground Shear Strain (GSS) In Mapping, Liquefaction Vulnerability Potential Using HVSR (Horizontal To Vertical Spectral Ratio) Method. Case Study: Sepaku Subdistrict, North Penajam Paser Regency, East Kalimantan

East Kalimantan has been designated as the new Capital Province of Indonesia. The criteria for determining the relocation of the nation's capital have passed various studies in the placement of a new state capital, one of which is by reviewing disaster aspects such as the risk of earthquakes, volcanoes and tsunamis. Based on historical records of earthquakes that have occurred, East Kalimantan is not a completely safe area from potential earthquakes or tsunamis. By using the HVSR method, the authors tried to determine the value of Ground Shear Strain (GSS) in the research area, namely Sepaku District, Penajam Paser Utara Regency. The result show that there were three zones namely Zona A, Zone B and Zone C had GSS values in Sepaku Subdistrict with landslides, soil compaction and liquefaction potential, namely 0.4 - 0.54. The result of HVSR interpretation curve are based on the phenomenon of the resulting ground shear strain value in accordance with the geological conditions in the research area of the Sepaku Subdistrict, Penajam Paser Utara Regency, East Kalimantan in characterizing the soil profile in the study area.

Constraining the depth extent of low-velocity zone along the Chenghai Fault by dense array ambient noise interferometry and horizontal-to-vertical spectral ratio

The depth extent of low velocity fault zones is still a matter of debate due to the lack of vertical resolution for most seismic methods. Autocorrelations and cross-correlations of seismic ambient noise provide us new constraints on this issue by obtaining the shallow subsurface reflections without the need of active sources or earthquakes. The horizontal-to-vertical spectral ratio (HVSR) method can also give constraints on this interface and does not need any extra cost of data collection. In this study, we used one-month continuous seismic data from a dense linear array deployed crossing the Chenghai Fault (CHF) in Yunnan, southwest China. Single station autocorrelation and adjacent station pair cross-correlation methods were used to retrieve body wave reflections from the interface of the shallow crustal low-velocity zone (LVZ). The HVSR method was also performed to delineate the shape of the LVZ. Results show an inverted trapezoidal LVZ which extends to ~1.0 km depth across the fault zone with a lateral extent of ~3 km. Previous studies based on ambient noise surface-wave tomography and teleseismic travel time analysis in the same area show similar characteristics, consistent with the interpretation of our results.

Site characterizations of BMKG seismic network, Indonesia, based on HVSR analysis

In the last few years BMKG (Agency for Meteorology Climatology and Geophysics) has increased the number of seismic stations significantly. Until mid-2020, when this study was conducted, the number of BMKG stations has reached 339 units. Development of the network is aimed to improve the data quality, speed of earthquake data processing and information dissemination to the public, accuracy of the hypocenter, as well as the magnitude. The objective of this study is to identify the site characteristics of the network. We analysed the noise recorded at BMKG stations throughout Indonesia using the HVSR (Horizontal to Vertical Spectral Ratio) method. The results of HVSR analysis were used to classify the site conditions of each station. We got the number of stations with the classification of hard rock, rock, hard soil, medium soil, and soft soil, 47, 57, 52, 30, and 84, respectively. This site condition represents the stations characteristics and affects the quality of seismic waveform data.

The horizontal-to-vertical spectral ratio and its applications

The horizontal-to-vertical spectral ratio (HVSR) has been extensively used in site characterization utilizing recordings from microtremor and earthquake in recent years. This method is proposed based on ground pulsation, and then it has been applied to both S-wave and ambient noise, accordingly, in practical application also different. The main applications of HVSR are site classification, site effect study, mineral exploration, and acquisition of underground average shear-wave velocity structure. In site response estimates, the use of microtremors has been introduced long ago in Japan, while it has long been very controversial in this research area, as there are several studies reporting difficulties in recognizing the source effects from the pure site effects in noise recordings, as well as discrepancies between noise and earthquake recordings. In practice, the most reliable way is the borehole data, and the theoretical site response results were compared with the HVSR using shear wave to describe site response. This paper summarizes the applications of the HVSR method and draws conclusions that HVSR has been well applied in many fields at present, and it is expected to have a wider application in more fields according to its advantages.

Polarization characteristics of Rayleigh waves to improve seismic site effects analysis by HVSR method

The horizontal-to-vertical spectral ratio (HVSR) method based on ambient noise recordings is a popular and efficient approach to estimate the seismic site effect parameters. This method is sufficiently robust in determining the predominant frequency (f0) of the site, whereas its estimate of the amplification factor lacks consistency because generation mechanisms, source locations and proportions of the wave types making up ambient noise vary locally and are not well understood, in addition to lack of clarity in the theory of the HVSR method. As the generation mechanisms, source locations and wave compositions of the single- and double-frequency (SF and DF respectively) microseisms (0.05–0.5 Hz) are better studied compared to the noise in other frequency bands, this study focuses on investigating the site effects on the fundamental Rayleigh (Rg) waves in the SF and DF microseisms recorded at 30 Transportable Array stations in the east coast of the United States. HVSRs are calculated first to obtain f0s as well as the unconsolidated sediment thicknesses (USTs). Then the polarization analyses are carried out to select the ellipsoids approximating the particle motions of the Rg waves and to estimate the HVSR of these Rg ellipsoids (HVSRR). The polarization parameters and HVSRR indicate that the Rg waves' particle motions form oblate ellipsoids in the UST range whose f0 values fall within the DF band. Comparison of the observed HVSR and HVSRR and the computed HVSR values based on the body wave and the Rayleigh wave theories shows that the observed HVSR in the DF band on bedrock and/or thin unconsolidated sediments are mainly due to the ellipticity of Rg waves, while those on thick unconsolidated sediments whose f0 values fall within the DF band result from the combined effects of the ellipticity of Rg waves and multiple reflection of S-waves. This outcome will help develop procedures to improve the consistency of HVSR's amplification estimates and hence its reliability in seismic site effects analysis.

Characterizing Fundamental Resonance Peaks on Flat-Lying Sediments Using Multiple Spectral Ratio Methods: An Example from the Atlantic Coastal Plain, Eastern United States

ABSTRACT Damaging ground motions from the 2011 Mw 5.8 Virginia earthquake were likely increased due to site amplification from the unconsolidated sediments of the Atlantic Coastal Plain (ACP), highlighting the need to understand site response on these widespread strata along the coastal regions of the eastern United States. The horizontal-to-vertical spectral ratio (HVSR) method, using either earthquake signals or ambient noise as input, offers an appealing method for measuring site response on laterally extensive sediments, because it requires a single seismometer rather than requiring a nearby bedrock site to compute a horizontal sediment-to-bedrock spectral ratio (SBSR). Although previous studies show mixed results when comparing the two methods, the majority of these studies investigated site responses in confined sedimentary basins that can generate substantial 3D effects or have relatively small reflection coefficients at their base. In contrast, the flat-lying ACP strata and the underlying bedrock reflector should cause 1D resonance effects to dominate site response, with amplification of the fundamental resonance peaks controlled by the strong impedance contrast between the base of the sediments and the underlying bedrock. We compare site-response estimates on the ACP strata derived using the HVSR and SBSR methods from teleseismic signals recorded by regional arrays and observe a close match in the frequencies of the fundamental resonance peak (f0) determined by both methods. We find that correcting the HVSR amplitude using source term information from a bedrock site and multiplying the peak by a factor of 1.2 results in amplitude peaks that, on average, match SBSR results within a factor of 2. We therefore conclude that the HVSR method may successfully estimate regional linear weak-motion site-response amplifications from the ACP, or similar geologic environments, when appropriate region-specific corrections to the amplitude ratios are used.

Mapping of weathered layer thickness and Seismic Vulnerability in Tegal using HVSR method

The Baribis-Kendeng-Tegal fault is part of the Baribis-Kendeng fault segment that extends from West to East in the north of Java Island. The existence of this fault segment affects seismic activity in Tegal. The research was conducted by microtremor measurement to determine the earthquake susceptibility index and the thickness of the weathered layer in Tegal City. The data collection method in the study used the HVSR (Horizontal-to-Vertical Spectral Ratio) method to determine the subsurface conditions of the study area. The results of the horizontal component analysis of the vertical component showed a spectrum peak at the dominant frequency and amplification factor. The thickness of the weathered layer was obtained from the comparison of the value of Vs 30 of the study area to its natural frequency value. The seismic vulnerability index was obtained by squaring the amplification and dividing it into its natural frequency. Measurements carried out at 37 points in Tegal City obtained varying thickness of the weathered layer, ranging from 4.29 meters to 88.14 meters and seismic vulnerability values ranging from 0.08 to 16.65.

Analysis of Horizontal to Vertical Spectral Ratio (HVSR) Curve Results from Processing using Interpolated Finite Impulse Response (IFIR) Filter

Microtremor is recorded with a device that is very sensitive to subsurface soil vibrations, that unwanted noise frequency waves are often also included. For this reason, processing needs to be done using filters that can attenuate the noise. In the previous research, signal decomposition was performed using the Interpolated Finite Impulse Response (IFIR) filter and resulted attenuation at the desired frequency value. Using the MATLAB program, microtremor data are processed with IFIR filtering, smoothing, calculation of: Horizontal to Vertical Spectral Ratio (H / V) curve, natural frequency (f0) and standard deviation (σ). The same process was carried out using Butterworth, Gaussian, and Cosine filters as a comparison with the Geopsy program. As a result, the amplitude characteristics obtained can provide a value of f 0 < 2.5 Hz. The value of f0 obtained using IFIR filter is also stable, proven by the standard deviation ranging from σ = 0.009 - 0.83 Hz, with an average of σ = 0.132 Hz. In conclusion, the IFIR filter can provide a well-attenuated H/V curve according to the specified cutoff frequency, which was 2 Hz. The IFIR filter is also considered capable of being applied to the HVSR method Microtremor data. Lowpass IFIR filter can attenuate high-frequency noise in microtremor data.