ли способ добавить вторичную легенду к диаграмме рассеяния, где размер разброса пропорционален некоторым данным?

Я написал следующий код, который генерирует график рассеяния. Цвет разброса представляет год (и берется из определяемого пользователем df), в то время как размер разброса представляет переменную 3 (также взятую из df, но это необработанные данные):

import pandas as pd 

colors = pd.DataFrame({'1985':'red','1990':'b','1995':'k','2000':'g','2005':'m','2010':'y'}, index=[0,1,2,3,4,5])

fig = plt.figure()
ax = fig.add_subplot(111)

for i in df.keys():
    df[i].plot(kind='scatter',x='variable1',y='variable2',ax=ax,label=i,s=df[i]['variable3']/100, c=colors[i])

ax.legend(loc='upper right')
ax.set_xlabel("Variable 1")
ax.set_ylabel("Variable 2")

Этот код (с моими данными) создает следующий график:

Таким образом, хотя цвета / годы четко и четко определены, размер разброса - нет.

Как добавить вторичную или дополнительную легенду, определяющую, что означает размер разброса?