Large model domains are making it so the smaller ones aren't always clickable - need to sort the features on the map so that the largest are the lowest layers. May be that sorting by area once loaded (with largest at top) will solve.
after loading:
all_df['area'] = all_df.area
all_df = all_df.sort_values(by=['area'],ascending=False)