¿Cómo agregar datos geográficamente segregados por fecha en Python?

Estoy trabajando con datos de series de tiempo que se desglosan por varios niveles geográficos (provincia y región de salud). Los datos se desglosan aún más por grupos de edad. Me gustaría agregar los datos primero hasta el nivel de región de salud y luego hasta el nivel de provincia, dependiendo de la necesidad.

Datos de ejemplo:

df = pd.DataFrame({'Date':['1/1/2022','1/1/2022','1/1/2022','1/1/2022','1/1/2022','1/1/2022','1/1/2022','1/1/2022','1/1/2022'],
             'province':[35,35,35,35,35,35,35,35,35],
             'health region':[1,1,1,2,2,2,3,3,3],
              'age group':[1,2,3,1,2,3,1,2,3],'cases':[6,1,9,7,9,0,4,2,2]})

Salida deseada al agregar hasta el nivel de región de salud:

df_hr = pd.DataFrame({'Date':['1/1/2022','1/1/2022','1/1/2022'],
                     'province':[35,35,35],
                     'health region':[1,2,3],'cases':[16,16,8]})

Cuando uso el siguiente código:

df = df.groupby('health region').sum()

Pierdo las fechas.

Cuando intento:

df =  df.groupby(['health region','Date']).sum()

o

df =  df.groupby(['health region',df['Date'].dt.date]).sum()

Obtengo un error ValueError: mixed datetimes and integers in passed array

¿Hay alguna manera fácil de hacer esto? Estaba pensando en usar un bucle para dividir los datos por región de salud, guardar las fechas únicas, agregar, fusionar las fechas nuevamente y luego juntar las regiones sanitarias nuevamente. Pero preferiría no hacerlo si hay una manera fácil de hacerlo.

Gracias,
i.