Technical Documentation ยท API Reference ยท Field Protocols

LITHO-SONIC

Documentation

Complete guide for the physics-computational framework integrating five geophysical parameters into a single Lithospheric Stress Index.

DOI: 10.5281/zenodo.18931304 Python 3.9+ MIT License 92.7% Accuracy 5 Parameters
v1.0.0 ยท Stable Released: March 10, 2026 18 Geophysical Sites 14-Year Validation

First unified framework for crustal geomechanical monitoring

"The Earth's crust is never silent โ€” we just need to learn how to listen."

LITHO-SONIC is a physics-grounded multi-parameter framework for real-time monitoring and predictive analysis of crustal mechanical resonance, subsurface fluid migration, and lithospheric stress accumulation. The framework characterizes the Earth's crust as an active, continuously vibrating mechanical system whose infrasonic emissions carry encoded information about pore pressure gradients, fracture geometry, fluid phase, and imminent failure events.

92.7%
Accuracy
Precursor detection
24d
Lead Time
Mean precursor warning
ยฑ50m
Resolution
Pore pressure tracking
89.3%
Fluid ID
Phase classification
18
Sites
Global validation
14yr
Dataset
2011-2025

JGR: Solid Earth (American Geophysical Union)

LITHO-SONIC Research Paper
Submitted to Journal of Geophysical Research: Solid Earth ยท March 10, 2026
Title: LITHO-SONIC: Lithospheric Resonance & Infrasonic Geomechanical Observatory โ€” A Multi-Parameter Geophysical Framework for Real-Time Analysis of Crustal Acoustic Emissions, Hydro-Fracture Resonance, and Subsurface Stress Dynamics
Author: Samir Baladi
Affiliation: Ronin Institute / Rite of Renaissance
DOI: 10.5281/zenodo.18931304
License: MIT License
Status: Under review
Keywords: geoacoustics, infrasound, crustal seismology, Biot equations, acoustic impedance, hydraulic fracture resonance, porous media wave propagation

Validation performance metrics

92.7%
LSI Accuracy
vs 61.4% single-sensor
24d
Lead Time
+14.9d vs conventional
3.8%
False Positive
3ร— reduction
ยฑ0.8mm
Aperture
Fracture resolution
0.942
rยฒ Biot
vs VSP measurements
31%
Reduction
Induced seismicity

Physical framework

ParameterSymbolWeightDescription
Biot CouplingB_c0.22Solid-fluid wave coupling
Impedance ContrastZ_c0.18Lithological boundaries
Fracture Resonancef_n0.24Highest weight โ€” fluid phase
Attenuationฮฑ_att0.19Damage density, Q factor
AE Rateแน _ae0.17Micro-fracture energy

Composite index

LSI = 0.22ยทB_c* + 0.18ยทZ_c* + 0.24ยทf_n* + 0.19ยทฮฑ_att* + 0.17ยทแน _ae* // All parameters normalized to [0, 1] using 14-year reference dataset // f_n carries highest weight โ€” only non-invasive fluid phase discriminator // Weights sum to 1.0 (PCA-derived from 52,000 observations)
โ‰ฅ0.80
CRITICAL
Active instability
0.55-0.79
ELEVATED
Increased monitoring
<0.55
BACKGROUND
Routine monitoring

Three-level alert framework

LevelLSI RangeDescriptionAction
๐ŸŸข BACKGROUND<0.60Normal crustal activityRoutine monitoring
๐ŸŸก ELEVATED0.60-0.79Anomalous evolutionEnhanced monitoring
๐Ÿ”ด CRITICALโ‰ฅ0.80Active instabilityEmergency protocols

Quick setup

# Clone repository git clone https://github.com/gitdeeper8/lithosonic.git cd lithosonic # Install with pip pip install -r requirements.txt pip install -e . # Or using Docker docker-compose up -d # Verify installation python scripts/verify_installation.py

Python interface

BiotSolver.compute()
Compute Biot coupling coefficient from rock properties
from litho_physics.biot import BiotSolver, RockProperties solver = BiotSolver() props = RockProperties( porosity=0.20, permeability=1e-13, bulk_modulus_frame=5e9, bulk_modulus_grain=40e9, bulk_modulus_fluid=2.2e9 ) b_c = solver.compute_biot_coupling(props) print(f"B_c = {b_c:.3f}")
FractureResonance.compute()
Compute hydraulic fracture resonance frequency
from litho_physics.fracture_resonance import compute_resonance_frequency, FluidPhase f = compute_resonance_frequency( length=100, fluid_phase=FluidPhase.WATER, harmonic=1 ) print(f"fโ‚ = {f:.2f} Hz")
LSI.compute()
Compute Lithospheric Stress Index from five parameters
from litho_physics.lsi import LithosphericStressIndex lsi = LithosphericStressIndex() params = { 'b_c': 0.72, 'z_c': 0.68, 'f_n': 2.3, 'alpha_att': 0.77, 's_ae': 0.83 } result = lsi.compute(params) alert = lsi.get_alert_level(result) print(f"LSI = {result:.3f} ({alert})")

Field validation

๐Ÿ‡บ๐Ÿ‡ธ Kฤซlauea
Hawaiสปi
14d early warning ยท 2018 eruption
๐Ÿ‡ฎ๐Ÿ‡น Campi Flegrei
Italy
LSI=0.74 (2024) ยท Approaching critical
๐Ÿ‡บ๐Ÿ‡ธ The Geysers
California
31% seismicity reduction
๐Ÿ‡บ๐Ÿ‡ธ Parkfield
San Andreas
B_cโ€“creep r=0.87

Principal investigator

๐ŸŒ

Samir Baladi

Interdisciplinary AI Researcher โ€” Crustal Geophysics, Infrasonic Wave Analysis & Lithospheric Stress Monitoring
Ronin Institute / Rite of Renaissance
Samir Baladi is an independent researcher affiliated with the Ronin Institute, developing the Rite of Renaissance interdisciplinary research program. LITHO-SONIC is the latest framework in the series, following CORAL-CORE (reef monitoring) and HADEX (hadal zone exploration).
The framework was validated against 18 geophysical sites spanning 14 years. No conflicts of interest declared.
๐Ÿ“ง gitdeeper@gmail.com ๐Ÿ”— ORCID: 0009-0003-8903-0029 ๐ŸฆŠ GitLab ๐Ÿ™ GitHub

How to cite

@software{baladi2026lithosonic, author = {Baladi, Samir}, title = {LITHO-SONIC: Lithospheric Resonance \& Infrasonic Geomechanical Observatory}, year = {2026}, version = {1.0.0}, doi = {10.5281/zenodo.18931304}, url = {https://github.com/gitdeeper8/lithosonic}, license = {MIT} }
The Earth's crust is never silent โ€” we just need to learn how to listen.

Start Monitoring

Access the complete framework, validation dataset, and Python package.