COMPORTAMIENTO DEL NIVEL DEL MAR EN LA COMUNIDAD “EL FARO” EN EL MUNICIPIO DE GUAYANILLA, PUERTO RICO, DESPUES DEL SISMO DE 6.4 MW OCURRIDO EL 7 DE ENERO DE 2020

Autor: Agrimensor Héctor M. Sanabria Valentín (PR)

hsanabria@hlcmgroup.com

Co Autor: C PhD. Rigoberto A. Moreno Vázquez (MX)

rigobertoamv@gmail.com

La comunidad “El Faro” se encuentra ubicada en la latitud N17° 49ʹ 47”, longitud W66° 47ʹ 06″ en el municipio de Guayanilla, P.R., es una comunidad costera. La costa de Guayanilla tiene una longitud de 29.224km lineales equivalentes a 18.159 millas y la componen cuatro barrios costeros a saber: Boca, Indios, Rufina y Playa.

Actualmente la población del municipio de Guayanilla y de acuerdo con el censo del 2010 es de 21,581 habitantes. Tiene un área total de 64.2 millas cuadradas de esas 21.9 millas cuadradas son cuerpos de agua, dejando 42.3millas cuadradas en área de terreno. Guayanilla su nombre para el tiempo de los Tainos era “Guaynia”, también es la cuna del cacique “Agüeybaná” cuyo nombre significa “El Gran Sol”.


Figura 1 Localización de la comunidad “El Faro”

En el área que comprende los municipios de Ponce, Peñuelas, Guayanilla, Yauco, Guánica, Sabana Grande y Lajas sufrieron percance y destrucción de estructuras por los sismos (terremotos) ocurridos en esa área desde el 28 de diciembre de 2019 hasta el 7 de enero de 2020 en donde ocurrió el mayor de ellos de una magnitud 6.4Mw. En esta región o área a un siguen ocurriendo sismos hoy en día.

De acuerdo con el mapa o panel (72000 C1985J) National Flood Hazard Layer FIRMette, (Figura 2) y el Advisory Map, (Figura 3) la comunidad “El Faro” se encuentra en una zona de inundable” A” que indica que los terrenos son costeros propensos a inundación por el alza de las aguas marina.


Figura 2 National Flood Hazard Layer FIRMette

Figura 3 Advisory Map

Para determinar lo que pasó en la comunidad El Faro, en la cual las aguas marinas penetraron tierra adentro después del sismo del día 7 de enero 2020 se decidió hacer un estudio utilizando la tecnología de Percepción Remota o también conocida como Teledetección (Remote Sensing) en inglés del programa Copérnico, Sentinel-2, se utilizaron también imágenes de la Plataforma Google Earth, la imagen de desplazamientos verticales de la (publicación número 3) y la tendencia del alza de  las aguas marinas de mar Caribe (sea level trend) en el mareógrafo de La Parguera, Lajas, P.R., y por último la información obtenida en el campo utilizando GNSS en tiempo real mediante la red de base virtuales de la compañía VRS, PR, INC.

La Percepción Remota (Remote Sensing) o Teledetección puede definirse como la ciencia y arte de obtener información de un objeto analizando los datos adquiridos mediante algún dispositivo que no está en contacto físico con dicho objeto.

Copérnico es el programa de observación de la Tierra más ambicioso de la historia, diseñado para proporcionar información precisa, actualizada y de fácil acceso para mejorar la gestión del medio ambiente, comprender y mitigar los efectos del cambio climático y garantizar la seguridad a la ciudadanía.

Copérnico es el nuevo nombre para el programa de Vigilancia Mundial del Medio ambiente y la Seguridad, antes conocido como “Global Monitoring For Enviroment Security” (GMES) por sus siglas en inglés.

El Programa Copérnico en un futuro cercano (2021) espera tener 6 misiones bajo el programa Sentinel, cada misión consta de al menos dos (2) satélites por misión para asegurar el cubrimiento establecido para cada misión, generando una robusta base de datos para los servicios Copérnico. Actualmente el programa Copérnico consta de las misiones Sentinel (1,2, y 3).

La teledetección es la técnica de la observación remota esta engloba dos procesos. Por un lado, la adquisición de información de la superficie terrestre o de la atmósfera captándola radiación electromagnética emitida o reflejada por éstas. Para ello utiliza normalmente sensores montados sobre satélites que obtienen imágenes en varias zonas del espectro electromagnético, que son visible, el ultravioleta, el infrarrojo o las microondas.

En segundo lugar, la información obtenida es trasmitida a centros terrestres, donde se almacena para posteriormente interpretarla y usarla. En función de la manera en la que cada cubierta terrestre (vegetación, agua, suelo urbano, suelo desnudo, etc.) refleje la luz solar en cada canal del espectro electromagnético, podremos, por ejemplo, clasificar los usos de suelos, o saber qué tipo de vegetación predomina en una zona determinada.

En este análisis se va a usar el Sentinel 2 este provee imágenes ópticas de resolución espacial a 10 m para servicios terrestres (por ejemplo, monitorización de la vegetación, suelo, zonas acuáticas interiores, ríos, lagos y regiones costeras).

Metodología Empleada En El Primer Análisis

El satélite Sentinel 2 del programa Copernico, las imágenes son proporcionadas por los “gemelos” Sentinel 2A y 2B, con una resolución de 10 metros en el visible y disponible para descarga gratuita. Con sus 13 bandas de trabajo se puede realizar filtros a color natural y falso color o componer índices espectrales a máxima resolución. La plataforma “Copernicus Open Access Hub” es el lugar para descargar las imágenes y herramientas gratuitas como el programa SNAP el cual permite procesar las bandas. Cada una de las imágenes tomadas muestra un barrido de hasta 290 kilómetros de anchura.

Cada satélite Sentinel 2 lleva un sensor “Multi-espectral Image” (MSI) por sus siglas en ingles con 13 bandas espectrales situadas entre la región espectral de visible hasta el infrarrojo de onda corta: con cuatro bandas de 10m de resolución espacial, seis bandas de 20m y tres bandas de 60m y una anchura de barrido de 290 km como indicáramos en el párrafo anterior. El S2 incorpora tres nuevas bandas en la región del red-edge que mejora la configuración del sensor para estudio de aguas y vegetación (Delegido et al;2011). Aunque está previsto que S2 proporcione productos de alto nivel (Nivel 2b/3), actualmente los datos S2 disponibles son un producto de radiancia en el techo de la atmósfera, formado por un conjunto de grillas (granules) contiguas de 100 km2,  correspondientes al nivel 1c. Para el procesado de los datos 1c se utiliza el software Sen2Cor (Sentinel 2 Corrections), basado en un conjunto de técnicas para la corrección atmosférica y un módulo de clasificación de escena (ESA,2016b). El resultado es un producto de reflectividades en superficie correspondientes al nivel 2a.

Para hacer el análisis espectral del satélite Sentinel 2 se seleccionaron el día 11 de diciembre de 2019 en donde no había ocurrido el sismo (terremoto) y el día 20 de enero de 2020 que es después del terremoto. Ver imágenes 4 y 5).


 Imagen 4 Zona de estudio Composición Infrarrojo

El análisis espectral para este estudio está basado en el infrarrojo, la cual se necesita en combinaciones RGB.

Rojo = Banda 8 Near infrared

Verde= Banda 4 Red

Azul= Banda 3 Green

Infrarrojo

Esta combinación de bandas tiene buena sensibilidad a la vegetación verde (la cual aparecerá representada en una tonalidad roja) debido a la alta reflectividad en el infrarrojo y la baja en el visible, representa de forma clara caminos y masas de agua. Además, muestra bosques coníferos con un rojo más oscuro mientras que los bosques caducifolios lo hacen con un rojo más claro.

Las tonalidades más habituales en una composición en falso color son:

  1. Rojo, indica una vegetación sana y bien desarrollada.
  2. Rosa, áreas vegetales menos densa o con vegetación menos desarrollada.
  3. Blanco, áreas con escasa o nula vegetación.
  4. Azul oscuro o negro indica la presencia de agua.
  5. Marrón, vegetación arbustiva muy variable.
  6. Beig-dorado, zonas de transición, prados secos asociados a matorral ralo.

Imagen 5 Zona de estudio Composición Infrarrojo

En esta imagen en el recuadro inferior derecho podemos contar tres (3) pixeles de 10 m cada uno lo que indica que son treinta (30) metros de penetración de las aguas marinas tierra adentro en La comunidad El Faro.

Metodología Empleada En El Segundo Análisis

En este análisis se utilizaron dos (6, 7) imágenes de Google Earth Pro la primera de ellas de fechada 10 de octubre de 2019 antes del sismo y la segunda de fecha agosto de 2020, la latitud y longitud de la orilla en la primera fecha son las siguientes:17°59’43.16000” N, 66°47’07.0000” W.

La segunda imagen con fecha agosto 2020 la latitud y longitud de la orilla después del terremoto son las siguientes: 17°59’42.44000”N, 66°47’.06.99000” W.


Imagen 6 Zona de estudio

Imagen 7 Zona de estudio

Imagen 7-A Identificación de las diferentes ubicaciones del nivel del mar

La distancia entre las dos imágenes identificadas con un PIN en amarillo y cuyas latitudes y longitudes se indicaron el párrafo anterior es de 30.4m. para calcular la misma se utilizó la herramienta geodésica “inverse” del NGS que se encuentra en la siguiente dirección: https://www.ngs.noaa.gov/cgi-bin/Inv_Fwd/inverse3.prl.

Metodología Empleada en El Tercer Análisis

En este vamos a utilizar la imagen (10) de mi publicación anterior de título: LOS EFECTOS DEL SISMO EN LA INFRAESTRUCTURA GEODESICA DE PUERTO RICO, OCURRIDOS EN LA ZONA SUROESTE DE LA ISLA.

Esta imagen nos indica los desplazamientos verticales en el área de estudio para la publicación mencionada, el punto de interés es el “Tidal” BM 975 8053 D, PID DO1339 este por motivos del sismo tuvo una Subducción de 0.261m, ver imagen 8.


Imagen 8 Identificación de desplazamientos

A este desplazamiento vertical hay que hacerle dos correcciones una de ellas es por el alza de las aguas marinas por motivos del cambio climático, el deshiele de los casquetes polares conocido en inglés como “Sea Level Trend”, esta tendencia en el mareógrafo de Isla Magüeyes es de 1.89mm/y. Los mareógrafos existentes en Puerto Rico se ilustran en la imagen número 9 y la tendencia del alza en las aguas marinas en la imagen número 10.


Imagen 9

Imagen 10

El mareógrafo de la Isla Magüeyes es el de más antigüedad en Puerto Rico, este data de los años 1955, del último ciclo metonico o ciclo cronal lunar enero 1983 a diciembre de 2001, lo que indica que son diecinueve (19) años. De esa fecha hasta el mes de junio de 2020 han transcurrido 19.5 años, si multiplicamos la tendencia en el alza de las aguas marinas que es 1.89mm/y por los 19.5 años nos da 0.0370m. El (NMML) “Nivel Medio del Mar Local”, en este mareógrafo es de 0.101m, pero el PRVD02 en dicho mareógrafo es 0.085m. Ver imagen 11


Imagen 11

El BM mareal 975 8053D, PID DO1339 tiene un desplazamiento vertical de 0.261m, que sumados a la tendencia del alza del agua marina (0.037) más la corrección de la diferencia entre el NMML y el PRVD02 que es de 0.016m, está la sumo a la tendencia del alza de las aguas marinas y es de 0.053m., aplicando las correcciones al desplazamiento vertical nos da 0.314m, asumiendo un metro (1m) de penetración tierra adentro por cada cm. De alza de aguas marinas, nos daría una penetración de 31.4m.

Metodología Empleada en cuarto (4) y Ultimo Análisis

La manera más exacta es la medida física o sea la muestra de campo. El día 4 y 7 de septiembre del 2020, el Sr. Juan Pablo Galarza Flores natural de Guayanilla, P.R., quien tiene un bachillerato en agrimensura de la UPPR., fue a la comunidad el Faro a medir físicamente en donde estaba la orilla antes del sismo y en donde se encuentra esta después del evento sísmico del 7 de enero de 2020.

El Sr. Juan Pablo Galarza Flores tiene un amigo que ha vivido toda su vida en dicha comunidad el Sr. Anthony Espada Figueroa quien sabe en donde estaba la orilla ya que en ella había un tubo en PVC. Se consiguieron fotos que muestran en donde estaba la orilla antes del suceso. Ver imagen número 12.


Imagen 12 Identificación física de la zona de estudio

Para localizar la orilla antes del sismo y después del sismo se utilizó una unidad receptora de GNSS marca Trimble Modelo R8S y la Red de Bases Virtuales de la compañía VRS System P.R., Inc., en la orilla antes del sismo se tomaron dos (2) archivos de cinco (5) segundo ya que el lugar no era propio para observar por más tiempo. Ver imagen número 13.


Imagen 13 Identificación física de la zona de estudio

En la próxima imagen es como está la orilla de la costa actualmente, ver imagen número 14


Imagen 14 Identificación física de la zona de estudio

En la costa actual se localizó un archivo de tres (3) minutos o 180 épocas. Ver imagen número 15.


Imagen 15 Identificación física de la zona de estudio

Las coordenadas para ambos puntos están en la tabla de la página siguiente.

COORDENADAS DE LA ORILLA ANTES Y DESPUES DEL SISMO

PuntoNorteEsteDistancia m.Descripción
1217964.184162754.388 Orilla Antes
2217992.712162737.82432.99Orilla Actual

Análisis de las tres primeras metodologías utilizadas para determinar cuánto penetraron las aguas marinas en la comunidad El Faro Guayanilla, P.R., después del sismo del 7 de enero del corriente año.
METODODISTANCIA mCOMENTARIO
Imagen espectral30.0Sentinel 2, 3pixeles de 10m
Google Earth Pro30.4Posiciones Estimadas
Desplazamiento vertical31.4BM., PID DO1339

El promedio de los tres métodos es 30.6m mientras que el método directo medido en el campo y de acuerdo con lo indicado por el Sr. Anthony Espada Figueroa en donde era la orilla antes del sismo por los escombros (debris) localizada físicamente y la orilla actual la distancia es 32.99m.

Conclusión:

Los tres métodos primeros utilizados son excelentes para tener una idea de cuanto penetraron las aguas marinas en la comunidad El Faro y en cualquier otro lugar que necesitemos hacer un estudio preliminar, definitivamente la medida física es el método para utilizarse y tener una medida más exacta.

Bibliografía:

  1. Gis & Beers  www.gisandbeers.com/lo-deberias-saber-imagenes-sentinel-2/
  2. Herramientas geodésicas del NGS. https://www.ngs.noaa.gov/cgi-bin/Inv_Fwd/inverse3.prl.
  3. Center for Operational Oceanographic Products and Services. https://www.tidesandcurrents.noaa.gov/
  4. Junta de Planificación de P.R., http://jp.pr.gov/

Contribución y Agradecimientos los colegas y otras personas que me ayudaron y contribuyeron para que este estudio de cuanto fue la penetración de las aguas marinas en la comunidad El Faro, Guayanilla, P.R., fuera posible.

  1. Luis A. Maldonado Pérez, PS, MGST
  2. Juan Pablo Galarza Flores, BSA.
  3. Anthony Espada Figueroa
  4. Michael P. Michalski, Oceanography, Center for Operational Oceanographic Products and Services, Silver Spring, Maryland.

LOS EFECTOS DEL SISMO EN LA INFRAESTRUCTURA GEODESICA DE PUERTO RICO, OCURRIDOS EN LA ZONA SUROESTE DE LA ISLA.

Autor:  Héctor M. Sanabria Valentín

hsanabria@hlcmgroup.com

Co Autor: Rigoberto Alejandro Moreno Vázquez

rigobertoamv@gmail.com

En Puerto Rico, la infraestructura geodésica ésta compuesta por una Red activa de bases GNSS que operan 24/7, estas se les conoce como CORS por sus siglas en inglés, Continuously Operating Reference Station. Los CORS con los que cuenta dicha infraestructura son los siguientes: PRHL, PRAR, PRJC, MAYZ, PRLT, PRMI, PRGY, P780, MIPR, PRN4, PRLP, PRFJ y CUPR, mismos que se encuentran expuestos en la siguiente imagen 01.


Imagen 01
Ubicación de los CORS en Puerto Rico.
Fuente: https://www.ngs.noaa.gov/CORS_Map/

Los otros dos componentes, son los controles pasivos o sea las marcas que están localizadas sobre la superficie terrestre, estas pueden tener posición horizontal solamente o pueden ser referencias de Nivel (RN) o BMS por sus siglas en inglés imagen 02.

La red de controles verticales en Puerto Rico está compuesta por un total de 678 RN distribuidas de la siguiente forma 618 en la isla grande, 34 en la isla municipio de Vieques, 17 en la Isla municipio de Culebra y 9 en la Isla protegida de Mona ver la imagen 03.


Imagen 02
Red de controle Verticales PRVD02

De la totalidad de la red de controles verticales (678BMS) del datum PRVD02, se observaron para el desarrollo del modelo geoidal hibrido un total de 121 BMS de estos se rechazaron 14 utilizándose para el desarrollo del modelo geoidal hibrido 2018 (107).


Imagen 03
Red de controles verticales en Puerto Rico observados por GNSS utilizados en el nuevo modelo geoidal 2018
Fuente: https://www.ngs.noaa.gov/opusmap/

QUE SON LOS TERREMOTOS O SISMOS

Son movimientos bruscos de la corteza terrestre. Los sismos o terremotos se producen por liberación de energía, idealmente en un punto. El desplazamiento relativo de dos bloques se da en un plano de falla. El punto, hipocentro (3D) o epicentro (2D en superficie) representa la posición central de ese plano de falla. Esta energía se transmite por el interior de la tierra de forma ondulatoria ondas P o primaria y ondas S o secundaria y también las ondas de cuerpo o superficie como las ondas Love y Rayleigh.  La interacción de estas ondas con la superficie terrestre produce las ondas L, causante de los daños sísmicos. El registro gráfico de los terremotos se realiza mediante aparatos o instrumentos llamados sismógrafos. La intensidad de los terremotos se mide de acuerdo con la escala de Mercalli Modificada (MM) en el Caribe. La antigua escala de magnitud conocida como Richter, mide el terremoto por su tamaño, basado en la energía liberada. Esta escala fue ideada por el japonés Wadati en 1931; pero la escala mantiene su nombre en honor al norteamericano Charles F. Richter, quien desarrolló el concepto en 1935 en el Estado de California. Se diferencia de la escala Mercalli Modificada, en que ésta interpreta la intensidad basada en las instalaciones humanas dañadas por el terremoto mientras que la Richter mide la energía liberada por el movimiento telúrico.


Imagen 04
Mostrando la Falla el Hipocentro, Epicentro y las ondas Sísmicas

Los terremotos son causados directa o indirectamente por movimientos de las placas tectónicas o por fallamientos tectónicos. Las placas tectónicas se mueven por corrientes de convección que son generadas en el manto de la Tierra. El magma contenido en el manto está en estado fundido o semifundido y se le considera líquido. En superficie su enfriamiento forma las placas tectónicas que si tienen estado sólido.  De una manera similar al famoso experimento de la “Gota de Brea” que se ha estado llevando a cabo en la Universidad de Queensland desde 1972. En el manto de la Tierra, las corrientes de convección son generadas por cambios en la temperatura y densidad. El material más caliente en la parte inferior del manto asciende y el más frio desciende generando cambios constantes por la circulación del material.

QUÉ SON LAS PLACAS TECTONICAS

Las placas tectónicas son enormes pedazos de corteza en la parte superior rígida del manto, que encajan en sus bordes y cubren la superficie de la Tierra. “Ver imagen 05”, existen sobre 12 placas grandes y varias pequeñas. Estas placas se mueven muy lentamente solo varios centímetros al año similar a la velocidad con la que crecen las uñas. Las placas tectónicas se mueven y le dan forma a la superficie de la Tierra. Las placas tectónicas se mueven en diferentes direcciones y a diferentes velocidades, una relativa a la otra en interacción entre sus bordes. Ese movimiento de las placas es una de las causas que hacen que cambien los valores de coordenadas en los marcos de referencia geodésicos.


Imagen 05
Placas Tectónicas

La Isla de Puerto Rico está localizada entre las placas de Norte América y la del Caribe. Hay evidencia de subducción oblicua y desplazamiento lateral entre las dos placas. La actividad sísmica se concentra en ocho (8) zonas ver imagen 06.

  1. La trinchera de Puerto Rico
  2. Las fallas de pendiente Norte y Sur de Puerto Rico
  3. Al noreste en la Zona del Sombrero
  4. Al oeste en el Cañón de la Mona
  5. Pasaje de Mona
  6. Al este, en las depresiones de Islas Vírgenes y Anegada
  7. Depresión de muertos al sur
  8. En el suroeste de Puerto Rico

Imagen 06
Actividad Sísmica

El área sur y suroeste de la Isla fue afectada el día 7 de enero del corriente año, por un terremoto de magnitud 6.4 Mw, que trajo como consecuencia el desplazamiento de las siguientes estaciones de monitoreo continuo de GNSS (CORS) por sus siglas en inglés: PRGY, PRMI, P780 y PRJC, las cuales ya la agencia federal “National Geodetic Survey”  (NGS) por sus siglas en inglés determinó sus nuevas posiciones y las mismas están publicadas en el página web de dicha agencia.

No solamente las estaciones de monitoreo continuo se afectaron, sino también los controles horizontales de ORDEN “B” que se establecieron en 1995 para el proyecto del Centro de Recaudaciones de Ingresos Municipales (CRIM), que están en el área impactada por el movimiento telúrico estas son (AB9841 Lajas2 y AB9842 Magas).  Las referencias de nivel de primer orden segunda clase (BMS) por sus siglas en inglés, que se establecieron entre los años de 2007 al 2011 para el proyecto del Datum vertical PRVD02 y que se encuentran en el área donde el sismo ocasiono daños a las obras civiles como a los CORS en el área de estudio. El estudio realizado comprende los municipios de Ponce, Peñuelas, Guayanilla, Yauco, Sabana Grande y La Parguera, Lajas, P.R., ver imagen 07.


Imagen 07
Municipios de la zona de estudio

Del área de estudio se seleccionaron once (11) referencias de control las cuales se ilustran en la imagen 07 con sus nombres. La zona de amortiguamiento en amarillo (Buffer Zone) contiene 95 BMS de los cuales dos son controles horizontales del proyecto del CRIM que se nivelaron en proyecto del PRVD02 mencionados en el párrafo anterior. De las once (11) referencia excluyendo las del orden B hay 7 que tenían posición horizontal publicada en https://www.ngs.noaa.gov/opusmap/. Como se ilustra en la imagen 03.

El método de observación fue el estático, las sesiones de observación tienen más de cuatro (4) horas estas se efectuaron con receptores de doble frecuencia  marca Trimble Modelos R6, R10, R4-3 y  exceptoR10-2 todos con trípode fijo con alturas variadas entre ( 1.80m, 2.000m y 2.050m ) exceptuando la observación de la RN DO1339 (975 8053D) que hubo que añadirle una extensión para evitar la posible multi trayectoria por una verja de alambre eslabonado y con una serpentina arriba, la altura de la antena en esta fue de 2.67m . La RN DO1332 (B1014) fue observada con una unidad Trimble modelo R10-2, la altura de la antena (HI) fue de 2.05m. Las imágenes que siguen a continuación (08 y 09) son con el propósito de ilustrar los equipos en las RN.      PID DO 1335 (975- 8053A) observada por el agrimensor Ángel Noel Colón Guzmán, la próxima imagen es en la RN., PID DO1273 (D1011) observada por el agrimensor Javier Osvaldo Escobales Medina.

Imagen 08
DO1335 (975-8053A)
Agrimensor Ángel Noel Colón Guzmán

Imagen 09
DO1275 (D1011)
Agrimensor Osvaldo Javier Escobales Medina, PS.
R. López De Azua & Asociados., PSC.

Todos los datos obtenidos en el campo se procesaron con el programa en línea OPUS versión   https://www.ngs.noaa.gov/OPUS/ , excepto los datos de la estación  Lajas2 (AB9841) el cual se procesó con el programa en línea Trimble RTX-PPP, de los resultados obtenidos, se cotejaron los indicadores estadísticos de la calidad de las observaciones en cuanto a: OBS. usada, ambigüedades resueltas, RMS y los Pick to Pick, encontrando que todas cumplían con estos criterios.

Se compararon las alturas ortométricas determinadas por nivelación geométrica diferencial del Datum vertical de Puerto Rico, (PRVD02) con las alturas ortométricas derivadas por GNSS utilizando el modelo geoidal 2018, para determinar los desplazamientos verticales (subsidencia, alzamiento y licuefacción) producidos por el terremoto ver mapa en la imagen 10.


Imagen 10
Mapa de representación de desplazamiento verticales (subsidencia).

Del control de Orden “B” AB9842 (Magas) se identificó subsidencia (Subsidence) de 0.148 cm y del control de orden “B” AB9841 (Lajas 2) se elevó (Uplift) 0.057cm y el DO1298 (C1012) que tuvo un (Uplift) de 0.032cm, los demás controles verticales demostraron subsidencia (Subsidence) entre 0.002 cm DO 1275 (D1011) hasta 26.3 cm DO1335 (975-8053A).

De los resultados obtenidos también se analizaron los desplazamientos horizontales y verticales observados, se analizaron aquellos que tenían posición horizontal previamente determinados y cuyos valores están publicados en la siguiente dirección https://www.ngs.noaa.gov/opusmap/ . Ver imagen 11


Imagen 11
Desplazamientos Horizontales

Los desplazamientos horizontales de los siguientes controles de ORDEN, B “Magas” se desplazó su posición original en un Azimuth de 231° 50ʹ 34″ y una distancia de 12.5 cm, “Lajas 2” se desplazó en un Azimuth de 346° 22ʹ 23″ y una distancia de 10.2 cm. En la imagen número 11 se ilustra el componente de la dirección y la magnitud (distancia) de todos los controles observados.

La placa del Caribe se mueve de manera natural sin actividad sísmica en dirección noreste, los efectos del sismo indican que los controles observados se desplazaron en sentido suroeste debido a que el terremoto ocurrido fue en la microplaca que va desde el oeste de la Isla de Puerto Rico en dirección sureste como ilustra la imagen 04, excepto D1011, C1013, M1913, R1013 y Lajas2 que se desplazaron en dirección noroeste.

Para validar el trabajo realizado en campo y entregarle al lector un dato adicional de la tecnología con que se cuenta hoy día se realizó un mapa de interferometría utilizando la información del satélite Sentinel-1 de la European Space Administration (ESA) dentro del programa Copérnico destinado a la monitorización terrestre, atmosfera y los océanos.

Su instrumento principal es un radar de apertura sintética (SAR) en la banda C, el cual es utilizado con frecuencia para la vigilancia del tráfico marítimo, el hielo marino, los derrames de petróleo, los terremotos, los desplazamientos de tierra y para la generación de cartografía.

El gran atractivo de esta tecnología de satélites SAR es que puede proveer imágenes durante el día y la noche y en cualquier condición meteorológica ya que es un sensor activo. Aunque se puede extraer información de las imágenes de SAR capturadas, por lo general estas imágenes se procesan posteriormente para su análisis. Para este estudio de interferometría se obtuvieron imagenes de la zona afectada del 2 de enero “antes del sismo” y del 26 de enero del presente año “después del sismo”.

Metodología Empleada:

Se busca corroborar la información obtenida en el campo mediante el uso de GPS/GNSS sobre referencias de nivel y que a su vez tuviesen posición horizontal para poder determinar desplazamientos horizontales y verticales, para esto se utilizó la técnica de interferometría diferencial (DinSar) para mapear la deformación sufrida en el área de estudio.

Se obtuvieron dos imágenes la primera (1) denominada (Master) del 2 de enero de 2020 antes del sismo y la segunda (2) denominada (Slave) del 26 de enero de 2020.

Se utilizó el programa SNAP proviene del acrónimo (Sentinel Application Platform) y responde a un programa gratuito ofrecido por la Agencia Espacial Europea para procesar y analizar las imágenes satelitales provenientes del satélite Sentinel-1, también se utilizó el programa QGIS versión 3.2. para la representación de los resultados Ver imagen 12.


Imagen 12
Resultado grafico de interferometría

El post proceso de las dos imágenes nos muestran unos contornos de la banda C utilizada para este estudio, que para el caso de Sentinel-1 trabaja en un rango de longitud de onda (λ) de 5.6cm, cada ciclo de fase equivale a 2.8cm de deformación.

En la imagen se aprecian ocho (8) simetrías del mismo color desde el punto más alto con dirección a la costa hasta el punto más bajo, si multiplicamos 2.8cm x 8 = 22.44 cm que es el desplazamiento en la línea de visión del satélite (line-of-sight), oblicuo.

Esto compara con el estudio efectuado por GPS/GNSS en el área de mayor impacto.

Conclusión:

Los profesionales de la agrimensura (topografía) deben tener mucho cuidado en la zona afectada por el terremoto ya que estos controles perdieron su posición vertical en su mayoría y aquellos que tenían posición horizontal también y no es aconsejable el utilizar los mismo para cualquier tipo de estudio sean estos:

Topografía, nivelación como tampoco aquellos que tenían posición horizontal por ajuste de redes (Magas y Lajas 2) y tampoco los que tenían posición publicada por OPUS SHARE.

Bibliografía:

  1. Red Sísmica de Puerto Rico
  2. Real Academia de Ciencia Exactas, Física y Naturales (RACEFM) España.
  3. Tesis de Maestría Sr. Gilberto Hermosillo Camacho (Universidad de Salamanca España)                               septiembre 2018
  4. Wikipedia.org
  5. http://www.gissandbeers.com

Contribución y agradecimientos: A los colegas agrimensores que al yo llamarlos y explicarles el propósito de mi intención de analizar los efectos en los controles horizontales y verticales en el área más afectada por el terremoto del día 7 de enero del corriente año, dijeron cuente con nosotros.

  1. Javier Osvaldo Escobales Medina, PS, (R. López de Azua & Asociados, PSC)
  2. Angel Noel Colón Guzmán, PS  (Global Survey, PSC)
  3. Luis S. Berrios Montes, PS  (Luis Berrios & Asociados)
  4. Carlos Lebron Cabrera, PS (JEB & Asociados)
  5. Victor M. Seda Figueroa  (V. Seda & Associates, PSC)
  6.  Luis E. Torres Román, PS
  7. Alex Hornedo Robles, PS, PE  (A Hornedo Robles & Asociados, PSC)
  8. Sr. Daniel Diaz Torres, MSPH
  9. Luis A. Maldonado Perez, PS, MGST.