fig, ax = plt.subplots()
ax.set_xlabel('Basin area [$\mathrm{km^2}$]')