diff --git a/modules/util/plot.py b/modules/util/plot.py
index f970a36cdfb3620f0e634471a7af23d4938e78ee..e9b959768f3c671b053c3d388b3c4959c17fb7af 100644
--- a/modules/util/plot.py
+++ b/modules/util/plot.py
@@ -292,10 +292,12 @@ def make_icing_image(h5f, probs, ice_lons, ice_lats, clvrx_str_time, satellite,
             lats = ice_lats_vld[ice_intensity == 1]
             ax.scatter(lons, lats, s=500.0, marker='o', color='white', transform=ccrs.PlateCarree())
             ax.scatter(lons, lats, s=200.0, marker='o', color='black', transform=ccrs.PlateCarree())
+            ax.scatter(lons, lats, s=1000.0, fillstyle=None, marker='o', color='black', transform=ccrs.PlateCarree())
             lons = ice_lons_vld[ice_intensity > 1]
             lats = ice_lats_vld[ice_intensity > 1]
             ax.scatter(lons, lats, s=500.0, marker='^', color='white', transform=ccrs.PlateCarree())
             ax.scatter(lons, lats, s=200.0, marker='^', color='black', transform=ccrs.PlateCarree())
+            ax.scatter(lons, lats, s=1000.0, fillstyle=None, marker='o', color='black', transform=ccrs.PlateCarree())
             lons = ice_lons_vld[ice_intensity == -1]
             lats = ice_lats_vld[ice_intensity == -1]
             # ax.scatter(lons, lats, s=500.0, marker='o', color='white', transform=ccrs.PlateCarree())