  We use a variational approach to study existence and regularity of solutions
for a Neumann $p$-Laplacian problem with a reaction term on metric spaces
equipped with a doubling measure and supporting a Poincar\'e inequality. Trace
theorems for functions with bounded variation are applied in the definition of
the variational functional and minimizers are shown to satisfy De Giorgi type
conditions.
