rot[x_, y_] = D[Sin[y + x^3], x] - D[Cos[x + y^3], y];
StreamDensityPlot[{{Cos[x + y^3], Sin[y + x^3]}, Log[Abs[rot[x, y]] + 1]}, {x, -2, 2}, {y, -2, 2}, StreamPoints -> Fine, ColorFunction -> "RedBlueTones", MaxRecursion -> 4, LightingAngle -> Automatic, StreamStyle -> Black]