The issue was caused by a bug in the GEOS version used in Shapely, as pointed out by RuthC. Installing conda with the latest versions of all packages solved the problem!