SlideShare a Scribd company logo
S-wave Velocity Structure of Mexico
City Obtained from Three-component
Microtremor Measurements and
Microtremor Array Measurements
1
Koichi Hayashi*, Geometrics
Atsushi Nozu, Port and Airport Research Institute,
Masanori Tanaka, Port and Airport Research Institute
Haruhiko Suzuki, OYO Corporation
Efraín Ovando Shelley, Universidad Nacional Autonoma
de Mexico
Outline
2
• Introduction
• Investigation site
• Data acquisition
Equipment
H/V spectrum
Dispersion curve
• Analysis results
• Comparison with United States
• Conclusions
Introduction
3
• The earthquake that struck Mexico on 19 September
1985 caused severe damage in Mexico City although
the city is located 400km away from the epicenter.
• The main reason for this damage is that the city is
located on a basin filled with very soft sediments.
• Distribution of these soft sediments has been
delineated by drillings and microtremor
measurements.
• A small number of attempts have been made to image
the S-wave velocity structure of the basin using
downhole seismic loggings.
• In order to delineate S-wave velocity structure of the
basin down to depth of approximately 200m, we have
performed three-component micro-tremor
measurements and microtremor array measurements.
1985 Mexico Earthquake
Mexico City
400km
Mw=8.3
4
1985 Mexico Earthquake
5
• Mw=8.3
• Most-often cited number of deaths is an
estimated 10,000 people but experts agreed
that it could be up to 40,000.
• Damage area corresponds to the western part
of the lake zone within 2 to 4 kilometers of the
Alameda Central.
• 6 to 15 story buildings are mainly damaged in
the city due to a frequency range of 0.25 to
0.5Hz (period of 2 to 4seconds)
> 2sec
Natural Period (H/V) of Mexico Basin
Lermo and Chavez-
Garcia (1994)
Central Mexico City
Airport
Damage Area
> 4sec
6
Investigation Site
7
• Investigation site is placed at the downtown of
Mexico City.
• 30km length survey line crosses the basin with a
west-southwest to east-northeast direction.
• 3 component microtremor measurements were
performed at more than 10 sites on the line.
• Microtremor array measurements were
performed at 6 sites on the line.
• Microtremor array measurements used 25 to
650m equilateral triangular arrays.
Chapultepec
Aragon
Texcoco No.7
Almeda
Texcoco No.8
Texcoco TXC
Investigation Site
Central Mexico City
Airport
8
Example of Array Configuration
19.422938
-99.182564
19.423838
99.184163
19.423818
99.180783
19.42138
99.182296
330m 19.462283
99.067562
19.463962
99.066875
19.4629
99.063903
19.464731
99.069911
19.461818
99.06931
19.460988
99.066317
19.459238
99.068806
650m
Chapultepec
9
Aragon
Data Acquisition
10
• Data acquisition was carried out during the daytime in
December 2008 and December 2009.
• Microtremor measurement systems (JU210) made by
Hakusan Corporation and data loggers (GPL-6A3P) made by
Mitsutoyo Corporation were mainly used for data
acquisition.
• Both systems use accelerometers for the sensors.
• In order to verify applicability of the accelerometers, servo-
type velocity meters made by Katsujima Corporation (SD-
110) and Tokyo Sokushin Corporation (VSE11F, VSE12F)
were also used in the 3 component microtremor
measurements.
• H/V spectra obtained through the accelerometers and the
velocity meters were compared.
• 30 min. to 1 hour of microtremors were recorded for each
three component measurement or array measurement.
Equipment(accelerometer)
Hakusan
11
Mitsutoyo
Comparison of H/V
It is clear that peak frequency of H/V decrease from west to east
(from edge to center of basin) 12
600
500
,-...
(/')
.'..E....., 4oo
>-.+-'
(.)
0 300(J)
>I
(J)
200
a..
100
0
0.1
• Texcoco Main No.7(2008)
• Texcoco North No.8(2008)
• Aragon(2009)
• Almeda(2009)
<>Chapultepec(2009)
 c
1
Frequency(Hz)
10 100
Sites where the peak frequency of H/V spectra is higher,
the phase velocity of the dispersion curve is also higher 13
,
•
9
<>
:
-
·l.llt>
.o.c.
Analysis (1)
14
• A joint inversion was applied to the observed
H/V spectra and dispersion curves, and S-wave
velocity models were analyzed for six sites.
• In the inversion, phase velocities of the
dispersion curves and peak frequencies of the
H/V spectra were used as the observation
data.
• Unknown parameters were layer thickness and
S-wave velocity.
• A Genetic algorithm was used for
optimization.
Analysis (2)
15
• Initial models were created by a simple
wavelength transformation in which
wavelength calculated from phase velocity
and frequency is divided by three and plotted
at depth.
• Theoretical H/V spectra and phase velocities
are generated by calculating the weighted
average of the fundamental mode and higher
modes (up to the 4th modes) based on
medium response.
• Rayleigh-Love ratio (R/L) is fixed as 0.7
S-wave Velocity Model(East)
Texcoco No.7
Aragon 0
50
100
150
200
250
300
350
400
450
500
0 0.2 0.4 0.6
Vs(km/ s)
Depth(m
)
S-wave Velocity Model(West)
Almeda
Chapultepec
17
Chapultepec
Aragon
Texcoco No.7
Almeda
Texcoco No.8
Texcoco TXC
Peak Frequency of H/V
and S-wave Velocity Model
0.7s
90m/s
1.5s
80m/s
3.9s
60m/s
3.3s
30m/s
1.1s
60m/s
>400m/s450m
<100m/s
≒200m/s
18
Comparison with United States
H/V Spectra Dispersion Curve S-wave Velocity Model
S-wave velocity of Mexico is extremely low
compare with San Jose and Redwood City.
Mexico City
San Jose (William St. Park), CA
Redwood City, CA
19
0
50
100
150
200
250
300
350
0.1 1 10 100
Phase-velocity(m/sec)
Frequency(Hz)
Mexico
San Jose
Redwood City
0.1
1
10
0.1 1 10
H/Vspectra
Period(sec)
Mexico
San Jose
Redwood City
0
20
40
60
80
100
120
140
160
0
S-wave velocity (m/sec)
100 200 300 400 500
Depth(m)
Mexico
San Jose
Redwood City
Conclusions
• We have performed the three-component
microtremor measurements and microtremor
array measurements in the Mexico basin and
estimated the S-wave velocity models down to a
depth of 200m.
• S-wave velocity in the middle of the Mexico basin
is lower than 150m/s to a depth of 70m and
much lower than typical alluvial plains in Japan
and United States.
• Peak frequencies of the H/V spectra in Mexico
City vary from 0.25 to 1Hz and it seems that
these peak frequencies are mainly due to the
low-velocity layer shallower than a depth of
100m.
20

More Related Content

PPTX
2D Microtremor Array
PPT
Ieee gold2010 sermi
PPT
IGARSS_xqwu_2011.ppt
PPT
Fukao Plenary.ppt
PPTX
Radiation Monitoring Data from Fukushima Area 05/13/2011
PPTX
Design and working of SkyTem
PPT
Application of Ground Penetrating Radar in Subsurface mapping
PDF
artigo publicado na IEEE Transactions on Arospace and Electronic Systems vol ...
2D Microtremor Array
Ieee gold2010 sermi
IGARSS_xqwu_2011.ppt
Fukao Plenary.ppt
Radiation Monitoring Data from Fukushima Area 05/13/2011
Design and working of SkyTem
Application of Ground Penetrating Radar in Subsurface mapping
artigo publicado na IEEE Transactions on Arospace and Electronic Systems vol ...

What's hot (20)

PPTX
ground penetrating radar
PDF
636908main paul presentation
PPT
NineYearsofAtmosphericRemoteSensingwithSCIAMACHY-InstrumentPerformance.ppt
PPTX
GROUND PENETRATING RADAR(GPR) ppt
PPTX
Oceansat
PPTX
LAICE Overview - NGC Presentation 12-02-15
PPTX
Ocean- sat for Oceanography
PPTX
SIXTEEN CHANNEL, NON-SCANNING AIRBORNE LIDAR SURFACE TOPOGRAPHY (LIST) SIMULATOR
PPT
FR2.L09 - PROCESSING OF MEMPHIS MILLIMETER WAVE MULTI-BASELINE INSAR DATA
PDF
Simulating tropical meteorology for air quality studies
PPT
TU3.T10.3.ppt
PDF
AGU2014-SA31B-4098
PDF
FAST実験6:新型大気蛍光望遠鏡による観測報告とピエールオージェ観測所への設置計画
PPTX
PDF
2007 t hz astronomy with sofia
PPS
Civilex Presentation Laser Scanning And Gpr March 2009 Ices
PDF
Comparing pulse doppler lidar with sodar and direct measurements for wind ass...
PDF
Setellite image sensors
PDF
Khara- Final-AMS- Poster
PDF
DUAL BAND GNSS ANTENNA PHASE CENTER CHARACTERIZATION FOR AUTOMOTIVE APPLICATIONS
ground penetrating radar
636908main paul presentation
NineYearsofAtmosphericRemoteSensingwithSCIAMACHY-InstrumentPerformance.ppt
GROUND PENETRATING RADAR(GPR) ppt
Oceansat
LAICE Overview - NGC Presentation 12-02-15
Ocean- sat for Oceanography
SIXTEEN CHANNEL, NON-SCANNING AIRBORNE LIDAR SURFACE TOPOGRAPHY (LIST) SIMULATOR
FR2.L09 - PROCESSING OF MEMPHIS MILLIMETER WAVE MULTI-BASELINE INSAR DATA
Simulating tropical meteorology for air quality studies
TU3.T10.3.ppt
AGU2014-SA31B-4098
FAST実験6:新型大気蛍光望遠鏡による観測報告とピエールオージェ観測所への設置計画
2007 t hz astronomy with sofia
Civilex Presentation Laser Scanning And Gpr March 2009 Ices
Comparing pulse doppler lidar with sodar and direct measurements for wind ass...
Setellite image sensors
Khara- Final-AMS- Poster
DUAL BAND GNSS ANTENNA PHASE CENTER CHARACTERIZATION FOR AUTOMOTIVE APPLICATIONS
Ad

Viewers also liked (6)

PPT
My Pathway in Science and Teaching
DOC
Effect of damper on sesmic demand ond diffrent soil with help of staad
PPTX
Seismic Reflection
PPT
Seismic survey
PPTX
New introduction to seismic method
PPTX
Introduction to Seismic Method
My Pathway in Science and Teaching
Effect of damper on sesmic demand ond diffrent soil with help of staad
Seismic Reflection
Seismic survey
New introduction to seismic method
Introduction to Seismic Method
Ad

Similar to S-wave Velocity Based On Microtremor Array (12)

PDF
10.1785@0120180100
PPTX
Seismic microzonation of hamirpur city using microtremor measurements
PPT
9oct 1 petronio-reflection seismic
PDF
Kent Odaklı Mikrobölgeleme: Sivas Örneği
DOC
Vrancea1 cor 6 6 2016 (2)
 
PPTX
Seismic Hazard Assessment
PPT
Caribbean3.ppt
PDF
Micro Earthquake Survey Design, Optimization with modeling
PDF
Deprem Tehlike Potansiyeli Araştırması: Antakya Örneği
PDF
Abimbola_NMO-1
PDF
Seismic Microzonation - Principles and Methodology
PDF
Alwyn Seismic Dataset Mini project- Refael (24081)
10.1785@0120180100
Seismic microzonation of hamirpur city using microtremor measurements
9oct 1 petronio-reflection seismic
Kent Odaklı Mikrobölgeleme: Sivas Örneği
Vrancea1 cor 6 6 2016 (2)
 
Seismic Hazard Assessment
Caribbean3.ppt
Micro Earthquake Survey Design, Optimization with modeling
Deprem Tehlike Potansiyeli Araştırması: Antakya Örneği
Abimbola_NMO-1
Seismic Microzonation - Principles and Methodology
Alwyn Seismic Dataset Mini project- Refael (24081)

More from Ali Osman Öncel (20)

PDF
Topografik ve sismolojik VS30 verilerinin karşılaştırılması; Balıkesir örneğ
PDF
APA Yazım Kuralları
PPTX
Gravimetri : Ders 14
PPTX
Gravimetri : Ders 13
PPTX
Gravimetri : Ders 12
PDF
Riskli Yapılar - Çevre ve Şehircilik
PPTX
Riskli Yapılar -Çevre ve Şehircilik
PPTX
Gravimetri : Ders 07
PPTX
Gravimetri : Ders 06
PPTX
Gravimetri: Ders 05
PPTX
Gravimetri : Ders 04
PPTX
Gravimetri : Ders 03
PPTX
Gravimetri Ders 02
PPTX
Gravimetri Ders 01
PDF
Kar Kar Geothermal Field Work
PDF
Beppu geothermal field
PDF
High Resolution Earth's Gravitational Field
PDF
Gravity Predictions for Earthquakes
PDF
Nakamura Technique for Soil Characterization
PDF
H/V User Guidelines
Topografik ve sismolojik VS30 verilerinin karşılaştırılması; Balıkesir örneğ
APA Yazım Kuralları
Gravimetri : Ders 14
Gravimetri : Ders 13
Gravimetri : Ders 12
Riskli Yapılar - Çevre ve Şehircilik
Riskli Yapılar -Çevre ve Şehircilik
Gravimetri : Ders 07
Gravimetri : Ders 06
Gravimetri: Ders 05
Gravimetri : Ders 04
Gravimetri : Ders 03
Gravimetri Ders 02
Gravimetri Ders 01
Kar Kar Geothermal Field Work
Beppu geothermal field
High Resolution Earth's Gravitational Field
Gravity Predictions for Earthquakes
Nakamura Technique for Soil Characterization
H/V User Guidelines

Recently uploaded (20)

PPTX
Introduction-to-Literarature-and-Literary-Studies-week-Prelim-coverage.pptx
PPTX
Orientation - ARALprogram of Deped to the Parents.pptx
PPTX
master seminar digital applications in india
PDF
LDMMIA Reiki Yoga Finals Review Spring Summer
PDF
Yogi Goddess Pres Conference Studio Updates
PDF
Anesthesia in Laparoscopic Surgery in India
PDF
Supply Chain Operations Speaking Notes -ICLT Program
PDF
Black Hat USA 2025 - Micro ICS Summit - ICS/OT Threat Landscape
PDF
2.FourierTransform-ShortQuestionswithAnswers.pdf
PPTX
UV-Visible spectroscopy..pptx UV-Visible Spectroscopy – Electronic Transition...
PDF
What if we spent less time fighting change, and more time building what’s rig...
PDF
RTP_AR_KS1_Tutor's Guide_English [FOR REPRODUCTION].pdf
PDF
Weekly quiz Compilation Jan -July 25.pdf
PPTX
school management -TNTEU- B.Ed., Semester II Unit 1.pptx
PDF
RMMM.pdf make it easy to upload and study
PDF
Chinmaya Tiranga quiz Grand Finale.pdf
PPTX
Final Presentation General Medicine 03-08-2024.pptx
PPTX
1st Inaugural Professorial Lecture held on 19th February 2020 (Governance and...
PDF
Updated Idioms and Phrasal Verbs in English subject
PDF
Complications of Minimal Access Surgery at WLH
Introduction-to-Literarature-and-Literary-Studies-week-Prelim-coverage.pptx
Orientation - ARALprogram of Deped to the Parents.pptx
master seminar digital applications in india
LDMMIA Reiki Yoga Finals Review Spring Summer
Yogi Goddess Pres Conference Studio Updates
Anesthesia in Laparoscopic Surgery in India
Supply Chain Operations Speaking Notes -ICLT Program
Black Hat USA 2025 - Micro ICS Summit - ICS/OT Threat Landscape
2.FourierTransform-ShortQuestionswithAnswers.pdf
UV-Visible spectroscopy..pptx UV-Visible Spectroscopy – Electronic Transition...
What if we spent less time fighting change, and more time building what’s rig...
RTP_AR_KS1_Tutor's Guide_English [FOR REPRODUCTION].pdf
Weekly quiz Compilation Jan -July 25.pdf
school management -TNTEU- B.Ed., Semester II Unit 1.pptx
RMMM.pdf make it easy to upload and study
Chinmaya Tiranga quiz Grand Finale.pdf
Final Presentation General Medicine 03-08-2024.pptx
1st Inaugural Professorial Lecture held on 19th February 2020 (Governance and...
Updated Idioms and Phrasal Verbs in English subject
Complications of Minimal Access Surgery at WLH

S-wave Velocity Based On Microtremor Array

  • 1. S-wave Velocity Structure of Mexico City Obtained from Three-component Microtremor Measurements and Microtremor Array Measurements 1 Koichi Hayashi*, Geometrics Atsushi Nozu, Port and Airport Research Institute, Masanori Tanaka, Port and Airport Research Institute Haruhiko Suzuki, OYO Corporation Efraín Ovando Shelley, Universidad Nacional Autonoma de Mexico
  • 2. Outline 2 • Introduction • Investigation site • Data acquisition Equipment H/V spectrum Dispersion curve • Analysis results • Comparison with United States • Conclusions
  • 3. Introduction 3 • The earthquake that struck Mexico on 19 September 1985 caused severe damage in Mexico City although the city is located 400km away from the epicenter. • The main reason for this damage is that the city is located on a basin filled with very soft sediments. • Distribution of these soft sediments has been delineated by drillings and microtremor measurements. • A small number of attempts have been made to image the S-wave velocity structure of the basin using downhole seismic loggings. • In order to delineate S-wave velocity structure of the basin down to depth of approximately 200m, we have performed three-component micro-tremor measurements and microtremor array measurements.
  • 4. 1985 Mexico Earthquake Mexico City 400km Mw=8.3 4
  • 5. 1985 Mexico Earthquake 5 • Mw=8.3 • Most-often cited number of deaths is an estimated 10,000 people but experts agreed that it could be up to 40,000. • Damage area corresponds to the western part of the lake zone within 2 to 4 kilometers of the Alameda Central. • 6 to 15 story buildings are mainly damaged in the city due to a frequency range of 0.25 to 0.5Hz (period of 2 to 4seconds)
  • 6. > 2sec Natural Period (H/V) of Mexico Basin Lermo and Chavez- Garcia (1994) Central Mexico City Airport Damage Area > 4sec 6
  • 7. Investigation Site 7 • Investigation site is placed at the downtown of Mexico City. • 30km length survey line crosses the basin with a west-southwest to east-northeast direction. • 3 component microtremor measurements were performed at more than 10 sites on the line. • Microtremor array measurements were performed at 6 sites on the line. • Microtremor array measurements used 25 to 650m equilateral triangular arrays.
  • 8. Chapultepec Aragon Texcoco No.7 Almeda Texcoco No.8 Texcoco TXC Investigation Site Central Mexico City Airport 8
  • 9. Example of Array Configuration 19.422938 -99.182564 19.423838 99.184163 19.423818 99.180783 19.42138 99.182296 330m 19.462283 99.067562 19.463962 99.066875 19.4629 99.063903 19.464731 99.069911 19.461818 99.06931 19.460988 99.066317 19.459238 99.068806 650m Chapultepec 9 Aragon
  • 10. Data Acquisition 10 • Data acquisition was carried out during the daytime in December 2008 and December 2009. • Microtremor measurement systems (JU210) made by Hakusan Corporation and data loggers (GPL-6A3P) made by Mitsutoyo Corporation were mainly used for data acquisition. • Both systems use accelerometers for the sensors. • In order to verify applicability of the accelerometers, servo- type velocity meters made by Katsujima Corporation (SD- 110) and Tokyo Sokushin Corporation (VSE11F, VSE12F) were also used in the 3 component microtremor measurements. • H/V spectra obtained through the accelerometers and the velocity meters were compared. • 30 min. to 1 hour of microtremors were recorded for each three component measurement or array measurement.
  • 12. Comparison of H/V It is clear that peak frequency of H/V decrease from west to east (from edge to center of basin) 12
  • 13. 600 500 ,-... (/') .'..E....., 4oo >-.+-' (.) 0 300(J) >I (J) 200 a.. 100 0 0.1 • Texcoco Main No.7(2008) • Texcoco North No.8(2008) • Aragon(2009) • Almeda(2009) <>Chapultepec(2009) c 1 Frequency(Hz) 10 100 Sites where the peak frequency of H/V spectra is higher, the phase velocity of the dispersion curve is also higher 13 , • 9 <> : - ·l.llt> .o.c.
  • 14. Analysis (1) 14 • A joint inversion was applied to the observed H/V spectra and dispersion curves, and S-wave velocity models were analyzed for six sites. • In the inversion, phase velocities of the dispersion curves and peak frequencies of the H/V spectra were used as the observation data. • Unknown parameters were layer thickness and S-wave velocity. • A Genetic algorithm was used for optimization.
  • 15. Analysis (2) 15 • Initial models were created by a simple wavelength transformation in which wavelength calculated from phase velocity and frequency is divided by three and plotted at depth. • Theoretical H/V spectra and phase velocities are generated by calculating the weighted average of the fundamental mode and higher modes (up to the 4th modes) based on medium response. • Rayleigh-Love ratio (R/L) is fixed as 0.7
  • 16. S-wave Velocity Model(East) Texcoco No.7 Aragon 0 50 100 150 200 250 300 350 400 450 500 0 0.2 0.4 0.6 Vs(km/ s) Depth(m )
  • 18. Chapultepec Aragon Texcoco No.7 Almeda Texcoco No.8 Texcoco TXC Peak Frequency of H/V and S-wave Velocity Model 0.7s 90m/s 1.5s 80m/s 3.9s 60m/s 3.3s 30m/s 1.1s 60m/s >400m/s450m <100m/s ≒200m/s 18
  • 19. Comparison with United States H/V Spectra Dispersion Curve S-wave Velocity Model S-wave velocity of Mexico is extremely low compare with San Jose and Redwood City. Mexico City San Jose (William St. Park), CA Redwood City, CA 19 0 50 100 150 200 250 300 350 0.1 1 10 100 Phase-velocity(m/sec) Frequency(Hz) Mexico San Jose Redwood City 0.1 1 10 0.1 1 10 H/Vspectra Period(sec) Mexico San Jose Redwood City 0 20 40 60 80 100 120 140 160 0 S-wave velocity (m/sec) 100 200 300 400 500 Depth(m) Mexico San Jose Redwood City
  • 20. Conclusions • We have performed the three-component microtremor measurements and microtremor array measurements in the Mexico basin and estimated the S-wave velocity models down to a depth of 200m. • S-wave velocity in the middle of the Mexico basin is lower than 150m/s to a depth of 70m and much lower than typical alluvial plains in Japan and United States. • Peak frequencies of the H/V spectra in Mexico City vary from 0.25 to 1Hz and it seems that these peak frequencies are mainly due to the low-velocity layer shallower than a depth of 100m. 20