Sometimes we get geom and sometimes we get geometry - we should select with dplyr::stars_with("geom") See https://github.com/BigelowLab/mpatools/blob/main/R/es50.R#L14
The MARINE variable has been numeric (0,1,2) in the past, but now we see TERRESTRIAL etc. Yikes! So, cast all to character and recode?