Time-varying primary magnetic fields generated outside Earth by the magnetospheric ring current induce electrical currents in Earth's interior, which give rise to secondary magnetic fields with a complementary geometry. Geomagnetic depth sounding involves the analysis of magnetic field data to compute frequency-dependent response functions which yield information about the electrical conductivity of Earth's interior. I explore methods and results of forward-modeling global electromagnetic induction under a variety of assumptions about Earth conductivity and the spatial structure of the primary field. I begin by developing computational tools to perform time- and frequency-domain simulations of global induction in models with arbitrary conductivity and primary field structure using FlexPDE, a general-purpose software package that employs the finite-element method to solve partial differential equations. The method is shown to produce solutions with better than 1% accuracy when the simulated fields and response functions are compared to analytic solutions for a variety of problems in electromagnetic induction, and to qualitatively reproduce fields and response functions measured by satellites and observatories. The technique is employed in combination with analytic methods to explore the effect on the response of Earth models to primary fields with asymmetric structure. Standard methods of producing response functions from scalar and vector magnetic data are compared, and scalar methods are found to generate responses with significantly greater spatial bias for models with non-zonal fields. I develop the mathematical formulation for including Earth-rotation in the forward models, and use it to calculate frequency-dependent estimates of the amount of non-zonal structure required to produce previously reported local-time bias in empirical satellite response functions. Because it is difficult to validate solutions to induction problems that lack analytic solutions, we participate in an ongoing project with other researchers who simulate the global induction problem with different methods. We compare the synthetic fields calculated with the FlexPDE method to those calculated with the integral equation method and with the time-domain spectral method for a variety of conductivity models