We advocate for the idea that there is a fundamentally new mechanism for axion production on Earth, as recently suggested in [1,2]. We specifically focus on production of axions within Earth, with low velocities such that they will be trapped in the gravitational field. Our computations are based on the so-called Axion Quark Nugget (AQN) dark matter model, which was originally invented to explain the similarity of the dark and visible cosmological matter densities. This occurs in the model irrespective of the axion mass ma or initial misalignment angle θ0. Annihilation of antimatter AQNs with visible matter inevitably produce axions when AQNs hit Earth. The emission rate of axions with velocities below escape velocity is very tiny compared to the overall emission, however these axions will be accumulated over the 4.5 billion year life time of the Earth, which greatly enhances the discovery potential. We perform numerical simulations with a realistically modeled incoming AQN velocity and mass distribution, and explore how AQNs interact as they travel through the interior of the Earth. We use this to estimate the axion flux on the surface of the Earth, the velocity-spectral features of trapped axions, the typical annihilation pattern of AQN, and the density profile of the axion halo around the Earth. Knowledge of these properties is necessary to make predictions for the observability of trapped axions using CAST, ADMX, MADMAX, CULTASK, ORPHEUS, ARIADNE, CASPEr, ABRACADABRA, QUAX, ORGAN, TOORAD, DM Radio.2 In particular, a first-order phase transition is not required as the axion domain wall plays the role of the squeezer. Another problem with [34,35] is that nuggets will likely evaporate on a Hubble time-scale even. For the AQN model this argument is not applicable because the vacuum-ground-state energies inside (color-superconducting phase) and outside (hadronic phase) the nugget are drastically different. Therefore, these two systems can arXiv:1905.00022v2 [astro-ph.CO]