diff --git a/docs/docs/assets/datasets/oisst.png b/docs/docs/assets/datasets/oisst.png
new file mode 100644
index 00000000..b81c9faa
Binary files /dev/null and b/docs/docs/assets/datasets/oisst.png differ
diff --git a/docs/docs/icechunk-python/virtual.md b/docs/docs/icechunk-python/virtual.md
index 9abc2f18..b9b30871 100644
--- a/docs/docs/icechunk-python/virtual.md
+++ b/docs/docs/icechunk-python/virtual.md
@@ -148,4 +148,12 @@ ds = xr.open_zarr(
 #    err      (time, zlev, lat, lon) float64 257MB dask.array<chunksize=(1, 1, 720, 1440), meta=np.ndarray>
 ```
 
-Success! We have created our full dataset with 31 timesteps spanning the month of august, all with virtual references to pre-existing data files in object store. This means we can now version control our dataset, allowing us to update it, and roll it back to a previous version without copying or moving any data from the original files. 
\ No newline at end of file
+Success! We have created our full dataset with 31 timesteps spanning the month of august, all with virtual references to pre-existing data files in object store. This means we can now version control our dataset, allowing us to update it, and roll it back to a previous version without copying or moving any data from the original files. 
+
+Finally, let's make a plot of the sea surface temperature!
+
+```python
+ds.sst.isel(time=26, zlev=0).plot(x='lon', y='lat', vmin=0)
+```
+
+![oisst](../assets/datasets/oisst.png)
\ No newline at end of file
diff --git a/docs/docs/sample-datasets.md b/docs/docs/sample-datasets.md
index 5021c171..688c2ed6 100644
--- a/docs/docs/sample-datasets.md
+++ b/docs/docs/sample-datasets.md
@@ -9,7 +9,7 @@
 
 > The NOAA 1/4° Daily Optimum Interpolation Sea Surface Temperature (OISST) is a long term Climate Data Record that incorporates observations from different platforms (satellites, ships, buoys and Argo floats) into a regular global grid
 
-Checkout an example dataset built using all virtual references pointing to daily Sea Surface Temperature data from 2020 to 2024 on NOAA's S3 bucket using python:
+Check out an example dataset built using all virtual references pointing to daily Sea Surface Temperature data from 2020 to 2024 on NOAA's S3 bucket using python:
 
 ```python
 import icechunk
@@ -23,4 +23,6 @@ storage = icechunk.StorageConfig.s3_anonymous(
 store = IcechunkStore.open_existing(storage=storage, mode="r", config=StoreConfig(
     virtual_ref_config=VirtualRefConfig.s3_anonymous(region='us-east-1'),
 ))
-```
\ No newline at end of file
+```
+
+![oisst](./assets/datasets/oisst.png)
\ No newline at end of file