diff --git a/tests/spherical_nsinker_bfbt.cc b/tests/spherical_nsinker_bfbt.cc
new file mode 100644
index 00000000000..19f0a384ea2
--- /dev/null
+++ b/tests/spherical_nsinker_bfbt.cc
@@ -0,0 +1,21 @@
+/*
+  Copyright (C) 2022 by the authors of the ASPECT code.
+
+  This file is part of ASPECT.
+
+  ASPECT is free software; you can redistribute it and/or modify
+  it under the terms of the GNU General Public License as published by
+  the Free Software Foundation; either version 2, or (at your option)
+  any later version.
+
+  ASPECT is distributed in the hope that it will be useful,
+  but WITHOUT ANY WARRANTY; without even the implied warranty of
+  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+  GNU General Public License for more details.
+
+  You should have received a copy of the GNU General Public License
+  along with ASPECT; see the file LICENSE.  If not see
+  <http://www.gnu.org/licenses/>.
+*/
+
+#include "../benchmarks/nsinker_spherical_shell/nsinker.cc"
diff --git a/tests/spherical_nsinker_bfbt.prm b/tests/spherical_nsinker_bfbt.prm
new file mode 100644
index 00000000000..43d20c5f782
--- /dev/null
+++ b/tests/spherical_nsinker_bfbt.prm
@@ -0,0 +1,32 @@
+set Dimension = 3
+
+include $ASPECT_SOURCE_DIR/benchmarks/nsinker_spherical_shell/gmg.prm
+
+subsection Solver parameters
+  subsection Stokes solver parameters
+    set Use full A block as preconditioner              = true
+    set Stokes solver type                              = block AMG
+    set Number of cheap Stokes solver steps             = 500
+    set Maximum number of expensive Stokes solver steps = 1000
+    set Linear solver tolerance                         = 1e-6
+    set GMRES solver restart length                     = 100
+    set Use weighted BFBT for Schur complement          = true
+    set Krylov method for cheap solver steps            = GMRES
+  end
+
+  subsection AMG parameters
+    set AMG aggregation threshold = 0.02
+  end
+end
+
+subsection Mesh refinement
+  set Initial global refinement                = 1
+  set Initial adaptive refinement              = 0
+end
+
+subsection Postprocess
+  subsection Memory statistics
+    # disable as it is not deterministic:
+    set Output peak virtual memory (VmPeak) = false
+  end
+end
diff --git a/tests/spherical_nsinker_bfbt/screen-output b/tests/spherical_nsinker_bfbt/screen-output
new file mode 100644
index 00000000000..8da06f3a37f
--- /dev/null
+++ b/tests/spherical_nsinker_bfbt/screen-output
@@ -0,0 +1,28 @@
+-----------------------------------------------------------------------------
+-----------------------------------------------------------------------------
+-----------------------------------------------------------------------------
+
+Loading shared library <./libspherical_nsinker_bfbt.debug.so>
+
+-----------------------------------------------------------------------------
+-----------------------------------------------------------------------------
+Number of active cells: 768 (on 2 levels)
+Number of degrees of freedom: 28,690 (20,790+970+6,930)
+
+*** Timestep 0:  t=0 seconds, dt=0 seconds
+   Rebuilding Stokes preconditioner...
+   Solving Stokes system (AMG-BFBT)... 65+0 iterations.
+
+   Postprocessing:
+     System matrix memory consumption:  35.11 MB
+     Writing graphical output:          output-spherical_nsinker_bfbt/solution/solution-00000
+
+Termination requested by criterion: end time
+
+
++---------------------------------------------+------------+------------+
++---------------------------------+-----------+------------+------------+
++---------------------------------+-----------+------------+------------+
+
+-----------------------------------------------------------------------------
+-----------------------------------------------------------------------------