INTRODUCCION

El anlisis del mecanismo focal de los terremotos tiene la finalidad de determinar la seria de procesos fsicos que tienen lugar en el foco ssmico para definir la serie de esfuerzos que actan en la regin focal estudiados teniendo en cuenta su aspecto cinemtico y dinmico. Al determinar el mecanismo de esfuerzos vamos a analizar parmetros como la orientacin de la fractura, la funcin temporal de las fuentes ssmicas y las dimensiones de la misma.

Existen diversos mtodos para determinar el mecanismo focal de los terremotos y su aplicacin se ha visto muy favorecida con el avance de la informtica y por la disponibilidad hoy en da, de informacin ssmica digital de buena calidad, permitiendo realizar estudios ms detallados sobre los procesos fsicos ocurridos en el foco ssmico.

FOCO EXTENSO. MODELOS COMPLEJOS

1.1 DIMENSIONES DE LA FUENTE La experiencia sismolgica pone de relieve que la fuente tiene unas dimensiones espacialesy que el proceso de ruptura no sucede simultneamente en toda la zona fracturada. Esto esparticularmente importante en el caso de la Ingeniera Ssmica puesto que los registros de altafrecuencia, y en particular los acelerogramas, no pueden explicarse sin tener en cuenta una fallaextensa sobre la que la ruptura se propaga con velocidad variable. Es necesario, por tanto,considerar modelos de ruptura ms complejos que se aproximen mejor a la realidad. Iniciaremos el estudio considerando el efecto que en el modelo de dislocacin introduce lapropagacin de la fractura a lo largo del plano de falla. Supongamos una falla de longitud L yanchura W, sobre la que la ruptura se propaga con velocidad constante r en un solo sentido(Figura 4.1).

Las ondas P y S en un punto alejado de la falla aparecen como impulsos rectangulares de anchuratemporal tO

donde c corresponde a la velocidad de las ondas P o S que se consideren. Es fcil ver que el primertrmino representa el tiempo transcurrido desde el inicio hasta el trmino de la fractura, mientrasque el segundo indica la diferencia entre los tiempos de recorrido hasta la estacincorrespondientes al punto de inicio y al punto final de la ruptura. Por ello, si representamos eldiagrama de los tiempos de recorrido, parecer que ha habido dos fuentes independientes situadasmuy prximas y separadas por un intervalo t. Esta nueva fase se conoce con el nombre de "fase deparada". Este sencillo ejemplo pone de manifiesto cmo la consideracin de caractersticas msreales en el proceso de ruptura se traduce en una mayor complicacin de los registros. Las "fases deparada", as como la "fase de ruptura" que se produce en el caso especial de que la ruptura llegue ala superficie, son ms fcilmente observables en el campo prximo. Esta ruptura en superficiegenera un impulso de aproximadamente un segundo de duracin que puede causar grandes daos alas estructuras. Este parece ser el origen de los daos sufridos por el Hospital Olive View por causadel terremoto de S. Fernando de 1971. Lgicamente, el patrn de radiacin, y por tanto las amplitudes de los

desplazamientosesperables, tambin se ven fuertemente afectados por la propagacin de la ruptura. La Figura 4.2,tomada de Hirasawa y Stauder (1965), representa los patrones de radiacin correspondientes a lasondas P y S originados por una fractura que se propaga con vr= 0.5 y vr= 0.9 en un solosentido. Las lneas continuas indican valores positivos y las discontinuas negativos. La fallaconsiderada es de desgarre con plano de fractura vertical. La comparacin con la Figura 3.2 ponede relieve la influencia de este fenmeno, llamado directividad, en la distribucin de amplitudes. Adems, estos resultados estn obtenidos considerando una ruptura que se propaga a lo largo de lafalla en una sola direccin. Los problemas se complican mucho al considerar el

movimiento de laruptura en lugares opuestos a lo largo de la falla. La propagacin puede generar efectos defocalizacin por los que las amplitudes y las frecuencias de las ondas radiadas aumentan en elsentido del movimiento y disminuyen en el sentido contrario. Es fcil imaginar la importancia queesto puede tener en Ingeniera Ssmica. Volveremos sobre esta cuestin en el Tema 8.

GRAFICO 4.2. Patrones de radiacin para una fractura mvil y 2 velocidades de ruptura diferente(Hirasawa y Stauder, 1965).

El modelo de fuente mvil descrito supone que el cambio en los esfuerzos o en los desplazamientos sobre la falla se representa con una funcin escaln unidad de Heaviside: H(t) = 0, t <0; H(t) = 1, t >0. Esta funcin, aunque de fcil manejo matemtico, tiene poca realidad fsica. Los modelos de Haskell (1964) y Brune (1970) intentan superar esta limitacin representando la evolucin temporal de la ruptura mediante funciones diferentes.

El modelo de Haskell considera que la dislocacin no se produce instantneamente sino que la evolucin temporal de los desplazamientos viene dada por una funcin del tipo: t = G(t X/vr). (4.2)

G(t) es una "funcin rampa" que vale cero para t 0 y aumenta linealmente hasta llegar a 1 cuando t alcanza un valor , llamado tiempo de formacin (rise time). Si ahora esta funcin se aplica a una falla rectangular como la descrita anteriormente, y a una velocidad r constante, -que suele tomarse como una fraccin (1/2 - 4/5) de la correspondiente a las ondas S en la regin focal-, la ruptura ocurre progresivamente a lo largo de la falla, de manera que toda ella resulta uniformemente afectada. Realmente, el modelo propuesto por Haskell (1964, 1966) es ms complejo, y slo en algunos segmentos de la falla el desplazamiento se representa por una onda coherente. El modelo de Brune supone que los esfuerzos tangenciales se generan al mismo tiempo sobre toda la falla (con lo que no existen efectos de propagacin) con una funcin escaln aplicada dentro de la fractura, que obliga a los lados opuestos de la falla a moverse en sentidos contrarios. La funcin de esfuerzo genera una onda pura de cizalla que se propaga en la direccin "y" perpendicular a la superficie de la falla. La expresin temporal del pulso reviste la forma:

La constante temporal se asemeja al "rise time" y equivale a r/ donde "r" representa elradio de la falla, que se supone circular. El correspondiente espectro de los desplazamientos en elcampo lejano, obtenido aplicando la Transformada de Fourier a los registros temporales, tiene laforma representada en la Figura 4.3. Este espectro contiene una parte plana en la zona de bajas frecuencias que se puede demostrar es proporcional a Mo. A altas frecuencias el espectro decae con 2 y las dos zonas delespectro se cortan en la frecuencia de esquina c. Esta frecuencia se relaciona

con la dimensin de la fuente. Para el modelo de Brune y el espectro de las ondas S de velocidad , el radio de la fractura vale:

Si se considera una fractura rectangular (modelo de Haskell y Savage), la expresin pasa a ser:

La anchura W suele estimarse a partir de la profundidad del foco o del espesor de la corteza o bien equipararse a una fraccin, (2/3 1/2), de L. Las dimensiones de la fuente tambin pueden estimarse a partir de la funcin de directividad (relacin entre las amplitudes espectrales de una onda superficial que deja la fuente con un ngulo y otra que lo hace con ngulo + ) (Ben Benahem 1961; 1962), o determinando el rea de las rplicas. La expresin para M es:

Siendo la amplitud espectral bajas frecuencias y R(,) el patrn de radiacin, que cuando no se conoce el mecanismo focal se aproxima a 1/2. Una vez calculado el momento, la cada de esfuerzos puede estimarse mediante la relacin:

 representa la disminucin uniforme del esfuerzo de cizalla que ha producido el deslizamiento ssmico. Como ejemplo de aplicacin, la Figura 4.4 representa dos espectros correspondientes al registro obtenido en la componente Z de la estacin de Toledo del sismo de Lorca (Espaa) del 6 de Junio de 1977. Este sismo tuvo una magnitud mbLg = 4,2, y coincidi con el funcionamiento en la zona de una red temporal de microsismicidad, instalada para estudiar la falla de Alhama de Murcia. Esta red, constituida por 5 estaciones porttiles MEQ, permiti localizar muchas de las rplicas que siguieron al temblor principal (Herriz, 1982; Mezcua et al., 1984). La Figura 4.5 esquematiza la situacin de la falla, la distribucin de las estaciones (tringulos) y la localizacin del sismo principal (estrella). Los sismos ocurridos con anterioridad a ste se representan con crculos blancos y las rplicas (eventos que suceden despus del principal y estn relacionados tectnicamente con l), con crculos negros. La falla de Alhama aparece simplificada con un trazo grueso. A su vez, la Figura 4.6 muestra la distribucin de profundidades focales de las rplicas mejor localizadas. La representacin responde a la proyeccin sobre el perfil AA' de la Figura anterior y las distancias han sido medidas a partir de la posicin de la falla en la superficie. El buzamiento de la falla y la zona activada en el proceso ssmico aparecen claramente esbozados. La extensin del rea as obtenida (56 Km2) coincide con la estimacin deducida aplicando el modelo de Brune a los espectros de las ondas P y S de los sismos ms importantes de la serie. La proporcionalidad de M con la amplitud del espectro en la zona de bajas frecuencias permiti en la segunda mitad de los sesenta medir el momento ssmico con una precisin mucho mayor que los dems parmetros focales porque, para esas fechas, las funciones de Green para las ondas superficiales de largo periodo ya eran bien reconocidas (Aki, 1966).

Modernamente se han desarrollado otros mtodos que utilizan registros Wood Anderson (Bolt y Herriz, 1983) o datos digitales de largo periodo (Dziewonski et al., 1987). De esta manera, el momento ssmico se ha convertido en un parmetro fundamental en Sismologa. Sin embargo, su utilidad en Ingeniera Ssmica ha sido inicialmente menor porque los movimientos fuertes dependen estrechamente de los detalles del proceso de ruptura y stos no pueden ser apreciados a travs de las ventanas de largo periodo utilizadas para determinar el momento mediante el anlisis espectral. La vinculacin de este parmetro con el proceso en el foco y su capacidad para estimar el tamao del terremoto, han hecho que actualmente sea ms utilizado en Ingeniera. Los modelos de Haskell (1964, 1966) y Brune (1970) ya citados son un ejemplo de los modelos cinemticos de la fuente, que describen el desplazamiento de la falla en cada punto y cada instante sin especificar los esfuerzos que los producen. Como se ha indicado, el modelo de Haskell asume que la dislocacin se propaga con una velocidad constante produciendo un rea de ruptura rectangular con un valor de deslizamiento igual para todos los puntos. Otra simplificacin discutible de este modelo es suponer que la ruptura se inicia simultneamente en una lnea y se detiene bruscamente, sin considerar lo que sucede al comienzo y al final de la ruptura. Por el contrario, el modelo de Brune s relaciona la dislocacin con la situacin de esfuerzos en el medio, y sienta las bases para obtener M, la dimensin de la fractura, y la cada de esfuerzos a partir del espectro de amplitudes de las ondas registradas en campo lejano. Ambos modelos se complementan, siendo el de Brune ms aplicado para terremotos de magnitud inferior a 5 ya que la limitada profundidad de la esquizosfera impide la consideracin de fallas circulares suficientemente extensas como para producir sismos de magnitudes mayores. Aunque el modelo de Brune asume una hiptesis fsicamente imposible, -velocidad de ruptura infinita-, conduce a resultados que se aproximan notablemente a los obtenidos con modelos ms complicados como los dinmicos (Madariaga, 1977). Estos se basan en la teora de la mecnica de fracturas y no imponen la funcin de deslizamiento sino que la

deducen a partir del estado de esfuerzos en la falla. El terremoto es producido por la propagacin rpida de una fractura de cizalla (shear crack) a lo largo del plano de falla, que genera fuertes concentraciones deesfuerzos y altas velocidades de deslizamiento cerca del frente de ruptura que acta, por ello, comoemisor mvil de energa. La radiacin de alta frecuencia se produce cuando la velocidad depropagacin de la ruptura cambia bruscamente al encontrar zonas resistentes (fases de parada). La distribucin irregular del deslizamiento sobre el plano de falla es hoy aceptada como un hecho suficientemente probado y constituye la hiptesis de partida para los modelos estocsticosque se describen en el apartado siguiente. MODELOS ESTOCSTICOS El anlisis de los acelerogramas ha puesto de relieve complejidades en el proceso de fractura que no son observables en los registros de campo lejano, en los que las altas frecuencias han sido filtradas por el medio. Estas observaciones han sido confirmadas recientemente con registros digitales de banda muy ancha obtenidos en campo lejano. Este tipo de registros ha permitido estudiar los espectros de aceleraciones de terremotos de gran magnitud, poniendo de manifiesto la existencia de altas frecuencias. Como ya se ha comentado, su origen se atribuye al encuentro del frente de ruptura con heterogeneidades del plano de falla que frenan o aceleran la propagacin (Madariaga, 1977; 1989). Los modelos que consideran una distribucin irregular de estas heterogeneidades se llaman estocsticos. Los ms importantes son el desarrollado por Hanks y Mc Guire (1981) y los que consideran la presencia de barreras y asperezas. El primero conduce a unos espectros que reflejan muy bien algunas observaciones y que sern comentados en el siguiente apartado. Los ltimos establecen que la disminucin de velocidad o incluso la detencin del movimiento de ruptura pueden achacarse a las barreras propuestas en el modelo de Das y Aki (1977), mientras que las aceleraciones seran atribuibles a las asperezas definidas en el modelo establecido por Kanamori y Stewart (1978). Este ltimo modelo fue introducido para explicar las formas de los

registros telessmicos de las ondas P correspondientes al terremoto de Guatemala de 1976. Das y Aki suponen un estado de esfuerzos homogneo en la falla en la que existen barreras originadas por cambios en la litologa, plegamientos bruscos, etc., que interfieren con el frente de ruptura. Entonces pueden darse tres situaciones: Primero, que la barrera sea dbil y la ruptura avance a travs de ella a menor velocidad. Segundo, que la barrera sea fuerte y detenga el movimiento hasta que la acumulacin de esfuerzos permita romperla originando un terremoto importante. La fractura continuara despus su avance. Tercero, que la barrera sea suficientemente fuerte y permanezca sin romperse hasta que el terremoto finalice. A estas tres situaciones se aade, si se considera el anlisis tridimensional, la posibilidad de que la fractura rodee la heterogeneidad dejando una barrera local en la falla. Las zonas no fracturadas pueden romperse originando las rplicas del sismo. Aki (1979) distingui dos tipos de barreras: las geomtricas y las de relajacin. Las primeras se refieren a cambios en la direccin de la fractura o a discontinuidades en la misma. Un ejemplo es la fractura de Mendocino situada al Norte de la falla de San Andrs que detuvo la ruptura de 1906. La falla que produjo el terremoto de El Asnam (Argelia) de 1980, es otro ejemplo de la presencia de barreras geomtricas. En este caso, la falla, de aproximadamente 36 km de longitud, presenta dos barreras que la dividen en tres tramos de tamao similar y que rompieron en tiempos sucesivos claramente diferenciados. A su vez, las barreras de relajacin, tambin llamadas homogneas, son las que no presentan una clara discontinuidad geomtrica. El borde Sur de la falla de San Andrs es un caso representativo. Para el modelo de asperezas el punto de partida es el opuesto: la falla tiene una distribucin heterognea de esfuerzos que se concentran en las zonas resistentes o asperezas. Las partes dbiles son capaces de acumular menos esfuerzos y se rompen dando lugar a los premonitores. Mientras tanto, las asperezas continan acumulando esfuerzos que al liberarse originan el terremoto principal. Son, por ello, zonas que facilitan la propagacin y la aceleran.

El resultado final es la relajacin en toda la falla. Lay y Kanamori (1981) estudiaron la relacin entre la distribucin de asperezas y el entorno tectnico para las zonas de subduccin del Pacfico y establecieron una clasificacin de estas zonas de acuerdo con el tamao de las asperezas. Las ms grandes corresponden al rea de Chile, donde el conjunto de la zona de ruptura forma una sola aspereza. El caso opuesto lo constituyen las Islas Marianas, zona en la que no hay asperezas. Ambos modelos plantean la introduccin de dos nuevos parmetros: la densidad de distribucin de barreras o asperezas y la dimensin mnima de la ruptura, es decir el tamao mnimo posible para las asperezas o para la distancia entre las barreras. Como se ver ms adelante, Aki (1984) relaciona esta dimensin con la frecuencia mxima de corte a alta frecuencia, max, caracterstica de los espectros de aceleracin obtenidos a partir de acelerogramas y, en varios casos, la evala en 200-300 m. Para longitudes menores la falla no es capaz de sustentar velocidades suficientes para generar ondas ssmicas. Por ltimo, conviene tener presente que los procesos contemplados en los modelos de asperezas y barreras, as como la complejidad geomtrica de la falla, pueden alargar la duracin de la fractura hasta el punto de que, en el caso de los terremotos de gran magnitud, pueda alcanzar 200 s. La presencia de premonitores y rplicas vinculados a un gran terremoto mueve a pensar que la realidad fsica es una mezcla de los dos modelos. Por todo ello se han desarrollado otros ms complicados que generalizan o combinan los modelos de asperezas y barreras (Papageorgiou y Aki, 1983; Rundle y Kanamori, 1987; Andrews, 1989; Chatterjee y Knopoff, 1990), o acentan el tratamiento estocstico del problema (Gusev, 1983; Lomnitz-Adler, 1985). Estos modelos explican algunos espectros de aceleracin ms complejos que se comentarn en seguida. La existencia de "terremotos caractersticos", sismos que se repiten con propiedades semejantes en una misma zona-, parece relacionarse tambin con la distribucin, ms o menos estable, de asperezas y barreras. La teora de los fractales desarrollada por Mandelbrot para describir estructuras no aleatorias o fragmentadas aleatoriamente, ha sido aplicada tambin con frecuencia a la geometra y distribucin de las fallas (King, 1983; Okubo y Aki, 1987; Avils et al., 1987).

ESPECTROS DE DESPLAZAMIENTO Y ACELERACIN Al comentar el modelo de Brune se introdujo el espectro de desplazamiento correspondiente a las ondas en el campo lejano (Figura 4.3). Ahora cabe preguntarse hasta qu punto este espectro se modifica al considerar los nuevos modelos introducidos, y cules son los parmetros que caracterizan los espectros resultantes. El uso de los espectros de amplitud, -mdulo del espectro de Fourier- se justifica porque el anlisis de las seales se realiza ms fcilmente en el dominio de la frecuencia. En la actualidad los espectros correspondientes al campo lejano se obtienen a menudo a partir de registros de estaciones digitales de banda muy ancha que registran frecuencias entre 0,001 y varios Hercios y permiten observar el espectro de grandes terremotos. Este es el caso, por ejemplo, de las redes digitales Iris y Geoscope. Comenzaremos recordando que dado el espectro de desplazamientos del modelo de Brune,es suficiente multiplicarlo por 2 para obtener el correspondiente espectro terico de la aceleracinen campo lejano (Figura 4.7b). La experiencia ha demostrado que la Figura 4.7a representa bien laenvolvente de los espectros reales y ha podido demostrarse tericamente que el decaimiento con2, caracterstico de este espectro para las altas frecuencias, es casi general entre los modelos defractura. Asimismo, la expresin 4.5 puede aplicarse para cualquier modelo de falla (Madariaga,1976). Sin embargo, el espectro de aceleraciones en campo lejano observado en la prctica no seasemeja al de la Figura 4.7b sino al de la Figura 4.7c, en el que puede verse una zona central plana,con valor o, acotada en las bajas frecuencias por la frecuencia de esquina ya definida, y en lasaltas por una frecuencia de corte Q. Este parmetro depende estrechamente de la atenuacin delmanto superior y vale alrededor de 1 Hz para las ondas P y 0,25 Hz para las ondas S (Madariaga,1989). Por su parte, o depende de la magnitud del terremoto y es una medida de su duracin. Los espectros de aceleracin en campo prximo tienen formas similares a las ya vistas parael campo lejano pero se encuentran desplazados hacia frecuencias ms altas (Figura 4.7d). Esta esla forma propuesta por Hanks y McGuire (1982) y que se conoce como

"modelo 2". Ahora lasfrecuencias esquina estn prximas a 1 Hz y el espectro se prolonga hasta alrededor de los 10 Hz,dependiendo este lmite de la regin concreta en la que se muevan las ondas registradas a menos de100 km. La frecuencia de corte pasa ahora a llamarse frecuencia mxima, max. Estedesplazamiento hacia las altas frecuencias se explica porque la litosfera, -zona en la que se sitanmayoritariamente estas ondas-, es menos disipativa que el manto para las altas frecuencias. Lasimilitud en la forma con los espectros anteriores pone de relieve que el modelo de Brune continasiendo aplicable para campo cercano, al menos en primera aproximacin. Ello es debido a que enesta zona la teora de rayos sigue siendo vlida si la longitud de onda considerada es menor que ladistancia entre la fuente y la estacin de registro. La frecuencia mxima, max,, adquiere un valor caracterstico en cada regin ssmica (8 Hzpara Friuli, Italia, y 6 Hz para el Sur de California) que permanece constante independientementede la magnitud, el momento ssmico, la trayectoria o el tiempo considerados. Anderson y Hough(1984) vinculan max con la atenuacin a la que se supone concentrada en los kilmetros mssuperficiales. Hanks (1982) y Haucksson et al. (1987) recalcan tambin el papel decisivo de lascapas ms superficiales en el control de este parmetro. Frente a estos autores que ponen el acentoen la influencia de la trayectoria, otros, como Campillo y Bouchon (1985); Gariel y Campillo(1989) o Aki (1984; 1987), atribuyen el origen de fmax al proceso en el foco. En particular, Aki(1984) lo relaciona con la dimensin mnima de ruptura introducida en el modelo de barreras. Roca(1990), utilizando datos de la red Smart-1, encontr que fmax depende significativamente de lamagnitud y no de la distancia. Este resultado indica que para esa zona de estudio, -Taiwan-, fmaxse manifiesta claramente como un parmetro caracterstico de la fuente. Adems, la confirma-cin experimental de que fmax disminuye con la magnitud tiene una gran importancia en laestimacin del riesgo ssmico ya que las aceleraciones mximas que se predicen para grandesterremotos cercanos a partir de sismos ms pequeos sern menores. En cualquier caso, el origende este parmetro contina siendo muy discutido entre los sismlogos.

Por ltimo, conviene sealar la existencia de espectro reales con la forma esquematizada enla Figura 4.7e. En ellos aparece una zona de transicin entre el tramo 2 y la parte plana. Ahoraexiste una frecuencia nueva llamada fpatch que, en los modelos de falla heterognea, representa lapresencia de zonas coherentes cuya longitud es inversamente proporcional a dicha frecuencia(Papageorgiou y Aki, 1983). La zona que une las frecuencias fc y fpatch es proporcional a -1. Paralos terremotos pequeos ambas se identifican en la frecuencia de esquina fc . Estos espectros no sonmuy habituales y el modelo ms frecuente responde a la forma esquematizada en la Figura 4.7d.