EVALUACIÓ EVALUACIÓN N DEL RIESGO BASADA EN LOS SIG DEBIDO A LAS ACTIVIDADES DE EXTRACCIÓN DE SAL S AL EN TUZLA (BOSNIA Y HERZEGOVINA) El hundimiento del suelo provocado por la extracción de sal a partir de los yacimientos situados bajo la ciudad de Tuzla (Bosnia y Herzegovina) es uno de los mayores peligros que acta en un !rea muy urbanizada desde "#$%& cuando comenzó la explotación de depósitos de sal por per'oraciones omo se demostró en este trabajo& el hundimiento indujo varios 'actores de riesgo tales como de'ormaciones severas del suelo& surgimiento de 'racturas pro'undas y super'iciales y un aumento r!pido de la tabla 're!tica& conectado con la extracción de salmuera& que ahora a'ecta a varios distritos *os 'actores antes mencionados han sido cuanti'icados por el uso de metodolog+as geom!ticas& incluyendo encuestas de campo y an!lisis de datos geogr!'icos on el 'in de estimar las tasas histór histórica icas s de hundim hundimien iento& to& los autore autores s proces procesaro aron n la gran gran cantid cantidad ad (y nunca nunca antes antes procesada) de datos topogr!'icos recogidos durante dos per+odos, -e "#$. a "##"& y de "##/ a /%%0& con datos muy escasos 1osteriormente& las encuestas tradicionales 'ueron detenidas total y de'initivamente El an!lisis revela un hundimiento acumulado de hasta "/ m durante todo el per+odo& causando da2os a edi'icios e in'raestructuras dentro de una zona zona que que incl incluy uye e gran gran part parte e de la ciud ciudad ad hist históri órica ca&& actu actual alme ment nte e casi casi tota totalm lmen ente te destruidas *as tasas de hundimiento modernas han sido monitoreadas con 314 est!tico& mientras que la presencia de 'racturas super'iciales monitoreadas con 314 cinem!tico *os 'actores relacionados con la presencia de 'racturas pro'undas y aumento de la tabla 're! 're!ti tica ca se han han eval evalua uado do medi median ante te t5cn t5cnic icas as de an!l an!lis isis is de curv curvat atur ura a y dato datos s piezom5tricos& respectivamente 1or ltimo& los 'actores de riesgo se han combinado en un mapa de riesgos utilizando las capacidades de !lgebra de mapas 364 (4istema de 6n'ormación 3eogr!'ica) y un simple an!lisis de decisiones multicriterio (7-8) 1ara ello& se ha obtenido obtenido un mapa de vulnerabilidad vulnerabilidad basado en la in'ormación in'ormación reportada en un par de im!genes de sat5lite de alta resolución recientemente detectadas El riesgo 'inal& surgido de la combinación de 'actores de riesgo nico y mapa de vulnerabilidad& resalta escenarios cr+ticos y amenazas insospechadas que est!n siendo consideradas por los responsables locales y los plani'icadores urbanos En particular& como se destaca en el mapa de riesgos 'inal& el aumento actual de la capa 're!tica& provocado por la disminución del bombeo de salmuera& est! poniendo seriamente una amenaza para una parte de la ciudad que no es la m!s involucrada en de'ormaciones del suelo
a litera literatur tura a relati relativa va a los proble problemas mas relaciona relacionados dos con la subsid subsidenc encia ia causad causados os por la disolución natural o inducida por el hombre de las evaporitas es rica en casos de estudio muy bien documentados ooper ("##.& /%%/)& 9euter y 4toyan ("##0) han investigado los sumideros en las regiones :!rsicas de 8lemania& 3uti5rrez et al (/%%"& /%%;) y 3uti5rrez y ooper (/%%/) in'ormaron una amplia gama de casos de estudio relacionados con con la disol disoluc ució ión n natu natura rall de los los evap evapor orit itos os Tri! Tri!si sico cos s y Tercia rciari rios os en Espa Espa2a 2a *as *as contri contribuc bucion iones es adicio adicional nales es de los invest investiga igador dores es espa2ol espa2oles es se propor proporcio cionan nan en una amplia zona de la cuenca del Ebro& donde la 'ormación de sumideros se encuentra con 'recuencia
1. INTRODUCCIÓN 8utor correspondiente -i!logo arquitectónico y urban+stico& 1olit5cnico de Bari&
& ?%"/$ Bari& 6talia Tel Tel @ A0# %;% $#.00##, ax@ A0# %;% $#.00>;
-irecciones
de
correo
electrónico@
'manciniCpolibait
(
7ancini)&
'rancescostecchi/Cuniboit ( ( 4tecchi)& giovannigabbianelliCuniboit (3 3abbianelli) %%"0D?#$/ F D ver in'ormación preliminar G /%%# Elsevier B< Todos los derechos reservados -oi@ "%"%". jenggeo/%%#%.%"; (3uti5rrez y 3uti5rrez& "##;& 4oriano y 4imón& "##$& Benito et al& "##;& /%%%& 4oriano y 4imón& /%%"& /%%/& 3uti5rrezD4antolalla et al& /%%$) En los Estados nidos los riesgos debidos debidos a la :arsti :arsti'ic 'icaci ación ón natura naturall han sido sido revisa revisados dos por Iohnso Iohnson n ("##?& ("##?& /%%$) /%%$) *as ocurrencias ocurrencias de subsidencia subsidencia relacionadas relacionadas con la solución de evaporita evaporita natural o inducida inducida tambi5n se han investigado en varias ciudades de Europa En la ciudad de 4tuttgart (8lemania) se han descrito pozos de drenaje en 3arle'' et al ("##?) y Toulemont ("#;>) estudiaron el hundimiento del suelo en las a'ueras de 1ar+s en rancia 8dem!s& los e'ectos relacionados con la disolución de yeso en 9ipon (Jorth Kor:shire& 3ran Breta2a) han sido discutidos por ooper ("#;;) y ooper y Laltham ("###)& mientras que 3uerrero et al (/%%>) y 4oriano y 4imón (/%%/) reportaron el peligro por disolución de evaporita en Maragoza (Espa2a) y 3uti5rrez y ooper (/%%/) investigaron los e'ectos de disolución en la histórica ciudad de alatayud (JE) *a subsidencia debida a la disolución de yeso tambi5n ha sido discutida por 1au:Ntys et al ("###) en las ciudades de Birai y 1asvalys (*ituania septentrional) y un estudio de caso
a litera literatur tura a relati relativa va a los proble problemas mas relaciona relacionados dos con la subsid subsidenc encia ia causad causados os por la disolución natural o inducida por el hombre de las evaporitas es rica en casos de estudio muy bien documentados ooper ("##.& /%%/)& 9euter y 4toyan ("##0) han investigado los sumideros en las regiones :!rsicas de 8lemania& 3uti5rrez et al (/%%"& /%%;) y 3uti5rrez y ooper (/%%/) in'ormaron una amplia gama de casos de estudio relacionados con con la disol disoluc ució ión n natu natura rall de los los evap evapor orit itos os Tri! Tri!si sico cos s y Tercia rciari rios os en Espa Espa2a 2a *as *as contri contribuc bucion iones es adicio adicional nales es de los invest investiga igador dores es espa2ol espa2oles es se propor proporcio cionan nan en una amplia zona de la cuenca del Ebro& donde la 'ormación de sumideros se encuentra con 'recuencia
1. INTRODUCCIÓN 8utor correspondiente -i!logo arquitectónico y urban+stico& 1olit5cnico de Bari& & ?%"/$ Bari& 6talia Tel Tel @ A0# %;% $#.00##, ax@ A0# %;% $#.00>;
-irecciones
de
correo
electrónico@
'manciniCpolibait
(
7ancini)&
'rancescostecchi/Cuniboit ( ( 4tecchi)& giovannigabbianelliCuniboit (3 3abbianelli) %%"0D?#$/ F D ver in'ormación preliminar G /%%# Elsevier B< Todos los derechos reservados -oi@ "%"%". jenggeo/%%#%.%"; (3uti5rrez y 3uti5rrez& "##;& 4oriano y 4imón& "##$& Benito et al& "##;& /%%%& 4oriano y 4imón& /%%"& /%%/& 3uti5rrezD4antolalla et al& /%%$) En los Estados nidos los riesgos debidos debidos a la :arsti :arsti'ic 'icaci ación ón natura naturall han sido sido revisa revisados dos por Iohnso Iohnson n ("##?& ("##?& /%%$) /%%$) *as ocurrencias ocurrencias de subsidencia subsidencia relacionadas relacionadas con la solución de evaporita evaporita natural o inducida inducida tambi5n se han investigado en varias ciudades de Europa En la ciudad de 4tuttgart (8lemania) se han descrito pozos de drenaje en 3arle'' et al ("##?) y Toulemont ("#;>) estudiaron el hundimiento del suelo en las a'ueras de 1ar+s en rancia 8dem!s& los e'ectos relacionados con la disolución de yeso en 9ipon (Jorth Kor:shire& 3ran Breta2a) han sido discutidos por ooper ("#;;) y ooper y Laltham ("###)& mientras que 3uerrero et al (/%%>) y 4oriano y 4imón (/%%/) reportaron el peligro por disolución de evaporita en Maragoza (Espa2a) y 3uti5rrez y ooper (/%%/) investigaron los e'ectos de disolución en la histórica ciudad de alatayud (JE) *a subsidencia debida a la disolución de yeso tambi5n ha sido discutida por 1au:Ntys et al ("###) en las ciudades de Birai y 1asvalys (*ituania septentrional) y un estudio de caso
ha sido denunciado por 8ndrajchou: y Olimchou: ("##0) en el oeste de crania 4on menos menos 'recue 'recuente ntes s los trabaj trabajos os relaci relaciona onados dos con las invest investiga igacio ciones nes topogr topogr!'i !'icas cas y geod5sicas dedicadas a un riguroso monitoreo de los 'enómenos de de'ormación del suelo debido a la disolución de sal en un !rea urbanizada *a de'ormación del suelo debido a la miner+a convencional o en solución en la mina de sal de
7!s recientemente& art+culos publicados en Oim et al& 4obre la evaluación de peligros relacionados con la disolución a trav5s de t5cnicas 463& han sido publicados (/%%.)& quien in'ormó sobre un caso relacionado con el peligro cerca de una mina de carbón subterr!nea abandonada en 4amcheo: ity (orea) y por ooper (/%%;)& quien revisó los riesgo riesgos s geológic geológicos os evapori evaporiteD teD:ar :arst st en 3ran 3ran Breta2 Breta2a a usando usando el en'oqu en'oque e 463 para el an!lisis de datos recolectados 8dem!s& 7erad et al (/%%>) sugirieron el uso de an!lisis de decisiones multicriterios para el manejo de riesgos de !reas sometidas a riesgos inducidos por la miner+a
2. EL CASO DE ESTUDIO Este art+culo presenta un nuevo caso relacionado con la subsidencia inducida por la extracción en solución de un depósito de sal situada justo debajo de la ciudad de Tuzla en Bosnia Bosnia y Herzeg Herzegovi ovina na (igur (igura a ") e ilustr ilustra a las posibi posibilid lidade ades s que o'rece o'recen n las t5cnic t5cnicas as geom!ticas para la de'inición de peligros y las capacidades 463 En la evaluación de riesgos mediante un an!lisis de decisiones multicriterio (en adelante 7-8)
-esde la edad romana de la ciudad de Tuzla se denuncia como Pla ciudad de los depósitos de sal y 4pringsP El nombre Tuzla 'ue atribuido en el siglo "$ de los otomanos despu5s de su conquista de la Bosnia medieval (-onia y ine& "##>) y deriva de PtuzP& la palabra turca para la sal El m5todo antiguo de extracción de salmuera era aprovechar *os manantiales naturales de salmuera& mientras que la t5cnica moderna consiste en unos "$% pozos con bombas conectadas a las salidas de salmuera subterr!neas naturales& a una pro'undidad media de >%%D$%% m Este tipo de explotación es a veces llamado Pbombeo de salmuera salvajeP& siendo imposible saber cómo est! ocurriendo la disolución *a operación minera cambió las condiciones hidrogeológicas que permiten el 'lujo hacia abajo de agua dulce causando la disolución adicional de la sal omo lo demuestra 7ancini et al (/%%#)& existe una relación lineal directa entre la cantidad de sal explotada y la tasa media de subsidencia& destacando que la tasa de explotación desempe2a un papel importante en el proceso de subsidios Esta evidencia surge del an!lisis de datos topogr!'icos recogidos a partir de "#$. cuando se iniciaron las encuestas tradicionales trigonom5tricas y de nivelación de esp+ritus (aunque las primeras encuestas geod5sicas datan de "#">) con periodicidad anual y sobre una red densa El seguimiento de la de'ormación 'inalizó en "##/ debido a la 3uerra ivil de los pa+ses
balc!nicos y en la d5cada de "##/ a /%%0 las mediciones 'ueron discontinuas y se redujeron en la cantidad de datos recopilados 8 partir de /%%0 el monitoreo geod5sico 'ue completamente detenido siendo el es'uerzo por actividades de campo no m!s asequible omo se discute en las siguientes secciones& el procesamiento de toda la serie histórica (que se reunió durante el per+odo "#$.D/%%0) revela una acumulación acumulativa de hasta "/ m Este impresionante 'enómeno de hundimiento no produjo 'osas& pero los da2os a los edi'icios& las aguas residuales y las l+neas de tr!'ico 'ueron sumamente relevantes En el per+odo mencionado se recogieron varias in'ormaciones adicionales& como la cantidad de agua de salmuera bombeada por los pozos& los niveles piezom5tricos& los par!metros litoestratigr!'icos y geot5cnicos y se creó una base de datos (-== 7ining 6nstitute Tuzla& /%%%) dentro de una cooperación cient+'ica entre el municipio de Tuzla y la ciudad de 9otterdam (1a+ses Bajos) on el 'in de investigar la evolución espacioDtemporal del hundimiento desde una perspectiva multidisciplinaria (geod5sica& geológica y geo'+sica)& se desarrolló un proyecto conjunto de colaboración entre el 6948 (entro 6nterdisciplinario de 9icerca per 4cienze 8mbientali& niversidad de Bolonia& 6talia) 4e establecieron los municipios de 9avenna (6talia) y Tuzla (Bosnia y Herzegovina) El principal =bjeto de la asociación 'ue el seguimiento de los +ndices modernos de subsidencia por las metodolog+as de geodesia espacial y las investigaciones de los geohazards relacionados con la hundimiento, -e'ormaciones& 'ormación de 'racturas super'iciales y pro'undas& elevación y pro'undización de las capas 're!ticas como consecuencia de las 'luctuaciones en la cantidad de agua salada explotada Estos& que representan individualmente 'actores de riesgo nico& inducen un riesgo total a la ciudad en relación con la vulnerabilidad de los edi'icios e in'raestructuras *os 'actores de riesgo que a'ectan al !rea han sido investigados por separado a trav5s de encuestas de campo y tecnolog+as modernas asistidas por computadora En la actualidad& el hundimiento est! siendo monitoreado por levantamientos 314 geod5sicos& lo que conduce a la evaluación de las tasas de hundimiento moderno (hasta /% cm a) racturas super'iciales han sido mapeadas y 'racturas pro'undas identi'icadas a trav5s de un an!lisis de curvatura matem!tica implementado en la super'icie de hundimiento En ambos casos& los 'actores de riesgo se de'inen en t5rminos de densidad de 'racturas *as 'luctuaciones del nivel de la mesa 're!tica representan uno de los 'enómenos relevantes del rodamiento de peligros omo consecuencia de la reciente reducción en el bombeo de salmuera& el monitoreo piezom5trico muestra una tendencia natural de subida de la capa 're!tica que actualmente est! amenazando la inundación de la zona m!s subsidiada 1ara
procesar un mapa de riesgo total en la ciudad de Tuzla& cada 'actor de riesgo ha sido representado por un mapa raster con cinco clases de intensidad y se ha implementado un an!lisis multicriterio para de'inir el peso relativo de los 'actores individuales inalmente& se obtuvo un mapa que reanuda el peligro total mediante procedimientos de !lgebra de mapas en un entorno 463 y se combina con el mapa de vulnerabilidad realizado sobre la base de la densidad de objetos obtenida mediante im!genes satelitales de alta resolución *os mapas de riesgo obtenidos por la multiplicación de los 'actores de riesgo y los mapas de vulnerabilidad () delimitan el !rea en riesgo donde se deben evitar nuevos asentamientos
3. INFORMACIÓN GENERAL Y ENTORNO GEOLÓGICO Tuzla est! situado en la parte noreste de Bosnia y Herzegovina& se asentó justo debajo de la cordillera 7ajevica a una altura media de /0# m sobre el nivel del mar *a zona central est! situada en la llanura del r+o Iala& que comparte los asentamientos en las partes norte y sur -espu5s de 4arajevo y Banja *u:a& Tuzla es actualmente la tercera ciudad m!s grande de Bosnia y Herzegovina con sus 0# distritos locales y ".%%%% habitantes *as 'ormaciones de Tuzla 'orman parte del relleno de la uenca del 1annonian del 4ur Esta cuenca neog5nica ha sido interpretada como una cuenca retroarc que comenzó por los procesos de ri'ting en el 1aratethys al comienzo del 7ioceno Temprano 4egn lo in'ormado por HrvatoviQ (/%%.)& consiste en calizas& margas& areniscas& y sal de la edad 8quitanoDBurdigaliana (Egeriano Tard+o& Egggenburgian& =tthnangian y Oarpatian siguiendo la división estratigr!'ica 1aratethyan entral) -espu5s de la terminación de la actividad de ri'ting& la hundimiento lento continuó hasta que una nueva regresión marina produjo una discon'ormidad regional *as condiciones marinas 'ueron restauradas en el 1annonian temprano segn lo revelado por el .%% m de espesor marl y las sucesiones arcillosas de la piedra caliza depositadas en estancado& anoxic& rico en carbono& !reas basales pro'undas -urante la 5poca pannoniana& los ambientes de deposición se trans'ormaron en cuencas de agua dulce y los turbiditos se superponen mediante gradación progresiva de la pendiente delta& delta 'rontal y delta plano Estas secuencias de areniscas& siltstones y margas intercaladas alcanzan los espesores m!ximos de m!s de "$%%D/%%% m en el !rea de inter5s (Tari y 1amiQ& "##;) En Oatzer ("#%0)& 4o:liQ ("#$#& "#.>& "#;/)& 4tevanoviQ ("#??)& iciQ (/%%/)& IovanoviQ ("#;%) y
sucesión neogena se encuentran calizas& margas& areniscas y sal de la 5poca 8quitaniaD Burdigaliana (HrvatoviQ& /%%.) omo se muestra en la ig / los depósitos de sal& relacionados con la evolución de la cuenca del sur de 1anonia& se encuentran bajo la ciudad de Tuzla y sus alrededores El cuerpo de la sal (halita) tiene una 'orma elipsoidal y cubre un !rea de aproximadamente / :m/ *a 'ig / se reanuda en la parte in'erior de la pro'undidad de los cuerpos salinos& tal como se revela mediante prospecciones de extracción y s+smica a lo largo de la sección transversal 88 R *a sal *as series contienen siltstone& arenisca 'inamente granulada& piedra caliza& dolom+tica& piedra caliza y tobas on respecto a los cuerpos de sal& la halita representa la evaporita principal& pero tambi5n hay anhidrita& yeso y otras evaporitas El espesor de la serie de bandas de evaporita es de >%%D$%% m& y la parte productiva superior es de aproximadamente /%% m omo puede verse en la ig /& la mayor parte de esta 'ormación se encuentra en el lado norte del r+o Iala donde los 'enómenos de de'ormación ocurren principalmente
4.
PELIGRO
HISTÓRICO
DE
SUBSIDENCIA
POR
LAS
MEDIDAS
TOPOGRFICAS DEL PASADO Y SUS RELACIONES CON LAS ACTIVIDADES DE EXPLOTACIÓN DE LA SAL. 8ntes del inicio de la explotación& las 'ormaciones de halita se encontraban en equilibrio con salmueras saturadas y selladas con siltones hmedos *a extracción no controlada indujo la salmuera a 'luir hacia los pozos de extracción& cambiando la hidrogeolog+a local e induciendo una disolución adicional *a extracción de sal 'ue inicialmente (durante unos seis meses) realizada por inyección de agua dulce para acelerar el proceso de disolución Este m5todo se detuvo inmediatamente despu5s de la aparición de un 'enómeno de pozo de $% m de di!metro en el !rea suburbana Este 'ue el nico sumidero documentado en la zona a lo largo de los ltimos cincuenta a2os& mientras que los mecanismos de hundimiento est!n bien descritos por una secuencia de etapas de hundimiento y 'racturamiento omo resultado& en los ltimos cincuenta a2os la hundimiento del suelo asciende a > millones de m0 en volumen y el nivel piezom5trico cayó "%% m (TatiQ& "#?#)
*os primeros datos disponibles sobre los levantamientos geod5sicos episódicos se remontan a "#">& cuando se llevaron a cabo trigonometr+a y nivelación de los esp+ritus con el 'in de detectar posibles de'ormaciones en el suelo que se produjeran y que probablemente ser+an causadas por la explotación de la sal 8dem!s& se realizó una investigación m!s sistem!tica mediante encuestas tradicionales desde "#$. con una 'recuencia anual y sobre una red topogr!'ica muy densa que inicialmente estaba compuesta por "$%% puntos de re'erencia distribuidos dentro de una extensión de $ :m/ que abarcaba el sector m!s urbanizado y sus alrededores& donde la salmuera *as actividades de los pozos *as encuestas topogr!'icas anuales& o incluso m!s cortas& se mantuvieron hasta "##/& cuando una interrupción casi completa de "% a2os se produjo como consecuencia de la guerra civil en la ex Kugoslavia *as mediciones muy pobres y menos 'iables disponibles en el per+odo "##/D/%%0 permiten una investigación menos completa de los 'enómenos de subsidencia en curso y su relación con las actividades mineras on el 'in de cuanti'icar las tasas históricas de hundimiento& los datos almacenados en la base de datos (comprendida entre "#$. y /%%0) se procesaron despu5s de la eliminación de los valores at+picos (principalmente debido a la presencia de errores groseros en el original Hoja de datos) y mediciones relacionadas con puntos de re'erencia inestables o da2ados -esa'ortunadamente& no se incluyó mucha in'ormación auxiliar en las documentaciones de la base de datos y se ha realizado un procedimiento de ajuste de datos con'iable de datos de elevación en una parte del conjunto de datos solamente En la ig 0 la cantidad total de subsidencia que ocurrió durante el per+odo "#$.D/%%0 se muestra adem!s de la cantidad acumulada de explotación de salmuera relacionada con pozos que han estado 'uncionando durante ese per+odo 0 cuanti'ica el hundimiento que ocurrió en la ciudad de Tuzla con una 'orma regular claramente visible del modelo de de'ormación En el per+odo mencionado& se ha detectado un valor acumulativo de hasta "/ m (parte central de la imagen) En esa zona algunos puntos de re'erencia muestran una tasa media anual de hundimiento de /$ cm a2o que produjo varios colapsos de edi'icios y cientos de propiedades gravemente da2adas El an!lisis de las series temporales representadas en la gr!'ica de la ig 0 destaca el aumento de las tasas de subsidencia promedio hasta el a2o "#;/ cuando& debido a la reducción de las actividades de explotación& el hundimiento se desaceleró 8 partir de "#;> una estabilidad en la cantidad de explotación de sal corresponde a una estabilidad en 'enómenos de hundimiento 4in embargo& como lo revela la gr!'ica& el agua salada bombeada por los pozos y las tasas de subsidencia est! 'uertemente e inversamente correlacionada 1or lo
tanto& se establece la principal causa de la de'ormación del suelo en Tuzla na descripción m!s detallada de tales procesos y una cuanti'icación de los volmenes involucrados se describen en 7ancini et al (/%%#)
7ancini
et
al
6ngenier+a
3eolog+a
"%#
(/%%#)
"?%D";/
*a 'ig / 7apa geológico del !rea de Tuzla con l+mite superpuesto del depósito de sal (l+nea continua oscura)
!. PELIGROS ACTUALES 1ara evaluar el riesgo de hundimiento en la ciudad de Tuzla hay que abordar y cuanti'icar los principales 'actores de riesgo de actuación adem!s de la vulnerabilidad de los sectores investigados& dependiente de la presencia de edi'icios y zonas residenciales a'ectadas por el suelo 6nestabilidades *a mayor parte de la in'ormación 4e obtuvo una evaluación de riesgos a partir de la nueva geodatabase realizada mediante la reorganización de los datos proporcionados por la investigación anterior y completada con la in'ormación proporcionada por las encuestas geod5sicas y geo'+sicas realizadas desde el inicio del proyecto cooperador *os datos implementados en una estructura de geodatabase hacen que el procesamiento en el entorno 463 sea m!s '!cil y e'icaz cuando los datos se georre'erencian cuidadosamente segn la escala geogr!'ica de la investigación *os principales 'actores que contribuyen a la hundimiento del suelo o provocados por los 'enómenos de de'ormación se discuten en las siguientes secciones@ tasas de hundimiento moderno& tendencia reciente en el aumento de la tabla 're!tica& presencia de 'racturas super'iciales o pro'undas ada uno de ellos constituye un 'actor de riesgo nico que debe combinarse con la vulnerabilidad del procedimiento de
evaluación
de
riesgos
*a 'ig 0 *+nea de contorno que representa la cantidad de hundimiento acumulada (m) como consecuencia de la explotación de la sal en el per+odo "#$.D/%%0 y la cantidad total (m0) de agua de salmuera bombeada por los pozos que operan durante el mismo per+odo *a gr!'ica in'erior muestra una comparación entre las tendencias de las tasas de hundimiento promedio anual (cm a2o)& como l+nea continua& y la cantidad anual de agua de salmuera bombeada (m0)& como l+nea punteada El per+odo de re'erencia de dicho an!lisis se limita al per+odo comprendido entre "#$. y "##/
!.1. TASAS DE SUBSIDENCIA MODERNAS
-espu5s de detener las encuestas tradicionales que ocurrieron en /%%0& se necesitó una observación m!s reciente sobre la de'ormación del suelo para recopilar nueva in'ormación sobre las tasas de subsidencia modernas& relacionarlas con las actividades mineras an realizadas y proporcionar un conocimiento til para la plani'icación 'utura de las actividades 1or lo tanto& a partir de /%%>& el posicionamiento 314 relativo est!tico 'ue seleccionado por el equipo topogr!'ico italiano que operaba en Bosnia y una red 314 compuesta de . estaciones de re'erencia y .% puntos de densi'icación incluidos en el !rea de subsistencia (4tecchi et al& /%%., 7ancini et al& /%%#) y encuestados& a trav5s del posicionamiento relativo est!tico y est!tico r!pido& cuatro veces en las actividades de campo de /%%>& /%%$& /%%. y /%%? *a 'ig *a 'igura > representa el dise2o de red con estaciones de re'erencia y de densi'icación *as actividades de medición llevadas a cabo durante las actividades de campo de /%%> y /%%$ se utilizaron para enmarcar la red 314 de Tuzla dentro del 7arco de 9e'erencia Terrestre 6nternacional (6T9& /%%%) y veri'icar la estabilidad de los monumentos utilizados como re'erencia on m!s detalle se han conectado a la red de Tuzla las estaciones permanentes de 7edicina (7E-6)& 7atera (78TE)& 4o'+a (4=6) y 1enc (1EJ) mediante la combinación de observaciones 314 y un ajuste de m+nimos cuadrados ) 1rocedimiento *os resultados relacionados con la red de re'erencia (puntos "D.) con'irmaron su estabilidad dentro de una precisión de "&/ cm& principalmente debido a la longitud de l+nea de base implicada *a comparación de las soluciones proporcionadas por el procesamiento de los datos recogidos durante las campa2as de 314 de /%%>& /%%$& /%%. y /%%? pone de relieve las tasas de subsidencia modernas resumidas en la Tabla "
*a 'ig > -istribución de las estaciones 314 de re'erencia y densi'icación alrededor del !rea de subsistencia de Tuzla
!.2. ELEVACIÓN DE LA MESA DE AGUA *as depresiones topogr!'icas que ocurren en las !reas c!rsticas de las evaporitas son 'recuentemente propensas a inundaciones ya sea por la concentración de escorrent+a super'icial o por las inundaciones de las aguas subterr!neas cuando la capa 're!tica se eleva por encima de su nivel del suelo (3uti5rrez et al& /%%;) *a comprensión de la hidrogeolog+a del !rea de estudio es un aspecto crucial del an!lisis del riesgo de hundimiento& siendo el 'lujo de aguas subterr!neas el principal agente geológico responsable de la :arsti'icación de las rocas evaporitas En
En particular& la posición de la capa 're!tica (nivel piezom5trico) y los cambios relativos a trav5s del tiempo y el espacio tienen que ser cuidadosamente investigados El bombeo de salmuera silvestre& cumplido a lo largo de casi "%% a2os& in'luyó 'uertemente en el nivel del agua subterr!nea 1ara el control del nivel piezom5trico se realizó un monitoreo por la acultad de 7inas& 3eolog+a e 6ngenier+a ivil (niversidad de Tuzla) que permitió la elaboración de mapas de la capa 're!tica de los a2os "#?#& "#;$& "##%& "##$& /%%% y /%%> (=ruS y 7andiQ& /%%$) Hoy en d+a es evidente que la parte m!s derrumbada de Tuzla es sometida a inundaciones durante episodios de lluvias 'uertes con complicaciones adicionales derivadas del aumento de la tabla 're!tica que a'ecta a la zona como consecuencia de la suspensión de reducción de bombeo de salmuera na valoración cuantitativa de la subida de la capa 're!tica se derivó de la numeralización en el 463 de los mapas iso're!ticos disponibles que permiten la investigación espacial de los desplazamientos de la capa 're!tica durante el per+odo "#;$D/%%> y para identi'icar las !reas m!s amenazadas por las inundaciones iudad& los niveles piezom5tricos se han correlacionado sucesivamente con el 7odelo de Terreno -igital (-T7) disponible Esto permitió la cuanti'icación de la distancia de la capa 're!tica por debajo de la super'icie del terreno& destacando las !reas donde el agua subterr!nea podr+a alcanzar la super'icie dependiendo de las tasas reales de subida 1ara proceder a esa simulación& se ha extrapolado un -T7 a partir del mapa de l+neas de contorno del a2o "#.$ sesgado por la sustracción de la cantidad acumulada de disminución dentro del per+odo *a comparación entre los niveles piezom5tricos& realizada por el c!lculo de !lgebra de mapa& representa una tabla de agua intensa bajando de la comparación entre los cambios de elevación de las estaciones de densi'icación detectadas durante el movimiento de agua salina y salina que enmascaran la s+smica y se2ales de resistividad respectivamente (comunicación b!sica& personal) -eterminar el estilo& la orientación y la densidad de las 'racturas en los estratos doblados o abovedados en el estudio del subsuelo& es un paso cr+tico en cualquier manejo exitoso del acu+'ero reservorio
4ubDsuper'icie y la delimitación de zonas de alto riesgo en la zona de Tuzla 1resentamos la e'ectividad del an!lisis de curvatura (*isle& "##>) aplicado a super'icies acumuladas de envolvente de hundimiento& como sustituto de la topogra'+a cl!sica y las t5cnicas geo'+sicas El an!lisis de curvatura permite identi'icar !reas en una super'icie donde la de'ormación est! m!s o menos localizada y& por lo tanto& ayuda en la identi'icación de zonas potencialmente de'ectuosas y 'racturadas 4e ensayaron varias t5cnicas de an!lisis de curvatura con un grado muy alto de concordancia 1or razones de brevedad sólo se discuten los resultados proporcionados por la curvatura gaussiana (ischer y Lil:erson& /%%%& Bergbauer y 1ollard& /%%0a& 1earce et al& /%%.) El an!lisis de curvatura gaussiana se ha aplicado al campo de de'ormación derivado del an!lisis de datos históricos& discutido en la 4ección 0& despu5s de la interpolación de los datos (4tecchi et al& /%%#) 1or lo tanto& la envolvente super'icial de la hundimiento total ha sido considerada como un registro de la partición de la tensión dentro de la ciudad de Tuzla y& dada su suavidad a peque2a escala& es una super'icie ideal para estudiar la curvatura *os resultados proporcionados por el an!lisis de curvatura se muestran en el a2o "#?# a "#;$& que corresponde a un per+odo de intensas actividades de bombeo de salmuera 1or el contrario& desde el a2o "#;$ hasta el a2o /%%> los datos ponen de relieve un aumento gradual de la capa 're!tica que alcanzó #% m en una amplia porción de la zona investigada *a 'ig $ representa tal tendencia en una doble perspectiva, El primero se presenta como un mapa en el que los niveles de gris son reveladores de los valores ascendentes (valores positivos) o de disminución (valores negativos)
!.3. FRACTURADO PROFUNDO MEDIANTE ANLISIS DE CURVATURA omo demostrado en el laboratorio con varias pruebas y modelos por Hongxing y Iac:son ("##;)& la disolución de sal es capaz de inducir la 'ormación de 'allas y 'racturas debido al 'enómeno de colapso en curso 1ara investigar la presencia y localización de las estructuras de'ormadas se utilizaron varias t5cnicas& tales como per'oraciones de pozos& resistividad geoel5ctrica y estudios de re'racción s+smica 4in embargo& los resultados no contribuyeron completamente a una mejor comprensión del patrón de 'allas y 'racturas debido a la existencia de una capa de baja velocidad *a identi'icación de !reas donde la de'ormación local produjo 'allas y 'racturas super'iciales se ha logrado mediante el mapeo 314 de 'allas y 'racturas visibles en el terreno -esa'ortunadamente& el sector m!s 'racturado y de'ectuoso corresponde a !reas donde los edi'icios han sido demolidos gradualmente por la de'ormación del suelo 1or consiguiente& estas !reas han sido pro'undamente modi'icadas por las actividades humanas de restauración y las 'allas super'iciales originales est!n ahora casi completamente borradas 1or el contrario& las 'racturas y 'allas siguen siendo evidentes
en las calles y paseo del centro& donde los edi'icios antiguos todav+a est!n presentes y algo da2ado omo se muestra en la ig . el patrón de 'racturas super'iciales detectadas (representadas por l+neas discontinuas) sigue las direcciones JLD4E& con'irmando el campo de tensión resaltado del an!lisis de curvatura gaussiana
". EVALUACIÓN DEL RIESGO TOTAL 4iguiendo el en'oque propuesto por )& existen tres componentes esenciales en la determinación del riesgo@ el riesgo& la probabilidad de experimentar cualquier peligro natural o tecnológico en un lugar o región, 6denti'icando y haciendo un inventario de personas o edi'icios o 'actores económicos que se ver+an a'ectados por el riesgo si se produjera y la vulnerabilidad de los elementos en riesgo& que se relaciona con la propensión de los edi'icios y propiedades a ser a'ectados por los peligros on p5rdidas económicas) n en'oque ligeramente di'erente de la evaluación de riesgos 'ue propuesto por rozier y 3lade (/%%$) que identi'icaron el nivel de riesgo como resultado de la intersección de peligros con los valores de los elementos en riesgo mediante un an!lisis de vulnerabilidad 4in embargo& en este trabajo se estimar! el riesgo anual potencial en un !rea dada usando la 'órmula ()@ 9 V H>E>< -onde 9 representa el riesgo& como v+ctimas o p5rdidas 'inancieras de propiedades dentro de un per+odo de re'erencia, H es el peligro, E la exposición o elementos en riesgo, < la vulnerabilidad& dada por la 'racción unitaria de la exposición *a evaluación del riesgo requiere una cuanti'icación de cada uno de los 'actores de riesgo individuales para combinarlos con los componentes E y < Este tipo de c!lculo tiene que estar relacionado con una subDregión bien de'inida del !rea investigada En la siguiente sección& los 'actores de riesgo nicos ser!n re'eridos a una cuadr+cula espaciada regular de $% x $% m con el 'in de utilizar datos r!ster y herramientas de !lgebra de mapas para producir un mapa tem!tico que represente el riesgo total debido a 'enómenos de subsidencia
*a 'ig $ 6ncremento de la cantidad de agua subterr!nea desde el a2o "#;$ hasta el /%%> (m)
*a 'ig . El mapa de curvatura gaussiana (" m/) obtenido por los datos de hundimiento de "#$.D/%%0
".1. CLASES DE INTENSIDAD GEOHAZARD En la sección anterior se trató la identi'icación y cuanti'icación de 'actores relacionados con los 'enómenos de colapso que se producen en Tuzla desde hace .% a2os Estos 'actores& que representan individualmente un solo peligro& inducen un riesgo para el municipio& en relación con la vulnerabilidad de los edi'icios e in'raestructuras on el 'in de proceder con una evaluación de riesgos actual& se escogió un corto intervalo de tiempo para el an!lisis de riesgo En particular& todas las tasas involucradas en el an!lisis se re'erir!n al intervalo de tiempo /%%>D/%%$ 1or lo tanto& un caudal geológico global se calcular! a partir de las tasas anuales de subsidencia& la intensidad de 'racturas super'iciales& la intensidad de 'racturas pro'undas y las tasas de subida anual de la capa 're!tica ada uno de los 'actores de peligro enumerados anteriormente se muestrearon en la red de c5lulas de $% m y se volvieron a clasi'icar en una rampa de $ clases que in'ormaban di'erentes niveles de intensidad (de baja a alta) *a Tabla / enumera una Pescala de gravedadP para los 'actores de riesgo tomados en cuenta En la Tabla /& las tasas anuales de hundimiento (m a2o) se obtuvieron a partir de mediciones est!ticas 314 El segundo 'actor& relacionado con el aumento de la tabla 're!tica segn se detecta a partir del registro de los niveles piezom5tricos& se basa en el supuesto de que las !reas m!s amenazadas se localizan donde& dadas las tasas actuales de aumento de la tabla 're!tica& per+odo Esto se ha logrado mediante la manipulación de datos en el 463& es decir& la relación entre el nivel de la capa 're!tica mapeada en /%%> (expresada como la distancia entre el l+mite piezom5trico y la super'icie del suelo) y la tasa anual de aumento del agua *as !reas con pendiente de moderada a alta& o localizadas en las colinas& se 'iltran del an!lisis considerando su baja o nula propensión a los riesgos de inundación En la Tabla / se expresa& por lo tanto& el 'actor de riesgo de Paumento de la tabla 're!ticaP como el tiempo esperado (a2o) para el a'loramiento de la capa 're!tica *a in'ormación necesaria se obtiene del -T7 introducido en la sección >/ El riesgo relacionado con la densidad de las 'racturas super'iciales y pro'undas (m ha) se ha in'erido a partir del an!lisis de curvatura gaussiana y la encuesta 314 di'erencial respectivamente on mayor detalle& la densidad resultante de 'racturas (tanto pro'undas
como super'iciales) 'ue muestreada por una herramienta espec+'ica en el entorno 463 capaz de encontrar la densidad de caracter+sticas lineales en el vecindario de una cierta c5lula r!ster mediante la selección de un radio de bsqueda apropiado *a densidad se calcula en unidad de longitud por unidad de !rea (hect!reas) Tabla /@ lases de intensidad para riesgos individuales
".2. ANLISIS DE DECISIONES MULTICRITERIO (MDA) En la ltima d5cada se investigaron varios m5todos basados en 7-8 con di'erentes modelos de decisión para evaluar di'erentes tipos de riesgos relacionados con la presencia de geohazards (arrara et al& "##$& 8leotti y hoUdhury& "###& 8yaleU y Kamagishi& /%%$& Oolat et al /%%.& 3uzzetti et al& /%%.& 7alczeUs:i& /%%.& -ouglas& /%%?& 3onWalves et al& /%%?& Hermans et al& /%%?& hang et al& /%%;& 4terlacchini et al& /%%;& Lang et al 8l& /%%;) En 7alczeUs:i ("###) se describe un amplio examen de los posibles en'oques de la 7-8& que investigó una amplia gama de metodolog+as para la ponderación de los criterios y su implementación en el entorno de los 463 4in embargo& la literatura no es tan rica de casos que impliquen un an!lisis de riesgo en !reas densamente urbanizadas debido a peligros derivados geológicos 1or ejemplo& se podr+a solicitar un an!lisis de riesgo en una zona asentada para seleccionar nuevos emplazamientos para hospitales& escuelas& '!bricas o en la gestión de los peligros naturales que puedan a'ectar a las in'raestructuras (autopistas& puentes& redes tecnológicas) *a decisión de este tipo debe tener en cuenta todos los 'actores de peligro
que actan en esa zona y debe llevarse a cabo un an!lisis geogr!'ico basado en criterios particulares en un entorno 463 En realidad& estas Ppreguntas geogr!'icasP implican una multitud de alternativas y& a menudo& criterios de evaluación contradictorios (7alczeUs:i& "###) 1or lo tanto& una regla de decisión tiene que ser seleccionada para dictar la mejor manera de ordenar peligros en el an!lisis de riesgo total El orden de pre'erencia de un 'actor de riesgo con respecto a otros se cuanti'ica generalmente estableciendo un sistema de ponderación (o una clasi'icación) entre los criterios (peligros desde nuestro punto de vista) El peso relativo del 'actor de riesgo nico se crea a continuación por un 7-8 muy simple y comn, *a comparación pairUise (Bana et al& /%%0& /%%.) 7ediante este en'oque& el peso normalizado de los 'actores considerados se obtiene por medio de una matriz de 'actores y la comparación de parejas (ver Tabla 0) Estas t5cnicas permiten crear una clasi'icación de los 'actores de riesgo en t5rminos de importancia& proporcion!ndoles un peso relativo (o normalizado) que representa el impacto o rendimiento de un 'actor con respecto a los otros 6ncluso si la comparación requiere una elección de expertos (una opinión personal por los tomadores de decisiones) en la etapa de juicio& 7-8 sensiblemente reduce los componentes subjetivos de la alternativa espec+'ica En la Tabla 0& el valor " indica que el 'actor en la 'ila de la matriz se considera m!s importante que el 'actor en la columna de la matriz (el valor % indica lo contrario) En la columna de Ppeso normalizadoP los valores se dan por la suma de los valores de la 'ila divididos por la suma de los valores de la columna PtotalP El 'actor 'icticio se introduce como escala de re'erencia (evitando un valor cero del 'actor menos signi'icativo) y& por de'inición& pierde todas las comparaciones *os criterios de peso se basan en la hipótesis de que la estabilidad y la integridad del edi'icio est!n m!s in'luenciadas respectivamente por 'racturas que por inundaciones de agua& subsidencia y 'inalmente por 'racturas in'eridas Esta elección se debe al hecho de que las 'racturas in'eridas son resultado de un c!lculo matem!tico y& por lo tanto& incluso si su impacto en la construcción podr+a ser considerable& su 'iabilidad no es lo su'icientemente consistente como para los 'actores observados& como el hundimiento y el aumento de la tabla 're!tica 8dem!s& se podr+a establecer hipot5ticamente una conexión entre la in'iltración de agua dulce (con disolución adicional de la sal) y las 'racturas inalmente& el puntaje de riesgo total se obtiene mediante un simple procedimiento de ponderación aditiva (48L) por las capacidades de
superposición de 463 (7alczeUs:i& "###) y se asociar! con la vulnerabilidad para continuar con el mapa 'inal de riesgo Tabla 0@ omparación por pares y pesos relativos (normalizados) del criterio de evaluación
".3. VULNERABILIDAD Y EXPOSICIÓN *a evaluación de la vulnerabilidad requiere un cuidadoso censo del elemento en riesgo En el caso de Tuzla& los elementos en riesgo son principalmente edi'icios e in'raestructuras -ebido a la 'alta de in'ormación& las in'raestructuras no se incluyen en el an!lisis 1ara lo que concierne a los habitantes de Tuzla& la probabilidad de muerte y lesiones es muy baja o ausente& porque los 'enómenos de colapso del suelo siempre ocurren como movimientos graduales que nunca provocan accidentes de construcción repentinos 1or lo tanto& el an!lisis de vulnerabilidad se basa en la presencia ausencia de Edi'icios en la zona actualmente amenazada& es decir asumir < V " para bienes y propiedades y < V % para la vida de las personas En realidad& un an!lisis detallado de la vulnerabilidad de los edi'icios debe hacerse teniendo en cuenta varios par!metros tales como@ edad& materiales& cimientos y presencia de grietas El gobierno local de Tuzla no proporcionó todav+a esa in'ormación& por lo que el mapa de vulnerabilidad 'ue creado por este en'oque El mapa de exposición se ha realizado teniendo en cuenta la densidad actual de edi'icios dentro de la zona amenazada -ebe considerarse que la hundimiento inducido por la miner+a ha causado en los ltimos $% a2os e'ectos dram!ticos sobre los edi'icios de Tuzla& lo que ha llevado a la evacuación de gran parte del !rea de hundimiento *a secuencia de tales cambios urbanos ha sido rastreada a partir del an!lisis de los mapas históricos de la ciudad (de "#$. y "#.$) y por las im!genes m!s recientes de muy alta resolución detectadas desde
6:onos y la plata'orma de sat5lite Xuic:Bird y con 'echa de /%%% y /%%$& respectivamente El an!lisis del cambio urbano& llevado a cabo en el 463 despu5s de la digitalización de los mapas y las im!genes de sat5lite panchromatic despu5s de ortorrecti'icación& destacó un saldo neto de "%/? edi'icios perdidos de "#.$ a /%%$ *as estimaciones o'iciales m!s antiguas sobre las p5rdidas causadas por 'enómenos de subsidencia in'ormó una p5rdida de /0%% apartamentos Hasta el a2o "##% y otros /%%% apartamentos irremediablemente da2ados con "$ %%% personas obligadas a trasladarse a las a'ueras de la ciudad (6breljiQ et al& /%%?) *os edi'icios restantes son aquellos considerados como PexpuestosP y actualmente amenazados por los 'actores de riesgo mencionados y ser!n considerados en el an!lisis de vulnerabilidadDexposición En particular& se ha calculado un 'actor de densidad del edi'icio mediante el mismo procedimiento utilizado en el an!lisis anterior en el que se crearon clases de intensidad En la ig 4e muestran ? mapas que in'orman de la densidad del edi'icio (expresada en metros de muros por hect!rea) y las clases de intensidad derivadas El mapa de la ig ? se podr+a considerar que incluye la vulnerabilidad y la exposición omponentes donde < es igual a " y E es proporcional a la densidad de los bienes
*a 'ig ? lases de intensidad de vulnerabilidad en la ciudad de Tuzla *as clases se agrupan con respecto a la densidad del edi'icio (m ha) por los intervalos %%>D""; (clase "), "&"# D /&?% (clase /), /&?" D >&$0 (clase 0), >&$> D .&.? (clase >), .&.;D""&"" (clase $) *os intervalos de clase se establecen sobre la base del principio de ruptura natural
#. ANLISIS DE RIESGO TOTAL *a 'órmula de representa los 'actores incluidos en el an!lisis& y U es el peso asociado a los peligros *os componentes E y < est!n de'inidos de 'orma nica *a 'órmula resultante se ha implementado '!cilmente en el entorno 364 mediante una herramienta de !lgebra de mapas para la manipulación de datos El an!lisis estad+stico de los datos tambi5n debe realizarse para atribuir un nivel de precisión al componente 9 -esa'ortunadamente& los errores para 'actores individuales son di'+ciles de modelar y el an!lisis con datos inciertos todav+a representa un problema abierto 4in embargo& los datos muy densos disponibles sobre el !rea de estudio y las metodolog+as rigurosas adoptadas en los procedimientos de recolección de datos asignan un alto grado de 'iabilidad a los resultados inalmente& los mapas de riesgo en la ig ; muestran las !reas m!s amenazadas por la suma de los
'actores de riesgo de actuación& en los que los plani'icadores y los administradores locales deben prestar especial atención a las acciones de plani'icación urbana
*a 'ig ; 7apa de riesgos total con cinco clases de gravedad
$. CONCLUSIONES El trabajo exploró una serie de oportunidades o'recidas por las herramientas de la 6ngenier+a 3eom!tica y 4istemas de 6n'ormación 3eogr!'ica en el an!lisis de riesgo de un !rea amenazada por 'actores de riesgo relacionados con la intensa explotación de depósitos de sal bajo la ciudad de Tuzla *os datos recogidos por las encuestas 314& registro piezom5trico y proporcionados por el an!lisis de curvatura se utilizaron para evaluar los 'actores que a'ectan a la parte de la ciudad donde los 'enómenos de hundimiento alcanzó "/ m como se detectó a partir del conjunto de datos topogr!'icos históricos *a combinación con la vulnerabilidad 'inalmente proporcionó el mapa de riesgo total& capaz de detectar sectores m!s propensos a da2os y p5rdidas económicas El an!lisis geogr!'ico 'ue apoyado por el entorno 364& 'acilitando la manipulación del conjunto de datos geore'erenciados y la producción de mapas e in'ormes tiles en las pol+ticas locales de gobierno& plani'icación urbana y gestión de temas relacionados con los posibles riesgos El mapa de riesgos podr+a representar un punto de partida en la
de'inición de las intervenciones de mitigación omo se discute en este trabajo& el levantamiento de la capa 're!tica podr+a ser considerado como el mayor riesgo y perseverancia en el monitoreo de la 'ractura& evitando la in'iltración de agua dulce y la disolución de la nueva sal& es un imperativo *a reducción destacada de las tasas de subsidencia& detectada a partir de las observaciones 314& es coherente con el abandono gradual de las actividades mineras -esa'ortunadamente& la tendencia en la subida de la capa 're!tica producida por el cese de la explotación intensa de la sal indujo al gobierno local a imponer recientemente un nuevo bombeo de salmuera para evitar la submersión de la zona m!s derrumbada 1or lo tanto& la disolución adicional de la sal seguida por el aumento de la subsidencia es probable suceder y el c+rculo vicioso parece ser inevitable
EXPRESIONES DE GRATITUD *os autores desean agradecer el apoyo 'inanciero de 9egione Emilia 9omagna y 1rovincia di 9avenna al proyecto P9ealizzaDzione a Tuzla (Bosnia y Erzegovina)& de un programa de asistencia t5cnica para la gestión del sistema territorial de control del control de monitoreo 1roblemi di subsidenza legati allo s'ruttamento dei sottostanti depositi di salgemma (oordinador@ 1ro' 3iovanni 3abDbianelli) Pque inició el proyecto *a investigación es ahora patrocinada por la -ivisión de -iplomacia 1blica de la =T8J en el marco de Piencia por la 1azP El 7unicipio de Tuzla (6nstituto de 7iner+a de Tuzla) proporcionó la serie histórica de los datos 8gradecemos especialmente a 8mer -zindo& 9usmir 4alihovic& Toni Ji:olic& 7auro 8ltizio& 7ario *aghi y 8ndrea 7inchio por el valioso apoyo durante las encuestas 314 *a 'ig " de este trabajo se produjo utilizando el paquete de gr!'icos 37T (Lessel y 4mith& "##;)
REFERENCIAS 8leotti& 1& hoUdhury& 9& "### Evaluación de riesgos de derrumbes@ resumen& revisión y nuevas perspectivas Bullettin de 6ngenier+a 3eolog+a y 7edio 8mbiente $; (")& /"D>>
8llUardt& 1& Bellahsen& J& 1ollard& --& /%%? urvatura y 'ractura basada en datos del sistema de posicionamiento global recopilados en el anticlinal de 4heep 7ountain& Lyoming 3eos'era 0 (.)& >%; D >/"
8yaleU& *& Kamagishi& H& /%%$ *a aplicación de la regresión log+stica basada en 463 para la cartogra'+a de susceptibilidad a deslizamientos en las monta2as Oa:udaDKahi:o& Iapón entral 3eomorphology .$& "$D0"
8ndrajchou:& <& Olimchou:& 8& "##0 El cambio ambiental y el impacto humano en el :arst en el oeste de crania En@ Lilliams& 1 (Ed)& Oarst Terrains& ambio 8mbiental e 6mpacto Humano& vol /$& p!ginas ">?D".%
Bana& E& osta& 8& 82tao -a 4ilva& 1& Junes orreia& & /%%0 Evaluación 7ulticriteraria de 7edidas de ontrol de 6nundaciones@ El aso de 9ibeira do *ivramento *a Escuela de Econom+a y iencias 1ol+ticas de *ondres& *ondres
Bana& E& osta& 8& =liveira& 4& >/ D >$%
Bergbauer& 4& 1ollard& --& /%%0a ómo calcular las curvaturas normales de las super'icies geológicas muestreadas Iournal o' 4tructural 3eology /$& /?? D /;#
Bergbauer& 4& 1ollard& --& /%%0b Erratum a Pómo calcular curvaturas normales de super'icies geológicas muestreadasP Iournal o' 4tructural 3eology /$& /".?
Benito& 3& 3uti5rrez& & 15rezD3onz!lez& 8& 7achado& 7I& "##; 9espuesta 'luvial a la hundimiento cuaternario por evaporización (9+o 3allego& uenca del Ebro& Espa2a) 3eomorphology //& />0 D /.0
Benito& 3& 3uti5rrez& & 15rezD3onz!lez& 8& 7achado& 7I& /%%% aracter+sticas geomor'ológicas y sedimentológicas en los sistemas 'luviales cuaternarios a'ectados por la hundimiento inducido por solución (uenca del Ebro& Espa2a) 3eomorphology 00& /%# D //>
arrara& 8& ardinali& 7& 3uzzetti& & 9eichenbach& 1& "##$ Tecnolog+a de 463 en la cartogra'+a del riesgo de deslizamiento En@ arrara& 8& 3uzzetti& (Eds)& 4istemas de in'ormación geogr!'ica en la evaluación de peligros naturales OluUer& -ordrecht& pp "0$ D "?$
hang& JB& 1arvathinathan& 3& Breeden& IB& /%%; ombinación de 463 con la toma de decisiones multicriterios di'usos para vertederos sentados en una región urbana de r!pido crecimiento Iournal o' Environmental 7anagement& ;? (")& "0#D"$0
iciQ& 4& /%%/ omposición geológica y terreno tectónico de Bosnia y Herzegovina En@ HrvatoviQ& H (Ed)& 3u+a 3eológica a trav5s de Bosnia y Herzegovina& 4arajevoD 3eological 4urvey de la ederación de Bosnia y Herzegovina& stanic:a& 4arajevo
ooper& 8H& "#;; 4ubsidencia resultante de la disolución del yeso de 1ermian en el !rea de 9ipon, 4u relevancia para la miner+a y la extracción de agua En@ Bell& 3& et al (Ed)& 3eolog+a de ingenier+a de los movimientos subterr!neos& vol $& p!gs 0;?D0#%
ooper& 8H& "##. 3ypsum :arst de 3ran Breta2a En@ Olimchou:& 8& *oUe& -& ooper& 8& 4auro& (Eds)& 3ypsum Oarst del 7undo@ 6nternational Iournal o' 4peleology& vol /$& 0 D >& p!gs "#$ D /%/
ooper& 8H& /%%/ Halter :arst geohazards (naturales y arti'iciales) en el 9eino nido Environmental 3eology >/& $%$ D $"/ ooper& 8H& /%%; El en'oque de los 463 para los peligros geológicos de evaporiteD:arst en 3ran Breta2a 3eolog+a 8mbiental $0& #;"D##/
ooper& 8H& Laltham& 8& "### 4ubsidencia causada por disolución de yeso en 9ipon& Jorth Kor:shire 9evista Trimestral de 6ngenier+a 3eolog+a 0/& 0%$D0"% rozier& 7I& 3lade& T& /%%$ 9iesgo y riesgo de deslizamiento de tierras@ problemas& concepto y en'oque En@ 3lade& T& 8nderson& 7& rozier& 7I (Eds)& 9iesgo de deslizamiento y riesgo Liley ;/> pp
-olo 6nstituto de 7inas de Tuzla& /%%% ontrol del problema de subsidencia en la ciudad de Tuzla& Bosnia y Herzegovina 9eporte t5cnico
-onia& 9I& ine& I<8& "##> BosniaDHercegovina& un tradicional traicionado Hurst \ o Ediciones *td 0"; pp -ouglas& I& /%%? 7odelización de la vulnerabilidad '+sica en la evaluación del riesgo de peligros naturales 9iesgos Jaturales y iencias del 4istema Terrestre ?& /;0D/;? ischer& 71& Lil:erson& 74& /%%% 1redicción de la orientación de las articulaciones desde la 'orma del pliegue@ resultados del modelado pseudoDtridimensional y an!lisis de curvatura 3eology /;& "$ D ";
