The dipper population of Taurus seen with K2

Dippers are typically low-mass, pre-main-sequence stars that display dips in their light curves. These dips have been attributed to dusty warps that form in the inner part of the disk. Our goal is to derive the properties of dipper stars in Taurus to assess the physical mechanisms that induce dipper light curves. We used the light curves of K2 C4 and C13 to select a dipper sample among 179 members and possible members of the Taurus star-forming region based on the light-curve morphology. We studied the periodicities by combining periodograms with wavelet analysis and derived the stellar parameters from the photometry. We also studied the morphology of the photometric dips. We find a dipper occurrence of ~30% in disk-bearing stars observed with K2 that were identified visually by us. This represents a lower limit to their true occurrence. About half of the dippers are aperiodic, and most of these are dominated by another type of variability. The chosen sample is of late spectral type (K/M), low mass and moderate mass accretion rates and has periods of a few days. We observed a transient dipper over a few rotation cycles and a dipper with a changing period. The structure of the dips can be complex and varies strongly over timescales of down to one stellar rotation. The corotation radii are located at a few stellar radii, and the temperatures at corotation allow dust survival. Many of the systems are seen at moderate to high inclination. We find that the angular extension of the dusty structure producing the dips is correlated with the stellar period. Magnetospheric accretion, which causes an accretion column and its base to occult the star, can explain most of the observed light curves. Although compatible with the model, many of the stellar inclination angles are moderate and do not exclude mechanisms other than the occultation by an inner disk warp to account for dipper light curves.


Introduction
Classical T Tauri stars (CTTSs) are still accreting from their circumstellar disk, and their activity is reflected in their complex photometric variability. Many physical mechanisms have been proposed to explain this variability, such as accretion hotspots, accretion bursts, or occultations by dusty structures in the disk (e.g., Cody et al. 2014;Alencar et al. 2010;Bouvier et al. 2003). At the next stage of protostellar evolution, weak-lined T Tauri stars (WTTSs) do not show accretion signatures in their spectral lines. Their magnetic field still produces cold spots, which give a sinusoidal modulation to the light curve.
In the past two decades, more attention has been paid to the low-mass star AA Tau. Its light curve showed a constant brightness and narrow quasiperiodic occultations that were less pronounced in the infrared (IR) than in the optical. These observations, combined with spectroscopic and spectropolarimetric data, led to the development of an occultation model. The magnetospheric accretion column would cause a dusty warp in the inner disk to obscure the star when it crosses the line of sight of the observer (Bouvier et al. , 2003. The possibility that CTTSs are still accreting from a circumstellar disk was proposed by Bertout et al. (1988) and had its roots in the work of Lynden-Bell & Pringle (1974). Although these authors already suggested that accretion might take place along magnetic field lines, the exact role played by the magnetic field in the accretion process was constrained only later. A rotating magnetosphere of a star with a strong magnetic field (approximately of some kG) disrupts the inner disk at a few stellar radii from the star and accretes the material along the closed field lines; the open field lines induce a strong disk wind (Camenzind 1990;Koenigl 1991;Shu et al. 1994). The most common magnetic field configuration for T Tauri stars is a strong dipole. Depending on the evolutionary status, other magnetic field topologies are possible, but the dipole in most cases remains the dominant component farther away from the star. Moreover, a misalignment of the dipole axis and the stellar rotation axis seems to be very common (Gregory et al. 2012). When there is such a misalignment, the rotation of the disk provokes a twist in the magnetic field lines, which are dragged in the rotation. The magnetic pressure is not the same in the upper and lower part of the disk, and the normal vector of the inner disk becomes tilted with respect to the rotation axis. This induces a disk warp (Lai 1999;Terquem & Papaloizou 2000). From the observational point of view, the strong emission lines observed in the spectra of T Tauri stars are related to accretion and not to winds (Hartmann et al. 1994). Radiative transfer modeling provided observational predictions for magnetospheric accretion (Muzerolle et al. 2001;Kurosawa et al. 2006;Lima et al. 2010).
Article number, page 1 of 29 arXiv:2106.02064v1 [astro-ph.SR] 3 Jun 2021 A&A proofs: manuscript no. 40646 Magnetohydrodynamical (MHD) simulations confirmed the possibility that magnetospheric accretion might warp the disk (e.g., Romanova et al. 2008Romanova et al. , 2013Romanova & Owocki 2015). In a stable configuration (stable accretion regime), two broad, stable, and symmetric matter streams form from the disk and are accreted as funnel flows onto the stellar surface, producing a hotspot on each hemisphere. In this scenario, the magnetosphere truncates the disk at a few stellar radii close to corotation, where the disk material rotates with ω ∼ Ω * . The base of the accretion column can host warm dust and be optically thick, which explains the photometric behavior of AA Tau.
After the the characterization of AA Tau, additional surveys identified a new class of AA Tau-like YSOs (young stellar objects), the so-called dippers. Dippers have been observed in NGC 2264, Upper Scorpius, ρ Ophiucus, and Orion (e.g., Alencar et al. 2010;Cody et al. 2014;McGinnis et al. 2015;Ansdell et al. 2016a;Cody & Hillenbrand 2018;Morales-Calderón et al. 2011). The occurrence rate of dippers is estimated to be up to 20%-30% of CTTSs. As in the case of AA Tau, the light curve exhibits a brightness continuum interrupted by narrow flux dips, which can be aperiodic or quasiperiodic. Quasiperiodic means that while the occurrence of the dips is periodic, their shape and amplitude are not constant over time. The dips are sharp and irregular, their periodicity is in the range of stellar rotation periods of K/M stars, and they last from 1-2 to 4-5 d. The amplitudes can range to up to 50 − 60% in flux.
The common consensus about the origin of dippers is that the occultation is produced by dusty structures in the inner disk. The correlation between the dip depth and the mid-infrared excess at 4.6 µm supports this hypothesis (Ansdell et al. 2016a) as IR excess at this wavelength is a tracer for warm dust at the corotation radius. This makes dippers a powerful tool for studying this region, which is difficult to resolve, yet physically complex. The dusty disk warp model ) requires dipper stars to be seen at high inclination, so that the line of sight crosses the base of the accretion column. The generalized magnetospheric accretion model (Bodman et al. 2017) allows also for moderate inclinations. The amplitude and the shape of the dips depends on parameters such as viewing angle, dust opacity, and the tilt between magnetic field and stellar rotation axis. Another possibility are vertical instabilities in the disk, which should be able to endure several rotation cycles. A star with an aligned magnetic field rather accretes via Rayleigh-Tailor (RT) instabilities (e.g., Romanova et al. 2008) in an unstable accretion regime. The socaused accretion tongues reach the stellar surface at different latitudes and induce a stochastic photometric variability. McGinnis et al. (2015) claimed that this mechanism proably lies at the origin of aperiodic dippers, which appear as stochastic occultations of the stellar photosphere, with smaller amplitudes than periodic dippers.
However, the high occurrence of dippers among YSOs is difficult to reconcile with the condition of a grazing viewing angle on these systems. It was found that the outer disk inclination of dippers can be low or even close to face-on (Ansdell et al. 2020). This either corroborates the possibility of an inner disk that is tilted with respect to the outer disk (e.g., Alencar et al. 2018) or requires other mechanisms that are compatible with lower inclination, such as disk winds (e.g., Bans & Königl 2012).
The Taurus star-forming region hosts a quite young stellar population 3 Myr (e.g., White & Ghez 2001;Kraus & Hillenbrand 2009) of ∼400 members (Kenyon et al. 2008) at a distance of ∼ 140 pc (Fleming et al. 2019;Galli et al. 2018). Rebull et al. (2020) presented the global sample of the Taurus region as covered by the K2 C13, with a few additions from C4. This work focuses on the dipper population of the same sample. In Sec. 3 the selection of the dipper sample, the periodicity study, and the stellar parameters are presented. In Sec. 5 the results and the possible scenarios to explain dippers are discussed. Final considerations and future work are presented in Sec. 6.

Observations and data reduction
The Kepler satellite (Borucki et al. 2010;Haas et al. 2010) was launched in 2009 with the main aim of detecting Earth-like exoplanets. It observed over 170,000 targets simultaneously in its long-cadence mode of 29.4 min. The failure of two reaction wheels led to a substantial change in observing strategy and to a renaming of the mission as K2 (Howell et al. 2014). The spacecraft used its two remaining reaction wheels in tandem with the solar radiation pressure (plus compensatory periodic thruster firings) to control movement in the third dimension. Instead of staring at one field, as the original Kepler mission had done, K2 was constrained to view fields in the ecliptic plane for about 70 days at a time (Howell et al. 2014). Within the framework of the K2 mission (2014 -2018), 19 fields were observed (Howell et al. 2014). Several discussions of data and reductions of Kepler and K2 data are available. In the context of this work, the two most important artifacts to keep in mind are the periodic thruster firings and the relatively low spatial resolution (pixels of ∼ 4 × 4 and the 95% encircled energy diameter of up to 7.5 px). The thruster firings happened every 0.245 days, and so apparent astrophysical periods near this must be scrutinized carefully to ensure that they are not affected by the spacecraft. The low spatial resolution means that source confusion can be an issue, and each source (and light curve) must be inspected individually to assess confusion. Customized apertures are sometimes needed to attempt to mitigate this. The effective observed range of K2 lies between 6 mag< K s <16 mag (e.g., Rebull et al. 2016). Saturation is not a relevant detection bias in Taurus, which does not host a massive star population. However, deeply embedded lowmass stars are likely to remain undetected with K2.
A few Taurus members were observed in K2 C4 (15 February -24 April 2015), but most of the Taurus members were observed during C13 from 8 March to 27 May 2017 1 . We are building on the analysis presented in Rebull et al. (2020). The K2 fields do not cover the entire Taurus cluster, but the stellar properties of the cluster are not thought to vary based on the spatial distribution (see, e.g., Fig. 1 in Rebull et al. 2020). We used light curves from several different data reduction pipelines (also see Rebull et al. 2017). First, the pre-search data conditioning (PDC) version from the Mikulski Archive for Space Telescopes (MAST). Second, a position-detrended version with moving apertures as in Cody & Hillenbrand (2018). Third, the 'self-flat-fielding' approach by Vanderburg & Johnson (2014) and the K2SFF pipeline as obtained from MAST. Last, the EVEREST2 pipeline, based on pixel-level decorrelation Luger et al. (2016), as downloaded from MAST.
A comparison of all these light-curve versions allows us to constrain the periodicity (or periodicities) and select the best available periodicity (see Rebull et al. 2020, and references therein) to represent that star's light curve. If the light curve does not present particular issues, the version with moving aperture as in Cody & Hillenbrand (2018)  tency. This version was also used for the dipper sample presented in this paper. As described in Rebull et al. (2020), we started from a more expansive definition of possible Taurus members with K2 light curves, and then weeded it down to 156 members and 23 possible members. The highest-confidence members are those listed in Luhman (2018). We searched for dippers in this set of 156+23 Taurus member light curves.

Identification of dippers
We identified dippers based on a visual inspection of the light curve (see lists in Table 1 and Table 2). A common feature in YSOs are sinusoidal patterns, which can be ascribed to cold spots on the stellar photosphere (e.g., Herbst et al. 1994). Irregular, aperiodic increases in the stellar brightness (bursts) are linked to accretion shocks for T Tauri stars . Planetary transits and eclipsing binaries do produce dips, but these are strictly periodic and in the first case, have a small amplitude in flux. The criteria for selecting dippers are thus the irregular and sharp shape of the dips, the duration of the dips of up to a few days, and an amplitude of at least ∼10% to ∼60% of the flux. A few examples of dipper light curves are displayed in Fig. 1. The light curves easiest to classify are those in which a stable brightness continuum is interrupted by irregularly shaped dips (e.g., JH 223 in Fig. 1). The bright-faint flux asymmetry is thus high (see, e.g., the M metric in Cody et al. 2014). The M metric is a useful tool to investigate a large sample of stars. However, it tends to classify as quasiperiodic symmetric (QPS) light curves with irregular dips that do not have a constant brightness continuum. We consider that from a physical point of view, this should not be a striking selection criterion and therefore preferred a manual classification. Dippers with an irregular brightness continuum are CFHT Tau 12 and ITG 34. More symmetric are HK Tau, GH Tau, GM Tau, and EPIC 247820821 (2MASS J04295950+2433078). The last two are more ambiguous, and we differentiated them from spots because of the sharp shape of their dips. We cannot exclude that the dips are superimposed on spots. GI Tau is contaminated by the neighbor GK Tau, and its light curve is complex. However, its own periodicity can be retrieved, and clear dips are visible in the light curve. DK Tau also exhibits a complex light curve. After ∼25 d of quiescent state, the star becomes a clear dipper. All dipper light curves are shown in Fig. B.1.
In addition to the main dipper sample, two other groups of light curves were taken into account: light curves that show a predominant behavior (e.g., spots, bursts, in Table 2) but also some dips, and low-quality dipper candidates that did not fulfill all of the criteria of the visual inspection and/or were too complex for discerning the different types of their variability (Table 3). In Fig. 2 a few example light curves are shown.
As a result, 22 objects were classified as dippers and 12 were classified as dippers dominated by another type of variability. Almost all of the latter exhibit some aperiodic dips superimposed Article number, page 3 of 29  Tau); a low-quality dipper candidate (IQ Tau). Although dips as described in Fig. 1 are present, the light curve is dominated by another type of variability.
on an otherwise variable light curve (Fig. B.2). Eleven additional stars are presented as low-quality candidates in Table 3 and Fig. B.3 and are not further considered in this study. They are rather listed as visual examples of the selection process. The final sample of 34 dippers represents 19% of high-confidence and possible Taurus members (156 + 23 as in Rebull et al. 2020) and 31% of the disked stars (94 + 7 as in Rebull et al. 2020) in this sample.
In addition to the dipper prototype AA Tau (which is now in a faint state and appears in Table 2), Rodriguez et al. (2017) identified five dippers in Taurus-Auriga with KELT light curves, four of which are confirmed in this sample. Rebull et al. (2020) counted 21 dippers in this sample. Sixteen of the dippers analyzed in this study are present in Rebull et al. (2020). Four out of the 5 further dippers in Rebull et al. (2020) appear here in Table 2. One further light curve (EPIC 247078342) indicated as a dipper in Rebull et al. (2020) is discarded in this study because its shape strongly depends on the reduction pipeline. Our classification is therefore not as strict as in Rebull et al. (2020). Further analysis of the Taurus members sample will be presented in Cody & Hillenbrand (in prep.). As mentioned above, the usage of the M metric or a manual classification can affect the classification for certain stars.

Wavelet analysis of dippers
Dippers can be quasiperiodic or aperiodic, their shape varying from cycle to cycle. Aperiodic narrow dips might appear in an otherwise periodic light curve, or a quasiperiodic dipper might be a transient phenomenon (e.g., McGinnis et al. 2015). Because the variability is irregular, the use of time-resolved period-search algorithms is of high interest.
Common tools used in astronomy such as the Lomb-Scargle periodogram (Lomb 1976;Scargle 1982) and Fourier analysis are not able to deliver information about transient phenomena because they merely resolve the time series in the frequency domain. The windowed Fourier transform (WFT) convolves the signal most commonly with a Gaussian window that is shifted along the signal before the power spectrum is computed. The resulting spectrograms to some extent allow retrieving a time resolution and were applied to dippers by Bodman et al. (2017). However, the width of the window and the time shift cannot be optimized for the entire frequency range, resulting in poor resolution.
When transient phenomena are to be characterized, the wavelet analysis represents a valid alternative because this method allows keeping track of the time variability of a periodic feature by partly losing frequency resolution (see for a summary Torrence & Compo 1998). It can be conceived as a time-resolved WFT with a variable window width that allows recognizing both high-and low-frequency features. A wavelet transform produces The logarithmic y-axis with the periods is the same as for the WPS. In the case of DK Tau (left), the two peaks in the periodogram at ∼8 and ∼10 d can be interpreted as a single period that changes during the observations by means of the WPS. GM Tau displays only one clear period in the WPS and the periodogram. The green stripe at ∼10 d is just a recurring pattern in the light curve: a deeper dip at t = 3030, 3040, 3050 d. This is not a relevant periodicity and does not appear in the periodogram. a power spectrum after the convolution of the signal and the wavelet, which in most cases is a sort of time-confined sinusoidal wave. We used here the complex Morlet wavelet, which is a complex exponential multiplied by a Gaussian envelope (Grossmann & Morlet 1984). The time and frequency resolution are reached by shifting the wavelet along the time series, then stretching or compressing the wavelet, and repeating the procedure again. This delivers a two-dimensional wavelet power spectrum (WPS), which has a power value for each combination of time and frequency (for more information, see Appendix A).
Wavelet analysis has been successfully applied to a number of phenomena in K2 light curves: stellar rotation (Mathur et al. 2014;García et al. 2014;Bravo et al. 2014), planetary transits, stellar variability, pulsation, and binaries (Bravo et al. 2014). The implementation in python runs with the package pywt.
Because the time series must be evenly sampled for wavelet analysis and some values in the light curve are flagged out during the data reduction, the light curves are linearly interpolated onto even time steps in a preliminary step.
An example of WPS is shown in Fig. 3, while the full atlas of the sample is presented in Appendix E. The x-and y-axis present the linear time and the logarithmic periods (equivalent to the scales, as explained in Appendix A). The contours correspond to the power of the wavelet spectrum, and the region omitted from the cone of influence (COI) is cross-hatched. The contours reach from blue (low power) to yellow (high power). To increase the readability of the plots, the contour plot is saturated at the 99th and third percentile of the power. Moreover, the periodogram and the WPS are cut for the plot at a lower limit of 0.5 d. The WPS is compared to the CLEAN (Roberts et al. 1987) periodogram (to the right of the WPS). For this sample of light curves, this periodogram shows a slightly higher frequency resolution than the Lomb-Scargle periodogram by better resolving close double-peaked periods. The WPS delivers a lower frequency resolution than the periodogram and close peaks are resolved only poorly; on the other hand, it becomes clear when a certain periodicity is present in the time series. High peaks in the periodogram that appear for timescales on the order of one pe-riod in the WPS are not real periodicities, but rather parts of the light curve that just highly correlate with the shape of the Morlet wavelet. When two major periodicities are present, the WPS allows us to distinguish which of the two is really predominant in the signal. In case of a changing period, the change can be timeresolved by the WPS, but will appear as several peaks in the periodogram, as is the case for DK Tau (Fig. 3). Depending on the needs of any given light curve, a higher time or frequency resolution can be obtained by changing the central frequency or the bandwidth of the Morlet wavelet (see Appendix A). This does not affect the physical interpretation of the results. Because of the lower frequency resolution of the wavelet transform, it is better to use the period as retrieved from the periodogram. It is possible to project the WPS onto the y-axis and thus average over time to obtain something similar to a periodogram. However, the lower frequency resolution of the wavelet transform does not provide any particular advantage for the period determination alone compared to a periodogram.

Period analysis
Dippers often appear to be quasiperiodic. Following the disk warp model (e.g., Bouvier et al. 2007), the dusty material obscuring the star must be located close to corotation, which is per definition the region in which circumstellar material rotates with the same angular velocity as the star. Stauffer et al. (2015) and Rebull et al. (2018) identified some light curves that showed both spots and dips. The overlapping periods of the two phenomena strongly supported the scenario of an inner disk that was rotationally locked to the star because stellar spots follow the stellar rotation. It is thus of interest to discuss the periodicities of dippers compared to stellar rotation.
The variations due to spots are smooth, largely sinusoidal variations. The quasiperiodic dippers produce rather sharp or complex dips, and these events are usually slightly different each time they occur. The distinction between variations originating in spots and dippers can be ambiguous in some unusual cases.  (a) and (c), folded at 2.99 d. The noise is high with respect to the dips, but no evident pattern is present.
Sometimes multiple real periodicities (not harmonics) can show up in the periodogram or in the WPS; in these cases, the physical interpretation of the periodicity can be complicated (see, e.g., discussion in Rebull et al. 2018Rebull et al. , 2020. In this sample, eight stars exhibit a single clear periodicity (JH 223, EPIC 246989752, CFHT Tau 12, V807 Tau, GK Tau, HK Tau, GM Tau, and IS Tau), three stars are aperiodic (HD 285893, FX Tau, and GO Tau), and for another three stars, the periodicity is unclear (Haro 6-37, St 34, and DK Tau). Two stars are transient quasiperiodic dippers (JH 112 A and DK Tau). The remaining seven dippers are quasiperiodic, but their light curve or periodogram is more complex than that of the other dippers: LkCa 15, EPIC 247885481, HP Tau, GH Tau, GI Tau, EPIC 247820821, and ITG 34. The periods are listed in Table 4. Those that are uncertain are flagged with a semicolon.
The periods coincide with those published by Rebull et al. (2020), with the exception of that of DK Tau (7.84 d vs. 7.69 d derived here; this small difference is due to the complex periodogram and the usage of two different algorithms). The WPS (Fig. 3) shows that the period increases during the K2 campaign. The periodogram exhibits two peaks at ∼8 and ∼10 d. By means of the WPS, they can be interpreted as a single changing period. This could explain why the period of DK Tau is never constrained precisely in the literature (see Appendix D for more details). The stars ITG 34 and GH Tau have a second reported period (Rebull et al. 2020), which is confirmed here. For GH Tau, a possible third period at 5.09 d is present. An interesting case is EPIC 247885481 (see Sec. 3.3 for a detailed discussion), where the main periodicity of 2.99 d is most probably caused by a spot, while the dips can be folded in phase with a period of 2.90 d. Rebull et al. (2020) reported 2.99 d for the main periodicity, while the dipper period is unresolved both in the periodogram here and in Rebull et al. (2020). In the case of GI Tau, the light curve is contaminated by the neighboring GK Tau, and the dominant peak in the periodogram is the period of GK Tau. Only the second period is therefore reported as being of GI Tau. EPIC 247820821 shows a second periodicity on the WPS at 7.00 d that does not appear in the periodogram. When folded in phase, the pattern appears to be related to a shift of the minimum of the dip. The periodogram of ITG 34 exhibits several double-peaked periods that cannot be harmonics of each other. It is difficult to recognize a dominant structure in the folded light curve. All the folded light curves are presented in Appendix F. The motions of the Kepler satellite might create spurious periods at 0.22, 1.75, and 1.97 d. None of them has been observed here.
A different period for spot and dips? The light curve for EPIC 247885481 (2MASS J05023985+2459337) shows narrow dips superimposed on a sinusoidal variability. The periodogram shows a clear peak at 2.99 d with σ = 0.08 d, but the folded light curve shows that the period refers to the sinusoidal variation (Fig. 4). The dips can be folded in phase around P − σ at ∼2.90 d. The peak in the periodogram does not appear to have any complex structure.
In order to better study the periodicity, the spot was fit with a sinusoidal wave with a period of 2.99 d (Fig. 4), which was then subtracted from the light curve (Fig. 5). The periodogram before and after the subtraction shows that the main contribution to the periodicity is given by the spot; no substructure appears after the spot subtraction. This can be explained by the fact that the amplitude of the eclipses is very small compared to the noise and that the dip is not clearly present in every period. The shape of the dips is nearly Gaussian, and the dip width is small compared to the period. The residual was fit with a periodic Gaussian pulse such as , with the period T as a free parameter. The retrieved period is thus 2.90 d for the dips, still very close to P − σ, and it cannot be confirmed that they are significantly different. Similar conclusions about other stars were reached by Stauffer et al. (2015). Nevertheless, the fact that the dips cannot be folded in phase with the period of the spot is a clear indication that they must be different. The difference between the two periods could be explained by the differential rotation of the stellar surface, assuming that the corotating material is aligned with the equator and the spot is at higher latitude.

Light-curve morphology
In the following section, the light curve properties of ideal and peculiar dippers are briefly presented and the determination of the dip properties, such as amplitude and width, are explained and discussed.
Ideal dippers As mentioned in Sec. 1 and in Sec. 3.1, the main characteristics of dippers are their dips. These are irregular in shape and can either be aperiodic or quasiperiodic. A good example for a quasiperiodic dipper is JH 223 (Fig. 1): The occurrence of the dips is periodic, their shape is irregular, and their amplitude can almost double in two neighboring dips. No other types of variability in form of bursts or other patterns affect the classification. Its period of 3.31 d and the dip duration of ∼2 d are in range of the rotation periods of CTTSs and of observed dip properties. Aperiodic dippers tend to present either very narrow dips whose shape is rather simple (e.g., Stauffer et al. 2015), which is the case of HD 285893 ( Fig. 1), or broad and complex dips.  a period of 2.21 d, and if directly linked to the stellar rotation, they are present for a timescale of about ten full rotations. Interestingly, the continuum brightness decreases only as long as the eclipses occur. DK Tau also changes from a fade state to a dipper with a higher brightness continuum, but its light curve is complex, and a strong overall variability has already been observed in former campaigns. The dipper status of quasiperiodic dippers is known to be transient over a few years (e.g., AA Tau), but it is rare to observe a quasiperiodic dipper on such a short timescale. The explanation of a disk warp would imply that the warp significantly changes its height and becomes not visible for the observer after a few rotations, or that the dust in it is completely dissipated on a short timescale.

Transient dippers
Double-and multi-peaked dippers In the sample studied, six periodic dippers exhibit clear double-peaked dips in their light curve: LkCa 15, JH 112 A, HP Tau, GI Tau, HK Tau, and IS Tau (Fig. 6). These double-or multiple-peaked dippers exhibit at least two well-detached dips in the folded light curve (see App. F). It should be noted that the dips themselves are in gen-eral not Gaussian and have a complex shape. The determination of the dip width for these sources depends on the desired information (the width of the primary peak or the total width of the multi-peaked dip) and is handled in Sec. 3.5. A more detailed description regarding the double dips of individual sources is provided in Appendix D. In the scenario in which a dusty warp occults the star, the constant presence (or absence) of a double dip delivers an indication about the stability and the shape of the warp. In general, it can be remarked that the dips, although periodic, do not occur exactly at the same phase. For HP Tau, the WPS shows that the periodicity slightly varies around a period of 4.33 d. The light curve displays clear shifts of the minima of the dip for different rotation cycles. Moreover, the variations in amplitude of the different dips do not seem correlated to each other. This suggests that either smaller and independent dusty warps occult the star, or that the shape of the dusty structure is highly unstable. A quantitative study of this variation will be object of future work.

Dip amplitude and dip width
Among all other sources of variability for YSOs, the long-term trends particularly affect the determination of the dip amplitude because the continuum brightness is unstable. To determine the dip amplitudes, the light curves were therefore detrended using a boxcar filler with a width chosen to be four times the period. This removed the long-term variability and did not interfere with the dips (Fig. 7). For nonperiodic dippers, a standard window of 7 d was chosen. The amplitudes were computed as the difference between the 90th and fifth percentile and peak-to-peak and are listed in Table 4. The light curves of Haro 6-37, JH 223, IS Tau, and DK Tau (for this star, the faint state of the first ∼30 d has not been included for this and the following analysis) did not need detrending because the continuum brightness is stable. The detrending was not applied for the study of the periodicity because it might affect the result and interfere with the physical phenomena that cause variability. It is only used in a prag- matic approach to study the dip morphology, without the effect of other stellar variability.
To define an eclipse width, it is useful to bin, that is, to average, the detrended and folded light curve (purple line in the upper panel of Fig. 8). The dip width (black horizontal line in the upper panel of Fig. 8) is then defined at the full width at half maximum (FWHM) of the resulting binned curve. The results are listed in Table 4 and shown in Appendix F. In general, the dip width is about half of the period. The largest dip widths are linked to double-or multi-peaked dippers, with the exception of JH 112 A and HK Tau, which both exhibit a relatively small dip width ( 0.5 P) and the shortest periods in this group (2.21 and 3.3 d, respectively). The discussion of the correlation between dip width and period is presented in Sec. 5.4.
We estimate an uncertainty on the dip width as follows. For each phase bin, the standard deviation σ(phase) of the flux between the different cycles was first computed (purple line in the upper panel of Fig. 8). For every point of the binned light curve, the standard deviation of the phase bin was then added in one case and subtracted in the other (purple line on center and bottom panel in Fig. 8). Finally, the dip width was computed again as FWHM of these modified binned light curves. Upper and lower error bars are defined as difference between the dip width of the modified binned light curve and the original one. Because the shape of the dips may vary strongly after this operation, these error bars can be large.
For double-peaked dippers, the line drawn at the FWHM can cross the binned light curve in more than two points: the error bars w + σ and w − σ were then chosen to represent the primary peak width and the total width of the dip, which is very large in these cases. This can result in strongly asymmetric error bars (see Sec. 5.4). In Table 4, some width values are flagged because they are not reliable. This can happen when the folded light curve is very noisy and the binned light curve does not show the dips as they can be seen by eye. These plots appear in Fig. F.1.
Another possible representation of the uncertainty on the dip width would be to compute the dip width at different thresholds of the same binned light curve, for example, at 0.3, 0.5 (FWHM), 0.7 of the flux maximum. However, the information about the width of the main dip and the total dip in double-peaked dips (see, e.g., IS Tau) would be lost.

Stellar parameters
Dipper stars are most commonly low-mass T Tauri stars. Their periods of a few days are consistent with the range of rotational periods found in general for young, low-mass stars. There is consensus that the dips in the light curve are caused by dust; in the dusty disk warp scenario, the inner disk warp is located at the corotation radius. Moreover, regardless of the position of the occulting dusty structure, the temperature must be low enough for dust to be able to be present, that is, not to sublimate.
In the following section, the derivation of the different stellar parameters for dipper stars is discussed. In order to verify the conditions and to compute the radius and temperature at corotation, the masses and radii of the star are required. The discussion about each parameter follows in the next paragraphs.

Effective temperatures and spectral types
Effective temperatures are listed in Table 1 and were derived according to the SpT-T eff conversions of Pecaut & Mamajek (2013). The main sources for the spectral types are Esplin & Luhman (2019) and Herczeg & Hillenbrand (2014), who agree for almost all objects. When spectral subclasses were not explicitely listed in the conversion tables, the value of T eff was linearly interpolated between the two closest subclasses. For spectral types later than M5, the temperature conversions by Herczeg & Hillenbrand (2014) were used.
The uncertainty on T eff follows the uncertainty on the spectral type: for spectral types up to K9, the uncertainty is one spectral subclass; for spectral types between M0 and M4, 0.4 subclasses; and for spectral types later than M4, 0.25 subclasses (Herczeg & Hillenbrand 2014). As a consequence, asymmetric error bars were computed from Pecaut & Mamajek (2013) for an earlier and later spectral type, respectively. If the later spectral type is missing at the end of the table, a symmetric error bar is produced. The systematic uncertainty linked to the choice of the models of up to 150 K is not taken into account here. A discussion of the consequences on the derived stellar parameters is given in Sec. 5.1.

VJHK photometry
The light curves of young, accreting stars might be affected by contamination (e.g., by hotspots) in the blue band. Photometry in the V and JHK bands is therefore preferred. The dipper list was cross-matched with the Two Micron All-Sky Survey (2MASS, Skrutskie et al. 2006) to obtain JHK photometry.
Because a photometric monitoring is the best way to determine the magnitude of variable stars, values reported in the photometric monitoring by Grankin et al. (2007) are preferred for the V band, where available. An inspection of the light curves presented for the objects studied here led to the choice of their V m as best estimate for the brightness continuum. However, many of the dipper stars are absent from that collection, and other sources have to be included as well. Three photometric catalogs are common to most sources and include V: GSC2.3.2 (Lasker et al. 2008), NOMAD-1 (Zacharias et al. 2005), and APASS (Henden et al. 2016). 2 For individual stars, other additional measurements were retrieved from the literature. All available values for each star were inspected to verify whether any strong inconsistencies existed between the collections. More entries are available in the GSC catalog than in any other collection. If the star was not already reported in Grankin et al. (2007) and there were no particular problems with the photometry, the entry of the GSC catalog was set as V value in Table 1, when available.
In the few cases where these single measurements were highly different (∆V > 1 mag) between the catalogs because of the intrinsic variability of young stars, special attention was given to those stars. One way to assess the plausibility of a V band measurement is to verify that it does not lead to a strongly negative extinction. In order to verify that the individually selected visual magnitudes from different sources were selfconsistent, the values of V were compared with the Gaia G (Gaia Collaboration et al. 2016, 2018) converted into V. No offset was found.

Visual extinction
Most stars in Taurus are subject to at least some reddening. Many stars have moderate or strong infrared excesses, which can affect even the near-IR bands. At the same time, the fainter stars are harder to detect at shorter wavelengths. We therefore calculated the visual extinction by comparing the (V-J) colors as observed to those for the corresponding spectral type in Pecaut & Mamajek (2013). Herczeg & Hillenbrand (2014) measured a quite high veiling at 7510 Å for some objects in this sample. For those with high veiling, extinction derived with photometry is not as accurate as with spectroscopy. For the stars with r > 0.1, the A V computed by Herczeg & Hillenbrand (2014) from spectroscopy is preferred.
We wished to minimize the risk of inconsistent values that is attached to very few, or just one, photometric measurement. For this reason, the available measurements of A V in the optical in the literature were collected and compared (see Table C.1). This allowed us to filter out several clear outliers that are flagged in Table C.1. For stars with lower veiling, the final A V is an average of the extinction computed here using (V−J) and the other values in the literature (a comparison of the values we derived and those in the literature is shown in Fig. C.1). The uncertainty was set to the rms of the different values, with a reasonable minimum of  (2019)   0.3 mag. For measurements of A V < 0, A v was set to 0 in order to be physical, as was the uncertainty if A V 0.

Luminosities
The stellar luminosity was computed from the bolometric magnitude M bol as in where M bol, = 4.25 mag. Because the brightness and bolometric corrections are mostly available for the J band, the bolometric magnitude was derived according to where m J is the J-band magnitude, p is the parallax (here from Gaia), A J is the extinction in the J band, and BC J is the bolometric correction for the J band. The bolometric corrections were retrieved from Pecaut & Mamajek (2013) and interpolated for intermediate spectral types. In case of spectral types later than M5 (not presented in the tables), the BC J was extrapolated. The computation of A J from (J − K) is affected by IR excess. Thus, we preferred to convert A v into A J as A J A v = 0.282 (Cardelli et al. 1989).  Notes. Columns: the photometric period P, the uncertainty estimated from the Gaussian fitting of the periodogram peak σ P , the dip width W measured as the FWHM of the folded and binned light curve in units of phase, the dip width estimated considering the binned light curve ±σ of the flux (W −σ , W +σ ), the dip amplitude with the 90 and 5 flux percentiles A, the peak-to-peak amplitude p2p. The stars are detrended before the width was determined, with the exception of those marked (a) . (b) Value discarded because it was not representative of the real dip width. (c) Double-peaked dip. (d) The peak of the periodogram is at 2.99 d and is linked to a cold spot. A more thorough analysis is needed to retrieve the dipper period at 2.90 d.
The Gaia parallax was not available for three objects (GH Tau, IS Tau, and FX Tau). We used the distance distributions in Taurus by Fleming et al. (2019) for these stars. FX Tau and GH Tau are located in the B18 cloud, which contains members of the near population at 127.4 ± 3.8 pc. A standard deviation of 7.9 pc was applied as uncertainty for unknown parallaxes. IS Tau lies on the filament L1495, which contains members of both populations. Its luminosity was computed for the two mean distances and is shown with an arrow in Fig. 9. IS Tau appears to more likely be a member of the far population when we assume that it has the same age as the other dippers. In the case of Haro 6-37, the 2MASS measurement in the J band is corrupted and no other measurements are available. The visual magnitude and the corresponding BC V were considered instead.
The procedure for deriving uncertainties on the bolometric correction is the same as for the effective temperature. A major issue for the determination of the luminosity uncertainty is the intrinsic photometric variability of dipper stars, which cannot be included here in absence of an extensive observation campaign. Nevertheless, for stars on their Hayashi tracks, the derivation of the mass depends much more on a precise spectral type than on luminosity (see Fig. 9). The uncertainty on the luminosity affects the derivation of age far more; the isochrones in Fig. 9 accordingly are approximations.

Masses and radii
In order to derive the mass of the dipper stars in this sample, evolutionary models that reach the lowest end of stellar formation are needed. We therefore preferred the tracks of Baraffe et al. (2015) (Fig. 9) as they include M dwarfs. Notes. Stellar inclinations reported as lower limits have sin i > 1. The aim here is to be self-consistent, therefore we derived all masses from the same model. The effective temperature is the dominant term for identifying the corresponding track for this sample. Thus, an identification of the corresponding tracks along the x-axis is the most efficient, with the only exception of the latest M dwarfs, for which the luminosity also plays an im-  Rebull et al. (2020) and appear as gray points on the HR diagram. The spectral types are converted into temperatures with the same conversion as we applied to dippers. For spectral types later than K4, the conversion provided by Herczeg & Hillenbrand (2014) exhibits a systematic temperature offset of about +60 K and up to +130 K, while for earlier spectral types, the temperatures proposed by Pecaut & Mamajek (2013) are higher. When several luminosities are provided for a star, the values were averaged in the HR diagram. Because dippers are associated with circumstellar extinction events and their V measurement might have been taken in the fainter state, it might be expected that they are less luminous than non-dipper stars. This does not seem to be the case here. Interestingly, many stars with a mass between 0.1 and 0.3 M lie well above the youngest   Baraffe et al. (2015). The brown dwarf limit lies at 0.073 M . Isochrones from top to bottom: 0.5, 1, 2, 5 Myr, and 1 Gyr. The red points are the Taurus members classified as dippers, while the gray points are Taurus members as in Herczeg & Hillenbrand (2014). HD 285893 (SpT F8) does not appear in this plot for the sake of readability because its derived T eff is much higher than that of the rest of the sample. The stellar ages scatter around 1 Myr, and with the exception of LkCa15 (marked with a black arrow) and HD 285893, they are on their Hayashi tracks and are still fully convective.
isochrone. On the one hand, this might be an indication of the difficulty of providing precise evolutionary models for stars with a very low mass. On the other hand, for faint and accreting brown dwarfs, the accretion luminosity might dominate the faint stellar luminosity. At the other extreme, with a T eff = 6100 K and L = 1.99 L , HD 285893 almost lies on the ZAMS on the HR diagram. A more detailed discussion is given in the individual notes in Appendix D. This 'old' dipper might belong to a recently identified class of debris disk dippers (Gaidos et al. 2019;Tajiri et al. 2020), for which aperiodic extinction events have been attributed to the disruption of planetesimals.
The stellar radii are derived according to the Stefan-Boltzmann law as and are presented along with the masses in Table 5. The uncertainties on T eff and the luminosity are propagated to the stellar radius.

Radius and temperature at corotation
In the accretion scenario, the dusty warp that obscures the star is located at corotation. The corotation radius defines the distance from the star at which the material of the Keplerian disk rotates with the same angular velocity as the star. Because in a stable orbit, the gravitational potential of the star is equal to the centrifugal force, the corotation radius is derived as where P is the stellar period. Considering the simple approximation of no energy transfer between a dust grain situated at R cor and the remaining disk, which reduces the effective irradiation by a factor 4 as in Bouvier et al. (1999), it is possible to derive the temperature at corotation as The corotation radii and temperatures computed in this way are listed in Table 5. The corotation radii extend to a few stellar radii, between 2.5 and 9 R * . The temperatures at corotation are ∼ 1000 K, and a few stars reach up to ∼ 1500 K. The error bars on this quantity are large because the uncertainties on the stellar radii are about 10%-20% for young stars.

Inclinations
In order to probe the capacity of the magnetospheric accretion model to account for dippers, the angle under which we observe the system is required because dippers cannot be seen close to face-on if the dusty part of the accretion column has to cross the observer's line of sight to produce dips in the light curve (e.g., McGinnis et al. 2015;Bodman et al. 2017). A more detailed discussion of this point is provided in Sec. 5.1. The stellar inclination is derived according to the formula v sin i = sin i 2πR * P .
The convention we used to define the inclination is 0 • for the star seen pole-on and 90 • for edge-on. The v sin i values were retrieved from the literature (see Table 5): the main sources were Nguyen et al. (2012) and Güdel et al. (2007), who retrieved their values for this sample from Rebull et al. (2004). The uncertainty on i grows with sin i, thus higher inclination angles have larger uncertainties. Some stars (Haro 6-37, EPIC 247820821, and DK Tau) exhibit a sin i > 1, which probably indicates that something is amiss with the stellar parameters. In these cases, the minimum inclination angle is indicated as a lower limit. The results are listed in Table 5. For Haro 6-37, the period is uncertain. It is possible that the observed period is not correct or not directly related with the rotation period. For DK Tau, different periods have been reported in the literature because the light curve is complex, and the WPS shows an increasing period during the K2 observations. Because the stellar rotation period is not assumed to vary much, other physical phenomena might affect the dips. For instance, Grankin et al. (2007) reported a long-term variability amplitude of ∆V ∼ 1.8 mag for this star, which propagates to any inclination angle between 30 • and 90 • if taken into account as the uncertainty on V.

Mass accretion rates and accretion regime
The mass accretion rates of the dippers were collected from the literature. When available, they are listed in the order of magnitude of 10 −8 M yr −1 (Table 5). This is in agreement with the observed mass accretion ranges for T Tauri stars (Gregory et al. 2006). Exceptions are Haro 6-37, HK Tau, and DK Tau, which are stronger accretors of about 10 −7 M yr −1 . LkCa 15 has a very low mass accretion rate for its mass, but it is also a transition object, and a lower mass accretion rate is expected at this stage. In general, quasiperiodic dippers are expected to be generated in a stable accretion regime (see Sec. 5.4) and aperiodic dippers instead in an unstable regime , which could be caused by a phase of enhanced mass accretion. The complex light curve of DK Tau, one of the strongest accretors, might be an example of unstable accretion. However, for the dipper sample presented in this paper, we found no correlation between mass accretion rate and periodicity.

Discussion
Several possible mechanisms for the origin of dippers have been presented in the literature, including dusty disk warps, dusty winds, and disk vortices. In the study of the protoypical dipper AA Tau, Bouvier et al. (1999) proposed that a magnetic dipole, tilted with respect to the stellar rotation axis, might contribute to an inner disk distortion. This can produce a dusty, optically thick disk wall that can account for the photometric variation, if observed at high inclination. This model was further explored by McGinnis et al. (2015), using as free parameters the stellar inclination and location, height, and azimuthal extension of the dusty wall. In the dipper sample of NGC 2264, which has a similar size as that of Taurus, inclinations down to ∼50 • could be measured, with large uncertainties. McGinnis et al. (2015) were able to fit the light curves with a dusty warp with a mean inclination of ∼70 • . Dippers close to edge-on are assumed to be undetectable in photometry because the disk would block the line of sight toward the stellar photosphere.
The first parameter we discuss for the compatibility with the model therefore is the inclination of the star and the disk (Sec. 5.1). For stable accretion, the magnetospheric truncation radius is close to corotation, and the photometric periods of dippers provide evidence for the dusty warp to be located at corotation. One other condition for the validity of the model is thus that the environment at corotation must be cold enough to avoid dust sublimation (Sec. 5.2). It would be helpful to derive the truncation radius to compare it with the corotation radius; unfortunately, not enough precise magnetic field measurements are available for this purpose. Bodman et al. (2017) provided a  rough estimate for the mass accretion rate and the magnetic field strength. However, this estimate has a very large uncertainty. This scenario was challenged by the discovery of dippers seen in millimeter wavelengths in a full range of inclination angles (Ansdell et al. 2016b(Ansdell et al. , 2020, on which the hypothesis was founded that the inner and outer disk might be misaligned. Binarity can be invoked as a cause of misalignement in a circumbinary disk (e.g., Facchini et al. 2013Facchini et al. , 2018Franchini et al. 2019). This opens the possibility of finding a dipper with a low outer disk inclination. Nonetheless, the inclination of the inner region must be reasonably high to detect dips. Two examples of wellstudied dippers with a misaligned inner disk are LkCa 15 and RX J1604.3-2130A (Alencar et al. 2018;Sicilia-Aguilar et al. 2020).
Another way of lifting dust away from the disk midplane is through dusty disk winds. These are spatially more extended than the base of accretion funnel flows, and might shield the stellar brightness at lower viewing angles. Magnetic field lines inclined with respect to the disk symmetry axis (θ ≥ 30 • ) are able to accelerate matter away from the disk. Dust, dragged away with gas, is able to survive in a disk-driven wind at temperatures of several 10 3 K because heating the dust through collisions with gas particles and sputtering is highly inefficient (Tambovtseva & Grinin 2008). A high mass-accretion rate and strong magnetic field (kGauss) are required to launch such a disk wind (Miyake et al. 2016;Labdon et al. 2019), and only the small dust grains are dragged by the wind, while the large-size grains remain in the disk midplane (Miyake et al. 2016). The resulting floating dust in the disk is able to produce optical fading events and NIR brightening. Bans & Königl (2012) modeled the NIR 3 µm bump as produced by absorbed stellar radiation in a dusty wind. However, it seems that such disk winds are observed under a rather high inclination ∼ 70 • (Vinković &Čemeljić 2020), much higher than what has been invoked to explain low-inclination dippers (e.g., Ansdell et al. 2020).
It is also possible to lift dust from the midplane through vortices in the disk that are caused, for example, by Rossby waves. However, the occultations caused in this way would have a very small amplitude and might thus account only for a small fraction of aperiodic dippers (Stauffer et al. 2015).

Inclination of star and disk
In this section, we investigate whether the stellar inclination agrees with the inclination of the outer disk and if this is compatible with an AA Tau-like star. For most of the objects considered here, a moderate to high inclination is clear despite the large error bars.
GH Tau and GM Tau are compatible with an inclination of ∼50 • , which is at the lowest end of the magnetospheric accretion scenario. This suggests that the magnetic field axis is highly tilted with respect to the stellar spin axis. In the case of HP Tau and V807 Tau, the proposed dusty warp scenario for the production of dippers fails to explain fading events seen at such a low inclination. We tested whether the photometric variability might significantly affect the estimation of the luminosity and propagate up to the derived stellar inclination. None of these stars shows an amplitude in the K2 light curve that would be large enough to explain the discrepancy. Nevertheless, for HP Tau, a long-term amplitude variation of ∆V ∼1 mag was reported by Grankin et al. (2007). This can be roughly converted into an upper limit of i ∼50 • , which is still at the lowest limit for the magnetospheric accretion model. For V807 Tau, which is a triple system, an overestimation of the brightness of the primary is possible. Schaefer et al. (2012) estimated that V807 A could be 0.5 mag fainter than the total system in the J band and fit an extinction of 0 to the system. From these parameters, a higher stellar inclination of 40 • can be derived, which is more plausible. This is not sufficient to explain the discrepancy, however. A stel-lar radius of ∼ 1.5 R would be required to reach an inclination angle of at least 50 • .
It is thus more probable that some other phenomenon contributes to the quasiperiodic extinction events for these stars that are inconsistent with a sufficiently high inclination. Dusty disk winds also need high inclination angles to be observable (Vinković &Čemeljić 2020), and the mass accretion rate of V807 Tau is too low ( Table 5). The small amplitude of ∼ 0.1 mag might indicate vortices caused by Rossby wave instabilities. However, given the clear periodicity of the light curve, this mechanism can be discarded. Moreover, the almost constant shape of the dips means that the structure occulting the star is probably stable, not very large, optically thick, and it might be that just a small fraction of it is seen by the observer under this low inclination.
One more possibility in the magnetospheric accretion scenario would be that the dust high in the accretion column increases its optical depth immediately before sublimation because evaporation and free-fall timescales are similar (Nagel & Bouvier 2020). This might imply that an optically thick part of the accretion column could cross the observer's line of sight, thus requiring a lower inclination to see dips.
We compared our dipper sample to inclinations of outer disks that were directly measured at millimeter wavelengths (Table 6). Akeson & Jensen (2014) derived the disk inclination in the image plane by fitting the clean continuum maps with 2D Gaussians. This results in larger error bars than with the methods used by other authors because the continuum map is reconstructed from the visibilities. In general, the inclination of the outer disk is slightly lower than the inclination derived for the inner region, with the exception of DK Tau, which seems strongly misaligned with the outer disk, although this value has to be considered with caution (see Sec. 4). Appenzeller & Bertout (2013) pointed out that the inclination of CTTSs derived from v sin i are not precise for i > 30 • because the uncertainty grows with i. The inclination might also be overestimated when not all line broadening effects are properly considered. For a sample of known CTTSs in different regions, i * and i disk are correlated, but the inclination derived from rotation is 19 • higher on average (Appenzeller & Bertout 2013). This trend is also observed in Taurus for HP Tau, GK Tau, GI Tau, and HK Tau. Nevertheless, the recent observations of misaligned inner disks support the possibility that this systematic also has a physical basis.
Another systematic uncertainty on the stellar inclination is produced by the use of conversion tables between spectral type and effective temperature, which are still not precise for young stars. The use of different conversions can lead to discrepancies of about 150 K, which for a late-type star with T eff ∼ 3000-3500 K means an additional uncertainty of ∼5%. Propagated to R * and sin i, this final uncertainty of ∼10% corresponds to about 5 • for i * < 30 • or 10 • for i * > 30 • . Thus, the stellar inclinations might be higher than derived by about 10 • .
No extreme cases of face-on dippers, as reported by Ansdell et al. (2016b), are observed in this sample. The stellar inclination was derived for 11 dippers (50%), and inclinations of the outer disk are available for 8 (36%) dippers. The dippers for which it was possible to derive a stellar inclination are compatible with the disk warp scenario (thus 9/11), with the exception of HP Tau, whose stellar parameters might be significantly affected by a strong overall photometric variability, and V807 Tau, for which a different scenario is required, even after decomposing the primary and secondary brightness of the system. In general,

Temperature at corotation and dust survival
Given the assumption that dusty material at corotation causes the dips, it is interesting to compare the temperature at cototation with dust sublimation temperatures in order to confirm or exclude certain characteristics of the dust grains. The sublimation temperature, T sub , of the grains depends on the gas pressure, latent heat, and molecular weight more than on the stellar parameters (Kobayashi et al. 2011). The full range of T cor in Table 5 extends from 800 K to 1600 K, with error bars of about 100 K. The only materials able to withstand these temperatures are olivine, pyroxene, and iron (Pollack et al. 1994, see their Table 3). Moreover, the sublimation temperature increases with the gas density; stars with T cor > 1400 K exclude gas densities < 10 −8 g cm −3 . The stars with the highest temperatures at corotation in Table 5 can only host dust for even higher gas densities > 10 −6 g cm −3 . Because it is currently difficult to resolve the inner disk rim for low-mass YSOs, models for this region have only been developed for intermediate-mass stars (e.g., Isella & Natta 2005;Tannirkulam et al. 2007). It is thus of interest to better constrain these properties for CTTSs in future studies to verify the conditions for dust survival. Setting a critical temperature for dust T subl = 1500 K, it is possible to derive from Eq. 5 and Eq. 4 the minimum stellar rotation period for a given mass and radius that permits a sufficiently cold environment for dust to survive at corotation as Fig. 10 shows that all periods of dipper stars are well above this limit. For the few stars close to the limit, the condition is satisfied if a T subl = 1600 K is assumed. Previous observations of dippers in ρ Oph and Upper Sco (Bodman et al. 2017) also indicate that it is possible for dust to survive at the corotation radius of dipper stars.

Dippers in Taurus and in other clusters
The star-forming regions that are known so far to host dippers include NGC 2264, which is 2-3 Myr old Cody et al. 2014;McGinnis et al. 2015;Stauffer et al. 2015), Orion, which is 3 Myr old (Morales-Calderón et al. 2011), ρ Oph and Upper Sco, which are 1 and 10 Myr old, respectively (Ansdell et al. 2016a;Hedges et al. 2018;Cody & Hillenbrand 2018;Rebull et al. 2018), and Taurus (Rebull et al. 2020). A feature that is shared by all of these regions is the late spectral type, K to M, of the identified dippers. A common explanation for this is the longer pre-main-sequence phase of late-type stars, which allows the disk to be present for a longer time, thus increasing the probability of observing a dipper. Moreover, the lower surface temperature enables dust to be present close to the stellar magnetosphere. In Taurus, as in ρ Oph and Upper Sco, the occurrence of M-type dippers is higher than that of K-type dippers because stars of lower mass always exist in a larger number.
For the fractional rate at which dippers are found among YSOs, the statistics differ from region to region. The main reasons for this are the different counting methods (i.e., fraction of disk-bearing stars vs. all members of the cluster), the different selection techniques, and the target selection process of the survey, which is not necessarily optimized for finding dippers. For ρ Oph and Upper Sco, dippers represent about 20% of disked YSOs (Hedges et al. 2018), which was increased to 30% by Cody & Hillenbrand (2018). These authors also pointed out that in the older Upper Sco, a higher dipper fraction is observed among disked stars. Alencar et al. (2010) claimed that 30% to 40% of disked stars in NGC 2264 might be dippers, depending on whether a sample of thick or anemic inner disks is considered. This amount was found to be 20% by Cody et al. (2014), who also studied NGC 2264. An exception is Orion, where only 5% of the stars are dippers. This might be due to the photometric accuracy of the survey and a biased target selection (Morales-Calderón et al. 2011), as well as the fact that they did not consider a disked-star sample. In Taurus, considering the additional 12 dippers in Table 2, the occurrence rate is 19% of all the members observed with K2, and 31% of the members hosting a disk. The ratio of quasiperiodic to aperiodic dippers is about 1:1 for NGC 2264, ρ Oph, and Upper Sco. In Orion, only one-third of the dippers is quasiperiodic. In contrast, most dippers appear to be quasiperiodic in Taurus. The sample analyzed in this study might be unusual in this sense (only 3 aperiodic and 3 uncertain dippers out of 22 in Table 4); it must be noted that for only 2 of the 12 additional dippers in Table 2 the reported periods are most probably linked to quasiperiodic dips. The predominant periodicity in the light curve for the remainder is caused by cold stellar spots, while the extinction events are stochastic. Only 17 out of the 34 presented dippers are strictly quasiperiodic, which brings the ratio to about 1:1. The periods marked as uncertain in Table 4 are not counted as quasiperiodic here.

Dip width and period
Previous studies investigated whether dip width, dip amplitude, and/or dip width are correlated for dippers. No such correlation has been found so far (e.g., Bodman et al. 2017), except for the class of short-period narrow dips presented by Stauffer et al. (2015). For dip width and amplitude, a dependence would be a constraint on the geometrical properties of the dusty warp as presented by Bouvier et al. (1999). The amplitude of the eclipses depends not only on the viewing angle, but also on the vertical   Table 4. A logarithmic fit a log x + b is performed just on the data points (red line), and for upper and lower limits of the dip width (orange interval). The double-dipped stars are represented as triangles. The size of the markers is proportional to the square root of the dip amplitude (between 0.05 and 1.3 mag). Two stars for which the dip width derived in this way is not reliable are not represented in this plot. The stellar inclination in degrees is coded as colors of the data points, when available. Although the exact coefficients of the fitting function cannot be constrained from this small sample, W correlates with ln P with Pearson's r = 0.77. No significant correlation with the dip amplitude can be measured. extent of the dusty material occulting the star. On the other hand, the dip width delivers an estimate of its azimuthal extent.
For our sample of Taurus dippers, the quasiperiodic dippers for which it was possible to derive a reliable width are shown in Fig. 11. The data suggest a logarithmic trend of the dip width and the period, with a high scatter due to the intrinsic variability of dippers. It is remarkable that many dippers are obscured 50% of the time, suggesting that the azimuthal extent of the dusty structure must be larger than 180 • around the star. According to the geometrical constraints of a dusty warp, the observed dip width should grow with the inclination. The inclination values listed in Table 6 are represented in Fig. 11 as a color code of the data points. The amplitude of the dips is proportional to the size of the data points as 2 A A min , where A min is the smallest amplitude of the sample. No correlation is observed within the inclination and the dip width, nor between inclination and dip amplitude. Toward the higher end of the widths, the dip is often double-or multi-peaked, and its shape and position vary from phase to phase. This is the case for HP Tau, LkCa 15, GI Tau, JH 112 A, HK Tau, and IS Tau (see Sec. 3.4). The large error bars on the dip width for these double-or multi-peaked stars reflect the difference between considering the width of the full occultation and the width of the main dip.
Although the dependence can be best described with a logarithmic curve, the size of the sample is too small to draw a solid conclusion on the nature of the best-fitting curve. It is also possible that two different regimes are present, such as a linear trend for short periods, followed by a plateau for longer periods. The two quantities dip width (W) and the natural logarithm of the period (ln P) are correlated with a Pearson correlation coefficient r = 0.77 and a p-value of 2.8·10 −4 . A logarithmic, rather than lin- ear, trend is expected because the star cannot be eclipsing close to 100% of the time, thus flattening the curve for longer periods. The two stars in the lower left corner are crucial for this correlation, however. EPIC246989752 and JH 112 A (see Fig. F.1 for the dip width) fall in the category of short-duration quasiperiodic dips (Stauffer et al. 2015), with the interesting case of JH 112 A being a transient dipper. The derived dip width describes the dip as observed in the folded light curve well. Moreover, the upper error bar is very generous because the binned light curve +σ is almost flat. We can therefore exclude that the dip width for these two stars is overestimated. Stauffer et al. (2015) identified a particular class of Gaussian-shaped, narrow, and shortperiodic dips. The definition of "narrow" is set where the ratio of the FWHM of the fit Gaussian and the period is lower than 0.25, while the narrow dips in NGC 2264 have a FWHM-toperiod ratio of 0.15 on average. In Taurus, EPIC 246989752, and JH 112 A would correspond to this category, with a FWHM-toperiod ratio of 0.22 and 0.20, respectively. Both exhibit a Gaussian shape. The short-period dippers in NGC 2264 tend to have narrow eclipses for later spectral types (Stauffer et al. 2015).
Here, this is true for EPIC 246989752, which is an M dwarf, while JH 112 A has a spectral type of K5.5. We could not find a correlation between dip width and spectral type in the dipper sample.
Another star that can be associated with the population defined by Stauffer et al. (2015) is EPIC 247885481 (Sec. 3.3). It does not appear in Fig. 11 because the width of the dips cannot be constrained reliably in a standardized way. The amplitude is very small compared to the noise, and the dips do not clearly appear in each phase (see Sec. 3.3). The star is an M dwarf and is therefore compatible with the group identified in NGC 2264.
A correlation between R cor and the dip width in units of phase is presented in Fig. 12. Although weaker (Pearson's r = 0.63), a behavior similar to that of Fig. 11 is evident. This is expected because the corotation radius as derived in Eq. 4 is a function of the stellar period. In this case, the correlation is evident even without considering the two data points in the lower left corner, which only argue for a logarithmic instead of a linear dependence.
It is of interest to verify whether this correlation of dip width with period is an implication of existing models and if the observed dip widths are compatible with the mechanisms that have been proposed to explain dippers. In order to generate a quasiperiodic dipper, the most stable configuration of the magnetic field seems to be an inclined dipole in the framework of magnetospheric accretion (Romanova et al. 2013). Two broad, stable, and symmetric matter streams form from the disk and are accreted as funnel flows onto the stellar surface, producing an hotspot in each hemisphere. The observer can see the star obscured either by the dusty warp or even by the funnel flow itself. In this case, the extent of the magnetosphere nearly coincides with the corotation radius, R T ≈ R cor . In case of R cor < R T , the inner region of the disk rotates with ω > Ω K (where Ω K is the Keplerian velocity at corotation) and the star is in the propeller regime. The magneto-centrifugal forces lead to the ejection of disk material into an outflow (Romanova & Owocki 2015). During phases of enhanced mass accretion, in which the dipole component of the magnetic field decreases, the disk compresses the magnetosphere and reduces R T (Romanova et al. 2013). For R T < R cor , the star is in the accretion regime and the magnetosphere rotates more slowly than the inner disk (ω < Ω K ). In case of a small or no misalignment between the stellar spin axis and the magnetic moment, the matter from the inner disk accumulates at R T and the accretion instead occurs through unstable accretion tongues (Romanova et al. 2013), which are taller and thinner than funnel flows that reach the stellar surface close to equator. They can be observable as irregular hotspots in the UV. As a result, the light curve is rather stochastic and bursting. An unstable accretion regime has been invoked to explain aperiodic dippers (e.g., McGinnis et al. 2015).
The extent of the corotation radius directly depends on the stellar rotation period. For slow rotators, the corotation radius is located at a larger distance from the star, and the inner disk therefore rotates faster than the star. This facilitates the accretion through instabilities (Romanova et al. 2013;.
An observer would expect broad dips to be quasiperiodic because the unstable accretion tongues are rather narrow and tall and would not be able to obscure the star for a significant fraction of time . The aperiodic dippers observed in Taurus indeed exhibit rather narrow dips in their light curves.
The quasiperiods observed in dippers might also be linked to a configuration in the accretion regime. In this case, trapped density waves in the inner disk can produce a warp that rotates more slowly than the star and beyond the corotation radius, with a period that might vary over time. In this case, the simulations show that the inner disk is tilted and also rotates more slowly (Romanova et al. 2013). This could explain the changing, rather long period observed in the WPS of the K2 light curve of DK Tau, which is associated with an outer disk seen at low inclination. However, this would not suffice to explain why the star is obscured in the first part of the light curve. The slowly rotating warp should also have a smaller height than a fast-rotating warp, and DK Tau has an amplitude of 1 mag.
A correlation of period and dip width thus suggests that slow rotators are more probably occulted by more azimuthally extended dusty structures. This aspect has not been directly considered in simulations to date. In the scenario of magnetospheric accretion, this might imply a broader funnel flow for slow rotators.
For T Tauri stars of late spectral type, the internal structure is either fully convective or has small radiative cores (Gregory et al. 2012). The HR diagram of Taurus dippers shows that with the exception of LkCa 15, all stars are on their Hayashi tracks and fully convective. In this case, the magnetic field is expected to be axisymmetric and the dipole component dominates, as is the case for AA Tau. Vidotto et al. (2014) found a correlation of magnetic flux and rotation period for accreting PMS stars. This implies that stars with a simple magnetic field (i.e., dipole) and with the strongest magnetic fields are also the slowest rotators and truncate the disk at larger distances R T . The star-disk interaction therefore affects the magnetic field, which in turn affects the stellar rotation. We speculate that the large-scale magnetic field topology might cause the correlation of dip width and period. A stronger magnetic field might have a stabilizing effect on a large warp.

Conclusion
We have studied a sample of 179 YSOs in Taurus, which are members or possible members of the cluster (Rebull et al. 2020). We identify a total of 34 dippers, 22 of which are dippers not dominated by another type of variability. This makes up ∼ 20% of the high-confidence and possible Taurus members and ∼ 30% of the disk-bearing Taurus members as observed with K2.
The observations in Taurus highlight again how ephemeral dippers are. A striking example is the dipper prototype AA Tau, which is now in a fainter state and is no longer classified as a quasiperiodic, prototypical dipper. Dipper light curves can persist over timescales of a few years or even a few days. This strongly suggests that the observed occurrence rates are a lower limit to their true occurrence among CTTSs.
The ratio of quasiperiodic to aperiodic dippers is 1:1, although most of the dippers studied in detail are quasiperiodic. The large majority of aperiodic dippers observed in Taurus is dominated by another type of variability, mainly that due to cold spots. As found in other surveys, dipper stars are of late spectral type K or M because the lower surface temperature allows dust to survive in the inner disk and create the dips. The dipper stars of the sample are fully convective, low-mass (< 1 M ) stars down to the brown dwarf limit, with the presence of a probable debris disk dipper.
Supporting previous dipper surveys, the observed periods are in the range of the rotation periods of low-mass CTTSs. In the sample, a transient period, a changing period, and a dipper period very close to the period of a stellar spot are also identified. The temperatures derived at corotation are also compatible with dust survival, and the mass accretion rates are typical for CTTSs.
Many of the dippers for which it was possible to derive a stellar inclination are seen at a rather moderate inclination angle, and the outer disk inclination is, when available, slightly lower than the stellar inclination. Whether the stellar inclination is systematically higher than the outer disk inclination cannot be constrained accurately in this small sample, but it would be consistent with recently observed stars with a tilted inner disk.
Magnetospheric accretion is able to explain most but not all Taurus dippers. However, dusty disk winds and Rossby wave instabilities do not seem convincing as explanations for the dippers seen at low inclination in this sample.
A fraction of the quasiperiodic dippers exhibits double or complex dips, whose structure varies with time and whose minima are shifted from one cycle to the next. These examples are worth exploring in a future more precise dynamical analysis.
Finally, a correlation is found between the dip width in units of phase (i.e., the equivalent of the angular extent of the occulting structure) and the rotational period of the dippers. We spec-ulate that this might indicate that the accretion columns of more slowly rotating stars have a larger base, possibly reflecting a different magnetic topology. This dependence needs to be investigated for dippers in other clusters in future work.
The Fourier transform (FT) is an extension of the Fourier series, which expresses periodic functions as sums of sinusoids. It decomposes the signal in the frequency domain, and its power spectrum represents each frequency that is present in the original signal. Fourier analysis requires (nearly) evenly sampled data as a condition. Fourier analysis is not able to provide information about the time localization of a given frequency because it has no time resolution in the frequency domain. This is because the sinusoids employed to transform the signal are infinite in time, thus the time information is lost, and only the frequencies are retrieved.
The windowed Fourier transform can help in solving this issue by considering a window (i.e., a slice) of the signal f and transforming only that section, then translating the window along the signal. Given the window g(t), the windowed signal is defined as (A.2) In physics, g(t) is commonly chosen to be a Gaussian. The window is multiplied with a chop of the signal, transformed separately, and then shifted along the signal. If g is centered in time and frequency,f t (ω) corresponds to the information of f around time t and frequency ω. The WFT underlies the uncertainty principle: the lower bound for the product of time and frequency resolution is For this method, the quality of the resulting analysis strongly depends on the chosen size τ of the window function. This determines the time and the frequency resolution a priori and might be inconvenient if both a good time and good frequency resolution are required.
The windowed Fourier transform is constrained by the fixed size of the window, while the Fourier transform delivers no time resolution at all. This issue can be solved by wavelets, by adapting the time-width of the window to the frequency that is to be investigated, where CWT is the continuous wavelet transform. The mother wavelet ψ(t) is in most cases a sinusoidal wave confined in time, stretched and compressed according to the scale s and shifted with a time step τ along the time series. The normalization factor 1 √ s ensures that the transformed signal has the same energy at each scale. The 2-dimensional wavelet power spectrum (WPS) is produced by a convolution between the wavelet and the signal, which is recomputed for each scale along the time series.
The varying size of the wavelet allows us to identify both high-and low-frequency features in the signal. Both the discrete (DWT) and continuous wavelet transform (CWT) exist; here only the CWT is presented. As for the WFT, the wavelet transform is subject to the uncertainty principle. Based on its shape, a certain wavelet has either a higher time or frequency resolution. Being confined in the time domain results in a broader peak in the frequency domain, which would be a δ-function in the case of the sinusoid of the FT. The complex Morlet wavelet is a good compromise in this respect; moreover, its form is similar to the sought-for signal.
The complex Morlet wavelet used in this work is defined as a superposition of a complex wave and a Gaussian, such as where the parameters F c and T b are the central frequency and the bandwidth parameter, respectively, and x is a unitless time parameter. This formulation differs from that in Torrence & Compo (1998), as they assumed a Gaussian with unit variance, which can be varied here, and a frequency ω 0 = 6. The advantage of freely setting the bandwidth and the central frequency is that we can customize the time and frequency resolution of the resulting wavelet plot, in order to highlight the position of a periodicity or to focus on the frequencies. A broader Gaussian in the time domain, thus a lower time resolution, is represented by a narrower peak in the frequency domain, translating into higher frequency resolution.
Increasing the bandwidth T b , that is, the width of the Gaussian, in the time domain decreases the bandwidth in the Fourier domain, that is, delivers a higher frequency resolution by losing time resolution. The central frequency should be close to the frequency of the signal of interest. In order to create an atlas of K2 light curves, the parameters were set to F c = 1.0 and T b = 1.5 as default. For individual cases of interest, these parameters may be changed. In the WPS, a certain scale s corresponds to a frequency F s as The set of scales of the wavelet transform is equivalent to the custom frequency resolution of a periodogram. They can be set linearly to extend from s 0 to s max , where s max corresponds to the longest period of interest according to Eq. A.6. For a faster and more efficient computation of the WPS, it is better to choose the set of scales according to a power law, s = s 0 2 jδ j , j = 0, 1, ..., J, δ j = log 2 s max J . (A.7) Here J is the total number of scales and s 0 should be chosen close to 2dt. The cone of influence (COI) defines the range within which the WPS is affected by edge effects, which grow with the scale, as the wavelet is more and more stretched in time. The power spectrum contained in this region should not be regarded as significant. Considering that for the Morlet wavelet the Fourier period λ is almost equivalent to the scale, as λ = 1.03s (Torrence & Compo 1998), the COI is plotted as √ 2p, where p is the period corresponding to a frequency. A central frequency different than 1.0 affects the COI as √ 2pF c .

Appendix B: Light curves of the dipper sample
Appendix C: Extinction leads to an increase in A V up to a factor 5 because of the infrared excess. The uncertainty on A V is set as the root-meansquare (rms) of the A V derived in this study compared with the literature. A minimum reasonable uncertainty of 0.3 mag is applied. Empirically, the uncertainty does not grow with extinction.

Appendix D: Notes on individual objects
In the following section, remarks about individual objects are presented. Three stars (St 34, HD 285893, and 2MASSJ05023985+2459337) are possible Taurus members (Rebull et al. 2020).
HD 285893 According to the HR diagram, the star lies almost on the ZAMS, which means that it is much older than the Taurus population. There are several possibilities to explain this. The star might not be a Taurus member (it is considered a possible member), the spectral type might not be correctly estimated, and/or the star might be seen edge-on, thus its luminosity would be strongly underestimated. The star shows a clear dipper behavior in the light curve, which is unexpected for a spectral type F8 because the higher temperature of the star depletes the inner disk of dust at distances much larger than corotation. The object is poorly studied in the literature and does not appear in more recent spectroscopic surveys. Its behavior appears to be related to the so-called 'old' dippers, which can host a debris disk and whose occultations might be caused by disrupting planetesimals (Gaidos et al. 2019;Tajiri et al. 2020). The dips are narrow, as is commonly the case for aperiodic dippers. The light curve continuum is unstable and aperiodic.
CFHT Tau 12 No binarity information has been reported for this star. The star is quasiperiodic (period agrees with Scholz et al. 2018), but there is no information about v sin i from which we could derive an inclination. The light curve never reaches a continuum, and no periodicity for this variation can be derived within the time window of the observations. The shape of the dips clearly distinguishes the dipper phenomenon from spots. The folded light curve is noisy.
2MASSJ04384725+1737260 No companion is known in the literature. The star never clearly reaches the continuum; the periodicity is pretty clear, but some dips seem to be aperiodic (see folded light curve in Fig. F.1). There might also be one or more spots.
2MASSJ04295950+2433078 No binarity is known (Davies et al. 2014). It has a high variability amplitude (0.66 mag). The sinusoidal pattern of the light curve looks rather spot-like, but the star is classified as a dipper based on the irregular shape of the dips. The WPS shows a significant periodicity at ∼ 7.5 d, which is much weaker in the periodogram. This might be due to the shape of the light curve, which is well-matched to the Morlet  Notes. Only the extinction measured by Luhman et al. (2017) with optical or CTTS methods is reported here. The veiling is measured at 7510 Å by Herczeg & Hillenbrand (2014). wavelet and is therefore highly correlated. Folded at 7 d, a periodic pattern of several dips appears. A harmonic is excluded.
2MASSJ05023985+2459337 No companion is known in the literature. The light curve exhibits both a spot and a dipper pattern. A dedicated discussion is presented in Sec. 3.3.
ITG 34 is known to be single (Davies et al. 2014;Akeson et al. 2019). The amplitude of the dips is large (0.39 mag). According to the light curve, the star either almost never reaches the continuum level or a few bursts occur during the K2 campaign. The periodogram shows two double-peaked periodicities. The main periodicity is also more stable over time according to the WPS. The second periodicity cannot be a harmonic of the first. The folded light curve is noisy, and no new features emerge when it is folded at a different period. The apparent beat cannot be explained with binarity.
JH 223 is the best example of a quasiperiodic dipper in the sample. The star is a binary with a separation of 2 . Both components appear to have a disk (Itoh et al. 2015). The star is reported in Kraus et al. (2012) as having a disk, but without evidence for accretion.
IS Tau is a binary with 0.2 separation (Schaefer et al. 2014). It is not known which (or if both) component hosts a disk. Watson et al. (2009) reported a very high mass fraction (80%-100%) of crystalline silicates in the inner disk. The light curve is characterized by a combination of broad and narrow eclipses. The strongest periodicity is given by the broad dips. No Gaia parallax is available. The star is also located on the L1495 filament of Taurus, which contains members of both the near and the far population (Fleming et al. 2019). In the HR diagram in Fig. 9, the possible luminosity of IS Tau is plotted according to mean distances of 130.6 and 160.2 pc (Fleming et al. 2019). The age is more consistent with that of the remaining sample if IS Tau is a member of the far population. The star shows a double dip; the primary has a larger, quite constant amplitude. For this star, the phase shift of the dips varies more than for other double-peaked dippers.
St 34 is a multiple system (Aab+B) with separation 1.18 , surrounded by a circumbinary, transitional disk (Rigliaco et al. 2015). The star is reported as SB2 by Akeson et al. (2019). The inner disk appears to be depleted, although some dust is still present (Hartmann et al. 2005 (White & Hillenbrand 2005) would lead to an age of 25 Myr, thus the authors must assume a distance of only 90 pc. In the HR diagram of our study (Fig. 9), the star appears to have an age of 1-2 Myr, given the limitations on the luminosity due to photometric variability. The Gaia parallax reports a distance of 143 pc, thus showing that St 34 cannot be a foreground star. Accordingly, the age-derivation based on the lithium depletion timescale apparently does not apply to this star. The light curve is irregular, and a period forest is present in both the periodogram and WPS. The presented period is uncertain. Some contamination might arise from multiplicity.
V807 Tau is a multiple system (A + Bab) and a separation of 0.3 for the wide components, 0.04 for the secondary (Schae-  The Gaia distance of V807 Tau is 113 pc, thus it is maybe underestimated. The system is multiple, which means that there might be an issue with the Gaia parallax. The star is located in the B18 region of Taurus (Fleming et al. 2019), which has a distance of 110-150 pc. This value has to be considered with caution, but it is not an outlier for this region. Rodriguez et al. (2017) derived a period of 0.809 d, which is inconsistent with K2 data. Schaefer et al. (2012) decomposed the brightness of the single components and obtained a primary star 0.5 mag fainter than the 2MASS J measurement. This leads to a much smaller stellar radius and higher, more plausible stellar inclination. The light curve presents an apparently sinusoidal long-term modulation; moreover, it is unclear if the star reaches its maximum brightness at all. The modulation does not appear to be related to binarity because the binary system has a period of 12 yr. There is only one clear periodicity, and the shape of the dips does not suggest stellar spots.
Haro 6-37 is a triple system. The close binary AB (0.3 Duchêne et al. 1999) is separated from C by 2.6 (White & Ghez 2001). The light curve presents unusually broad dips for a dipper with high amplitude (0.5 mag). The periodicity is unclear in the WPS and periodogram. Our luminosity is highly different from that of Herczeg & Hillenbrand (2014) (here 0.8L , there 0.07L ). The derived inclination suggests that the period is either incorrect or that the occultations have another origin.   Table C.1. Red points show stars with high veiling. For these, the spectroscopic measurement by Herczeg & Hillenbrand (2014) is preferred. The dashed line represents A V (V − J) = A V (lit). The stars with the largest scatter are strongly veiled, and the extinction measured with photometry is less reliable. For GO Tau (four black symbols at the bottom, close to A V (V − J) = 0), the A V (V − J) value is discarded because it is inconsistent with the literature. Another outlier is the A V of V807 Tau as measured by KH95; this value is discarded for the same reason.