Rainfall vs Soil Moisture

import matplotlib.pyplot as plt
import numpy as np
from matplotlib.colors import LinearSegmentedColormap

np.random.seed(42)

# Environmental data
rainfall = np.random.exponential(20, 4000)  # mm
moisture = 15 + rainfall * 1.2 + np.random.normal(0, 8, 4000)
moisture = np.clip(moisture, 5, 100)

fig, ax = plt.subplots(figsize=(10, 8), facecolor='#020B14')
ax.set_facecolor('#020B14')

# Custom colormap: dark to blue to cyan
colors = ['#020B14', '#0d1a3d', '#276CF5', '#27D3F5', '#27F5B0']
cmap = LinearSegmentedColormap.from_list('water', colors, N=256)

h = ax.hist2d(rainfall, moisture, bins=45, cmap=cmap, cmin=1)
cbar = plt.colorbar(h[3], ax=ax, pad=0.02)
cbar.set_label('Measurements', color='white', fontsize=11)
cbar.ax.yaxis.set_tick_params(color='white')
plt.setp(plt.getp(cbar.ax.axes, 'yticklabels'), color='white')

ax.set_xlabel('Rainfall (mm)', fontsize=11, color='white', fontweight='500')
ax.set_ylabel('Soil Moisture (%)', fontsize=11, color='white', fontweight='500')
ax.set_title('Rainfall vs Soil Moisture', fontsize=14, color='white', fontweight='bold', pad=15)

ax.tick_params(colors='white', labelsize=9)
for spine in ax.spines.values():
    spine.set_color('#333333')

plt.tight_layout()
plt.show()