The response equations governing nonequilibrium solvation in the continuum approach are derived and implemented for second order in response theory. Applications to a simple benchmark system consisting of a solvated water molecule are presented. We consider the role of nonequilibrium solvent configurations on the first hyperpolarizability. Effects of solvation parameters, such as optical and static dielectric constants, and frequency dispersion are discussed.