EstimaciónEstimation
El número acumulado de terremotos de magnitud igual o superior a M decae exponencialmente con la magnitud. Por encima de la magnitud de completitud Mc (el umbral a partir del cual la red detecta todos los eventos) el parámetro b se estima por máxima verosimilitud:
b = log10(e) / [ ⟨M⟩ − (Mc − ΔM/2) ]
(Aki, 1965; corrección por discretización de Utsu, 1966, con ΔM = 0.1). La incertidumbre se acota con
σb = 2.3 · b² · √[ Σ(Mi − ⟨M⟩)² / (N(N − 1)) ]
y a = log10 N(≥Mc) + b·Mc mide la productividad sísmica de la ventana. Por defecto se usa Mc = 4.5, la completitud aproximada de la red global del USGS. Como alternativa, Mc puede estimarse por máxima curvatura (moda de la distribución no acumulada) con la corrección de +0.2 de Woessner & Wiemer (2005), pero ese método solo es válido en catálogos espacialmente homogéneos.
Advertencia metodológica: el feed global mezcla redes con completitudes muy distintas (California y Alaska detectan M < 1; gran parte del planeta solo M > 4) y escalas de magnitud heterogéneas (ML, Md, mb, Mww). El b resultante describe este catálogo compuesto, no una región tectónica. Un análisis riguroso exigiría homogeneizar a Mw y restringir espacialmente.
The cumulative number of earthquakes with magnitude greater than or equal to M decays exponentially with magnitude. Above the magnitude of completeness Mc (the threshold above which the network detects every event) the b parameter is estimated by maximum likelihood:
b = log10(e) / [ ⟨M⟩ − (Mc − ΔM/2) ]
(Aki, 1965; binning correction by Utsu, 1966, with ΔM = 0.1). Its uncertainty is bounded by
σb = 2.3 · b² · √[ Σ(Mi − ⟨M⟩)² / (N(N − 1)) ]
and a = log10 N(≥Mc) + b·Mc measures the seismic productivity of the window. By default Mc = 4.5 is used, the approximate completeness of the USGS global network. Alternatively, Mc can be estimated by maximum curvature (mode of the non-cumulative distribution) with the +0.2 correction of Woessner & Wiemer (2005), but that method is only valid for spatially homogeneous catalogues.
Methodological caveat: the global feed mixes networks with very different completeness (California and Alaska detect M < 1; much of the planet only M > 4) and heterogeneous magnitude scales (ML, Md, mb, Mww). The resulting b describes this composite catalogue, not a tectonic region. A rigorous analysis would require homogenising to Mw and restricting spatially.
Significado físico de bPhysical meaning of b
Para la sismicidad tectónica global b ≈ 1: cada unidad de magnitud se producen unas diez veces menos eventos. Como la energía radiada escala como log10E ≈ 1.5M + 4.8, los grandes terremotos, aunque raros, dominan el balance energético.
b como esfuerzómetro. En experimentos de emisión acústica en roca, b disminuye al aumentar el esfuerzo diferencial y antes de la rotura (Scholz, 1968). En la corteza, Schorlemmer et al. (2005) muestran una dependencia sistemática con el régimen de esfuerzos: fallas inversas (b bajo, mayor esfuerzo diferencial) < desgarre < normales (b alto).
b alto (> 1.2): medios muy heterogéneos o fracturados, altas presiones de fluido y gradientes térmicos — enjambres volcánicos y geotérmicos, sismicidad inducida. b bajo (< 0.8): zonas de alto esfuerzo, asperezas bloqueadas.
Geometría fractal. Bajo hipótesis de autosimilitud, la dimensión fractal de la distribución de rupturas cumple D ≈ 2b (Aki, 1981): b codifica cómo se reparte la fracturación entre escalas.
Se han propuesto variaciones temporales de b como indicador de carga en asperezas, pero su poder predictivo prospectivo no está establecido: la estimación es muy sensible a Mc, al tamaño de muestra y a la homogeneidad del catálogo.
For global tectonic seismicity b ≈ 1: each unit of magnitude brings about ten times fewer events. Since radiated energy scales as log10E ≈ 1.5M + 4.8, large earthquakes, although rare, dominate the energy budget.
b as a stress meter. In acoustic emission experiments on rock, b decreases as differential stress increases and before failure (Scholz, 1968). In the crust, Schorlemmer et al. (2005) show a systematic dependence on the faulting regime: thrust faults (low b, higher differential stress) < strike-slip < normal (high b).
High b (> 1.2): highly heterogeneous or fractured media, high pore-fluid pressure and thermal gradients — volcanic and geothermal swarms, induced seismicity. Low b (< 0.8): highly stressed regions, locked asperities.
Fractal geometry. Under self-similarity, the fractal dimension of the rupture distribution satisfies D ≈ 2b (Aki, 1981): b encodes how fracturing is distributed across scales.
Temporal b-value variations have been proposed as an indicator of asperity loading, but their prospective predictive power is not established: the estimate is very sensitive to Mc, sample size and catalogue homogeneity.