diff --git a/ProcessLib/LiquidFlow/CMakeLists.txt b/ProcessLib/LiquidFlow/CMakeLists.txt
index 8890ed1cd9b9df31d1ea736e15100d2eb80980f5..0cad455746411746646e36a9449853d699d6ab89 100644
--- a/ProcessLib/LiquidFlow/CMakeLists.txt
+++ b/ProcessLib/LiquidFlow/CMakeLists.txt
@@ -77,6 +77,21 @@ AddTest(
     hex.vtu isotropic_gravity_driven3D_pcs_0_ts_1_t_1.000000.vtu analytic_pressure pressure
 )
 
+# Coupling
+AddTest(
+    NAME StaggeredTH_ThermalDensityDrivenFlow2D
+    PATH StaggeredCoupledProcesses/TH/ThermalDensityDrivenFlow2D
+    EXECUTABLE ogs
+    EXECUTABLE_ARGS thermal_gravity_driven_flow2d.prj
+    WRAPPER time
+    TESTER vtkdiff
+    REQUIREMENTS NOT OGS_USE_MPI
+    ABSTOL 1e-20 RELTOL 1e-10
+    DIFF_DATA
+    mesh2D.vtu gravity_driven_pcs_1_ts_5_t_300.000000.vtu pressure_ref pressure
+    mesh2D.vtu gravity_driven_pcs_1_ts_5_t_300.000000.vtu temperature_ref temperature
+)
+
 # PETSc/MPI
 AddTest(
     NAME LiquidFlow_LineDirichletNeumannBC
@@ -154,3 +169,18 @@ AddTest(
     DIFF_DATA
     hex.vtu isotropic_gravity_driven3D_pcs_0_ts_1_t_1_000000_0.vtu analytic_pressure pressure
 )
+
+# Coupling
+AddTest(
+    NAME StaggeredTH_ThermalDensityDrivenFlow2D
+    PATH StaggeredCoupledProcesses/TH/ThermalDensityDrivenFlow2D
+    EXECUTABLE_ARGS thermal_gravity_driven_flow2d.prj
+    WRAPPER mpirun
+    WRAPPER_ARGS -np 1
+    TESTER vtkdiff
+    REQUIREMENTS OGS_USE_MPI
+    ABSTOL 1e-10 RELTOL 1e-10
+    DIFF_DATA
+    mesh2D.vtu gravity_driven_pcs_1_ts_5_t_300_000000_0.vtu pressure_ref pressure
+    mesh2D.vtu gravity_driven_pcs_1_ts_5_t_300_000000_0.vtu temperature_ref temperature
+)
diff --git a/Tests/Data b/Tests/Data
index 2204d266bd9a25175d947be6a1c2a136c54d763f..2f53e1700c2bafb806c8fd2d7e9adfe85fa9d389 160000
--- a/Tests/Data
+++ b/Tests/Data
@@ -1 +1 @@
-Subproject commit 2204d266bd9a25175d947be6a1c2a136c54d763f
+Subproject commit 2f53e1700c2bafb806c8fd2d7e9adfe85fa9d389