diff --git a/CHANGELOG.md b/CHANGELOG.md
index 34626ee5..d632bd4a 100644
--- a/CHANGELOG.md
+++ b/CHANGELOG.md
@@ -1,3 +1,18 @@
+# 2.1.0
+* Added `pylops.signalprocessing.DCT`, `pylops.signalprocessing.NonStationaryConvolve1D`,
+  `pylops.signalprocessing.NonStationaryConvolve2D`, `pylops.signalprocessing.NonStationaryFilters1D`, and
+  `pylops.signalprocessing.NonStationaryFilters2D` operators
+* Added `pylops.waveeqprocessing.BlendingContinuous`, `pylops.waveeqprocessing.BlendingGroup`, and
+  `pylops.waveeqprocessing.BlendingHalf` operators
+* Added `kind='datamodel'` to `pylops.optimization.cls_sparsity.IRLS`
+* Improved inner working of `pylops.waveeqprocessing.Kirchhoff` operator significantly
+  reducing the memory usage related to storing traveltime, angle, and amplitude tables.
+* Improved handling of `haxes` in `pylops.signalprocessing.Radon2D` and `pylops.signalprocessing.Radon3D` operators
+* Added possibility to feed ND-arrays to `pylops.TorchOperator`
+* Removed `pylops.LinearOperator` inheritance and added `__call__` method to `pylops.TorchOperator`
+* Removed `scipy.sparse.linalg.LinearOperator` and added `abc.ABC` inheritance to `pylops.LinearOperator`
+* All operators are now classes of `pylops.LinearOperator` type
+
 # 2.0.0
 PyLops has undergone significant changes in this release, including new ``LinearOperator``s, more features, new examples and bugfixes.
 To aid users in navigating the breaking changes, we provide the following document [MIGRATION_V1_V2.md](https://github.com/PyLops/pylops/blob/dev/MIGRATION_V1_V2.md).
diff --git a/docs/source/changelog.rst b/docs/source/changelog.rst
index 33073794..a56dd9aa 100644
--- a/docs/source/changelog.rst
+++ b/docs/source/changelog.rst
@@ -3,6 +3,26 @@
 Changelog
 =========
 
+Version 2.1.0
+-------------
+
+*Released on: 17/03/2023*
+
+* Added :class:`pylops.signalprocessing.DCT`, :class:`pylops.signalprocessing.NonStationaryConvolve1D`,
+  :class:`pylops.signalprocessing.NonStationaryConvolve2D`, :class:`pylops.signalprocessing.NonStationaryFilters1D`, and
+  :class:`pylops.signalprocessing.NonStationaryFilters2D` operators
+* Added :class:`pylops.waveeqprocessing.BlendingContinuous`, :class:`pylops.waveeqprocessing.BlendingGroup`, and
+  :class:`pylops.waveeqprocessing.BlendingHalf` operators
+* Added `kind='datamodel'` to :class:`pylops.optimization.cls_sparsity.IRLS`
+* Improved inner working of :class:`pylops.waveeqprocessing.Kirchhoff` operator significantly
+  reducing the memory usage related to storing traveltime, angle, and amplitude tables.
+* Improved handling of `haxes` in :class:`pylops.signalprocessing.Radon2D` and :class:`pylops.signalprocessing.Radon3D` operators
+* Added possibility to feed ND-arrays to :class:`pylops.TorchOperator`
+* Removed :class:`pylops.LinearOperator` inheritance and added `__call__` method to :class:`pylops.TorchOperator`
+* Removed `scipy.sparse.linalg.LinearOperator` and added :class:`abc.ABC` inheritance to :class:`pylops.LinearOperator`
+* All operators are now classes of `:class:`pylops.LinearOperator` type
+
+
 Version 2.0.0
 -------------
 
@@ -98,7 +118,7 @@ Version 1.18.0
 * Extended :py:func:`pylops.Laplacian` to N-dimensional arrays
 * Added `forward` kind to :py:class:`pylops.SecondDerivative` and
   :py:func:`pylops.Laplacian`
-* Added `chirp-sliding` kind to :py:func:`pylops.waveeqprocessing.seismicinterpolation.SeismicInterpolation`
+* Added `chirp-sliding` kind to :py:class:`pylops.waveeqprocessing.seismicinterpolation.SeismicInterpolation`
 * Fixed bug due to the new internal structure of `LinearOperator` submodule introduced in `scipy1.8.0`