diff --git a/yt/visualization/_handlers.py b/yt/visualization/_handlers.py
index 659b1eb6d3f..ec52acf6c05 100644
--- a/yt/visualization/_handlers.py
+++ b/yt/visualization/_handlers.py
@@ -6,7 +6,7 @@
 import matplotlib as mpl
 import numpy as np
 import unyt as un
-from matplotlib.colors import Colormap, LogNorm, Normalize, SymLogNorm
+from matplotlib.colors import Colormap, LogNorm, Normalize, SymLogNorm, TwoSlopeNorm
 from unyt import unyt_quantity
 
 from yt._typing import Quantity, Unit
@@ -61,13 +61,14 @@ class NormHandler:
         "_dynamic_range",
         "_norm_type",
         "_linthresh",
-        "_norm_type",
+        "_vcenter",
         "_norm",
         "prefer_log",
     )
     _constraint_attrs: list[str] = [
         "vmin",
         "vmax",
+        "vcenter",
         "dynamic_range",
         "norm_type",
         "linthresh",
@@ -84,6 +85,7 @@ def __init__(
         norm_type: Optional[type[Normalize]] = None,
         norm: Optional[Normalize] = None,
         linthresh: Optional[float] = None,
+        vcenter: Optional[float] = None,
     ):
         self.data_source = weakref.proxy(data_source)
         self.ds = data_source.ds  # should already be a weakref proxy
@@ -95,6 +97,7 @@ def __init__(
         self._dynamic_range = dynamic_range
         self._norm_type = norm_type
         self._linthresh = linthresh
+        self._vcenter = vcenter
         self.prefer_log = True
 
         if self.norm is not None and self.has_constraints:
@@ -196,6 +199,22 @@ def vmax(self, newval: Optional[Union[Quantity, float, Literal["max"]]]) -> None
         else:
             self._set_quan_attr("_vmax", newval)
 
+    @property
+    def vcenter(self) -> Optional[Union[un.unyt_quantity, Literal["zero"], float]]:
+        return self._vcenter
+
+    @vcenter.setter
+    def vcenter(
+        self, newval: Optional[Union[un.unyt_quantity, Literal["zero"], float]]
+    ) -> None:
+        self._reset_norm()
+        if newval == "zero":
+            self._vcenter = 0.0
+        else:
+            self._set_quan_attr("_vcenter", newval)
+        if newval is not None:
+            self.norm_type = TwoSlopeNorm
+
     @property
     def dynamic_range(self) -> Optional[float]:
         return self._dynamic_range
@@ -364,6 +383,12 @@ def get_norm(self, data: np.ndarray, *args, **kw) -> Normalize:
 
             kw.setdefault("linthresh", linthresh)
             kw.setdefault("base", 10)
+        elif norm_type is TwoSlopeNorm:
+            if self.vcenter is not None:
+                vcenter = self.to_float(self.vcenter)
+            else:
+                vcenter = 0.0
+            kw.setdefault("vcenter", vcenter)
 
         return norm_type(*args, **kw)