Ingeniero Civil / En este trabajo de memoria de título se analiza y estudia el efecto que tiene compatibilizar registros de terremotos reales con un espectro de diseño o Target, para diferentes tipos de suelos. Entre los objetivos que se destacan del trabajo son la implementación en el software Matlab de un programa de compatibilización espectral, realizar una base de datos de registros compatibilizados y la generación de curvas de fragilidad a partir de dicha base datos. En primer lugar, se revisó antecedentes generales acerca de un método de compatibilización espectral, el cual modifica el registro que se quiere compatibilizar en el tiempo agregando funciones wavelet, para posteriormente implementar dicho método en el software Matlab. Luego se hicieron diferentes estudios con el programa ya implementado, utilizando registros de terremotos chilenos compatibilizados a un Target, los cuales fueron comparados mediante la obtención de sus parámetros sísmicos tales como la aceleración máxima del suelo (PGA), Potencial Destructivo, etc. Después de efectuar diversos estudios con el programa implementado en Matlab, se hizo una base de datos de registros de terremotos chilenos compatibilizados a diferentes espectros utilizados como Target (considerando diferentes niveles de PGA), a partir de la cual finalmente se generaron curvas de fragilidad mediante la obtención de la respuesta de desplazamiento no lineal de un modelo estructural simplificado (para distintos rangos de PGA), el cual representa las características de los principales edificios chilenos de hormigón armado que se vieron dañados sustancialmente por el terremoto del Maule el 27 Febrero de 2010 en la Región Metropolitana, de manera de determinar distintos estados de daño de la estructura. Los registros compatibilizados de la base de datos, y registros de terremotos relevantes de Chile, fueron analizados y comparados a partir de sus gráficos en escala logarítmica de Potencial Destructivo en función del PGA, de forma de visualizar la dispersión y tendencia de los datos. Mediante las comparaciones realizadas, se observó que la tendencia y propiedades de los registros compatibilizados dependen principalmente del espectro utilizado como Target y del registro inicial que se utiliza para realizar la compatibilización. Finalmente, a partir de las curvas de fragilidad, se analizó el estado de daño del modelo estructural, obteniéndose un incremento de probabilidad de daño a medida que el tipo de suelo utilizado para generar el Target exigía mayor demanda sísmica a la estructura.

Se presenta un estudio de la respuesta sísmica de la cuenca de Santiago consistente en dos partes principales. Una primera, que corresponde a un esfuerzo de compilación y análisis de los principales rasgos físicos que caracterizan la cuenca de Santiago; donde se discuten ampliamente, aspectos que van desde la génesis de la cuenca hasta la conformación actual de la red drenaje, pasando por la geología de los depósitos y rocas presentes, la forma del basamento rocoso de acuerdo a interpretación de gravimetría, la distribución de velocidad de ondas de acuerdo a recopilación de datos de sísmica de refracción, la geomorfología dada por modelos digitales de elevación, el daño ocurrido por causa del terremoto de marzo de 1985 en términos de intensidades, mediciones de vibraciones ambientales y el control sismotectónico reconocido. Se conceptualiza un modelo paramétrico de la cuenca consistente en una capa homogénea, delimitada por la topografía del terreno y por la batimetría del lecho rocoso. Una segunda parte, consiste en la modelización matemática de la respuesta sísmica de la cuenca cuando es solicitada por un frente incidente, vertical, compresional, plano; dado por una ondícula de Ricker. Para ello, se han revisado los conceptos fundamentales de la teoría de la elasticidad pertinentes a la modelización. Posteriormente, se ha implementado numéricamente el modelo matemático obtenido, por medio del método de elementos espectrales. Con este método, se obtiene la variación en el tiempo del campo de deformación, el que es utilizado para explorar la respuesta sísmica de la cuenca. Se revisaron tres tipos. La respuesta en términos de la razón espectral de Nakamura, de sitio versus roca para la componente vertical y de sitio versus roca para la componente horizontal. De estas respuestas calculadas, se ha encontrado que la componente geométrica es preservada en la razón de Nakamura y en la razón sitio versus roca para la componente horizontal. La respuesta en términos de la razón sitio/roca para la componente vertical no muestra mayor correlación con la geometría de la cuenca. Se ha revisado la correlación existente entre las respuestas calculadas y factores que se creen son importantes con respecto a la amplificación de sitio, como lo es por ejemplo la profundidad del basamento rocoso y la profundidad del nivel freático. En el caso de la profundidad del lecho, se ha encontrado correlación entre amplificación y la forma del basamento rocoso. Para el caso del nivel freático, no ha sido posible hallar una relación contundente. Finalmente, se analiza la correlación entre la amplificación estimada y aquella encontrada a partir de sismos de magnitud moderada registrados en la cuenca, no encontrándose mayor correlación, se discute una posible explicación a este hecho.

Geólogo / Isla de Pascua es una isla oceánica ubicada sobre la Placa de Nazca, en el borde occidental de la Dorsal Asísmica de Pascua. Está generada a partir de un hotspot del cual se discute su posición. Esta isla se encuentra a aproximadamente 350 km al este del Eastern Pacific Rise (EPR) y la Microplaca de Pascua. La isla posee un alineamiento de estructuras en dirección NE-SW. Sin embargo, no se sabe la causa de este alineamiento. El presente trabajo pretende encontrar, mediante el estudio de datos sismológicos y modelación de esfuerzos, las causas tectónicas para el alineamiento de estructuras presentes en la isla, y crear un modelo general de formación de ésta. Además, se pretende estudiar la presencia de movimiento magmático actual bajo la isla. Para esto, se analizaron datos sismológicos de la estación VA02 de dos años, en busca de sismos volcánicos y tectónicos. También, se analizaron datos de mecanismo de foco de sismos catalogados por el Centroid Moment Tensor Project (CMT) y United States Geological Survey (USGS) para encontrar posibles fallas y caracterizar el movimiento de la litósfera. Finalmente, se utilizó la información obtenida de dichos datos para generar un modelo tectónico y analizar los esfuerzos traspasados desde la EPR a la Placa de Nazca, utilizando el software COMSOL Multiphysics . Al observar los datos de mecanismo focal de la Microplaca de Pascua, se puede establecer una primera aproximación para caracterizar la posición, rumbo y movimiento de las dorsales oceánicas de dicha placa, los cuales difieren de Bird (2003). Analizando los mecanismos focales en la región cercana a la isla, se muestra una zona de extensión de dirección NW-SE sobre Isla de Pascua, la cual concuerda con lo descrito anteriormente por Vezzoli & Acocella (2009). Esta extensión se debería a fallas creadas durante la formación de la corteza al pasar por una zona de esfuerzos extensionales, y que permanecen activas en la actualidad. Estas fallas serían aprovechadas por esfuerzos creados sobre la corteza a raíz de un levantamiento termal producido por la Pluma de Pascua para crear un rifting activo en la zona, el cual daría lugar a la extrusión de magma para la formación de Isla de Pascua. Actualmente, no existe una pluma mantélica en la zona somera cercana a la isla (Zhao, 2007), por lo que la isla estaría pasando por un proceso de subsidencia termal y un consiguiente rifting pasivo, lo que explicaría la presencia de sismos de tipo extensional en la región en la actualidad.

El Edificio de la Cámara Chilena de la Construcción (CChC) cuenta desde 1997 a la fecha con 12 sensores uniaxiales de aceleración, 4 de ellos ubicados en su base y el resto distribuido en otros tres pisos de la estructura, que permiten registrar tanto microvibraciones ambientales como movimientos sísmicos. Desde su instalación, el sistema de sensores ha registrado la respuesta dinámica del edificio para 54 temblores de mediana intensidad. El objetivo de este trabajo de titulo, consiste en la continuación de una serie de estudios sobre la variabilidad de las propiedades modales en el Edificio de la Cámara Chilena de la Construcción, ante distintas variables ambientales. El trabajo se orienta particularmente en el estudio del comportamiento de la estructura ante eventos sísmicos no asociados a daño, por lo cual, basa su análisis en los 54 registros de movimiento antes mencionados. El proceso de identificación de los parámetros modales del edificio, a partir de los registros sísmicos, se realiza mediante un método de identificación basado en múltiples señales de entrada y salida (MIMO: Multiple Input - Multiple Output), el cual requiere tanto de los registros de la excitación externa aplicada sobre la estructura (Input) como los de respuesta dinámica en distintos puntos de ella (Output). Para el edificio de la Cámara Chilena de la Construcción, la excitación externa aplicada en eventos sísmicos se asume conocida, e igual a los registros de los sensores ubicados en su base. Dado que parte de esta memoria consiste en el desarrollo del algoritmo computacional de identificación, basado en el método MIMO, su eficacia es comprobada a partir de modelos analíticos de la estructura en estudio, en los cuales el error en los parámetros dinámicos identificados, salvo el caso de un modo torsional particular, no supera el 1%. Para considerar el comportamiento no lineal de la estructura en estudio a lo largo de los eventos sísmicos, el proceso de identificación es realizado por ventanas de avance en los registros de movimiento, de manera de observar la variación en el tiempo de los parámetros modales. De igual manera, el uso de ventanas permite obtener una mayor cantidad de datos, dando mayor robustez estadística, para el análisis de la relación entre los distintos parámetros modales y la amplitud del movimiento de la estructura Para el análisis de la relación entre los parámetros modales identificados en la estructura y la amplitud de su movimiento, se emplean distintas definiciones de esta última variable, donde la amplitud modal de movimiento, que corresponde a la componente de la respuesta asociada a los rangos de frecuencia característica de cada modo, es la que muestra mejor correlación con las variables obtenidas. En el caso de las frecuencias modales identificadas, se observa que sus valores decaen con la amplitud del movimiento de la estructura, llegando a rangos de variabilidad entre 4 y 8%, esto para valores de amplitud modal (según es definida en esta memoria) que oscilan entre 0 y 25 cm/s2 que, en cuanto a los registros estudiados, equivale a un rango de aceleración máxima basal entre 0 y 65 cm/s2. En lo que respecta a los amortiguamientos identificados, sus valores se incrementan con la amplitud de la vibración, llegando, en el mismo rango de amplitudes detallado para las frecuencias, a órdenes de variación del 60%, en el caso de los dos primeros modos traslacionales de la estructura. Para el resto de los modos, la variabilidad observada en su amortiguamiento, según sea el caso, supera el 180%.

La sismicidad presente en el Perú y la orogenia andina poseen un origen común, la subducción la Placa de Nazca debajo de la Placa Sudamericana, esta convergencia tiene velocidades del orden de 8 cm/año dando forma a la actual topografía caracterizada por la presencia de la Cordillera de los Andes. Durante el periodo íntersísmico la sismicidad posee una distribución espacial bastante compleja, nosotros observamos que la sismicidad se correlaciona con la topografía, donde la placa subducida (slab) presenta una sismicidad con fallas inversas en la zona de antearco volviéndose normales en la zona de trasarco, disminuyendo esta sismicidad entre ambas zonas, específicamente en los Altos Andes (2000m). El peso litosférico de la alta topografía modifica el estado de los esfuerzos tectónicos generando un incremento en el esfuerzo vertical litosférico, el cual produce la compensación del esfuerzo tectónico generado por la subducción, produciendo de esta manera una disminución en la sismicidad. / Seismicity in Peru and the Andean orogeny possess a common origin, the subduction of the Nazca plate beneath the South American plate. This convergence has velocities of 8 cm/yr., giving form to the current topography characterized by the presence of the Andes Cordillera. During interseismic phases the seismicity has a complex spatial distribution, with seismic activity correlating with topography, where the subducting plate (slab) presents seismicity of inverse fault type in the forearc zone, normal seismicity in the backarc zone, diminishing this seismicity between the two zones, specifically in the high Andes (2000m). The lithospheric weight of high topography modifies the state of tectonic stresses, generating an increment in the vertical lithospheric stress, which produces the compensation of tectonic stress generated by subduction, thus producing a decrease in the seismicity.

Earthquakes are natural phenomena that can cause severe damage to civilian infrastructure and prolonged disruption to society. Depending on their magnitude, epicenter location, infrastructure characteristics, and many other features, earthquakes may generate large amounts of debris and waste. The large amounts of debris generated after the disaster become one of the main problems for a population facing a health issue and rapid reconstruction of the city. A proper characterization and quantification of debris and subsequent waste management and reconstruction plan is essential for the restoration of an area affected by an earthquake. This study presents a methodological approach to characterize and quantify the debris produced as a consequence of earthquakes, as well as the flow of materials required for the reconstruction of the area affected. The proposed methodology includes an infrastructure characterization stage, a probabilistic estimation of damage by characterizing the vulnerability functions using CAPRA-GIS tool, and material flow analysis (MFA) for the characterization and quantification of debris associated with the event of an earthquake and new materials for the reconstruction stage. A case study was developed to test this methodological approach. The residential sector of Tacna, a city with high seismic risk located on the southern coast of Peru, was selected. Moreover, five different construction systems (i.e. Reinforced Masonry Bearing Walls with Concrete Diaphragms, Adobe, Wood, Concrete Shear Walls and Straw) presented in the residential sector of Tacna were characterized. Also three possible earthquake scenarios (i.e. 8.6 Mw., 7.5 Mw. and 6.2 Mw.) were analyzed, each one with three different end-of-life management situations. Simultaneously, the origin and quantities of new materials needed for the reconstruction of civil infrastructure were determined. The flow of new materials considered productivity rates in the construction and manufacturing sectors. The results show that in the presence of the greater earthquake (8.6 Mw.), adobe and straw homes suffered greater damage, with damage percentages of 63 and 48, yielding 27000 tonnes and 1390 tonnes of debris, respectively. Also, 204,000 tonnes of concrete, 7,400 tonnes of steel and 461,400 tonnes of clay brick were included as debris generated in this scenario. Furthermore, for all scenarios, the MFA estimates a regional import of materials (e.g. cement, steel, brick and wood) for the reconstruction phase. The following methodology is applicable to developed and undeveloped countries with different housing types, their respective vulnerability functions and constant earthquake recurrence. / Tesis

La presente investigación evalúa el aporte de ejercicios de preparación para desastres al nivel distrital en el Perú (simulacros) al mejoramiento de planes y políticas distritales, provinciales, regionales y nacionales para la preparación, reacción y respuesta ante los desastres. A través de un estudio de caso de un simulacro para la respuesta ante tsunamis en Cerro Azul, Cañete, de 2011 a 2012, este trabajo analiza hasta qué punto el marco para los ejercicios de respuesta frente los desastres en el Perú concuerda con un enfoque hacia las lecciones aprendidas para mejorar el nivel de preparación ante los desastres. Según los hallazgos, puede que la eficacia del marco de los ejercicios se limite por:  La rigidez de la directiva política que gobierna los ejercicios de respuesta ante desastres  La falta de claridad, y por lo tanto, la falta de comprensión, del propósito de simulacros entre los que lideran, planean y participan en los ejercicios  Una cultura organizacional a cada nivel de la gestión de desastres (GRD) que no ha integrado una cultura de aprendizaje por completo  Inconsistencias en el marco para la evaluación del simulacro y su implementación  Capacidad, recursos o compromiso insuficientes como para planear a implementar, evaluar y realizar seguimiento para los simulacros  Las políticas y los procedimientos para la institucionalización de las lecciones aprendidas a través de los simulacros que no son completos o no existen La investigación determina que un enfoque más preciso hacia la mejora continua a través de las lecciones aprendidas puede aumentar la utilidad del simulacro y por lo general, fortalecer la preparación frente los desastres al nivel nacional, regional y local. Con base en los hallazgos, la autora les brinda recomendaciones estratégicas a todos los niveles de gobierno para fortalecer la eficacia de ejercicios de preparación para los desastres. La importancia de este trabajo radica en que el Perú es propenso a desastres como terremotos y tsunamis y los sistemas de preparación ante desastres y reducción de los riesgos constituyen medidas principales de reducir los daños y pérdidas asociados a los mismos. / Tesis

Ingeniera Civil / Producto de la alta sismicidad que caracteriza al país se hace necesario desarrollar herramientas para el estudio de peligro sísmico. Éstas proveen información significativa en lo que respecta a la estimación de demanda sísmica sobre las estructuras y del movimiento del suelo. Así, en el presente trabajo de título se tiene como objetivo desarrollar curvas de atenuación para dos parámetros sísmicos de registros de movimiento fuerte en Chile, la Intensidad de Arias y la duración del movimiento fuerte, utilizando el método de ajuste Bayesiano. Se lleva a cabo una revisión bibliográfica donde se escogen algunos funcionales aplicables para la curva de atenuación de los parámetros en estudio. Se desarrolla una base de datos de sismos chilenos con la información necesaria de todos los eventos y estaciones utilizadas, sumando un total de 1117 registros correspondientes a 311 eventos con magnitudes desde 4,5 Mw hasta 8,8 Mw y distancias a la ruptura que varían de 22,1 [km] hasta 1026,01 [km]. El ajuste de las constantes se realiza utilizando un método de actualización Bayesiano, el cual permite actualizar la estimación de parámetros con incertidumbre. Para escoger el funcional que mejor se ajusta a cada parámetro se verifican los residuos, se comparan las desviaciones estándar, se comprueba que los resultados tengan sentido físico, entre otros. Es así como se obtienen los modelos predictivos de movimiento del suelo para la Intensidad de Arias, la duración significativa y la duración entre intervalos, en función de la magnitud de momento (Mw ), la fuente sismogénica, la profundidad del sismo (h), la distancia más cercana a la falla (R) y la velocidad de onda de corte promedio en los primeros 30 metros (Vs30 ). La curva obtenida para la Intensidad de Arias presenta una desviación estándar del error de 1,19, la cual se justifica por la variabilidad natural que presenta el parámetro. En cuanto a la duración significativa se utiliza un intervalo de 5%-95% de energía liberada, obteniendo una curva con una desviación estándar de 0,47. Para la duración entre intervalos se utiliza un límite de 0,05g, donde se obtiene una desviación estándar del error de 1,09, que se justifica en parte por la cantidad más limitada de registros utilizados, esto debido a que el ajuste se realiza únicamente con registros con duraciones entre intervalos mayores a cero, incorporando el resto por medio de un método de clasificación basado en Support Vector Machine. Los residuos obtenidos para las tres curvas no presentan ninguna tendencia en relación a Mw, Vs30 y R, además para la Intensidad de Arias los residuos presentan una distribución normal. Para la duración significativa los residuos muestran una distribución cercana a la normal y para la duración entre intervalos no presentan una distribución normal.

Doctor en Ciencias, Mención Geología / Se invierte la distribución de deslizamiento de cinco terremotos, 3 ocurridos en Chile (Tocopilla 2007, Mw 7.8; Michilla 2007, Mw 6.7; Maule 2010; Mw 8.8) y dos de Japón (Iwate 2008, Mw 6.8 y Tohoku 2011, Mw 9.0). Se realizan inversiones cinemáticas para ellos con la excepción de Michilla 2007 e inversiones dinámicas para Michilla 2007 e Iwate 2008. Las inversiones son hechas proponiendo a priori una distribución geométrica del área de ruptura formada por una o dos elipses y distribución gaussiana de deslizamiento. La búsqueda de la mejor solución se realiza utilizando el algoritmo de vecindad. Acelerogramas y GPS continuos (cGPS) fueron invertidos. Para Tocopilla 2007 se obtiene una distribución de slip caracterizado por 2 elipses. Para el terremoto del Maule dos elipses fueron propuestas encontrándose que el máximo deslizamiento se ubica en la zona norte de una ruptura de casi 500 km; para este terremoto además se identificaron las asperezas que controlaron el movimiento en el rango de frecuencias intermedias (0.02 Hz a 0.1 Hz) en la zona norte de la ruptura. El terremoto de Tohoku 2011 pudo ser caracterizado por la ruptura de una elipses y luego se realizó una búsqueda de la mejor solución utilizando un método de Monte Carlo fijando algunos parámetros y liberando solo 3 de ellos: velocidad de ruptura, deslizamiento máximo y el tamaño de la elipse (manteniendo fija la razón de aspecto entre sus ejes), encontrando que entre estos tres parámetros existen fuertes acoplamientos, confirmando que la solución no es única. También para este terremoto se realiza desarrolla una preliminar inversión utilizando una discretización clásica de rectángulos, encontrándose resultados similares a la inversión por elipse. Finalmente realizamos la inversión de dos terremotos intraplaca de profundidad intermedia de magnitud cercana a Mw 6.8. Para estos terremotos nosotros realizamos las primeras inversiones dinámicas liberando todos los parámetros. Aquí se confirma, realizando inversiones del tipo Monte Carlo, que las inversiones no son únicas y que la ruptura de los terremotos queda controlada por los parámetros de la ley de fricción, pudiendo tomar diferentes valores pero agrupándose en valores específicos de momento sísmico y kappa (kappa es un parámetro que relaciona la energía liberada con la energía disponible para que el terremoto se propague).

Este estudio recoge los principales resultados del trabajo sobre modelización de la anelasticidad en la Península Ibérica a partir de ondas de Coda correspondientes a una selección de 935 registros sísmicos digitales registrados en 27 estaciones equipadas con instrumental de corto periodo y pertenecientes al Instituto Geográfico Nacional (IGN). Tras la selección, se ha aplicado sistemáticamente la corrección por respuesta instrumental sobre todas las señales a analizar. Los sismos disponibles para este estudio ocurrieron en el intervalo 1997-1999 y se sitúan en las siguientes tres áreas en las que se concentra la sismicidad natural de la Península Ibérica: El Noroeste Peninsular (Galicia); El Noreste Peninsular (Cataluña, Aragón, Navarra y, parcialmente, el sur de Francia); El Sur Peninsular (sur de Valencia, Murcia, Andalucía y, parcialmente, el sureste de Portugal).Mediante métodos basados en la dispersión simple se ha calculado el factor Q de coda y su dependencia frecuencial en todas las estaciones en el rango frecuencial comprendido entre 1.5 Hz y 9.0 Hz. La variación lateral de Q de coda se resume en los mapas de iso-Q correspondientes a la banda frecuencial 6.0-8.0 Hz. Esta banda frecuencial se sitúa en la parte central del ancho de banda de las estaciones empleadas. Los mapas de iso-Q a 7 Hz han permitido detectar anomalías de Qc relacionadas con características geotectónicas de cada región. En general se obtienen valores de Qc relativamente bajos en las tres regiones estudiadas. Estos valores son ligeramente más altos en la zona noroeste y más bajos en la zona sur. También se han ajustado a la ley de dependencia frecuencial Q = Q0 ּ f ν, los valores obtenidos para el factor Qc a distintas frecuencias, obteniéndose, en general, valores de Qc0 bajos y exponentes ν altos. A raíz de los resultados obtenidos en este estudio, se advierte que la ley de dependencia frecuencial utilizada, y normalmente empleada para analizar la variación de los parámetros de atenuación inelástica, intrínseca y dispersiva, resulta no ser adecuada cuando se analiza una banda frecuencial amplia.Un segundo aspecto analizado en este trabajo es la variación temporal de los parámetros de atenuación. La ocurrencia de dos sismos de intensidad VII, uno en el Sur y otro en el Noroeste Peninsular, ha permitido analizar esta variación temporal utilizando sismos anteriores y posteriores a las crisis sísmicas de Mula en el Sur (Febrero de 1999) y de Lugo en el Noroeste (Mayo de 1997). Como consecuencia de este análisis, los resultados apuntan variaciones significativas de Qc antes y después del sismo principal.La aplicación sistemática del método de ventanas temporales múltiples que tiene en cuenta la dispersión múltiple y permite separar la atenuación intrínseca y dispersiva, ha permitido confirmar aspectos tradicionalmente aceptados y avanzar en otros en los que existe una controversia importante como son, la relación entre la Q de coda y la Q total, intrínseca y dispersiva y la dependencia frecuencial de todas ellas. Los resultados obtenidos con este análisis permiten afirmar que la atenuación intrínseca domina sobre la dispersiva para altas frecuencias en todas las regiones analizadas. En cambio, para frecuencias bajas, inferiores a 3Hz, la atenuación dispersiva domina en la región sur. / In the present study we analyze the anelastic structure of the Iberian Peninsula by using seismic coda waves of local and regional earthquakes. We use digital data collected at the seismic network of the Spanish Geographic Institute (Instituto Geográfico Nacional, IGN). Our data reflect the natural seismicity of the Iberian Peninsula and lead to consider and analyze three separated regions with different geotectonic features: The Northwestern (NO), the Northeastern (NE) and the Southern (S) Iberian Peninsula.We have analyzed 935 high quality seismograms recorded at 27 stations, equipped with digital short period and some broadband instruments, in order to observe the Qc quality factor and its frequency dependence in the frequency range between 1 and 10 Hz. The best resolution has been obtained for frequencies between 6.0 and 9.0 Hz. This frequency band, which in most cases is the center frequency of the instrument response, has been used to plot iso-Q7Hz maps for the three studied regions. The iso-Q maps at 7 Hz show a good agreement with the main geotectonic characteristics of the mapped areas. Relatively low Qc values are found in the three regions, but Qc values are higher in the NO region and lower in the S region.The frequency dependence of the quality factor has been also analyzed by using the well known law Q = Q0  f. We obtain low Qc0 values and high . It is noted that probably the law mentioned before, which is mostly used in attenuation studies, is not adequate when a relatively broad frequency band is involved.We have used two seismic events occurred in the NE region and in the S region to analyze the time dependence of the quality factor. The seismic intensity of both main earthquakes was VII in the MSK scale. We have analyzed events occurred before and after the main shock occurred during the seismic crisis of Mula in the Southeastern Iberian Peninsula (February, 1999), and the same procedure was used with events associated with the seismic crisis of Lugo, in the Northwestern Iberian Peninsula (May, 1997). Our results seems indicate temporal differences before and after the main shocks in both seismic crisis.The application of the multiple lapse time window analysis (MLTWA), which takes into account the multiple scattering and leads to the separation of the intrinsic and scattering attenuation, confirms most of the results published but also gives some indications to go further in some subjects, as the relation between coda Q (Qc), intrinsic Q (Qi), scattering Q (Qs) and total Q (Qt) and their frequency dependence. The final results indicate that the intrinsic attenuation is more important that the scattering attenuation at high frequencies in all three analyzed regions. In the southern region and for low frequencies, lower than 3 Hz, the scattering loss is the main attenuation mechanism.