The problem with my derivation is I am not sure if it correct and I do not have a way to check it.
Try dimensional analysis to check your solution. Is your solution consistent with dimensional analysis? Estimate a solution to a real example, based upon you knowledge of geophysics and geology. Does your solution even come close to a reasonable estimate? If you are unfamiliar with dimensional analysis and estimation methods, I recommend the online book (MIT course) Street Fighting Mathematics for some possible ways to stimulate ideas of ways to check your derivation.
For example, as you have written your solution, the $b$ and the $\alpha$ are constant over time ($t$) and location ($x)$. First substitute your solution into the differential equation, and carry out the differentiation to confirm that it is in fact a solution. If it can be verified as a solution, turn your attention to the dimensions of b and alpha. The dimension of b must be time and the dimension of alpha must be length. You are modeling the system as a viscoelastic solid, so the the initial condition ($t=0$) for your displacements $w(0,x)$ are just the elastic solution, and the time-dependent term is unity. Use the elastic solution to solve for alpha (the characteristic wavelength) because this does not change over time. Name the purely elastic displacement '$e$' and it has dimension of length.
You know you want to include viscosity somewhere, so let's try the kinematic viscosity $v$ (dimensions $M*M/s$). So try $b = -(L * L)/v$, or $b = -(L * 1)/v$, because this is a 2D problem and can be one unit deep. You'll need the negative sign to make the displacement grow faster for a smaller viscosity.
What could $L$ be? I would try $L=e$, or $L=e/2$. Because you really don't know the kinematic viscosity of the lithosphere, the factor is not critical.
Substitute $b$ and $\alpha$ into your solution. Does the result seem to make physical and geological sense?