You can use Interpolation, after converting your triples {x,y,z} into the structure {{x,y},z}:
f = Interpolation[root1 /. {x_Real, y_, z_} :> {{x, y}, z}]; Plot3D[f[x, y], {x, .1, .4}, {y, .1, 9}]
