diff --git a/docs/changes/39.bugfix.rst b/docs/changes/39.bugfix.rst
new file mode 100644
index 0000000..9a98182
--- /dev/null
+++ b/docs/changes/39.bugfix.rst
@@ -0,0 +1,3 @@
+- Fix gridding in `pyvisgen.simulation.observation.Observation` methods `create_rd_grid` and `create_lm_grid`
+  methods resulting in rotated images
+- Fix `pyvisgen.simulation.observation.ValidBaselineSubset` dataclass field order
diff --git a/docs/changes/39.feature.rst b/docs/changes/39.feature.rst
new file mode 100644
index 0000000..d2cda72
--- /dev/null
+++ b/docs/changes/39.feature.rst
@@ -0,0 +1,4 @@
+- Add class `Polarisation` to `pyvisgen.simulation.visibility` that is called in `vis_loop`
+  - Add linear, circular, and no polarisation options
+- Update `pyvisgen.simulation.visibility.Visibilities` dataclass to also store polarisation degree tensors
+- Add keyword arguments for polarisation simulation to `pyvisgen.simulation.observation.Observation` class
diff --git a/docs/changes/39.maintenance.rst b/docs/changes/39.maintenance.rst
new file mode 100644
index 0000000..83a2365
--- /dev/null
+++ b/docs/changes/39.maintenance.rst
@@ -0,0 +1,6 @@
+- Change pyvisgen.simulation.visibility.Visibilities dataclass component names from stokes components (I , Q, U, and V) 
+  to visibilities constructed from the stokes components (`V_11`, `V_22`, `V_12`, `V_21`)
+- Change indices for stokes components according to AIPS Memo 114
+  - Indices will be set automatically depending on simulated polarisation
+- Update comment strings in FITS files
+- Update docstrings accordingly in `pyvisgen.simulation.visibility.vis_loop` and `pyvisgen.simulation.observation.Observation`