We re-examine current methodology for estimating seismic b-values.
We propose a new work-flow for calculating completeness magnitudes.
A pragmatic minimum complete catalogue size of 200 events is proposed.
The error in b-value from the selection of the completeness magnitude is quantified.
Temporal b-value variations at El Hierro and Mt. Etna are analysed.