Differences in numpy.polyfit and stats.linregress on geospatial netcdf data

Question

I have monthly averaged precipitation data for the last 40 years across the entire globe. I want to find the trend in precip over that duration. There are two approaches (th