Geostatistics dataset: finding the closest ten observations at a new coordinate

I have a geostatistics dataset Y(x,y) where x and y are the coordinates and I want to find the 10 nearest observations to the new location (x* ,y*).

Is there any available package in R that gives me the 10 nearest observations. my data look like:

 x        y             Y(s)
 0.52916369 0.7428021  3.11069813
 0.98035458 0.4520996 -1.62532766
 0.09133259 0.3221473  2.28771148
 0.32759395 0.2892624  0.06164611