Help with NMDS plot customization

[jahid10]

Hi, Im relatively new in R and Vegan package. I have performed NMDS on my taxonomic abundance data which gave me a 2 dimensional scatter plots. However, I want to have an oval shade on points where they are clusturing, like the image below. How Can I modify the code below to make my graph have this representation (given that I already plotted up to the scatter graph

[jarioksa]

Check vegan function ordiellipse(..., kind = "ehull") that draws enclosing ellipses.