13th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2004 Paper No. 2840

RECENT ADVANCES IN SEISMIC MICROZONATION STUDIES IN : THE MANIZALES CITY CASE

Luis YAMÍN1, Omar D. CARDONA2, Mauricio GALLEGO1, Camilo PHILLIPS1, Santiago ARÁMBULA1

SUMMARY

The article presents the recent advances on the seismic microzonation studies that have been carried out for some Colombian cities, making emphasis in the microzonation study for the city of Manizales, which is located on volcanic ash soils. An analysis of the Colombian territory's seismological variables was carried out and the attenuation laws for different types of seismological sources were defined based on regional recent strong motion records. The non-linear dynamic properties of the city's characteristic soils were studied using in situ and laboratory dynamic tests. A non-linear transfer function evaluation model for the volcanic ash deposits was developed in order to obtain Fourier Amplitude Spectra (FAS) at ground surface level. Using the FAS at surface level, uniform hazard spectra were constructed at different locations according to the classic hazard theory. Simultaneously, sensibility analyses for the different variables involved were carried out as well as Montecarlo simulations of non -stationary processes. Bidimensional analyses were also performed in order to evaluate the amplification effects generated by abrupt topographic characteristics of the ground on which the city has been built. Artificial earthquakes were generated for the time response analysis of the soil deposits and representative structures using the seismological theory of the Omega Squared model following a Green empirical function procedure. Recommendations for the installation of a seismic instrumentation network were made in order to improve the proposed models. The spectra establish the basis for the city's seismic design code and city's earthquake loss scenarios for emergency response planning. Finally, taking into account disaster risk management and insurance public policies, a seismic risk evaluation was performed on a representative group of government's building infrastructure in order to establish the optimum seismic risk transfer and retention strategies for the city.

INTRODUCTION

The three major cities in the Colombian Coffee Growing Region are, in their order, Pereira (380.000 inhabitants), Manizales (370.000 inhabitants) and Armenia (270.000 inhabitants). Considering the high seismic activity characteristic of the zone, these three cities are undertaking, from some years ago, seismic microzonation studies aimed to estimate the expected seismic response of the several distinctive zones for purposes of structural design and rehabilitation, to estimate damage and loss future scenarios and to

1 Universidad de Los Andes, CITEC, Bogotá, Colombia, [email protected] 2 Universidad de Los Andes, CEDERI, Bogotá, Colombia, [email protected] develop prevention and emergency management plans. The following general methodology has been used to carry on the above mentioned microzonation studies:

ƒ Conduction of basic studies, including the following: regional and local geology, neotechtonics, historical and instrumental seismicity, accelerographic local network set up, geophysical characterization (seismic refraction and microtremors), geotechnical characterization by compiling previous studies, geomorphology, topography, and cartography. ƒ Assessment of local seismic hazard and definition of design earthquakes, based on existing studies at national level [1, 2, 3, 4, 5] and the existing instrumental and historical seismic information. ƒ Detailed geotechnical investigation by means of borings in the several different characteristic zones of the city, field tests to measure wave velocities and laboratory and in situ tests for soil dynamic behavior characterization. ƒ Assessment of dynamic response of the different characteristic deposits by using unidimensional analytical models. ƒ Evaluation of dynamic response related to geometric effects such as fills, slopes and hill zones, by using two-dimensional analytical models. ƒ Definition of earthquake-resistant uniform hazard design spectra. ƒ Instrumentation (using accelerographic networks) in order to calibrate the analytical results. ƒ Seismic risk evaluation ƒ Seismic risk management strategies

SEISMIC HAZARD ANALYSIS

An analysis of the Colombian territory's seismological variables was carried out. Figure 1 presents the analytical model developed to evaluate the seismic hazard for Colombia. The most important seismic sources are shown. Each one of them is characterized using the seismic information available from the national accelerographic network.

Caicedo Copacabana a Caracoli Bello San Vicente Landazuri Itagui Peñol Guatape Urrao n SANTANDER Ebejico San Carlosi Granada t Pto Nare M El Peñon Betulia Medellín s Amaga Cocorna Bolivar La Union SanLuis e Salgar venecia La Ceja l Pto. Boyaca ur Pto Triunfo Montebello La Belleza Fredonia a Bolivar San Francisco indo Santa Barbara BOYACA ANTIOQUIA Abejorral P Albania Puente Nacional Jericó s Sonson Otanche Argelia Pauna aSaboya jardin l Tunungua

r Andes a Yacopi Quibdó La Dorada g Caldas Aguadas l Pauna n LLoro a Nariño a La Victoria S i

Caramanta Samana o Maripi t l Pacora r r e Bagadó u Simijaca Supia e P a CHOCÓ Riosucio Pensilvania Muzo Susa Marmato Salamina Caparrapi B S Mistrato Paime Fuquene m CALDAS Victoria Ubate

La Palma Cucunaba a Marquetalia San Pablo Quinchia Filadelfia Supatasa Tadó o Topaipi Pto Rico Marulanda Villagoméz Arnzazu Choconta h Utica Nimaima Manzanares Tausa RISARALDA Anserma R Mariquita Neira Macheta

i Fresno Apia CUNDINAMARCA Honda Guaduas Pacho

a Viterbo Nocaima Suesca Manizales Falan Condoto Supata Herveo Santuario Palestina Libano Armero Sasaima Nemocon Sesquile S San Francisco La Celia Villamaria Chaguani Alban Tabio Guatavita

o San José del Palmar Tocancipa Lerida Bituima Sopó Manta Novita El Aguila Balboa Guayabal Itsmina Tenjo Chia l a Dosquebradas Murillo VilletaSan Juan Santa Rosa de Cabal Facatativa Guasca a Zipacon Cota Gacheta Ansermanuevoc Pereira Santa Rosa de Cabal Ambalema Madrid El Cairo Quipile Funza n La Calera Argelia u Santa Isabel Venadillo Anolaima Bojaca Puli Sipi Cartago Alcala Mosquera o Toro Tena Salento Beltran Versalles Anzoategui Choachi Filandia San Antonio La Victoria La Mesa Soacha B Santa Isabel Jerusalen El Colegio Fomeque LaCa Union Ubaque Armenia Coello Granada El Dovio Calarca Apulo Sibate .

La Victoria e C Anapoima . Chipaque QUINDIO Silvania Cajamarca Ibague Guataqui D Caqueza Meta

n Cordoba s Nariño á a Zarzal t Une t Agua de Dios Fusagasuga o Fosca Roldanillo Caicedonia Pijao g n Girardot Quetame a e Tibacuy Pasca o

p Melgar B Genova u a Rovira Arbelaez r Sevillai g Gutierrez Trujillo Bugalagrande Flandes Venecia Guayabetal r o a Pandi b Espinal San Bernardo a

Ga f Roncesvalles Tulua I Icononzo Acacias Villavicenc

a San Luis Venecia

Guamon La Victoriac I Riofrio Carmen de Apicala e

n TOLIMA San Pedro Guamal u San Antonio l Villarica La Victoria Ortega Saldaña Cabrera

a t

a META Buga e Purificación Restrepo d Buenaventura C r Chaparral El Castillo Dagua Guacari Coyaima Prado Vijes l g Ginebra m a VALLE DEL t Lejanias a El Cerrito Colombia

CAUCA e n Rioblanco Alpujarra Dolores Yumbo

Palmira d M HUILA Ataco La Uribero

Cali i Natagaima

a F La Victoria La Victoria

Figure 1: Analytical model for seismic hazard analysis

Attenuation laws for different types of seismological sources were defined based on regional recent strong motion records. Using a probabilistic integration scheme, exceedence rates for firm soil deposits are estimated for the city. Figure 2 illustrates the different acceleration exceedence rates for firm soil for each one of the seismic sources of the zone. In addition, the maximum acceleration exceedance rates considering the integrated effects of all relevant seismic sources is presented as well as the relative participation factor for each one of them in the global hazard curve for the city.

Tasas de Excedencia en Roca para Manizales 1 Participación de las Fuentes para un Tr = 475 años

70%

0 60%

50%

40%

0 30% Participación 20%

10% 0 0% 110100Aceleración (gal)

Total Subducción Benioff I Benioff I Cauca Romeral Murindo Palestina Benioff II Benioff Bennioff I Benioff III

Benioff III Benioff Prof Cauca Frontal Subducción Benioff Prof Benioff Fuente Ibague Palestina Romeral Figure 2: Exceedance rates for relevant seismic sources and relative participation factors

Following the general random vibration theory, different seismic intensity maps such us acceleration, velocity, displacement or any other related parameter, are plotted for different return periods. Figure 3 presents the maximum probable acceleration in firm soil for a 475 year return period, corresponding to the design return period established by local structural design codes. Simultaneously, uniform hazard spectra are presented for the city, for different periods of analysis return.

Aceleración Máxima del suelo; Tr=475 años, (gal)

Caicedo Espectros de Amenaza Uniforme; Manizales Copacabana Caracoli Bello San Vicente Landazuri Itagui Peñol Guatape Urrao Ebejico San Carlos SANTANDER Granada Pto Nare El Peñon 700 Betulia Medellín Amaga Cocorna Bolivar La Union SanLuis Salgar venecia La Ceja Pto. Boyaca 6.0 Pto Triunfo Fredonia Montebello La Belleza Bolivar San Francisco Santa Barbara BOYACA Puente Nacional Tr=50 años ANTIOQUIA Abejorral Albania Jericó 600

Sonson Otanche Argelia Pauna Saboya jardin Tunungua Andes La Dorada r Yacopi Caldas Quibdó Aguadas a Pauna g LLoro Nariño l La Victoria a

S Caramanta Tr=100 años Samana Maripi

o Pacora t r

e Bagadó u Simijaca Supia P 5.5 CHOCÓ Riosucio Pensilvania Muzo Caparrapi Susa Marmato Salamina 500 Fuquene Mistrato Paime CALDAS Victoria Ubate La Palma Filadelfia Marquetalia Cucunaba San Pablo Quinchia Supatasa Tr=200 años Tadó Pto Rico Topaipi Villagoméz Marulanda Choconta Arnzazu Manzanares Utica Nimaima Mariquita Tausa RISARALDA Anserma Neira Fresno CUNDINAMARCA Macheta Apia Honda Guaduas Viterbo Manizales Nocaima Pacho Suesca Falan 400 Condoto Supata Herveo Tr=475 años Santuario Palestina Libano Armero Sasaima Nemocon Sesquile San Francisco La Celia Villamaria Chaguani 5.0 Alban Guatavita Tabio Tocancipa San José del Palmar Lerida Bituima Sopó Balboa Guayabal Manta Itsmina Novita El Aguila Villeta Tenjo Chia Dosquebradas Murillo San Juan Facatativa Guasca Santa Rosa de Cabal Cota Zipacon Gacheta Ansermanuevo Pereira Santa Rosa de Cabal Madrid Tr=1000 años Ambalema El Cairo Quipile Funza 300 Latitud Venadillo La Calera Argelia Santa Isabel Anolaima Bojaca Puli Sipi Cartago Alcala Mosquera Toro Tena Salento Beltran Versalles Anzoategui Choachi Filandia San Antonio La Victoria La Mesa Soacha Santa Isabel Jerusalen El Colegio Fomeque La Union Armenia Coello Ubaque El Dovio Calarca Apulo Granada Sibate La Victoria

. 4.5 Anapoima Chipaque QUINDIO C Cajamarca Ibague Guataqui Silvania . Caqueza Meta D 200 Cordoba Nariño á Une Zarzal t Fosca El Calvario Caicedonia Pijao Agua de Dios Fusagasuga o Roldanillo g Girardot Tibacuy Pasca Quetame o

Melgar B Genova Rovira Arbelaez Sevilla Gutierrez Trujillo Bugalagrande Flandes Venecia Guayabetal Pandi

Espinal San Bernardo Pseudoaceleración (gal) Tulua Roncesvalles Icononzo Acacias San Luis Venecia 100 La Victoria Guamo Carmen de Apicala Riofrio San Pedro TOLIMA Guamal 4.0 Cubarral San Antonio Villarica La Victoria Ortega Saldaña Cabrera

Buga META Purificación Buenaventura Restrepo Chaparral El Castillo Dagua Guacari Coyaima Prado Vijes 0 Ginebra VALLE DEL Lejanias CAUCA El Cerrito Colombia Rioblanco Alpujarra Dolores Yumbo Palmira HUILA Mesetas Ataco 00.511.522.53 3.5 Cali Natagaima La Victoria La Victoria

-77.0 -76.5 -76.0 -75.5 -75.0 -74.5 -74.0 Periodo (Segundos) Longitud Figure 3: Maximum probable acceleration for firm soil in the area of Manizales and uniform hazard spectra for different return periods

Based on the hazard information available, artificial earthquakes were generated for the time response analysis of the soil deposits and representative structures using the seismological theory of the Omega Squared model following a Green empirical function procedure. Some representative records are presented in Figure 4. Figure 5 presents the response spectra of the artificial earthquakes generated as compared to the design spectrum defined by the local code.

0.30 0.20

0.15 0.20 0.10 0.10 0.05

0.00 0.00

-0.05 Aceleración (g) -0.10 Aceleración (g) -0.10 -0.20 -0.15

-0.30 -0.20 0 1020304050 0 1020304050 Tiempo (seg) Tiempo (seg)

0.20 0.20

0.15 0.15

0.10 0.10

0.05 0.05

0.00 0.00

-0.05 -0.05 Aceleración (g) Aceleración (g) Aceleración -0.10 -0.10

-0.15 -0.15

-0.20 -0.20 0 20406080100 0 20406080100 Tiempo (seg) Tiempo (seg) Figure 4: Representative artificial records

Espectros de Aceleración 1.00 Sismo Romeral Deconvolución Sismo Romeral Sintético 0.75 Sismo Benioff Sintético Sismo Benioff Calima NSR - 98 S=1.00 0.50 NSR - 98 S=2.00 Sa (g)Sa

0.25

0.00 0.00.51.01.52.02.53.0 Periodo (seg) Figure 5: Acceleration response spectra for artificial earthquakes

DYNAMIC SOIL RESPONSE

From the stand point of dynamic response analysis, it is necessary to characterize surface deposits which due to its dynamic behavior features (wave velocity, stiffness or density) show a contrasting condition with deeper deposits. The City of Manizales is located at the western flank of the Colombian Central mountain range, wherein there are found more than 23 of the identified active volcanoes of the country, with high activity in the past. For such reason, surface deposits of the entire zone are characterized by the presence of volcanic ashes or pyroclastic fall deposits of variable thickness, ranging from 3 m to 35 m. These materials were deposited over Tertiary and Cretaceous units or over more recent deposits made up by debris flows, alluvial or colluvial deposits. Also plenty of hydraulic, mechanical, or sanitary fills are found at surface level with variety of geometry and high heterogeneous distribution. The fills were used to form flat zones in initially wavy topography with the purpose of developing housing construction activities. In order to study the dynamic properties of these soils, a comprehensive boring plan was carried out as appears illustrated in Figure 6. A total of 24 boreholes were carried out at different depths which were supplemented with more than 200 previous boreholes.

Figure 6: Distribution of boreholes

Geotechnical information as well as microtemors measurements carried out in different locations of the city allowed to develop a thickness map of superficial volcanic ash soils deposits.

Based on the information obtained, recommendations for the installation of a seismic instrumentation network were made in order to improve the proposed models. The location of the accelerographs should match the representative explored locations of the city selected from the points indicated in Figure 6. The installation of at least fifteen (15) strong motion accelerographs at surface level and at least two (2) at base rock level was recommended in order to determine the characteristics of typical earthquake records and the empiric transfer functions for each one of the points of analysis within the study area.

The non-linear dynamic properties of these characteristic soils deposits were studied using in situ and laboratory cyclic tests on unaltered samples. Cyclic triaxial tests (standard ASTM D5311-92), resonant column tests (standard ASTM D4015-92), and wave velocity tests on unaltered samples confined in triaxial chamber ("Bender Element", standard ASTM D2845-95) were made. In addition, seismic cone in situ test, cyclic in situ presionmeter test and down hole shear wave velocity tests were performed. Findings from the tests conducted show in general a high variability of ash dynamic behavior.

Typical results from in situ cyclic presiometer test as well as correlation between different tests is presented in Figure 7.

1 25 0.9 Antigua alcaldia 14.8 m El cable 5m 0.8 20 el cable 8.5m 0.7 el cable 12.8 m 0.6

15 0.5 0.4 0.3 10 0.2 0.1 5 0 0.000001 0.00001 0.0001 0.001 0.01 0.1 1

0 γ 0 2 4 6 8 10 12 14 16 u (mm) Bender Element Triaxial cíclico Columna Resonante Presiometro Figure 7: Typical results from in situ cyclic presiometer test and correlation between different tests

UNIFORM HAZARD SPECTRA

A non-linear transfer function evaluation model for the volcanic ash deposits was developed in order to obtain Fourier Amplitude Spectra (FAS) at ground surface level. Using the FAS at surface level, uniform hazard spectra were constructed at different locations according to the classic hazard theory. Simultaneously, sensibility analyses for the different variables involved were carried out as well as Montecarlo simulations of non-stationary processes. The spectra establish the basis for the city's seismic design code. Figure 8 illustrates uniform hazard spectra at different locations of the city as compared with deterministic analysis using classical non-linear equivalent one-dimensional response analysis with four different accelerographic artificial records.

Espectros de Aceleración en Superficie Espectros de Aceleración en Superficie Sondeo 1 Sondeo 5 .0 .0 Romeral Deconvolución Romeral Deconvolución

Romeral Sintético Romeral Sintético

Benioff Sintético Benioff Sintético .5 .5 Benioff Calima Benioff Calima

Espectro de Peligro Espectro de Peligro Uniforme Tret = 475 años Uniforme Tret = 475 años .0 .0

.5 .5

.0 .0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 Periodo (seg) Periodo (seg) Espectros de Aceleración en Superficie Espectros de Aceleración en Superficie Sondeo 13 Sondeo 4 .0 .0 Romeral Deconvolución Romeral Deconvolución

Romeral S intético Romeral Sintético

Benioff Sintético Benioff Sintético .5 .5 Benioff Calima Benioff Calima

Espectro de Peligro Espectro de Peligro Uniforme Tret = 475 años Uniforme Tret = 475 años .0 .0

.5 .5

.0 .0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 Periodo (seg) Periodo (seg) Figure 8: Uniform hazard and deterministic spectra for different locations TOPOGRAPHIC EFFECTS

The effects of the highly abrupt superficial topography of the city in the surface dynamic response, is studied. A two-dimensional model is resolved by using QUAD-4M computer program [6, 7]. Figure 9 shows the scheme of the model developed and results regarding maximum superficial acceleration for different locations in the modeled slope shown. Figure 10 presents acceleration response spectra and amplification factors for different locations. Amplification factors are evidenced for the period range under analysis in relation with the equivalent one-dimensional response. The analysis is made for several different points along the slope.

Aceleración Máxima del Terreno Modelo 5 - h=30 m y θ=30º 400

300

200

100

0 0 40 80 120 160 200 Longitud (m) Romeral Deconvolución Benioff Sintético Benniof.ace Benioff Calima

Figure 9: Two- dimensional model and maximum acceleration at different locations

Amplificación Modelo 30m 45 grados Espectros de Aceleración Modelo 30m 45 grados Romeral por Deconvolución Romeral por Deconvolución Punto 1 Shake 91 2.00 1.20 Punto 2 Punto 1 Punto 3 Punto 2 Punto 4 Punto 3 1.75 1.00 Punto 5 Punto 4 Punto 6 Punto 5 1.50 Punto 7 Punto 6 0.80 Punto 7 Punto 8 Punto 8 Punto 9 1.25 Punto 9 Punto 10 Punto 10 Punto 11 0.60 Punto 11 1.00 Punto 12 Punto 12

Punto 13 (g) Aceleración 0.40 Punto 13 Punto 14 Punto 14

Factor de Amplificación de Factor 0.75 Punto 15 Punto 15 Punto 16 Punto 16 0.20 0.50 Punto 17 Punto 17 Punto 18 Punto 18 Punto 19 Punto 19 0.25 Punto 20 0.00 Punto 20 0.00 0.50 1.00 1.50 2.00 2.50 3.00 Punto 21 0.00 0.50 1.00 1.50 2.00 2.50 3.00 Punto 21 Punto 22 Punto 22 Periodo (seg) Periodo (seg) Punto 23 Punto 23

Figure 10: Acceleration response spectra and amplification factors at different locations

SEISMIC MICROZONATION GIS

Considering the high variability of soils deposits, both from the geometric and physical, mechanical and dynamic properties point of view, the results of seismic microzonation are presented using a geographical information system. The system allows the management of all the city’s geographical information including general cartography, rivers, roads, places of interest, location of essential and special buildings, geology, geotechnical and in general all the related information available. The system allows visualizing the different seismic intensity parameters (acceleration, velocity or displacement) for any return period and for any structural period of analysis. Simultaneously, for each specific location the system calculates, by means of interpolation, the acceleration spectrum for different return periods and for the damping coefficient selected.

One of the main advantages of the system is that it specifies the design spectra for each location, eliminating the use of traditional seismic microzones which generally bear a considerable overdesign due to the need of specifying enveloping spectra for the specified uniform microzones. Figure 11 presents the main window of the Geographical Information System to calculate seismic hazard for the city and the typical results given by the system for a particular location.

Figure 11: System main window and typical results

SEISMIC RISK EVALUATION

A seismic risk evaluation was performed on a representative group of government’s constructions in order to establish the basis for more detailed analyses related to general seismic risk management for the city.

To carry out seismic risk assessment and to evaluate expected losses caused by earthquakes, a modified version of the RN-COL v 2.1 system developed by Evaluación de Riesgos Naturales, ERN – Colombia, was used. This system includes the most comprehensive information available on the country’s seismic conditions. The system was integrated to the geographical information system for hazard evaluation in the city, including subsoil amplification effects. Additionally, the system includes the information that characterizes building typologies more often used in the city. The geographical location of the 40 main city buildings analyzed, corresponding to the municipality’s public buildings is presented in Figure 12.

Figure 12: Geographical location of analyzed public risks

15 14 13 12 11 10 9 8 7 6 5 4 Number of Buildings 3 2 1 0 0 - 50 50 - 100 100 - 200 200 - 500 500 - 1000 1000 - 1500 1500 - 2000 2000 - 5000 Insured Value (Thousands $USD)

Figure 13: presents a frequency curve of insurable values of the different buildings analyzed.

The results of the analysis are presented in Figures 14 to 16. Figure 14 shows a histogram of expected losses for the buildings analyzed. These results indicate that expected values of the individual loss of the buildings present a great variability with values ranging from 8% to 96% with respect to insurable values. Figure 15 presents a concentration curve of the group of buildings. This curve indicates that about 5% of the most critical buildings concentrate nearly 85% of risk for the city administration. On the other hand, Figure 16 depicts probable maximum loss variation, PML, for the group of buildings as a function of the return period of analysis. It is concluded that the city’s public buildings PML varies between 31% and 41% for a range of return periods between 500 and 1500 years. These results set the basis for a risk transfer and retention strategy for the city administration.

100 90 80 70 60 50 40 30

Expected Loss(%) 20 10 0 12345678910111213141516171819202122232425262728293031323334353637383940 Building Number

Figure 14: Expected losses for the buildings analyzed

00 11100% 90% 80% 70% 60% 50% 40% Concentration Index 30% IC = 0.91 Cumulative prime prime Cumulative 20% 10% 0% 0% 20% 40% 60% 80% 100% Number of buildings

Figure 15: Concentration curve of the group of buildings

60 49.7 50 40.9 43.1 37.8 40 31.2

30 20.1

PML (%) PML 20 8.9 10 0 0 1000 2000 3000 4000 5000 Return period (years)

Figure 16: Probable Maximum Loss variation

CONCLUSIONS

Based on the previous findings, the following conclusions can be drawn:

ƒ Uniform hazard spectra resulting from an analysis using the random vibrations theory together with a non-linear transfer function evaluation model developed in order to obtain Fourier Amplitude Spectra (FAS) at ground surface level, indicates very similar results as compared with acceleration spectra resulting from the classical non-linear equivalent one-dimensional response analysis with a group of accelerographic artificial records. ƒ The characterization of soils dynamic properties for volcanic ash deposits requires different laboratory and field tests in order to adequately characterize their behavior for the entire range of deformations. ƒ The results given by different laboratory tests such as wave velocity, cyclic triaxial and resonant column present an acceptable agreement with equivalent results obtained from field tests such as down-hole, cyclic presiometer and seismic cone test. ƒ Abrupt topography can generate geometric amplification effects that can be in the order of 1.2 to 1.3 in the range of structural periods between 0.2 and 0.5 seg (rigid buildings) depending on the relative location of the building within the topographical profile. ƒ The results of hazard analyses in microzonation studies are interpreted in a more convenient way by using a geographical information system (GIS) which can evaluate the design spectra at each exact location by means of continuous interpolation, thus eliminating traditional seismic microzones which generally bear a considerable over design. ƒ The results of seismic risk evaluation for public buildings of the city of Manizales indicate that expected values of the individual loss of the buildings present a great variability with values ranging from 8% to 96% with respect to insurable values. ƒ City’s public buildings PML varies between 31% and 41% for a range of return periods between 500 and 1500 years. These results set the basis for a risk management strategy for the city.

REFERENCES

1. AIS, “Estudio General de Amenaza Sísmica de Colombia. Santa Fe de Bogotá”, Colombia. October 1996, 254pp. 2. CARDER, “Exploración geotécnica, Investigación de laboratorio y Microzonificación sísmica de las áreas urbanas y suburbanas de las ciudades de Pereira, Dosquebradas y Santa Rosa de Cabal”. Pereira, Colombia. 1999, pp30-37. 3. OMPAD, “Estudios de amenaza sísmica y evaluación sismogeotécnica preliminar para la zonificación sísmica de Manizales”. Colombia, March 1996, pp39-65. 4. OPES, “Proyecto Microzonificación Sísmica de Santa Fe de Bogotá. Santa Fe de Bogotá, Colombia”. November 1996, pp 57-63. 5. Universidad de los Andes, “Microzonificación Sísmica de la Ciudad de Manizales”, Alcaldía de Manizales , Colombia, 2002. 6. Hudson, M., Idriss, I. and Beikae, M, “QUAD-4M: A computer program evaluate the seismic response of soil structures using finite element procedures and incorporating a compliant base”. Davis: University of California. May 1994, 68pp. 7. Idriss, I., and Sun, J., “SHAKE-91: A computer program for conducting equivalent linear seismic response analyses of horizontally layered soil deposits”. Davis: University of California. November 1992, 35pp.