Paleomagnetic poles, or paleopoles, quantify the past position of rocks relative to the geomagnetic pole and constrain tectonic reconstructions and apparent polar wander paths (APWPs; e.g., Besse & Courtillot, 2002;Torsvik et al., 2012). The calculation of paleopoles relies on the assumption that the time-averaged geomagnetic field approximates a geocentric axial dipole (GAD), but is complicated by short-term deviations from this field (e.g., Cromwell et al., 2018; Oliveira et al., 2021) known as paleosecular variation (PSV). To obtain a paleopole, paleomagnetists therefore average virtual geomagnetic poles (VGPs), whereby every VGP is then assumed a "spot reading": an instantaneous reading of the past geomagnetic field collected from a rock unit ("site") that represents an increment of geological time, such as a lava flow (Butler, 1992;Tauxe et al., 2010). However, not every VGP represents an accurate spot reading because artifacts may be introduced by measuring errors or remagnetization (e.g., Bilardello et al., 2018;Butler, 1992;Irving, 1961Irving, , 1964. Therefore, the paleomagnetic community commonly uses a set of data filters to acquire a set of reliable spot readings and to objectively eliminate outliers and unreliable data. However, these filters vary between authors and were not determined by studies aiming to constrain paleopoles.The studies that established the data filters, investigated PSV and geomagnetic field behavior by determining the between-site scatter of a set of VGPs (e.g.,