1 Resumen Ejecutivo ›
Este documento define la estrategia de integracion con los productos de la GEM Foundation (Global Earthquake Model) para la plataforma PSHA Intelligence. Tras un analisis exhaustivo de licencias, capacidades tecnicas y riesgos legales, se adopta una arquitectura hibrida.
Arquitectura hibrida: Motor PSHA propio (NumPy/SciPy) para produccion + OpenQuake en Docker para validacion cruzada.
Rationale: Mantiene independencia de licencias AGPL, permite benchmarking riguroso y evita dependencia de un unico motor de calculo.
Hallazgos Clave
- Solo GAF-DB (Global Active Faults Database) tiene licencia CC BY-SA 4.0, compatible con uso comercial SaaS sin restricciones.
- Todos los demas productos GEM (Hazard Map, Mosaic, Risk Profiles, Exposure, Vulnerability) son CC BY-NC-SA 4.0 — uso comercial prohibido sin licencia adicional.
- OpenQuake Engine es AGPL-3.0 — cualquier
import openquakeactiva copyleft. Solo la comunicacion via REST API (Docker container aislado) es segura. - GEM tiene una politica de 18 meses tras la cual algunos datasets migran de CC BY-NC-SA a CC BY-SA — monitorizar activamente.
2 Matriz de Licenciamiento GEM ›
Analisis completo de los 7 productos principales de GEM Foundation y su viabilidad para uso comercial en la plataforma SaaS.
| Producto GEM | Licencia | Uso Comercial SaaS | Uso como Input | Accion Requerida |
|---|---|---|---|---|
| GAF-DB (Fallas Activas) | CC BY-SA 4.0 | SI | SI | Atribucion en reportes PDF y UI |
| Global Seismic Hazard Map | CC BY-NC-SA 4.0 | NO | NO | Solo validacion interna |
| OpenQuake Engine | AGPL-3.0 | Solo REST API | N/A | Docker container, sin import |
| Mosaic (34 modelos regionales) | CC BY-NC-SA 4.0 | NO | NO | Contactar GEM para licencia comercial |
| Risk Profiles (190 paises) | CC BY-NC-SA 4.0 | NO | NO | PNG solamente, sin datos subyacentes |
| Global Vulnerability Model | CC BY-NC-SA 4.0 | NO | NO | No necesario para Tiers 1-3 |
| Global Exposure Model | CC BY-NC-SA 4.0 | NO | NO | Contactar GEM para Tier 0 |
product@globalquakemodel.org para negociar acceso a Mosaic, Exposure y Risk Profiles. Precio estimado: negociacion caso a caso, tipicamente suscripcion anual.
Resumen de Estrategia por Producto
3 Estrategia de Integracion OpenQuake ›
OpenQuake Engine (AGPL-3.0) es el estandar de facto para PSHA en la industria. Sin embargo, su licencia AGPL impone requisitos de copyleft que afectan directamente a cualquier software que lo importe. A continuacion se analiza la matriz de riesgo AGPL.
Matriz de Riesgo AGPL
| Metodo de Integracion | Ejemplo | Riesgo Copyleft | Veredicto |
|---|---|---|---|
import openquake |
Importar modulos OQ en nuestro codigo Python | COPYLEFT TRIGGERED | NO USAR |
subprocess.call() |
Invocar binario OQ como subproceso | INCIERTO | EVITAR |
| REST API (Docker) | HTTP requests a oq-engine:8800 |
NO TRIGGERED | VIABLE |
| Validacion interna | Comparar outputs sin distribuir OQ | NO TRIGGERED | SEGURO |
from openquake.hazardlib import *, import openquake.engine, o cualquier referencia directa. El copyleft AGPL se extiende a todo el software que enlaza con codigo AGPL.
Motor PSHA propio para produccion + OpenQuake Docker container (REST :8800) exclusivamente para validacion cruzada.
Rationale: Separacion total de procesos. Nuestro codigo nunca toca el codigo AGPL. OpenQuake corre en su propio container aislado, comunicacion unicamente via HTTP.
Arquitectura Hibrida
Flujo de Validacion Cruzada
services:
oq-engine:
image: openquake/engine:latest
ports:
- "8800:8800"
environment:
- OQ_CONFIG_FILE=/etc/openquake/openquake.cfg
volumes:
- oq-data:/var/lib/openquake
networks:
- validation-net # Red aislada, sin acceso al pipeline
deploy:
resources:
limits:
cpus: '2'
memory: 4G
# Nuestro pipeline NUNCA esta en esta red
# Comunicacion SOLO via HTTP desde red externa
volumes:
oq-data:
networks:
validation-net:
internal: false
4 Formato NRML v0.5 ›
NRML (Natural Risks Markup Language) v0.5 es el formato XML estandar de OpenQuake para definir modelos de fuentes sismicas, arboles logicos y resultados de PSHA. Nuestra plataforma implementa soporte completo de lectura y escritura (round-trip) para interoperabilidad con OpenQuake y otros motores.
Namespace y Estructura
<?xml version="1.0" encoding="UTF-8"?>
<nrml xmlns="http://openquake.org/xmlns/nrml/0.5">
<sourceModel name="PSHA Platform - SSM">
<!-- Source definitions -->
</sourceModel>
</nrml>
Tipos de Fuente Soportados
| Tipo de Fuente | Elemento NRML | Descripcion | Uso Tipico |
|---|---|---|---|
| Falla Simple | simpleFaultSource |
Falla definida por traza superficial + dip + profundidad | Fallas corticales bien caracterizadas |
| Falla Compleja | complexFaultSource |
Falla definida por bordes superior e inferior 3D | Zonas de subduccion |
| Fuente de Area | areaSource |
Poligono con sismicidad distribuida uniformemente | Zonas de sismicidad difusa (background) |
| Fuente Puntual | pointSource |
Punto con MFD y geometria de ruptura | Fuentes individuales o grid de puntos |
Arboles Logicos
El formato NRML define dos tipos de arboles logicos que nuestra plataforma soporta completamente:
- SSM Logic Tree (
ssmLT.xml) — Define alternativas para el modelo de fuentes sismicas: variaciones de geometria, parametros de recurrencia, magnitud maxima. - GMPE Logic Tree (
gmmLT.xml) — Define las ecuaciones de movimiento del suelo aplicables a cada tipo tectonico: cortical, subduccion, estable.
<nrml xmlns="http://openquake.org/xmlns/nrml/0.5">
<logicTree logicTreeID="lt1">
<logicTreeBranchSet uncertaintyType="sourceModel"
branchSetID="bs1">
<logicTreeBranch branchID="b1">
<uncertaintyModel>ssm_lower.xml</uncertaintyModel>
<uncertaintyWeight>0.2</uncertaintyWeight>
</logicTreeBranch>
<logicTreeBranch branchID="b2">
<uncertaintyModel>ssm_mean.xml</uncertaintyModel>
<uncertaintyWeight>0.6</uncertaintyWeight>
</logicTreeBranch>
<logicTreeBranch branchID="b3">
<uncertaintyModel>ssm_upper.xml</uncertaintyModel>
<uncertaintyWeight>0.2</uncertaintyWeight>
</logicTreeBranch>
</logicTreeBranchSet>
</logicTree>
</nrml>
Implementacion en el Pipeline
Nuestro parser y serializer NRML garantizan round-trip fidelity: un archivo NRML parseado y re-serializado produce un output semanticamente identico al original.
| Modulo | Ruta | Funcion |
|---|---|---|
| Parser | pipeline/extract/nrml_parser.py |
Lee NRML v0.5 → objetos Python (dataclasses) |
| Serializer | pipeline/transform/nrml_serializer.py |
Objetos Python → NRML v0.5 XML valido |
5 Biblioteca de GMPEs (Prioridad) ›
La seleccion de GMPEs (Ground Motion Prediction Equations) sigue las recomendaciones de los informes de referencia de Dynamis y Langan. Se implementan en tres fases, priorizando las ecuaciones con mayor aplicabilidad global.
Fase 1: NGA-West 2 (Estandar Universal)
Las cuatro GMPEs del programa NGA-West 2 son el estandar de la industria para regiones corticales activas. Se asigna peso igual de 0.25 a cada una en el arbol logico.
| GMPE | Referencia | Peso | Rango Mw | Rango Rrup (km) |
|---|---|---|---|---|
| ASK14 | Abrahamson, Silva & Kamai (2014) | 0.25 | 3.0 - 8.5 | 0 - 300 |
| BSSA14 | Boore, Stewart, Seyhan & Atkinson (2014) | 0.25 | 3.0 - 8.5 | 0 - 400 |
| CB14 | Campbell & Bozorgnia (2014) | 0.25 | 3.3 - 8.5 | 0 - 300 |
| CY14 | Chiou & Youngs (2014) | 0.25 | 3.5 - 8.5 | 0 - 300 |
Fase 2: NGA-Sub (Subduccion)
Para zonas de subduccion, las GMPEs del programa NGA-Subduction complementan a NGA-West 2.
| GMPE | Referencia | Aplicabilidad |
|---|---|---|
| AG2015 / AG2020 | Abrahamson & Gulerce (2015, 2020) | Interface + Intraslab |
| Parker2020 | Parker et al. (2020) | Interface + Intraslab |
| Kuehn2020 | Kuehn et al. (2020) | Interface + Intraslab |
Fase 3: Modelos Regionales
| GMPE | Region | Contexto Tectonico |
|---|---|---|
| Akkar et al. (2014) | Europa y Mediteraneo | Cortical activa |
| Zhao et al. (2006) | Japon / Pacifico | Subduccion + cortical |
calculate_mean(), calculate_sigma()). Unit tests contra tablas publicadas por los autores. Validacion cruzada contra OpenQuake hazardlib (via REST API).
6 GEM Hazard Mosaic — Datos de Referencia ›
El GEM Global Seismic Hazard Mosaic es un mapa global de peligrosidad sismica construido a partir de 34 modelos regionales armonizados. Se distribuye como GeoTIFF y es una referencia clave para validacion.
Especificaciones Tecnicas del GeoTIFF
| Parametro | Valor |
|---|---|
| Dimensiones | 7200 x 3000 pixels |
| Tipo de dato | float64 |
| Resolucion | 0.05 grados (~5.5 km en el ecuador) |
| Variable | PGA (g) para periodo de retorno de 475 anos |
| Cobertura | Global (180W-180E, 75S-75N) |
| Licencia | CC BY-NC-SA 4.0 |
Uso como Benchmark
El Mosaic se utiliza exclusivamente como lookup de PGA de referencia para validar los resultados de nuestro motor propio. No se incluye en ningun output entregado al cliente.
- Validacion: Comparar PGA calculado por nuestro motor vs. PGA del Mosaic para la misma coordenada.
- Tolerancia: Discrepancias > 30% activan revision manual (diferencias esperadas por diferencias metodologicas).
- No produccion: Nunca se entrega el valor del Mosaic al cliente como resultado.
Paleta de Color para Mapbox
Se adopta la escala de 11 bins de PGA definida por GEM para mantener consistencia visual con el estandar de la industria.
| Bin | Rango PGA (g) | Color | Nivel de Peligrosidad |
|---|---|---|---|
| 1 | < 0.02 | #FFFFFF | Muy baja |
| 2 | 0.02 - 0.04 | #D4EDFC | Baja |
| 3 | 0.04 - 0.08 | #ABDDA4 | Baja-Moderada |
| 4 | 0.08 - 0.12 | #66C2A5 | Moderada |
| 5 | 0.12 - 0.20 | #FDAE61 | Moderada-Alta |
| 6 | 0.20 - 0.28 | #F46D43 | Alta |
| 7 | 0.28 - 0.40 | #D53E4F | Alta |
| 8 | 0.40 - 0.56 | #B71C3D | Muy Alta |
| 9 | 0.56 - 0.80 | #9E0142 | Muy Alta |
| 10 | 0.80 - 1.20 | #67001F | Extrema |
| 11 | > 1.20 | #330010 | Extrema |
7 Producto Tier 0 — Explorador Regional ›
Tier 0 es un nuevo concepto de producto que ofrece un explorador regional de riesgo sismico como puerta de entrada al funnel comercial. Combina datos de contexto (exposicion, riesgo) con nuestros propios calculos PSHA.
Fases de Implementacion
Fase A: Contexto de Exposicion
Datos del GEM Global Exposure Model para contextualizar el entorno construido en la zona de estudio.
- Stock de edificaciones por tipologia constructiva
- Costos de reposicion estimados
- Poblacion expuesta
- Requiere licencia GEM o fuentes alternativas (censo nacional, OSM)
Fase B: Contexto de Riesgo
Metricas de riesgo sismico regional para dar perspectiva al cliente.
- Average Annual Loss (AAL) regional
- Distribucion historica de danos
- Requiere licencia GEM o computo propio
Fase C: Metricas PSHA Propias
Resultados generados por nuestro motor, sin dependencia de licencias externas.
- PGA y Sa para periodos de retorno estandar
- Mapa de fallas activas del entorno (GAF-DB)
- Catalogo sismico historico filtrado (USGS/ISC)
- Sin restricciones de licencia
Enfoque hibrido: datos de exposicion para contexto (con licencia GEM o alternativas open-source) + metricas PSHA propias para valor tecnico diferenciado.
Estimacion: 2-3 semanas para MVP. El Tier 0 sirve como demo interactiva y herramienta de captacion comercial.
Flujo del Producto
8 Arbol Logico (Estructura Estandar) ›
La estructura del arbol logico sigue las recomendaciones del informe Dynamis para capturar la incertidumbre epistemica en todas las etapas del analisis PSHA.
Fuentes Tipo Falla
Cada falla individual genera un sub-arbol con las siguientes ramas de incertidumbre:
| Parametro | Ramas | Pesos | Descripcion |
|---|---|---|---|
| Dip (buzamiento) | 3 | [0.2, 0.6, 0.2] | Bajo, central, alto |
| Modelo de recurrencia | 2 | [0.5, 0.5] | Gutenberg-Richter / Characteristic |
| Magnitud maxima (Mmax) | 3 | [0.2, 0.6, 0.2] | -0.2, central, +0.2 unidades Mw |
| Tasa de deslizamiento (slip rate) | 3 | [0.2, 0.6, 0.2] | Bajo, central, alto |
| GMPE | 4 | [0.25, 0.25, 0.25, 0.25] | ASK14, BSSA14, CB14, CY14 |
Fuentes Tipo Area (Background)
Las zonas de sismicidad distribuida tienen un arbol logico mas compacto:
| Parametro | Ramas | Pesos | Descripcion |
|---|---|---|---|
| Valor b (G-R) | 3 | [0.2, 0.6, 0.2] | b - sigma, b, b + sigma |
| Magnitud maxima (Mmax) | 3 | [0.2, 0.6, 0.2] | Mmax observada, +0.3, +0.5 |
| GMPE | 4 | [0.25, 0.25, 0.25, 0.25] | NGA-West 2 (4 modelos) |
9 Acciones Prioritarias ›
Lista priorizada de acciones derivadas de este analisis de integracion GEM.
product@globalquakemodel.org solicitando informacion sobre licencias comerciales para Risk Profiles, Exposure Model y Mosaic. Necesario para Tier 0 y posible expansion futura.docker-compose.validation.yml con container OQ aislado en red separada. Verificar que el REST API responde en :8800. Crear script de validacion automatizada.10 Roadmap de Implementacion ›
Cronograma de implementacion alineado con la estrategia de integracion GEM y el desarrollo del motor PSHA propio.
- Desplegar OpenQuake Docker container con REST API en
:8800 - Implementar export NRML v0.5 desde nuestro pipeline (
nrml_serializer.py) - Crear script de validacion automatizada que envie jobs a OQ y compare resultados
- Establecer pipeline CI con tests de integracion contra OQ
- Implementar integral Cornell-McGuire vectorizada (NumPy)
- Implementar 4 GMPEs NGA-West 2 (ASK14, BSSA14, CB14, CY14) + 1 adicional
- Implementar arbol logico basico (SSM + GMPE)
- Validacion cruzada contra OpenQuake: tolerancia ≤ 2-3% en hazard curves
- Primeros outputs Tier 1 y Tier 2 generados con motor propio
- Expandir biblioteca a 15-20 GMPEs (NGA-Sub, regionales)
- Implementar disaggregacion (M-R-epsilon)
- Arboles logicos completos con ramas de incertidumbre parametrica
- UHS (Uniform Hazard Spectra) para multiples periodos de retorno
- Optimizacion de rendimiento: paralelismo, caching de GMPEs
- Tier 3 PSHA completo con todas las funcionalidades
- Suite de benchmarking automatizada: nuestro motor vs. OpenQuake vs. R-CRISIS
- Objetivo de precision: ≤ 2-3% vs. OQ en todos los escenarios de test
- Documentacion tecnica de validacion para revision por pares
- Negociacion de licencia GEM completada (si aplica para Tier 0)
El motor PSHA propio se considera validado cuando produce hazard curves con discrepancia ≤ 2-3% respecto a OpenQuake para un conjunto de al menos 10 escenarios de referencia que cubran diferentes contextos tectonicos (cortical, subduccion, estable).