In this study, we examine the scaling of earthquake waiting time distributions in Northern Chile. For this purpose, we analyzed 7-yr high spatial resolution and low complete- ness magnitude IPOC seismic catalog and 43-yr USGS catalog. A unified earthquake moment-space-time scaling relation is empirically evaluated by calculating waiting times for different ranges of magnitude and epicentral area linear dimension and estimating the associated scaling coefficients, β analog to b-value and γ, the correlation fractal dimension. We find a scaling function that can be characterized with 3 distinct regions, regions whose behaviour depends on whether seismicity is in the coastal area or at intermediate depth. Moreover, high resolution localizations from the IPOC catalog allows us to further observe differences: the lower plane seismicity of the double seismic zone behaves just like that at intermediate depth , while the crustal upper plate and the upper plane from the double seismic zone seismicity behaves more like interface seismicity. Thus, the earthquake waiting time distribution primarily depends on whether seismicity is located in a crust-on-crust or in a crust-mantle system: having respectively high/low correlated behaviour in the short scale region, non-exponential/exponential decay in the transition middle region and, in all cases, long-term clustering with a slower than exponential decay in the long scale.