diff --git a/CHANGELOG.md b/CHANGELOG.md
new file mode 100644
index 0000000..36a353c
--- /dev/null
+++ b/CHANGELOG.md
@@ -0,0 +1,105 @@
+# Changelog
+All notable changes to this project will be documented in this file.
+
+The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/),
+and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html).
+
+## Unreleased
+### Changed
+- `CXXFLAGS` are now passed to the linker
+
+## v2.2.1 -- 2021-03-02
+### Added
+- Output registered fields and dumpers in `print(model)`
+- Added `operator[]` to the C++ model class (for fields)
+- Added `traction` and `displacement` properties to Python model bindings
+- Added `operators` property to Python model bindings, which provides a
+  dict-like access to registered operators
+- Added `shape` and `spectrum` to properties to Python surface generator
+  bindings
+- Surface generator constructor accepts surface global shape as argument
+- Choice of FFTW thread model
+### Changed
+- Tests use `/tmp` for temporary files
+- Updated dependency versions (Thrust, Pybind11)
+### Deprecated
+- Most `get___()` and `set___()` in Python bindings have been deprecated. They
+  will generate a `DeprecationWarning`.
+### Removed
+- All legacy code
+
+## v2.2.0 -- 2020-12-31
+### Added
+- More accurate function for computation of contact area
+- Function to compute deviatoric of tensor fields
+- MPI implementation
+- Convenience `hdf5toVTK` function
+- Readonly properties `shape`, `global_shape`, `boundary_shape` on model to give
+  shape information
+### Changed
+- Preprocessor defined macros are prefixed with `TAMAAS_`
+- Moved `tamaas.to_voigt` to `tamaas.compute.to_voigt`
+### Fixed
+- Warning about deprecated constructors with recent GCC versions
+- Wrong computation of grid strides
+- Wrong computation of grid sizes in views
+
+## v2.1.4 -- 2020-08-07
+### Added
+- Possibility to generate a static `libTamaas`
+- C++ implementation of DFSANE solver
+- Allowing compilation without OpenMP
+### Changed
+- NetCDF dumper writes frames to a single file
+### Fixed
+- Compatibility with SCons+Python 3
+
+## v2.1.3 -- 2020-07-27
+### Added
+- Version number to `TamaasInfo`
+### Changed
+- Prepending root directory when generating archive
+
+## v2.1.2 -- 2020-07-24
+This release changes some core internals related to discrete Fourier transforms
+for future MPI support.
+### Added
+- Caching `CXXFLAGS` in SCons build
+- SCons shortcut to create code archive
+- Test of the elastic-plastic contact solver
+- Paraview data dumper (`.pvd` files)
+- Compression for UVW dumper
+- `__contains__` and `__iter__` Python bindings of model
+- Warning message of possible overflow in Kelvin
+### Changed
+- Simplified `tamaas_info.cpp`, particularly the diff part
+- Using a new class `FFTEngine` to manage discrete Fourier transforms. Plans are
+  re-used as much as possible with different data with the same shape. This is
+  in view of future MPI developments
+- Redirecting I/O streams in solve functions so they can be used from Python
+  (e.g. in Jupyter notebooks)
+- Calling `initialize()` and `finalize()` is no longer necessary
+### Fixed
+- Convergence issue with non-linear solvers
+- Memory error in volume potentials
+
+## v2.1.1 -- 2020-04-22
+### Added
+- SCons shortcut to run tests
+### Fixed
+- Correct `RPATH` for shared libraries
+- Issues with SCons commands introduced in v2.1.0
+- Tests with Python 2.7
+
+## v2.1.0 -- 2020-04-17
+### Added
+- SCons shortcuts to build/install Tamaas and its components
+- Selection of integration method for Kelvin operator
+- Compilation option to remove the legacy part of Tamaas
+- NetCDF dumper
+### Fixed
+- Link bug with clang
+- NaNs in Kato saturated solver
+
+## v2.0.0 -- 2019-11-11
+First public release. Contains relatively mature elastic-plastic contact code.