diff --git a/3_visualize/src/plot_utils.R b/3_visualize/src/plot_utils.R
index 626d1aab1d56642b979e04e663e198011c6ff98a..e534ea190e85fcb1c2e2ef5be6be69f06938091d 100644
--- a/3_visualize/src/plot_utils.R
+++ b/3_visualize/src/plot_utils.R
@@ -129,7 +129,7 @@ plot_raster <- function(in_raster, conus_sf, conus_inner, outfile_path,
   showtext_auto(enable = TRUE)
   
   # doesn't work with cowplot::get_plot_component()
-  legend_label <- expression(bold("Population per 1 km"^2)) #bquote('Population per 1'~km^2)
+  legend_label <- expression(paste(bold("Population per 1 "), bold(km^2), bold(", 2020")))
   
   (pop_map <- ggplot() +
       #geom_spatraster(data = pop_usa_proj)+