diff --git a/src/analysis/forecasted_storm_impacts.py b/src/analysis/forecasted_storm_impacts.py
index eca4b80..51fb920 100644
--- a/src/analysis/forecasted_storm_impacts.py
+++ b/src/analysis/forecasted_storm_impacts.py
@@ -29,7 +29,9 @@ def forecasted_impacts(df_profile_features, df_forecasted_twl):
 
     # Join with df_profile features to find dune toe and crest elevations
     df_forecasted_impacts = df_forecasted_impacts.merge(
-        df_profile_features.query("profile_type=='prestorm'")[["dune_toe_z", "dune_crest_z"]], how="left", left_index=True, right_index=True
+        df_profile_features.query("profile_type=='prestorm'").reset_index('profile_type')[["dune_toe_z",
+                                                                                         "dune_crest_z"]], how="left",
+        left_on=['site_id'], right_on=['site_id']
     )
 
     # Compare R_high and R_low wirth dune crest and toe elevations