  Substantial progress has been made recently in the generation of master-field
ensembles. This has to be paired with efficient techniques to compute
observables on gauge field configurations with a large volume. Here we present
the results of the computation of hadronic observables, including hadron masses
and meson decay constants, on large-volume and master-field ensembles with
physical volumes of up to $(18\,\mathrm{fm})^4$ and $m_\pi L$ up to $25$,
simulated using $N_{\mathrm{f}}=2+1$ stabilized Wilson fermions. We obtain
sub-percent determinations from single gauge configurations with the combined
use of position-space techniques, volume averages and master-field error
estimation.
