We simulate shaking in Tacoma, Washington, and surrounding areas from 6.5 and 7.0 earthquakes on the Tacoma fault. Ground motions are directly modeled up to 2.5 Hz using kinematic, finite‐fault sources; a 3D seismic velocity model considering regional geology; and a model mesh with 30 m sampling at the ground surface. In addition, we explore how adjustments to the seismic velocity model affect predicted shaking over a range of periods. These adjustments include the addition of a region‐specific geotechnical gradient, surface topography, and a fault damage zone. We find that the simulated shaking tends to be near estimates from empirical ground‐motion models (GMMs). However, long‐period (T = 5.0 s) shaking within the Tacoma basin is typically underpredicted by the GMMs. The fit between simulated and GMM‐derived short‐period (T = 0.5 s) shaking is significantly improved with the addition of the geotechnical gradient. From comparing different 6.5 earthquake scenarios, we also find that the response of the Tacoma basin is sensitive to the azimuth of incoming seismic waves. In adding surface topography to the simulation, we find that average ground motion is similar to that produced from the nontopography model. However, shaking is often amplified at topographic highs and deamplified at topographic lows, and the wavefield undergoes extensive scattering. Adding a fault damage zone has the effect of amplifying short‐period shaking adjacent to the fault, while reducing far‐field shaking. Intermediate‐period shaking is amplified within the Tacoma basin, likely due to enhanced surface‐wave generation attributable to the fault damage zone waveguide. When applied in the same model, the topography and fault damage zone adjustments often enhance or reduce the effects of one another, adding further complexity to the wavefield. These results emphasize the importance of improving near‐surface velocity model resolution as waveform simulations progress toward higher frequencies.