Organonitrates (ON) are key components of secondary organic aerosols (SOA) with potential environmental and climate effects. However, ON formation from limonene, a major monoterpene, and its sensitivity to oxidation pathways remain insufficiently explored due to the absence of models with explicit chemical mechanisms. This study advances the representation of limonene-derived ON formation across three oxidation pathways (O<sub>3</sub>, OH, NO<sub>3</sub>) into both a chemical box model and a global model. Box model sensitivity experiments revealed that competition among major oxidation pathways, coupled with the high yield of limonene-derived ON from O<sub>3</sub>-initiated oxidation, leads to increased limonene-derived ON production when the O3-initiated pathway is enhanced, whereas strengthening the OH- or NO<sub>3</sub>-initiated pathways reduces ON formation. Compared to the box model, the global simulation exhibits stronger nonlinear responses and great spatiotemporal variability in limonene-derived ON formation across different oxidation pathways. This is primarily driven by the complex distribution of precursors and oxidants, as well as changing in dominate chemical pathways under various meteorological conditions. In the presence of the other two pathways, increasing the O<sub>3</sub>- or NO<sub>3</sub>-initiated pathway reduces the global limonene-derived ON burden by 19.9 % and 17.3 %, respectively, whereas enhancing the OH-initiated pathway increases it by 44.7 %. The refined limonene-derived ON chemistry improves the global model's ability to simulate ON formation, validated against observations. This study also establishes a methodological framework based on explicit chemical mechanisms, providing a foundation for simulating SOA formation processes.