diff --git a/mlair/data_handler/data_handler_wrf_chem.py b/mlair/data_handler/data_handler_wrf_chem.py
index 56c90c6e52ffc9dd0596f00b237500a04941382e..e117effad33bcfd28effeb1eb85c78a81c128def 100644
--- a/mlair/data_handler/data_handler_wrf_chem.py
+++ b/mlair/data_handler/data_handler_wrf_chem.py
@@ -711,7 +711,12 @@ class SingleGridColumnWrfChemDataLoader(BaseWrfChemDataLoader):
     def __toarstats_aggregation(self, data, local_time_zone, meta, target_sampling, var):
         with TimeTracking(name=f"{self.station}: apply toarstats `{self.statistics_per_var[var]}({var})`"):
             spatial_dims = list(remove_items(data[var].dims, self.physical_t_coord_name))
-            df = data[var].to_dataframe()[[var]].reset_index(level=spatial_dims)
+            try:
+                df = data[var].to_dataframe()[[var]].reset_index(level=spatial_dims)
+            except Exception as e:
+                logging.info(f"can't reset_index in __toarstats_aggregation: {e}")
+                df = data[var].to_dataframe()[[var]].reset_index(level=spatial_dims)
+
             df = df[[var] + spatial_dims]
             df = self.set_external_time_zone_and_convert_to_local_time_zone(df, local_time_zone)
             df = df.groupby(spatial_dims)