You can build on the axis of the joint with jointplot as follows:
a=sns.jointplot(x="volatility", y="return", data=data, size=10) a.ax_joint.plot([0.2],[0.2],'ro')
As an aside, you can also access the x and y subtitles with
a.ax_marg_x a.ax_marg_y
source share