We analyze the prospects of employing a distributed global network of precision measurement devices as a dark matter and exotic physics observatory. In particular, we consider the atomic clocks of the Global Positioning System (GPS), consisting of a constellation of 32 medium-Earth orbit satellites equipped with either Cs or Rb microwave clocks and a number of Earth-based receiver stations, some of which employ highly-stable H-maser atomic clocks. High-accuracy timing data is available for almost two decades. By analyzing the satellite and terrestrial atomic clock data, it is possible to search for transient signatures of exotic physics, such as "clumpy" dark matter and dark energy, effectively transforming the GPS constellation into a 50,000km aperture sensor array. Here we characterize the noise of the GPS satellite atomic clocks, describe the search method based on Bayesian statistics, and test the method using simulated clock data. We present the projected discovery reach using our method, and demonstrate that it can surpass the existing constrains by several order of magnitude for certain models. Our method is not limited in scope to GPS or atomic clock networks, and can also be applied to other networks of precision measurement devices.

Cosmological observations indicate that 85% of all matter in the Universe is dark matter (DM), yet its microscopic composition remains a mystery. One hypothesis is that DM arises from ultralight quantum fields that form macroscopic objects such as topological defects. Here we use GPS as a ~ 50,000 km aperture DM detector to search for such defects in the form of domain walls. GPS navigation relies on precision timing signals furnished by atomic clocks hosted on board GPS satellites. As the Earth moves through the galactic DM halo, interactions with topological defects could cause atomic clock glitches that propagate through the GPS satellite constellation at galactic velocities ~ 300 km/s. Mining 16 years of archival GPS data, we find no evidence for DM in the form of domain walls at our current sensitivity level. This allows us to improve the limits on certain quadratic scalar couplings of domain wall DM to standard model particles by several orders of magnitude.