李錫堤等人﹙2001﹚以區域地質資料及少數的地質鑽探資料，參考1997 UBC的分類方法，針對當時已有的708處臺灣強地動觀測網﹙TSMIP﹚自由場強震站做地盤分類，提出臺灣地區的地盤分類圖。自該場址分類提出後，氣象局至94年已完成其中266個測站之調查工作，包括鑽探調查與波速度的量測。因此本研究利用這一些實測資料，以及更多的鑽探調查資料，進一步評估及更新各站之場址分類。
地表下三十公尺以內的地層特性，為主要影響震波放大效應的關鍵，故近年來場址效應研究也多以地表下三十公尺平均剪力波速度﹙Vs30﹚為評估場址條件的重要指標。本研究主要針對臺灣Vs30的分佈進行研究；首先搜集並檢核強震測站鑽探資料，以重新檢查過的高品質剪力波速度資料進行與SPT-N值關係研究，建立剪力波速度推估模式。藉由此模式對地調所工程地質資料庫中的各個鑽探資料進行不同深度剪力波速度之推估，並運用地質統計方法進行空間資料分析，採Kriging with varying local means方法對各個鑽孔所計算的Vs30進行內插，完成全臺灣之Vs30分佈圖，提出以Vs30為基礎之全臺灣強震站場址分類結果。

Based on the available geologic geomorphologic data and some borehole data, and refered the 1997 Uniform Building Code (UBC) provisions, Lee et al. (2001) classified 708 free-field strong-motion station sites of Taiwan Strong-Motion Instrumentation Program (TSMIP). They also published a geologic site class map for whole Taiwan. Till year 2005, the Central Weather Bureau has finished core drilling and PS logging for 266 strong-motion sites. This study utilize the 266 drilling and logging data and more other new data to evaluate and revise the site classification, and update the site class map of Taiwan.
Because the property of the upper 30 meters soil profile is a key point to the amplification of ground-motion, most of the researches on site effect consider the average shear-wave velocity in upper 30 meters (Vs30) as an important factor. The purpose of this study was mainly on mapping the distribution of Vs30 in Taiwan. First, I collect and examine the drilling and PS logging data for the strong-motion station, and recheck these data. Second, by using these rechecked high quality shear-wave velocity data and SPT-N value, I build a model to estimate the shear-wave velocity from SPT-N value. Third, by using this model, I estimate the shear-wave velocity of each drilling data from the Central Geological Survey Engineering Geology Database. Finally, I use a geostatistical analysis method called “Kriging with verying local means method” to generate a Vs30 distribution map of Taiwan with the estimated Vs30 value from each drilling data. Furthermore, a site classification result of Taiwan’s strong-motion stations is updated base on the Vs30 map.

Anderson, J. G., Y. Lee, Y. Zeng, and S. Day﹙1996﹚Control of strong motion by the upper 30 meters, Bulletin of the Seismological Society of America, 86, 1749–1759.
Atkinson, G. M.﹙1993﹚Notes on ground motion parameters for Eastern North America: duration and H/V ration, Bulletin of the Seismological Society of America, 83, 1749–1759.
Boore, D. M., and W. B. Joyner﹙1989﹚Peak acceleration from the 17 October 1989 Loma Priera earthquake, Seismological Reasearch Letters, 60, 151–166.
Boore, D. M., W. B. Joyner, and T. E. Fumal﹙1993﹚. Estimation of response spectra and peak accelerations from Western North American earthquakes: an interim report, U. S. Geological survey Open-File Report, 93–509.
Borcherdt, R. D.﹙1970﹚Effects of local geology on ground motion near San Francisco Bay, Bulletin of the Seismological Society of America, 60, 29-61.
Borcherdt, R. D.﹙1994a﹚Estimates of site-dependent response spectra for design﹙methodology and justification﹚, Earthquake Spectra, 10, 617-653
Borcherdt, R. D.﹙1994b﹚An integrated methodology for estimates of site-dependent response spectra, seismic coefficients for site dependent building code provisions, and predictive GIS maps of strong ground shaking, in Proceedings of Seminar on New Developments in Earthquake Ground Motion Estimation and Implications for Engineering Design Practice, Applied Technology Council ATC 35-1, 10-1–10-44.
Building Seismic Safety Council﹙BSSC﹚﹙1998﹚1997 Edition NEHPR Recommended Provisions for Seismic Regulations for New Buildings and Other Structures, FEMA 302/303, Part 1﹙Provisions﹚and Part 2﹙Commentary﹚, developed for the Federal Emergency Management Agency, Washing, DC., 337p.
Building Seismic Safety Council﹙BSSC﹚﹙2004﹚2003 Edition NEHPR Recommended Provisions for Seismic Regulations for New Buildings and Other Structures, FEMA 450, Part 1﹙Provisions﹚, Prepared by the Building Seismic Safety Council for the Federal Emergency Management Agency, Washing, DC., 288p.
Castro, R. R., M. Mucciarelli, F. Pacor, and C. Petrungaro﹙1997﹚. S-wave site-response estimates using horizontal-to-vertical spectral rations, Bulletin of the Seismological Society of America, 87, 256–260.
Chen, M. H., K. L. Wen, C. H. Loh, and R. L. Nigbor﹙2004﹚Experience of Suspension P-S Logging Method and empirical formula of shear wave velocities in Taiwan, submitted to Journal of Soil Foundation.
Chiou, B., and K. L. Wen﹙2006﹚Estimation of VS30 at CWB’s Free-Field Sites, Chiou and Youngs PEER-NGA Empirical Ground Motion Motion Model for the Average Horizontal Component of Peak Acceleration and Pseudo-Spectral Acceleration for Spectral Periods of 0.01 to 10 Seconds, Appendix C, C-1–C-5.
Darragh, R. B. and A. F. Shakal﹙1991﹚The site response of two rock and soil station pairs to strong and week ground motion, Bulletin of the Seismological Society of America, 81, 1885–1899
Diehl J.﹙2001﹚GEOVision P-S Log Notes & Procedures。
Imai, T and R. Tonouchi﹙1982﹚Correction of N value with S-wave velocity, Proceeding of Second European Symposium on Penetration Test, Amsterdam,24-27 May, 67-72.
Imai, T.,﹙1977﹚P- and S-wave velocities of the ground in Japan, 9th ICSMFE, 257-260.
International Conference of Building Officials﹙ICBO﹚ ﹙1991﹚Uniform Building Code, Whittier, California., 1050 p.
International Conference of Building Officials﹙ICBO﹚ ﹙1997﹚Uniform Building Code, Whittier, California., 492 p.
Japan Road Association﹙JRA﹚ ﹙1996﹚ Japan Specifications for Highway Bridges, Part V: Seismic Design, Tokyo, Janpan.
Jinan, Z.﹙1987﹚Correlation between seismic wave velocity and the number of blow of SPT and depth, Selected Papers from the Chinese Journal of Geotechnical Engineering, ASCE, 92-100.
Joyner, W. B., and T. E. Fumal﹙1985﹚Predictive mapping of earthquake ground motion in Evaluating Earthquake Hazards in the Los Angeles Region: An Earth-Science Perspective, J. E. Ziony﹙Editor﹚, U. S. Geological survey professional paper, 1360, 203–220.
Lee, C. T., C. T. Cheng, C. W. Liao, and Y. B. Tsai﹙2001﹚Site classification of Taiwan Free-Field Strong-Motion Stations, Bulletin of the Seismological Society of America, 91, 1283–1297.
Lee, S.H.H﹙1988﹚Regression model for shear wave velocity of sandy and silty soils of Taipei basin, Proc. of the CCNAA-AIT Joint SeminarWorkshop on Research and Application for Multiple Hazards Mitigation, pp.489-494, Taipei.
Lew, M. and K. W. Campbell﹙1985﹚Relations between shear wave velocity and depth of overburden, Proc. Measurement and Use of Shear Wave Velocity for Evaluation Dynamic Properties, ASCE Convention in Denver, Colorado, 63-77.
Liao, S., and Whitman, R. V.﹙1985﹚Overburden Correction Factors for SPT in Sand, Journal of Geotechnical Engineering, ASCE, Vol. 114, No. GT3, 373–377.
Matsuoka, M., K. Wakamatsu, K. Fujimoto, and S. Midorikawa﹙2006﹚Structral Engineering, JSCE, 23, No.1, 57–68.
Mohraz, B.﹙1976﹚A study of earthquake response spectra for different geological conditions, Bulletin of the Seismological Society of America, , 66, 915–935.
Nakamura, Y.﹙1989﹚A method for dynamic characteristics estimation of subsurface using microtremor on the ground surface, Quarterly Report of Railway Technical Reaserch﹙RTRI﹚, 30, No 1.
Ohta, Y. and N. Goto﹙1978﹚Empirical shear wave velocity equations in terms of characteristic soil index, Earthquake Engineering and Structural Dynamics, 6, 167-187
Park, S., and S. Elrick﹙1998﹚Predictions of shear-wave velocities in southern California using surface geology, Bulletin of the Seismological Society of America, 88, 677–685.
Parkin, T. B.﹙1993﹚Spatial variability of microbial process in soil, Journal of Environmental Quality, 22, 409–417.
Phung, V., G. M. Atkinson, and D. T. Lau﹙2006﹚Methodology for site classification estimation using strong ground motion data from the Chi-Chi, Taiwan, Earthquake, Earthquake Spectra, 22, 511–531.
Schmertmann, J. H.﹙1978﹚Use the STP to measure dynamic soil properties？Yes, but…!, ASTM STP 654, 341-344.
Seed, H. B., C. Ugas, and J. Lysmer﹙1976﹚Site-dependent spectra for earthquake-resistant design, Bulletin of the Seismological Society of America, 66, 221–243.
Sokolov, V. Y., C. H. Loh, and W. Y. Jean﹙2007﹚Application of horizontal-to-vertical (H/V) Fourier spectral ratio for analysis of site effect on rock (NEHRP-class B) sites in Taiwan, Soil Dynamics and Earthquake Engineering, 27, 314-323
Theodulidis, N. P. Y. Bard, R. Archuleta, and M. Bouchon﹙1996﹚Horizontal-to-vertical spectral ratio and geological conditions: the case of Garner Balley downhole array in southern California, Bulletin of the Seismological Society of America, 86,306–319.
Wald, D. J., and T. Allen﹙2006﹚Topographic slope as a proxy for seismic site conditions and amplification, Submitted for publication in Bulletin of the Seismological Society of America.
Yamazaki, F., and M. A, Ansary﹙1997﹚Horizontal-to-vertical spectrum ratio of earthquake ground motion for site characterization, Earthquake Engineering & Structural Dynamics, 26, 671–689.
Yamazaki, F., K. Wakamatsu, J. Onishi, K. T. Shabestari﹙2000﹚Relationship between geomorphological land classification ad site amplification ratio based on JMA strong motion records, Soil Dynamics and Earthquake Engineering, 19, 41–53.
Zare, M., P. Y. Bard, and M. Ghagory-Ashtiany﹙1999﹚Site characterizations for the Iranian strong motion network, Soil Dynamics and Earthquake Engineering, 18, 101–123.
李咸亨、吳志明﹙1990﹚台北盆地土壤之動態性質研究﹙III﹚~下井探測法與剪力波速度