Empirical Correlation between Standard Penetration Resistance (SPT-N) and Shear Wave Velocity (V s ) for Soils in Metro Manila, Philippines

: The measurement of the shear wave velocities (V s ) of soils is an important aspect of geotechnical and earthquake engineering, due to its direct relation to the shear modulus (G), which in turn inﬂuences the stress–strain behavior of geomaterials. V s can be directly measured or estimated using a variety of onsite tests or in a laboratory. Methods such as downhole PS logging require boreholes and may not be logistically and economically feasible in all situations. Many researchers have estimated V s from other geotechnical parameters, such as standard penetration test resistance (SPT-N), by means of empirical correlations. This paper aimed to contribute to this subject by developing an empirical relationship between V s and SPT-N. Data from twenty sites in Metro Manila were obtained from geotechnical investigation reports. V s proﬁles of the same sites were also acquired using the refraction microtremor method. New empirical relationships were developed for all, sandy, and clayey soil types, using a non-linear regression method that is applicable for Metro Manila soils. Statistical evaluation and comparison of the proposed correlations with other previous works suggested the viability of the empirical model.


Introduction
The measurement of the shear wave velocities (V s ) of soils is an important aspect in the application of geotechnical and earthquake engineering.This geophysical parameter is related to the shear modulus (G) through the elastic theory, in the equation G = ρV s 2 , where ρ is the density of material.The stress-strain behavior of materials is then dictated by the shear modulus, as the strain is directly affected by this parameter [1].This relationship is also important in ground response analysis, as it is used to evaluate the dynamic response of soils during earthquakes.The seismic shear wave velocity of the upper 30-m soil layer (V s30 ) largely influences the ground motion amplification and is therefore considered an important parameter in earthquake engineering [2].Additionally, the shear wave velocity of geomaterials is also used in studies on soil stratigraphy, liquefaction analysis, and site-specific subsurface modeling [3][4][5].
Shear wave velocities can be measured or estimated using a variety of laboratory or field tests.Bender element and resonant columns are some of the laboratory techniques that measure V s .The accuracy of these techniques is greatly dependent on the degree of disturbance of the soil samples [6].In situ field measurements include invasive methods such as down-hole, cross-hole, and up-hole tests.These techniques require the drilling of boreholes to a particular depth and, thus, may not be economically feasible in all situations.On the other hand, non-invasive geophysical methods, such as spectral analysis of surface waves (SASW), multichannel analysis of surface waves (MASW), refraction microtremor Appl.Sci.2022, 12, 8067 3 of 24 of the development of the Central Valley, rapid sedimentation from highlands persisted until the Late Cretaceous, intercalated with tuff and other volcanic products [37].During the Late Tertiary and Quaternary periods, complex tectonic events led to the basic geologic structure of Metro Manila at present.This geologic evolution has produced three distinct terranes in the metropolis: the Coastal Lowland, Central Plateau, and the Marikina Valley.East of the metropolis, the uplands of Rizal province with basement-ophiolite complexes can be found.To the west, this region is bounded by Manila Bay.

The Study Area
Metropolitan Manila (Metro Manila) is the capital region of the Philippines (Figure 1), which is also the center of the economy and the most populated areas.The region is composed of several cities and municipalities, sprawled over an area of around 620 square kilometers.It lies in the southern extension of the Luzon Central Valley.In the early stages of the development of the Central Valley, rapid sedimentation from highlands persisted until the Late Cretaceous, intercalated with tuff and other volcanic products [37].During the Late Tertiary and Quaternary periods, complex tectonic events led to the basic geologic structure of Metro Manila at present.This geologic evolution has produced three distinct terranes in the metropolis: the Coastal Lowland, Central Plateau, and the Marikina Valley.East of the metropolis, the uplands of Rizal province with basementophiolite complexes can be found.To the west, this region is bounded by Manila Bay.BMG (1983) [38] and geomorphologic map obtained from MMIERS (2004) [39].Locations of the testing sites are also depicted in black squares.
The Coastal Lowland is the flat and low plain section of Metro Manila, facing the Manila Bay.This section can be further subdivided into deposits from sand bars, backmarshes, tidal flats, the Pasig River delta, and some reclaimed land.This Quaternary Alluvium deposit has a thickness that can reach up to 40 m [40].Sixty percent of the Metro Manila area is underlain by the Central Plateau, which is the geomorphic manifestation of the bedrock of the region, namely the Diliman Tuff member of the Pleistocene Guadalupe Formation.This lithologic member is made of flat-lying sequences of medium to thinbedded, fine-grained vitric tuff and welded pyroclastic breccias.Minor fine to mediumgrained tuffaceous sandstone was also reported to be part of this unit [41].The Marikina Valley is the flat section on the eastern bounds of the metropolis, with lower elevations than the Central Plateau.This is produced by a pull-apart basin associated with the West Valley Fault (WVF) and East Valley Fault (EVF), collectively called the Marikina Valley Fault System [42].Quaternary alluvium deposits in this section are associated with the floodplain of the Marikina River and the delta along the Laguna de Bay.Surface deposits on this sector are reported to reach up to 50 m in thickness [40].
Geologic field mapping and different subsurface investigations have revealed that most of the Metro Manila soils and sediments are composed of thick sequences of pyroclastic and epiclastic deposits.Pyroclastic deposits, primarily tuff materials, are derived from eruptions of nearby volcanoes in the region.The epiclastic deposits come from pyroclastic materials that were reworked by fluvial and lacustrine processes.Tan (1983) [43] A B  [38] and geomorphologic map obtained from MMIERS (2004) [39].Locations of the testing sites are also depicted in black squares.
The Coastal Lowland is the flat and low plain section of Metro Manila, facing the Manila Bay.This section can be further subdivided into deposits from sand bars, backmarshes, tidal flats, the Pasig River delta, and some reclaimed land.This Quaternary Alluvium deposit has a thickness that can reach up to 40 m [40].Sixty percent of the Metro Manila area is underlain by the Central Plateau, which is the geomorphic manifestation of the bedrock of the region, namely the Diliman Tuff member of the Pleistocene Guadalupe Formation.This lithologic member is made of flat-lying sequences of medium to thin-bedded, fine-grained vitric tuff and welded pyroclastic breccias.Minor fine to medium-grained tuffaceous sandstone was also reported to be part of this unit [41].The Marikina Valley is the flat section on the eastern bounds of the metropolis, with lower elevations than the Central Plateau.This is produced by a pull-apart basin associated with the West Valley Fault (WVF) and East Valley Fault (EVF), collectively called the Marikina Valley Fault System [42].Quaternary alluvium deposits in this section are associated with the floodplain of the Marikina River and the delta along the Laguna de Bay.Surface deposits on this sector are reported to reach up to 50 m in thickness [40].
Geologic field mapping and different subsurface investigations have revealed that most of the Metro Manila soils and sediments are composed of thick sequences of pyroclastic and epiclastic deposits.Pyroclastic deposits, primarily tuff materials, are derived from eruptions of nearby volcanoes in the region.The epiclastic deposits come from pyroclastic materials that were reworked by fluvial and lacustrine processes.Tan (1983) [43] reported that these epiclastic fluvial deposits are generally made of a sequence of loose to firm silty fine sand at the surface, underlain by very soft clayey silt or silty clay, further underlain by stiff to very stiff clay, and finally hard clay and silt, sometimes intercalated with very dense sand and/or gravel that grade into bedrock.The thickness and lateral continuity of these layers is highly variable, such that during borehole investigations on a particular site, it is normal to encounter different thicknesses and sequences of soil and sediment deposits.

Geotechnical Investigations
The most common type of geotechnical site investigation in the country is through the standard penetration test (SPT).The National Structural Code of the Philippines [44] requires at least one borehole for two-story structures and greater, with the number of boreholes determined depending on the footprint area of the building.This in situ geotechnical test involves the driving of a hollow tube into the soil, while noting the number of blows needed to push the tube sampler down a specific vertical distance.A 63.5 kg hammer is repeatedly dropped from a height of 76 cm to advance the tube, in three successive increments of 150 mm.The first 150 mm is considered as a seating penetration and is disregarded, while the total number of blows needed to penetrate the second and third 150-mm iterations is the determined N-value of that layer.This procedure is described in the ASTM D1586 [45] guideline.The number of blow counts or N-value can be affected by a variety of factors.The most important factors include the amount of energy delivered to the tube in each impact of the hammer, the type of hammer, diameter of the borehole, and the length of the rod.
Standard penetration test reports of twenty sites in Metro Manila (Figure 1) were obtained from the Department of Education, and Department of Public Works and Highways of the Philippines.These SPT investigations were typically carried out as a prerequisite before the construction of school buildings and other structures.The sites selected to be included in this study were identified as having sufficient space, where geophysical surveys can be conducted and geotechnical data are available.Moreover, the spatial distribution was intended so that there are representatives of the various major geomorphic units identified in the coastal portion of the area.These sites are located along the western coast of Metro Manila, atop the Coastal Lowlands, and are expected to have relatively thick soil deposits.Specifically, the sites were built on top of sandbars, backmarshes, and lowland deposits.
Most of the boreholes have depths that reach up to 15 m, the typical depth requirement for low-story buildings in the country.In some cases, the depth can reach up to 30 m, particularly for sites with thicker soil profiles.On the other hand, there are also SPT investigations where the penetration only reaches up to five to ten meters, for sites with a shallower bedrock depth.In such cases, a coring procedure is employed to penetrate the target depth.The obtained soil samples from the SPT investigation were subjected to typical laboratory tests, such as grain size analysis, moisture content determination, Atterberg Limit test, soil classification, and unconfined compression test.A summary of the measured N-values and soil types identified for each depth in the twenty sites is presented in Table 1.
Table 1.Summary of the SPT from the twenty sites, indicating the location, water table depth, and Unified Soil Classification System (USCS) soil types at varying depths.The USCS scheme is adapted as the standard practice for soil classification for engineering purposes [46].SM-silty sand, GM-silty gravel, ML-silt, CL-low plasticity clay, MH-high plasticity silt, SC-clayey sand, SP-poorly-graded sand, GP-poorly-graded gravel, CH-high plasticity clay, SW-well-graded sand, and GC-clayey gravel.

Geophysical Investigations
In the country, the most popular method among geotechnical engineers for assessing the strength of soil at a site is through standard penetration testing.Cross-or downhole profiling is also common to obtain in situ measurements of shear wave velocity with depth.These techniques may be precise and can provide measurement at, typically, one-meter resolutions.However, these methods can also rather be difficult and expensive to carry out in urban areas, as they require borehole drilling.Hence, non-invasive geophysical techniques are also utilized to overcome the problems posed by the practicality of borehole drilling.Seismic exploration techniques, e.g., multichannel analysis of surface waves (MASW), spectral analysis of surface waves (SASW), and refraction microtremors (hereafter "ReMi") are relatively more economical, faster, and can be performed in a variety of sites, including urban areas.This study utilized ReMi surveys to derive one-dimensional shear wave velocity profiles of the surveyed sites in Metro Manila.
The ReMi technique was first developed by Louie (2001) [47], to demonstrate the application of ambient noise in determining the shallow shear wave velocity structure of a site.The idea behind this technique is that the vertical component of ambient noise, primarily dominated by Rayleigh waves that are caused by human activities, can be picked up by a geophone array.
In this study, ReMi surveys were performed at the twenty sites in Metro Manila where the SPT investigations were performed.The surveys are done by coupling a straight linear array of twelve vertical geophones into the ground.These geophones have a frequency of 4.5 Hz and were arranged to have a spacing that ranged from four to eight meters.The array was then connected to a DAQLink II multichannel recorder to gather thirty iterations of thirty seconds unfiltered noise.The sampling rate of the recording was set at 2 ms (500 Hz).In the surveyed sites, typical sources of noise include pedestrians, vehicle traffic, construction, and wind.To provide high frequency noise for the recordings, a 10 lbs.sledgehammer was introduced, to hit a 6 inches × 6 inches × 1 inch steel plate placed at each end of the array, and which was utilized as an active noise source.
Data processing was done using ReMi Vspect and ReMi Disper modules.The Rayleigh waves were identified from other wave arrivals using slowness-frequency (p-f) transforms of the noise recording, using the ReMi Vspect module.The fundamental mode phase velocity was chosen along the minimum velocity envelope of energy within the p-f domain.Louie (2001) [47] demonstrated that selecting along this trend is the best procedure for choosing a dispersion curve to obtain the best estimate of the true phase velocities of shear waves (Figure 2A).The ReMi Disper module allows the modelling of a dispersion curve with multiple layers and shear wave velocities, to match the dispersion curve of the field data.The modeler varies the layer thickness and velocities until the resulting dispersion curve matches the previously selected dispersion points (Figure 2B,C).Vertical variations of shear wave velocities on the upper 20 m of the sites are illustrated in Figure 3.For most of the sites, the shear wave velocity of the upper 20 m generally ranged from ~100 m/s to ~220 m/s, corresponding to SPT-N values typically ranging from 2-10.For sites with much stiffer strata, the V s at some depths can reach up to ~250 m/s to ~500 m/s.This was in areas with higher recorded N-values, that ranged from 15 to 50.In terms of the average V s of the upper 30 m of the sites (V s30 ), most of the western portion of Metro Manila was classified into Site Class D, based on the scheme of the National Structural Code of the Philippines [48].This indicates that the area is generally underlain by stiff soil material, with a V s30 ranging from 180-360 m/s.This was also observed in the

Methodology
A regression analysis is a statistical technique employed for analyzing empirical relationships between a dependent variable and one or more independent variables.This technique can be primarily applied to infer a relationship between the dependent and independent variable in a fixed dataset.It can also help in the visualization of variation of the values of the dependent variable, when one of the independent variables is adjusted while keeping the other independent variables constant.A particular type of regression analysis, the non-linear regression technique involves the modelling of the observed dependent variables using a mathematical function that combines one or more independent variables in a non-linear function.Power functions, exponential functions, logarithmic functions, and trigonometric functions are several typical examples of non-linear regression functions.
In this research, a power function was employed to correlate the estimated shear wave velocities, to measure the SPT-N values.This non-linear function is widely utilized in rapid evaluation of V s and associated geotechnical parameters.The power function has a basic form of where x and y are two variables, with x as the independent variable and where y is the response variable, while a and b are parameters of the power function [49].a is a constant controlling the amplitude of the function, and b controls the relationship curve.This research examined the correlation based upon corrected and uncorrected SPT-N values and the V s , irrespective of the depth, overburden pressure, geologic age of the deposit, and fines content.
Applying the natural logarithm function to Equation ( 1) transforms the non-linear function into a linear regression equation with the form where Y = ln y, b 0 = ln a, b 1 = b, and X = ln x.The parameters of regression b 0 and b 1 are determined using the least squares method (3) with x i substituted by ln x, y i by ln y, and n as the number of data points.When the values of b 0 and b 1 are determined, then the original parameters of the power function at Equation ( 1), a and b, can be also be calculated.
To measure the quality of the non-linear regression analysis, various statistical parameters need to be determined as well.The residual (ε i ) of the regression analysis is the difference from the observed value (y i ) of the dependent variable in comparison to its predicted value (y i * ), using the regression function (Equation ( 5)).
The total variability (S YY ) of the dependent variable (y i ) is determined from two components: the (1) variability in the observations of the y i , which is accounted for by the regression function (SS R ); and the (2) residual variation, which is left unexplained by the regression function (SS E ) [49].The method of determining S YY is shown in Equations ( 6) and (7).In these equations, y is the average of the observed values of the dependent variable.
The coefficient of regression (R 2 ) is the fraction of the total variability in the dependent variable that is accounted for by the regression function.This has a value that ranges from 0 to 1 and is expressed as To define the strength and nature of relationship between the dependent and independent variables, the coefficient of correlation (r) is also determined [49,50].This value is expressed as Another method to examine the quality of the obtained correlation is to determine Spearman's rank correlation coefficient (r s ) [50,51].This non-parametric measure of correlation utilizes ranks to measure the correlation between an independent and response variable.The ordered datasets are replaced by rankings and the correlation coefficient (r s ) is calculated on the ranks, measuring the strength of association between the two ranked variables.This coefficient indicates how closely two sets of rankings agree with each other.r s is obtained using the equation where d i is the difference between the ranks and n is the number of members for each parameter.For the purposes of this work, the smallest value of each parameter (i.e., SPT N-values, shear wave velocities) was assigned with rank 1, and the rest of the dataset were ranked subsequently, in increasing order.Identical values in a dataset were assigned a rank equal to the average of their positions in the ascending order of values.Aside from calculating different statistical parameters, graphical analysis of residuals can also be used to evaluate the proposed correlations.Histograms of standardized residuals can identify potential outliers and can determine whether the variance is normally distributed.A probability-probability (p-p) plot compares the empirical cumulative distribution function with a theoretical cumulative distribution function.p-p plots can determine if a model follows the assumed normality, as the points must follow the 1:1 diagonal line to be normally distributed data.

Results and Discussion
Several empirical relations of shear wave velocity (V s ) and penetration resistance (N) have been developed by a variety of researchers.In some references, corrected N-values were used in the regression analysis.These empirical correlations are summarized in Table 2.These correlations have a basic power law form V s = a × N b , where a and b are coefficients that govern the shape of the empirical function.From Table 2, it is observable that for all previously proposed empirical correlations, a typically ranges from 19 to 121 and b ranges from 0.24 to 0.85.It can be observed that the correlation by Kanai (1966) [52] can be used as the lower bound, while the work of Athanasopoulos (1995) [23] and Jafari et al. (1997) [64] provide a suitable upper bound for the proposed empirical relation of this study.The other twentyone correlations were constrained by these three functions.Furthermore, the proposed curve of the present study lies closest to that proposed by Iyisan (1996) [26].
This study also explored correlating the shear wave velocity to corrected SPT-N values by applying several correction factors (Figure 4B).These correction factors are for the rod length (C R ), sampling method (C S ), borehole diameter (C B ), energy ratio (C E ), and overburden pressure (C N ).There are different equations for the correction of SPT-N values (N 1 ) 60 , but the equation applied in this work is as follows: For overburden pressure correction, we adapted the equation proposed by Kayen et al. (1992) [69], which is as follows: σ vo is the effective overburden pressure at the depth in kPa, and P a is the reference pressure equal to 100 kPa.Rod length correction (C R ) was applied using the factors provided by Youd et al. (2001) [70].A sampling method (C S ) correction factor of 1 was used, as most SPT drillings in the Philippines utilize standard samplers.Borehole diameter (C B ) and energy ratio (C E ) factors were applied following the scheme of Skempton (1986) [71].After the application of these correction factors, non-linear regression analysis was employed again, to determine the empirical function to correlate V s and corrected SPT-N values.
Appl.Sci.2022, 12, x FOR PEER REVIEW 15 of 26 twenty-one correlations were constrained by these three functions.Furthermore, the proposed curve of the present study lies closest to that proposed by Iyisan (1996) [26].This study also explored correlating the shear wave velocity to corrected SPT-N values by applying several correction factors (Figure 4B).These correction factors are for the rod length (CR), sampling method (CS), borehole diameter (CB), energy ratio (CE), and overburden pressure (CN).There are different equations for the correction of SPT-N values (N1)60, but the equation applied in this work is as follows: ( 1 ) 60 =  60 ×   (13) For overburden pressure correction, we adapted the equation proposed by Kayen et al. (1992) [69], which is as follows: σ'vo is the effective overburden pressure at the depth in kPa, and Pa is the reference pressure equal to 100 kPa.Rod length correction (CR) was applied using the factors provided by Youd et al. (2001) [70].A sampling method (CS) correction factor of 1 was used, as most SPT drillings in the Philippines utilize standard samplers.Borehole diameter (CB) and energy ratio (CE) factors were applied following the scheme of Skempton (1986) [71].After the application of these correction factors, non-linear regression analysis was employed again, to determine the empirical function to correlate Vs and corrected SPT-N values.
(A)  This relation has a regression coefficient R 2 = 0.5685, indicating a moderate tion, which is slightly lower than the regression coefficient when uncorrected SP ues are used in the analysis.This suggests that uncorrected SPT-N values correl wave velocity better than corrected values.

Sandy Soils
A similar non-linear regression analysis was applied for sandy soil materia the 265 data points, 139 were identified to be of sandy soil.The proposed correl sandy soil types is The determined correlation from the data is V s = 61.77(N 1 ) 60 0.4528 (15) This relation has a regression coefficient R 2 = 0.5685, indicating a moderate correlation, which is slightly lower than the regression coefficient when uncorrected SPT-N values are used in the analysis.This suggests that uncorrected SPT-N values correlate shear wave velocity better than corrected values.

Sandy Soils
A similar non-linear regression analysis was applied for sandy soil materials.Out of the 265 data points, 139 were identified to be of sandy soil.The proposed correlation for sandy soil types is V s = 45.07N0.5534 (16) This proposed correlation has a coefficient of regression value of R 2 = 0.7922.This value is relatively higher than the coefficient for all soil types and can be regarded as a strong correlation [68].Figure 5 shows the plotted V s and N-value for sandy soil and the corresponding proposed empirical relation.The empirical correlations for sandy soils proposed by other researchers were also compiled.Thirteen correlations follow the power law form V s = a × N b while utilizing the uncorrected SPT-N values in the analysis, while Pitilakis et al. (1999) [22] and Hasancebi and Ulusay (2007) [28] utilized corrected N-values for their statistical analysis.From the compiled curves, the a coefficient ranges from 31.7 to 162, while the b coefficient ranges from 0.17 to 0.54.On the other hand, Fumal and Tinsley (1985) [56] and Kayabali (1996) [27] proposed an empirical relation with the form V s = a(N + b) c .Pitilakis et al. (1999) [22] and Hasancebi and Ulusay (2007) [28] utilized corrected N-values for their statistical analysis.From the compiled curves, the a coefficient ranges from 31.7 to 162, while the b coefficient ranges from 0.17 to 0.54.On the other hand, Fumal and Tinsley (1985) [56] and Kayabali (1996) [27] proposed an empirical relation with the form Vs = a(N + b) c .
When the proposed empirical relation is compared to the works of other researchers, it can be observed that the relations proposed by Okamoto et al. (1989) [59] and Shibata (1970) [10] can act as the upper and lower limits for empirical relationships for sandy soil types, respectively.The other eleven published empirical relations lie between these constraints, including the one proposed by this study.As the curve of this present study lies in the middle of all the other curves, this indicates a very good agreement with respect to the previous works.

Clayey Soils
A non-linear regression analysis was also applied to clayey soil types.A total 120 data pairs of Vs-SPT-N values were identified and utilized for the statistical analysis.The proposed correlation for clayey soil in this study is This proposed correlation has a coefficient of regression value of R 2 = 0.73.This value is relatively lower than the coefficient for all soil types but can still be regarded as having a strong correlation [68].Figure 6  It can be observed from the compiled empirical curves for clayey soils that the functions by Pitilakis et al. (1999) [22] and Athanasopoulos (1995) [23] can serve as the upper bound of the proposed relations, with the curve of Pitilakis et al. (1999) [22] for N values less than 22 and N values greater than 22 for Athanasopoulos (1995) [23].For the lower bound, the curve by Dikmen (2009) [30] can be applied for N-values less than 24.Further constraints by Imai (1977) [14] can be applied as a lower bound for N-values greater than 24.The other empirical relations lie within these boundaries, including the one proposed  [10,11,14,18,22,24,28,30,35,55,59,60,62].
When the proposed empirical relation is compared to the works of other researchers, it can be observed that the relations proposed by Okamoto et al. (1989) [59] and Shibata (1970) [10] can act as the upper and lower limits for empirical relationships for sandy soil types, respectively.The other eleven published empirical relations lie between these constraints, including the one proposed by this study.As the curve of this present study lies in the middle of all the other curves, this indicates a very good agreement with respect to the previous works.

Clayey Soils
A non-linear regression analysis was also applied to clayey soil types.A total 120 data pairs of V s -SPT-N values were identified and utilized for the statistical analysis.The proposed correlation for clayey soil in this study is V s = 70.26N 0.4220 (17) This proposed correlation has a coefficient of regression value of R 2 = 0.73.This value is relatively lower than the coefficient for all soil types but can still be regarded as having a strong correlation [68].Figure 6 illustrates the plotted data pairs of V s and N-values for clayey soil and the corresponding empirical relations.Empirical relations proposed by other authors were also compiled and compared.Eight correlation functions used the basic power law function while incorporating uncorrected SPT-N values in the analysis.In this compilation, coefficient a ranges from 27 to 114.43, and coefficient b ranges from 0.217 to 0.73.
It can be observed from the compiled empirical curves for clayey soils that the functions by Pitilakis et al. (1999) [22] and Athanasopoulos (1995) [23] can serve as the upper bound of the proposed relations, with the curve of Pitilakis et al. (1999) [22] for N values less than 22 and N values greater than 22 for Athanasopoulos (1995) [23].For the lower bound, the curve by Dikmen (2009) [30] can be applied for N-values less than 24.Further constraints by Imai (1977) [14]  It can be observed from the compiled empirical curves for clayey soils that the functions by Pitilakis et al. (1999) [22] and Athanasopoulos (1995) [23] can serve as the upper bound of the proposed relations, with the curve of Pitilakis et al. (1999) [22] for N values less than 22 and N values greater than 22 for Athanasopoulos (1995) [23].For the lower bound, the curve by Dikmen (2009) [30] can be applied for N-values less than 24.Further constraints by Imai (1977) [14]

Statistical Evaluation
Several methods were used to evaluate the proposed Vs-SPT-N values empirical lationships for all, sandy, and clayey soil types (Table 3).From the regression analy evaluated regression lines have an R 2 that ranges from 0.73 to 0.80, suggesting a moder to strong fit of the proposed empirical function to the obtained dataset.Evaluation us the coefficient of correlation (r) and the Spearman rank coefficient (rs) suggested that th is moderate to strong correlation between the Vs and SPT-N values from the dataset.An alternative method applied to evaluate the validity of the proposed empirical lations was through an analysis of the residuals.Graphical methods for analyzing the siduals applied in this research included the development of residual histograms a probability-probability (p-p) plots (Figure 7).Histogram plots of the residuals determin whether the variance of the datasets was normally distributed.An ideal histogram of residuals, which follows a bell-shaped distribution and which is distributed around ze suggests that the normality assumption is likely true.In this analysis, histograms w produced from standardized residual values, which is solely the raw residual divided its standard error.The method of standardizing residuals is applied to transform da such that its mean centers at zero and the standard deviation equals one.This techniq also identified potential outliers, such that residuals with standard deviations less than or greater than 3, could be identified.In this analysis, histograms from the datasets su gested that the residuals generally followed a normal distribution, and outliers w standard deviations greater than |3| were not detected.[14,22,23,28,30,35,60,61,66].

Statistical Evaluation
Several methods were used to evaluate the proposed V s -SPT-N values empirical relationships for all, sandy, and clayey soil types (Table 3).From the regression analysis, evaluated regression lines have an R 2 that ranges from 0.73 to 0.80, suggesting a moderate to strong fit of the proposed empirical function to the obtained dataset.Evaluation using the coefficient of correlation (r) and the Spearman rank coefficient (r s ) suggested that there is moderate to strong correlation between the V s and SPT-N values from the dataset.An alternative method applied to evaluate the validity of the proposed empirical relations was through an analysis of the residuals.Graphical methods for analyzing the residuals applied in this research included the development of residual histograms and probability-probability (p-p) plots (Figure 7).Histogram plots of the residuals determined whether the variance of the datasets was normally distributed.An ideal histogram of the residuals, which follows a bell-shaped distribution and which is distributed around zero, suggests that the normality assumption is likely true.In this analysis, histograms were produced from standardized residual values, which is solely the raw residual divided by its standard error.The method of standardizing residuals is applied to transform data, such that its mean centers at zero and the standard deviation equals one.This technique also identified potential outliers, such that residuals with standard deviations less than −3 or greater than 3, could be identified.In this analysis, histograms from the datasets suggested that the residuals generally followed a normal distribution, and outliers with standard deviations greater than |3| were not detected.
the empirical cumulative distribution function with the theoretical cumulative distribution function.Moreover, this plot can also show whether the residuals or errors follow the assumed normality.Points on the p-p plot must follow the 1:1 diagonal line for normally distributed data.
Finally, the measured Vs and estimated Vs from each empirical correlation were compared, in order to assess the performance of the regression models (Figure 8).Data pairs mostly lay in the region between the lines 1:1.50 and 1:0.65.

Application for Site-Specific Characterization
One of the main applications of these empirical correlations is the estimation of the average shear wave velocity of the upper 30 m of a site (V s30 ).This metric is generally used for the classification of sites into the general characteristic of the subsurface, for earthquake engineering design.In the Philippines, V s30 is the most important parameter needed for the site classification scheme stipulated by the National Structural Code of the Philippines (2015) [44].This scheme classifies a site into six categories used for earthquake engineering design (Table 4).Soil requiring site-specific evaluation 1 Soil profile type S E also includes any soil profile with more than 3.0 m of soft clay defined as soil with plasticity index > 20, moisture content ≥ 40%, and undrained shear strength < 24 kPa.
The magnitude of V s30 can be computed in accordance with the expression where h i and v i are the layer thickness (m) and shear wave velocity (m/s) of the ith layer, respectively, in a total of n layers existing in the top 30-m of a site.The magnitude of the V s is estimated from SPT-N values using the proposed empirical relations applicable for all soil types.Sample boreholes with penetration depths reaching 30 m were examined, and the proposed empirical correlation was applied to the datasets (Table 5).The calculated V s30 s highlight the site classifications obtained from SPT data.Site Class D classification is dominant in Metro Manila, particularly on the western sector of the region along the coast.This suggests that the sites are generally underlain by stiff soil deposits.This is further supported by the Quaternary Alluvium deposits mapped along the western coast, proximal to Manila Bay.An additional utility of these correlations is in the refinement of the corresponding N-values for a specific range of shear wave velocities.Based on the outcomes of this study, a simplified correlation can be established for basic soil strength classification (i.e., soft soils, stiff soils, and very dense soils) (Table 6).This demonstrates the utility of the proposed empirical correlations for geotechnical description at site-specific scales, as well as seismic microzonation efforts in the study area.Specific ranges of V s have a corresponding range of N-values, where geophysical measurements of V s can complement point borehole profiles and establish the lateral continuity of subsurface material between drilling points.In situations where SPT boreholes might not be feasible, due to economic or logistic reasons, refraction microtremor or other geophysical tests that can reliably measure shear wave velocities and can serve as a proxy to determine the relative strength of soils, especially for smaller construction projects.This was demonstrated by Pancha and Apperley (2021) [73], through utilizing refraction microtremors in two-dimensional shear wave velocity mapping, thereby extending subsurface characterization and mapping lateral variations between spot boreholes.
Alternatively, shear wave velocity surveys might not be feasible in some urban areas, due to a lack of availability of open space to lay the linear geophone array, which is a technical requirement of some geophysical tests.These sites instead might have robust borehole survey results acquired before the building's construction.Compilation of these SPT datasets can be utilized and converted to shear wave values, to produce or refine seismic microzonation maps.The derived empirical correlations of V s and reliable static field data such as the SPT N-value can be an alternative way to determine these seismic parameters, such as the V s30 .
However, an important caveat, before the application of the proposed empirical relations, is to ensure the reliability of obtained SPT data to estimate the corresponding V s values.On the other hand, the processing of geophysical data to acquire V s must also be performed with utmost caution, to minimize uncertainty before applying it for basic soil strength classification.Furthermore, these empirical correlations must be applied with

Figure 1 .
Figure 1.(A) Geological and (B) geomorphologic map of Metro Manila, Philippines, and the locations of the refraction microtremor surveys.Geologic map modified from BMG (1983)[38] and geomorphologic map obtained from MMIERS (2004)[39].Locations of the testing sites are also depicted in black squares.

Figure 1 .
Figure 1.(A) Geological and (B) geomorphologic map of Metro Manila, Philippines, and the locations of the refraction microtremor surveys.Geologic map modified from BMG (1983)[38] and geomorphologic map obtained from MMIERS (2004)[39].Locations of the testing sites are also depicted in black squares.

Figure 2 .
Figure 2. (A) Sample of a generated p-f image from the ReMi surveys and the identified dispersion curve and picks (black box) in the ReMi Vspect module.(B) Fitting of dispersion curve model into the dispersion picks in the ReMi Disper module.(C) Corresponding one-dimensional shear wave velocity profile from the dispersion curve model.Vertical variations of shear wave velocities on the upper 20 m of the sites are illustrated in Figure3.For most of the sites, the shear wave velocity of the upper 20 m generally ranged from ~100 m/s to ~220 m/s, corresponding to SPT-N values typically ranging from 2-10.For sites with much stiffer strata, the Vs at some depths can reach up to ~250 m/s to ~500 m/s.This was in areas with higher recorded N-values, that ranged from 15 to 50.In terms of the average Vs of the upper 30 m of the sites (Vs30), most of the western portion of Metro Manila was classified into Site Class D, based on the scheme of the National Structural Code of the Philippines[48].This indicates that the area is generally underlain by stiff soil material, with a Vs30 ranging from 180-360 m/s.This was also observed in the data obtained from the refraction microtremor surveys, as the Vs30 of the sites mostly fell in the Site D range.

Figure 2 .
Figure 2. (A) Sample of a generated p-f image from the ReMi surveys and the identified dispersion curve and picks (black box) in the ReMi Vspect module.(B) Fitting of dispersion curve model into the dispersion picks in the ReMi Disper module.(C) Corresponding one-dimensional shear wave velocity profile from the dispersion curve model.

Figure 3 .
Figure 3. Variation of shear wave velocity (Vs) and penetration resistance (SPT-N) with depth for all twenty sites.Legend: blue lines for shear wave velocities; red lines for SPT-N values.

Figure 3 .
Figure 3. Variation of shear wave velocity (V s ) and penetration resistance (SPT-N) with depth for all twenty sites.Legend: blue lines for shear wave velocities; red lines for SPT-N values.

Table 2 .
Existing correlation between shear wave velocity and standard penetration resistance.

Table 3 .
Summary of the proposed correlations from this study.

Table 3 .
Summary of the proposed correlations from this study.

Table 4 .
Site classification scheme according to the National Structural Code of the Philippines, employing the V s30 metric to determine the general soil description and its typical corresponding N-values.

Table 5 .
Typical average shear wave velocity to a 30-m depth (V s30 ) of all soils at various locations in Metro Manila.

Table 6 .
Simplified correlation of V s and SPT-N for different soil strength classifications in this study.