diff --git a/web/content/docs/benchmarks/thermo-richards-mechanics/DECOVALEX2023-TaskC.md b/web/content/docs/benchmarks/thermo-richards-mechanics/DECOVALEX2023-TaskC.md
new file mode 100644
index 0000000000000000000000000000000000000000..bd9a0010784b079292441f59cb8965b44047b0ed
--- /dev/null
+++ b/web/content/docs/benchmarks/thermo-richards-mechanics/DECOVALEX2023-TaskC.md
@@ -0,0 +1,54 @@
++++
+project = "ThermoRichardsMechanics/TaskCDECOVALEX2023/Decovalex-0.prj"
+
+title = "A test based on DECOVALEX2023 Task C for water vapour diffusion model for non-isothermal Richards flow"
+date = "2022-02-07T13:19:49+01:00"
+author = "Wenqing Wang and Sonja Kaiser"
+
+weight = 70
+
+[menu]
+  [menu.benchmarks]
+    parent = "thermo-richards-mechanics"
+
++++
+
+{{< data-link >}}
+
+This test is prepared exactly according to the specifications of Step 0c of
+ [Task C of
+ the DECOVALEX 2023 project](https://decovalex.org/D-2023/task-c.html),
+ which is aimed to simulate the coupled THM processes
+ in the full-scale emplacement experiment (FE experiment) at the Mont Terri
+ Underground Rock Laboratory [[1]](#1).
+
+The description of the Step 0c can be found in the task specifications.
+
+Running the entire test takes long time (three hours on a PC with i7-8565U CPU).
+ As a test for code, it is limited to run only 20 steps.
+
+The test is used to mainly verify the implementation of water vapour diffusion model,
+ heat latency model, and equation wise residuum compensation with non equilibrium
+ initial state, and etc.
+
+ The results obtained are comparable with that obtained by other teams in
+ Task C. The following figure shows the distribution of temperature in the domain,
+ water saturation in the vicinity of the heater, and displacement magnitude
+ in the domain after nearly 3 years' heating:
+{{< img src="../decovalex_2023_c.png" >}}
+
+The following two figures show the temporal variations of temperature and water
+ saturation, respectively, at a node near the heater:
+
+<img src="../decovalex_2023_c_T_t.png" alt="drawing" width="450"/>
+<img src="../decovalex_2023_c_S_t.png" alt="drawing" width="450"/>
+
+As shown the water saturation variation curve, the de-saturation -
+ re-saturation process is well captured by the numerical simulation.
+
+## References
+<a id="1">[1]</a>
+Müller, H.R., Garitte, B., Vogt, T. et al. Implementation of the full-scale
+ emplacement (FE) experiment at the Mont Terri rock laboratory.
+ Swiss J Geosci 110, 287–306 (2017). DOI:
+[10.1007/s00015-016-0251-2](https://sjg.springeropen.com/articles/10.1007/s00015-016-0251-2)
diff --git a/web/content/docs/benchmarks/thermo-richards-mechanics/decovalex_2023_c.png b/web/content/docs/benchmarks/thermo-richards-mechanics/decovalex_2023_c.png
new file mode 100644
index 0000000000000000000000000000000000000000..d87a8c8c9bb729a29ee0a4c3c8dee8c4df9cb1a8
Binary files /dev/null and b/web/content/docs/benchmarks/thermo-richards-mechanics/decovalex_2023_c.png differ
diff --git a/web/content/docs/benchmarks/thermo-richards-mechanics/decovalex_2023_c_S_t.png b/web/content/docs/benchmarks/thermo-richards-mechanics/decovalex_2023_c_S_t.png
new file mode 100644
index 0000000000000000000000000000000000000000..01c12cb8946e07be587801db212323c4bae999f7
Binary files /dev/null and b/web/content/docs/benchmarks/thermo-richards-mechanics/decovalex_2023_c_S_t.png differ
diff --git a/web/content/docs/benchmarks/thermo-richards-mechanics/decovalex_2023_c_T_t.png b/web/content/docs/benchmarks/thermo-richards-mechanics/decovalex_2023_c_T_t.png
new file mode 100644
index 0000000000000000000000000000000000000000..5a80ee143b1960cabcf7ec61fde6290b59685233
Binary files /dev/null and b/web/content/docs/benchmarks/thermo-richards-mechanics/decovalex_2023_c_T_t.png differ