I wanted to push the limits and see when NN subtraction performance starts to break by changing the number of seismometers and the size of the array. For aLIGO, 10 seismometers in a doubly-wound spiral around the test mass with outer radius 8m is definitely ok. Only if I simulate a seismic field that is stronger by a factor 20 than the 90 percentile curve observed at LHO does it start to get problematic. The subtraction residuals in this case look like

The 20 seismometer spiral is still good, but the 10 seismometer spiral does not work anymore. It gets even worse when you consider arrays with circular shape (and one seismometer at the center near the test mass):

This result is in agreement with previous results that circular arrays have trouble in general to subtract NN from locally generated seismic waves or seismic transients (wavelets).

I should emphasize that the basic assumption is that I know what the minimum seismic wavelength is. Currently I associate the minimum wavelength with a Rayleigh overtone, but scattering could make a difference. It is possible that there are scattered waves with significantly smaller wavelength.