This thesis covers recent advances in the adjoint variable method for the sensitivity estimations through time-domain electromagnetic simulations. It considers both frequency-independent and frequency-dependent response functions, and at the same time, provides a novel adjoint treatment for addressing dispersive sensitivity parameters in the material constitutive relation. With this proposed adjoint technique, response sensitivities with respect to all N sensitivity parameters can be computed through at most one extra simulations regardless of the value of N. This thesis also extends the existing adjoint technique to estimate all N^2 second-order sensitivity entries in the response Hessian matrix through N additional simulations. All adjoint sensitivity techniques presented in this thesis are numerically validated through various practical examples. Comparison shows that our produced adjoint results agree with those produced through central finite-difference approximations or through exact analytical approaches. / Dissertation / Doctor of Engineering (DEng)