diff --git a/geomagio/netcdf/NetCDFFactory.py b/geomagio/netcdf/NetCDFFactory.py
index 920c9d8c276384e61aa5f49f77879574c831fa23..590e1eadbc633a0f4156c889f09058190d050ee0 100644
--- a/geomagio/netcdf/NetCDFFactory.py
+++ b/geomagio/netcdf/NetCDFFactory.py
@@ -198,7 +198,7 @@ class NetCDFFactory(TimeseriesFactory):
             # Create an empty NetCDF structure with minimal CF compliance
             with netCDF4.Dataset(fh, "w", format="NETCDF4") as nc_dataset:
                 nc_dataset.title = "Geomagnetic Time Series Data"
-                nc_dataset.Conventions = "CF-1.4"
+                nc_dataset.Conventions = "CF-1.6"
             return
 
         timeseries.merge()
@@ -214,7 +214,7 @@ class NetCDFFactory(TimeseriesFactory):
             print("No matching channels found after filtering.")
             with netCDF4.Dataset(fh, "w", format="NETCDF4") as nc_dataset:
                 nc_dataset.title = "Geomagnetic Time Series Data"
-                nc_dataset.Conventions = "CF-1.4"
+                nc_dataset.Conventions = "CF-1.6"
             return
 
         timeseries.sort(keys=["starttime"])
@@ -237,7 +237,7 @@ class NetCDFFactory(TimeseriesFactory):
         # Create NetCDF dataset with CF conventions
         with netCDF4.Dataset(fh, "w", format="NETCDF4") as nc_dataset:
             # Define global attributes with CF compliance
-            nc_dataset.Conventions = "CF-1.4"
+            nc_dataset.Conventions = "CF-1.6"
             nc_dataset.title = "Geomagnetic Time Series Data"
             nc_dataset.institution = getattr(stats, "agency_name", "")
             # nc_dataset.source = getattr(stats, "source", "")