diff --git a/src/astral/sun.py b/src/astral/sun.py
index cee4440..ded32af 100644
--- a/src/astral/sun.py
+++ b/src/astral/sun.py
@@ -376,9 +376,6 @@ def time_of_transit(
     timeDiff = 4.0 * delta
     timeUTC = 720 + timeDiff - eq_of_time(t)
 
-    if timeUTC < 0:
-        raise ValueError(f"Sun never transits at a zenith of {zenith} on {date}")
-
     td = minutes_to_timedelta(timeUTC)
     dt = datetime.datetime(date.year, date.month, date.day) + td
     dt = pytz.utc.localize(dt)  # pylint: disable=E1120
diff --git a/src/test/test_wellington.py b/src/test/test_wellington.py
new file mode 100644
index 0000000..1dab1b6
--- /dev/null
+++ b/src/test/test_wellington.py
@@ -0,0 +1,17 @@
+import datetime
+
+import pytz
+from astral.geocoder import database, lookup
+from astral.sun import sun
+
+from almost_equal import datetime_almost_equal
+
+
+def test_Wellington():
+    wellington = lookup("Wellington", database())
+    dt = datetime.date(2020, 2, 11)
+    tz = pytz.timezone(wellington.timezone)
+    s = sun(wellington.observer, dt, tzinfo=tz)
+    assert datetime_almost_equal(
+        s["sunrise"], tz.localize(datetime.datetime(2020, 2, 11, 6, 38, 42))
+    )