Dask of june breaks xclim
The most recent dask release deprecated (with no warning) passing both `meta` and `output_dtypes` to `apply_gfunc` which is itself called by `xarray.apply_ufunc` when `dask='parallelized'`.

The fix is simply to remove one of them, dask recommends using `meta`. AFAIK, this will be retrocompatible.
diff --git a/xclim/indices/fwi.py b/xclim/indices/fwi.py
index 92e166d1..f4435a64 100644
--- a/xclim/indices/fwi.py
+++ b/xclim/indices/fwi.py
@@ -1118,7 +1118,6 @@ def fire_weather_ufunc(
         input_core_dims=input_core_dims,
         output_core_dims=output_core_dims,
         dask="parallelized",
-        output_dtypes=output_dtypes,
         dask_gufunc_kwargs={
             "meta": tuple(np.array((), dtype=dtype) for dtype in output_dtypes)
         },
diff --git a/xclim/sdba/utils.py b/xclim/sdba/utils.py
index 78879806..ee12c190 100644
--- a/xclim/sdba/utils.py
+++ b/xclim/sdba/utils.py
@@ -656,7 +656,6 @@ def get_clusters(data: xr.DataArray, u1, u2, dim: str = "time"):
         input_core_dims=[[dim], [], []],
         output_core_dims=[["cluster"], ["cluster"], ["cluster"], ["cluster"], []],
         kwargs={"N": N},
-        output_dtypes=[int, int, int, data.dtype, int],
         dask="parallelized",
         vectorize=True,
         dask_gufunc_kwargs={
['xclim/indices/fwi.py', 'xclim/sdba/utils.py']