The development and initial application of a mechanistic model (TOUGHREACT‐N) designed to characterize soil nitrogen (N) cycling and losses are described. The model couples advective and diffusive nutrient transport, multiple microbial biomass dynamics, and equilibrium and kinetic chemical reactions. TOUGHREACT‐N was calibrated and tested against field measurements to assess pathways of N loss as either gas emission or solute leachate following fertilization and irrigation in a Central Valley, California, agricultural field as functions of fertilizer application rate and depth, and irrigation water volume. Our results, relative to the period before plants emerge, show that an increase in fertilizer rate produced a nonlinear response in terms of N losses. An increase of irrigation volume produced NO₂⁻ and NO₃⁻ leaching, whereas an increase in fertilization depth mainly increased leaching of all N solutes. In addition, nitrifying bacteria largely increased in mass with increasing fertilizer rate. Increases in water application caused nitrifiers and denitrifiers to decrease and increase their mass, respectively, while nitrifiers and denitrifiers reversed their spatial stratification when fertilizer was applied below 15 cm depth. Coupling aqueous advection and diffusion, and gaseous diffusion with biological processes, closely captured actual conditions and, in the system explored here, significantly clarified interpretation of field measurements.