Seismic Activity Horizon

import matplotlib.pyplot as plt
import numpy as np

COLORS = {
    'bands': ['#4d0d0d', '#8a1a1a', '#EF4444', '#F87171'],
    'background': '#0a0a0f',
    'text': '#ffffff',
}

np.random.seed(1212)
hours = np.arange(0, 720)  # 30 days
# Seismic: mostly quiet with occasional activity clusters
base = np.random.exponential(0.5, len(hours))
aftershocks = np.zeros(len(hours))
main_events = np.random.choice(len(hours), 5, replace=False)
for evt in main_events:
    decay = np.exp(-np.abs(np.arange(len(hours)) - evt) / 20)
    aftershocks += np.random.uniform(2, 5) * decay
activity = base + aftershocks
activity = np.clip(activity, 0, 8)
activity_norm = activity / 8

fig, ax = plt.subplots(figsize=(14, 3), facecolor=COLORS['background'])
ax.set_facecolor(COLORS['background'])

band = 0.25

ax.fill_between(hours, 0, np.clip(activity_norm, 0, band), color=COLORS['bands'][0])
ax.fill_between(hours, 0, np.clip(activity_norm - band, 0, band), color=COLORS['bands'][1])
ax.fill_between(hours, 0, np.clip(activity_norm - 2*band, 0, band), color=COLORS['bands'][2])
ax.fill_between(hours, 0, np.clip(activity_norm - 3*band, 0, band), color=COLORS['bands'][3])

ax.set_xlim(0, 719)
ax.set_ylim(0, 0.3)

ax.set_title('Seismic Activity (Magnitude)', color=COLORS['text'], fontsize=12, fontweight='bold', pad=10)
ax.set_xlabel('Day', color=COLORS['text'], fontsize=10)
ax.set_ylabel('Magnitude', color=COLORS['text'], fontsize=10)
ax.set_xticks([0, 168, 336, 504, 672])
ax.set_xticklabels(['Week 1', 'Week 2', 'Week 3', 'Week 4', ''])

for spine in ax.spines.values():
    spine.set_visible(False)
ax.tick_params(colors=COLORS['text'], labelsize=9)

plt.tight_layout()
plt.show()