code: def calculate_equinoxes(self, year, timezone='UTC'):
          
  
          tz = pytz.timezone(timezone)
  
          d1 = ephem.next_equinox(str(year))
          d = ephem.Date(str(d1))
          equinox1 = d.datetime() + tz.utcoffset(d.datetime())
  
          d2 = ephem.next_equinox(d1)
          d = ephem.Date(str(d2))
          equinox2 = d.datetime() + tz.utcoffset(d.datetime())
  
          return (equinox1.date(), equinox2.date())
docs: calculate equinox with time zone