The Specified Dynamics version of the Whole Atmosphere Community Climate Model (SD-WACCM) and the Goddard Space Flight Center two-dimensional (GSFC 2-D) models are used to investigate the effect of galactic cosmic rays (GCRs) on the atmosphere over the 1960–2010 time period. The Nowcast of Atmospheric Ionizing Radiation for Aviation Safety (NAIRAS) computation of the GCRcaused ionization rates are used in these simulations. GCRcaused maximum NOx increases of 4–15 % are computed in the Southern polar troposphere with associated ozone increases of 1–2 %. NOx increases of ∼ 1–6 % are calculated for the lower stratosphere with associated ozone decreases of 0.2–1 %. The primary impact of GCRs on ozone was due to their production of NOx . The impact of GCRs varies with the atmospheric chlorine loading, sulfate aerosol loading, and solar cycle variation. Because of the interference between the NOx and ClOx ozone loss cycles (e.g., the ClO + NO2 + M → ClONO2 + M reaction) and the change in the importance of ClOx in the ozone budget, GCRs cause larger atmospheric impacts with less chlorine loading. GCRs also cause larger atmospheric impacts with less sulfate aerosol loading and for years closer to solar minimum. GCR-caused decreases of annual average global total ozone (AAGTO) were computed to be 0.2 % or less with GCR-caused column ozone increases between 1000 and 100 hPa of 0.08 % or less and GCR-caused column ozone decreases between 100 and 1 hPa of 0.23 % or less. Although these computed ozone impacts are small, GCRs provide a natural influence on ozone and need to be quantified over long time periods. This result serves as a lower limit because of the use of the ionization model NAIRAS/HZETRN which underestimates the ion production by neglecting electromagnetic and muon branches of the cosmic ray induced cascade. This will be corrected in future works.