如何使用xarray来计算每个网格单元在不同时间纬度上的趋势？

大家好，本人第一次接触xarray还不太熟练，想请教一下各位如何使用python来计算每个网格单元在不同时间纬度上的趋势？
目前的数据组中有一下四组数据：

Dimensions:    (latitude: 721, longitude: 1440, time: 72)
Coordinates:
  * time       (time) datetime64[ns] 2000-01-01 2001-01-01 ... 2021-01-01
  * longitude  (longitude) float32 0.0 0.25 0.5 0.75 ... 359.25 359.5 359.75
  * latitude   (latitude) float32 90.0 89.75 89.5 89.25 ... -89.5 -89.75 -90.0
Data variables:
    pH          (time, latitude, longitude) float32 49654.793 49654.793 ... 49654.793

我目前没有搞明白xr.DataArray应该如何使用。我输入了一下code

import xarray as xr
da = xr.DataArray(data=pH, coords={"lat": (["x","y"], lat),
                          "lon": (["x","y"], lon),
                          "time":(["x","y"], time)},
                  dims=["x","y"])

但是系统显示无法将time中的年份转换成variable，想请教一下各位该如何修改这一部分？十分感谢！