It cannot happen in a continuous time 2D system, simply due to uniqueness of ODE property. At least (2+1)-D is needed, i.e. this phenomenon can be seen in 2D maps derived from taking time-T sections of a time-dependent 2D ODE system.
A good resource is this document by Wiggins: http://docenti.lett.unisi.it/files/4/2/13/1/Wiggins_CHAOS.pdf
Also take a look at Wiggin's "Global bifurcations and chaos".
Other keywords to search for : "Melnikov method", "conley-moser conditions".