# The Disk Substructures at High Angular Resolution Project (DSHARP). II. Characteristics of Annular Substructures

Jane Huang, Sean M. Andrews, Cornelis P. Dullemond, Andrea Isella, Laura

M. Pérez, Viviana V. Guzmán, Karin I. Öberg, Zhaohuan Zhu, Shangjia

Zhang, Xue-Ning Bai, et al.

► To cite this version:

Jane Huang, Sean M. Andrews, Cornelis P. Dullemond, Andrea Isella, Laura M. Pérez, et al. The Disk Substructures at High Angular Resolution Project (DSHARP). II. Characteristics of Annular Substructures. The Astrophysical Journal Letters, 2018, 869, 10.3847/2041-8213/aaf740. insu-03635023

HAL Id: insu-03635023

https://insu.hal.science/insu-03635023v1

Submitted on 1 Sep 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# The Disk Substructures at High Angular Resolution Project (DSHARP). II. Characteristics of Annular Substructures

Jane Huang¹, Sean M. Andrews¹, Cornelis P. Dullemond², Andrea Isella³, Laura M. Pérez⁴, Viviana V. Guzmán⁵,⁶, Karin I. Öberg¹, Zhaohuan Zhu⁷, Shangjia Zhang⁷, Xue-Ning Bai⁸, Myriam Benisty⁹,¹⁰, Tilman Birnstiel¹¹, John M. Carpenter⁵, A. Meredith Hughes¹², Luca Ricci¹³, Erik Weaver³, and David J. Wilner¹

¹Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA; jane.huang@cfa.harvard.edu

²Zentrum für Astronomie, Heidelberg University, Albert Ueberle Str. 2, D-69120 Heidelberg, Germany

³Department of Physics and Astronomy, Rice University, 6100 Main Street, Houston, TX 77005, USA

⁴Departamento de Astronomía, Universidad de Chile, Camino El Observatorio 1515, Las Condes, Santiago, Chile

⁵Joint ALMA Observatory, Avenida Alonso de Córdova 3107, Vitacura, Santiago, Chile

⁶Instituto de Astrofísica, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, 7820436 Macul, Santiago, Chile

⁷Department of Physics and Astronomy, University of Nevada, Las Vegas, 4505 S. Maryland Pkwy, Las Vegas, NV 89154, USA

⁸Institute for Advanced Study and Tsinghua Center for Astrophysics, Tsinghua University, Beijing 100084, People's Republic of China

⁹Unidad Mixta Internacional Franco-Chilena de Astronomía, CNRS/INSU UMI 3386, Departamento de Astronomía, Universidad de Chile, Camino El Observatorio 1515, Las Condes, Santiago, Chile

¹⁰Univ. Grenoble Alpes, CNRS, IPAG, F-38000 Grenoble, France

¹¹University Observatory, Faculty of Physics, Ludwig-Maximilians-Universität München, Scheinerstr. 1, D-81679 Munich, Germany

¹²Department of Astronomy, Van Vleck Observatory, Wesleyan University, 96 Foss Hill Drive, Middletown, CT 06459, USA

¹³Department of Physics and Astronomy, California State University Northridge, 18111 Nordhoff Street, Northridge, CA 91130, USA

Received 2018 October 8; revised 2018 November 19; accepted 2018 November 28; published 2018 December 26

## Abstract

The Disk Substructures at High Angular Resolution Project (DSHARP) used ALMA to map the 1.25 mm continuum of protoplanetary disks at a spatial resolution of ~5 au. We present a systematic analysis of annular substructures in the 18 single-disk systems targeted in this survey. No dominant architecture emerges from this sample; instead, remarkably diverse morphologies are observed. Annular substructures can occur at virtually any radius where millimeter continuum emission is detected and range in widths from a few astronomical units to tens of astronomical units. Intensity ratios between gaps and adjacent rings range from near-unity to just a few percent. In a minority of cases, annular substructures coexist with other types of substructures, including spiral arms (3/18) and crescent-like azimuthal asymmetries (2/18). No clear trend is observed between the positions of the substructures and stellar host properties. In particular, the absence of an obvious association with stellar host luminosity (and hence the disk thermal structure) suggests that substructures do not occur preferentially near major molecular snowlines. Annular substructures like those observed in DSHARP have long been hypothesized to be due to planet-disk interactions. A few disks exhibit characteristics particularly suggestive of this scenario, including substructures in possible mean-motion resonance and “double gap” features reminiscent of hydrodynamical simulations of multiple gaps opened by a planet in a low-viscosity disk.

**Key words:** planets and satellites: formation – protoplanetary disks – techniques: high angular resolution

**Supporting material:** machine-readable table

## 1. Introduction

Exoplanet surveys have uncovered remarkable diversity in the masses, radii, and composition of exoplanets, as well as in the overall architectures of exoplanetary systems (e.g., Howard et al. 2010; Borucki et al. 2011; Winn & Fabrycky 2015; Fulton et al. 2017). These wide variations in properties are thought to be closely tied to differences in where and when planets formed in the protoplanetary disks orbiting young stars (e.g., Öberg et al. 2011b; Mordasini et al. 2012; Lee & Chiang 2016). Recent discoveries of complex features in protoplanetary disks, such as annular substructures, spiral arms, and high-contrast azimuthal asymmetries, have fueled interest in using observed disk morphologies to infer the properties of an unseen population of newly forming planets (e.g., Muto et al. 2012; van der Marel et al. 2013; ALMA Partnership et al. 2015; Isella et al. 2016).

Of these various types of small-scale features, annular substructures have been reported most often. They have predominantly been detected in millimeter continuum emission tracing larger (~millimeter-sized) dust grains in the disk midplane (e.g., Andrews et al. 2016; Fedele et al. 2018; Long

et al. 2018a), but have also been observed in scattered light tracing sub-micron-sized dust grains in the disk atmosphere (e.g., van Boekel et al. 2017; Avenhaus et al. 2018; Muro-Arena et al. 2018) and in molecular emission (e.g., Isella et al. 2016; Fedele et al. 2017; Teague et al. 2017). The detections of these annular substructures are particularly intriguing with respect to planet formation theory because numerical simulations have long predicted that a sufficiently massive protoplanet will trigger density waves that eventually shock and create disk gaps (e.g., Papaloizou & Lin 1984; Bryden et al. 1999; Nelson et al. 2000).

However, it has been unclear whether annular substructures (or any other kind of substructure, for that matter) necessarily arise from planet-disk interactions, and a number of other hypotheses have been proposed to explain their origins. One set of hypotheses suggests that condensation of certain molecular species may affect how dust grains fragment, grow, and drift, in turn causing dust to accumulate in proximity to various snowlines (e.g., Banzatti et al. 2015; Zhang et al. 2015; Okuzumi et al. 2016; Pinilla et al. 2017). Another group of proposed mechanisms explores the effects of coupling between magnetic fields and disk material, such as zonal
flows arising from magnetorotational instability (MRI) turbulence (e.g., Johansen et al. 2009; Bai & Stone 2014; Simon & Armitage 2014), surface density enhancements at the edges of dead zones (e.g., Flock et al. 2015; Lyra et al. 2015; Ruge et al. 2016), and the spontaneous concentration of magnetic flux (e.g., Bai & Stone 2014; Béthune et al. 2017; Suriano et al. 2018). Other proposed mechanisms have focused on instabilities arising from interactions between gas and dust, including secular gravitational instabilities resulting from friction between gas and dust (e.g., Takahashi & Inutsuka 2014) or dust-driven viscous ring-instability resulting from dust drift altering the viscosity profile of the disk and amplifying density perturbations (Dullemond & Penzlin 2018).

While much of our current knowledge of disk substructure comes from serendipitously discovered case studies, establishing the origins of annular substructures in disks can be facilitated by a systematic analysis of their properties. To that end, we undertook the Disk Substructures at High Angular Resolution Project (DSHARP), the first high angular resolution ALMA survey of disks (Andrews et al. 2018a). Substructures are detected in the 1.25 mm continuum emission of all 20 sources, lending weight to the notion that millimeter continuum substructures are common and perhaps even ubiquitous. In this work, we present an analysis of the annular substructures observed in single-disk systems. Analysis of the two multiple-disk systems in DSHARP is presented in Kurtovic et al. (2018). Section 2 provides a brief overview of the disk sample. In Section 3, we measure the locations, widths, and contrasts of the substructures. In Section 4, we examine trends in their positions and sizes. Section 5 comments on possible origins for the substructures and discusses implications for their prevalence in the general disk population. Section 6 summarizes the main findings.

## 2. Sample Overview

DSHARP targeted 20 systems, 18 of which are single-disk systems (i.e., all sources except AS 205 and HT Lup). In the context of this Letter, “the DSHARP sample” refers specifically to these 18 disks. Angular resolutions range from ~30 to 60 mas, probing scales of ~5 to 8 au. Stellar properties, including distances, masses, accretion rates, ages, and luminosities, are listed in Table 1 of Andrews et al. (2018a). The analysis in this work is based on the “fiducial” images presented in Andrews et al. (2018a), which also provides details about the observational setup and the calibration and imaging procedure. Table 4 of Andrews et al. (2018a) provides basic information about these images, including synthesized beam sizes, peak intensities, continuum flux densities, and image rms noise levels.

Some portions of the analysis fold in publicly available continuum observations of HL Tau and TW Hya taken at comparable spatial resolution and sensitivity. The TW Hya analysis uses the 1.0 mm image presented in Huang et al. (2018a) and the HL Tau analysis uses the 1.0 mm image presented in ALMA Partnership et al. (2015). Although 1.3 mm continuum observations are available for both sources (ALMA Partnership et al. 2015; Tsukagoshi et al. 2016), the 1.0 mm images are used instead because of their better angular resolution and uv coverage.

## 3. Disk Features

### 3.1. Radial Locations of Annular Substructures

We use the term “annular substructures” to refer collectively to the “bright” and “dark” annuli visible in the DSHARP 1.25 mm continuum images. Colloquially, the “bright” and “dark” features are often referred to as “rings” and “gaps,” respectively (e.g., ALMA Partnership et al. 2015; Isella et al. 2016; Fedele et al. 2017). Figure 1 shows the 1.25 mm continuum images of the DSHARP sources with labeled annular substructures. The three-part procedure used to identify and label these substructures is explained below.

First, many of the emission features manifest as well-separated series of ellipses. The positions of such distinct substructures are measured in a manner similar to that employed for the HL Tau disk in ALMA Partnership et al. (2015). We perform a preliminary deprojection of each disk image using estimates of the disk orientation from either a Gaussian fit to the image or from measurements of previous ALMA observations of certain axisymmetric sources (Fedele et al. 2018 for AS 209, Dipierro et al. 2018 for Elias 24, and Andrews et al. 2018b for Sz 114 and GW Lup). For each substructure of interest, the radial intensity profile is measured along a series of evenly spaced azimuthal angles in the deprojected image. The azimuthal angle spacing is set to be on the order of a spatial resolution element—that is, if the FWHM of the major axis of the synthesized beam is $b_{maj}$ and the feature occurs roughly at radius $r$ in units of arcseconds (estimated from visual inspection), then the azimuthal angle spacing in degrees is $\Delta\theta \sim 360^{\circ} \times b_{maj}/(2\pi r)$ (with some minor adjustment so that $\Delta\theta$ is a divisor of 360). In each of the $n$ intensity profile cuts, we search for either a local maximum (if the substructure of interest is a bright ring) or minimum (if the substructure is a gap), generally using a search window of ~6 to 20 au around $r$ depending on the apparent width of the feature. If an extremum cannot be identified in a given intensity profile cut due to insufficient angular resolution arising from the disk projection or to low signal-to-noise ratios (S/Ns), that azimuthal angle is excluded. In addition, the azimuthal angles where the spiral arms cross the gap ($-117^{\circ} < \theta < -76^{\circ}$ and $68^{\circ} < \theta < 123^{\circ}$) are excluded for the Elias 27 disk, and the azimuthal angles spanned by the bright crescent on the southeast side of the HD 143006 disk ($80^{\circ} < \theta < 144^{\circ}$) are excluded when fitting its outermost ring. The azimuthal angle conventions are explained in Appendix A.

The $n$ locations of the sampled extrema are mapped back to pixel coordinates for the original image. An ellipse is then fit to the $n$ coordinates extracted for a given annular substructure. The model ellipse is defined by five parameters: the offset of the ellipse center from the phase center ($\Delta x$ and $\Delta y$), the semimajor axis of the ellipse ($r_0$), the length ratio of the semiminor to semimajor axis (i.e., the cosine of the disk inclination, $\cos i$), and the position angle (P.A.) (the angle east of north by which the major axis is rotated). The absolute offsets are not physically significant because phase centers were sometimes adjusted to align execution blocks (see details in Andrews et al. 2018a), but they are necessary for subsequent calculations of the radial profile. The quality of the fit of the model ellipse to the measured coordinates of the annular substructure is assessed by computing the orthogonal distances $d_i$ between the measured coordinates and the model ellipse. The orthogonal distance between a point $(x, y)$ and a curve $f(x, y)$
<figure><img src="image_3.png" /><figcaption>**Figure 1.** ALMA 1.25 mm continuum images of the DSHARP sample, ordered by decreasing stellar luminosity (see Table 1 of Andrews et al. 2018a) from left to right and top to bottom. An arcsinh stretch is applied to the color scale of each disk to increase the visibility of substructures in the outer regions. Axes are labeled with angular offsets from the disk center. Annular substructures are marked with dotted arcs (gaps) or solid arcs (bright emission rings). The synthesized beam is shown in the lower left corner of each panel.</figcaption></figure>

is the minimum Euclidean distance between the point and the curve. Formulae for calculating orthogonal distances from an ellipse are provided in Zhang (1997).

The log-likelihood takes the form

$$ \ln L = -\frac{1}{2} \sum_{i=1}^{n} \left[ \frac{d_i^2}{\sigma^2} + \ln(2\pi\sigma^2) \right], \quad (1) $$

where $\sigma$ is the standard deviation of the orthogonal distances. It is not straightforward to estimate $\sigma$ directly from the image; although it is tied to the angular resolution, it also depends in some fashion on the S/N of the substructure emission. Thus, in addition to fitting for the five parameters defining the position and orientation of the ellipse, we allow $\ln \sigma^2$ to be a free parameter. Uniform priors are specified for $\Delta x$, $\Delta y$, and $r_0$ over the range
<figure><img src="image_4.png" /><figcaption>Figure 1. (Continued.)</figcaption></figure>

($-4''5, 4''5$) based on the typical image size. A uniform prior is also specified for $\ln \sigma^2$ over the range ($-15, -5$), under the assumption that the synthesized beam roughly sets the largest possible value of $\sigma$ and the pixel size roughly sets the lowest. In practice, the best-fit values for $\ln \sigma^2$ generally end up being between $-10$ and $-8$, with uncertainties of a few tenths. Gaussian priors for $\cos i$ and P.A. are specified for the aforementioned axisymmetric sources with previous ALMA measurements. Otherwise, a uniform prior for $\cos i$ in the range (0, 1) and P.A.

in the range (0°, 180°) are adopted (i.e., the smaller angle east of north is adopted as the P.A., and the direction of the angular momentum vector is not considered because it cannot be determined unambiguously for many of the disks). For HD 143006 and WaOph 6, the bounds for the P.A. are set to (90°, 270°) because the major axis of the disk appears to be relatively close to 180° (e.g., Cox et al. 2017; Benisty et al. 2018).

The posteriors are explored with the affine invariant MCMC sampler implemented in `emcee` (Goodman &
<figure><img src="image_5.png" /><figcaption>Figure 1. (Continued.)</figcaption></figure>

Weare 2010; Foreman-Mackey et al. 2013). Each ensemble consists of 30 walkers and proceeds for 20,000 steps, with the first 500 steps discarded as burn-in. Convergence of the MCMC chains is assessed by checking that the number of steps significantly exceeds the autocorrelation time (typically on the order of $10^2$).

The posterior medians for each well-resolved substructure are listed in Table 1, with error bars calculated from the 16th and 84th percentiles. If the annular substructure's coordinates were extracted by searching for local minima in the radial

profile cuts, it is labeled with the prefix “D” (for “Dark”), followed by the semimajor axis of the best-fit ellipse rounded to the nearest whole number of astronomical units. Similarly, if the annular substructure’s coordinates were extracted by searching for local maxima, it is labeled with the prefix “B” (for “Bright”), again followed by the semimajor axis of the best-fit ellipse rounded to the nearest whole number of astronomical units. This naming convention is adapted from that used for HL Tau in ALMA Partnership et al. (2015). However, instead of naming the substructures based on the
<table><caption>**Table 1**
Properties of Annular Substructures</caption><thead><tr><th>Source<br>(1)</th><th>Feature<br>(2)</th><th>Δx<br>(mas)<br>(3)</th><th>Δy<br>(mas)<br>(4)</th><th>r<sub>0</sub><br>(mas)<br>(5)</th><th>r<sub>0</sub><br>(au)<br>(6)</th><th>Incl.<br>(°)<br>(7)</th><th>P.A.<br>(°)<br>(8)</th><th>Method<br>(9)</th><th>Width<br>(au)<br>(10)</th><th>Depth<br>(11)</th></tr></thead><tbody><tr><td rowspan="15">AS 209</td><td>D9</td><td>3 ± 1</td><td>-3.9 ± 0.9</td><td>71.8 ± 0.8</td><td>8.69 ± 0.11</td><td>35.6<sup>+0.7</sup><sub>-0.8</sub></td><td>85.7 ± 0.7</td><td>E</td><td>4.7 ± 0.2</td><td>0.45 ± 0.02</td></tr><tr><td>B14</td><td>2.0 ± 1.3</td><td>-4.3<sup>+1.2</sup><sub>-1.1</sub></td><td>117.6 ± 1.1</td><td>14.2 ± 0.14</td><td>34.9<sup>+0.7</sup><sub>-0.8</sub></td><td>85.8 ± 0.7</td><td>E</td><td>8.9 ± 0.2</td><td>...</td></tr><tr><td>D24</td><td>2.7 ± 1.4</td><td>-3.8 ± 1.2</td><td>197.0 ± 1.3</td><td>23.84 ± 0.15</td><td>35.6 ± 0.7</td><td>85.8 ± 0.7</td><td>E</td><td>3.4 ± 0.3<sup>a</sup></td><td>0.895 ± 0.013</td></tr><tr><td>B28</td><td>3 ± 2</td><td>-4 ± 2</td><td>230 ± 2</td><td>27.8 ± 0.3</td><td>35.1 ± 0.7</td><td>86.1 ± 0.7</td><td>E</td><td>4.7 ± 0.3</td><td>...</td></tr><tr><td>D35</td><td>4.7 ± 1.7</td><td>-3.5 ± 1.5</td><td>289.6 ± 1.6</td><td>35.04 ± 0.18</td><td>35.2 ± 0.6</td><td>85.4 ± 0.7</td><td>E</td><td>3.0 ± 0.3<sup>a</sup></td><td>0.97 ± 0.01</td></tr><tr><td>B39</td><td>4 ± 2</td><td>-3 ± 2</td><td>320 ± 2</td><td>38.7<sup>+0.3</sup><sub>-0.2</sub></td><td>34.9 ± 0.7</td><td>85.7 ± 0.7</td><td>E</td><td>3.4 ± 0.3<sup>a</sup></td><td>...</td></tr><tr><td>D61</td><td>-2 ± 3</td><td>-4 ± 3</td><td>503 ± 3</td><td>60.8 ± 0.3</td><td>35.8 ± 0.6</td><td>86.1 ± 0.6</td><td>E</td><td>15.5 ± 0.3</td><td>0.014 ± 0.011</td></tr><tr><td>B74</td><td>-0.8 ± 1.1</td><td>-4.4 ± 1.0</td><td>612.8 ± 1.3</td><td>74.15 ± 0.12</td><td>34.92 ± 0.3</td><td>85.7 ± 0.4</td><td>E</td><td>9.3 ± 0.3</td><td>...</td></tr><tr><td>D90</td><td>3 ± 3</td><td>-3 ± 3</td><td>742.6 ± 3</td><td>89.9 ± 0.3</td><td>34.4 ± 0.5</td><td>86.0 ± 0.6</td><td>E</td><td>4.5 ± 0.6</td><td>0.48 ± 0.13</td></tr><tr><td>B97</td><td>1 ± 5</td><td>-3 ± 4</td><td>799 ± 4</td><td>96.7 ± 0.5</td><td>34.5 ± 0.6</td><td>85.2 ± 0.7</td><td>E</td><td>8.1 ± 0.6</td><td>...</td></tr><tr><td>D105</td><td>1 ± 3</td><td>0 ± 3</td><td>872 ± 3</td><td>105.5 ± 0.4</td><td>35.0 ± 0.5</td><td>85.6 ± 0.6</td><td>E</td><td>14.7 ± 0.4</td><td>0.016 ± 0.009</td></tr><tr><td>B120</td><td>0.2 ± 1.1</td><td>2.1 ± 1.0</td><td>992.9 ± 1.3</td><td>120.14 ± 0.12</td><td>34.92 ± 0.19</td><td>85.6 ± 0.3</td><td>E</td><td>11.2 ± 0.4</td><td>...</td></tr><tr><td>D137</td><td>...</td><td>...</td><td>1132 ± 8</td><td>137 ± 1</td><td>...</td><td>...</td><td>R</td><td>4.2 ± 1.4<sup>a</sup></td><td>0.90 ± 0.07</td></tr><tr><td>B141</td><td>...</td><td>...</td><td>1165 ± 8</td><td>141 ± 1</td><td>...</td><td>...</td><td>R</td><td>2.8 ± 1.4<sup>a</sup></td><td>...</td></tr><tr><td rowspan="5">DoAr 25</td><td>D74</td><td>...</td><td>...</td><td>~540</td><td>~74</td><td>...</td><td>...</td><td>V</td><td>...</td><td>...</td></tr><tr><td>B86</td><td>...</td><td>...</td><td>~620</td><td>~86</td><td>...</td><td>...</td><td>V</td><td>...</td><td>...</td></tr><tr><td>D98</td><td>40 ± 3</td><td>-495 ± 2</td><td>710 ± 3</td><td>98.0 ± 0.3</td><td>68.2 ± 0.3</td><td>110.9 ± 0.3</td><td>E</td><td>15.5 ± 0.5</td><td>0.22 ± 0.04</td></tr><tr><td>B111</td><td>34 ± 4</td><td>-493 ± 3</td><td>806 ± 4</td><td>111.3 ± 0.4</td><td>66.1<sup>+0.3</sup><sub>-0.4</sub></td><td>110.2 ± 0.4</td><td>E</td><td>14.3 ± 0.5</td><td>...</td></tr><tr><td>D125</td><td>...</td><td>...</td><td>906 ± 7</td><td>125 ± 1</td><td>...</td><td>...</td><td>R</td><td>10.0 ± 1.4</td><td>0.63 ± 0.05</td></tr><tr><td rowspan="2">DoAr 33</td><td>D9</td><td>...</td><td>...</td><td>~60</td><td>~9</td><td>...</td><td>...</td><td>V</td><td>...</td><td>...</td></tr><tr><td>B17</td><td>...</td><td>...</td><td>~120</td><td>~17</td><td>...</td><td>...</td><td>V</td><td>...</td><td>...</td></tr><tr><td rowspan="4">Elias 20</td><td>D25</td><td>-54.4 ± 1.5</td><td>-491.0 ± 1.3</td><td>181.6 ± 1.5</td><td>25.07 ± 0.17</td><td>49 ± 1</td><td>153.2 ± 1.3</td><td>E</td><td>3.5 ± 1.0</td><td>0.75 ± 0.03</td></tr><tr><td>B29</td><td>...</td><td>...</td><td>210 ± 7</td><td>29 ± 1</td><td>...</td><td>...</td><td>R</td><td>5.2 ± 1.0</td><td>...</td></tr><tr><td>D33</td><td>...</td><td>...</td><td>239 ± 7</td><td>33 ± 1</td><td>...</td><td>...</td><td>R</td><td>2.5 ± 1.4<sup>a</sup></td><td>0.95 ± 0.03</td></tr><tr><td>B36</td><td>...</td><td>...</td><td>261 ± 7</td><td>36 ± 1</td><td>...</td><td>...</td><td>R</td><td>1.9 ± 1.4<sup>a</sup></td><td>...</td></tr><tr><td rowspan="4">Elias 24</td><td>D57</td><td>110.0 ± 1.8</td><td>-385.8 ± 1.9</td><td>418 ± 2</td><td>56.8 ± 0.3</td><td>27.7 ± 0.9</td><td>51 ± 2</td><td>E</td><td>22.8 ± 0.3</td><td>0.03 ± 0.01</td></tr><tr><td>B77</td><td>111 ± 1</td><td>-387 ± 1</td><td>564.0 ± 1.2</td><td>76.71 ± 0.13</td><td>29.2 ± 0.4</td><td>45.2 ± 0.8</td><td>E</td><td>12.2 ± 0.3</td><td>...</td></tr><tr><td>D89</td><td>...</td><td>...</td><td>654 ± 7</td><td>89 ± 1</td><td>...</td><td>...</td><td>R</td><td>...</td><td>...</td></tr><tr><td>B123<sup>b</sup></td><td>...</td><td>...</td><td>904 ± 7</td><td>123 ± 1</td><td>...</td><td>...</td><td>R</td><td>...</td><td>...</td></tr><tr><td rowspan="2">Elias 27</td><td>D69</td><td>-5 ± 5</td><td>-8 ± 3</td><td>596<sup>+8</sup><sub>-7</sub></td><td>69.1 ± 0.4</td><td>56.2 ± 0.8</td><td>118.8 ± 0.7</td><td>E</td><td>14.3 ± 1.1</td><td>0.73 ± 0.02</td></tr><tr><td>B86</td><td>...</td><td>...</td><td>741 ± 9</td><td>86 ± 1</td><td>...</td><td>...</td><td>R</td><td>21.2 ± 1.1</td><td>...</td></tr><tr><td rowspan="4">GW Lup</td><td>D74</td><td>-2.6 ± 1.5</td><td>0.8 ± 1.5</td><td>479.2 ± 1.8</td><td>74.3 ± 0.2</td><td>38.7<sup>+0.4</sup><sub>-0.5</sub></td><td>37.7 ± 0.8</td><td>E</td><td>12.1 ± 0.4</td><td>0.31 ± 0.03</td></tr><tr><td>B85</td><td>-2.3 ± 0.2</td><td>1.1<sup>+1.8</sup><sub>-1.7</sub></td><td>551 ± 2</td><td>85.4 ± 0.3</td><td>38.7 ± 0.4</td><td>37.5±0.7</td><td>E</td><td>11.3 ± 0.4</td><td>...</td></tr><tr><td>D103</td><td>...</td><td>...</td><td>665 ± 6</td><td>103 ± 1</td><td>...</td><td>...</td><td>R</td><td>4.3 ± 1.4<sup>a</sup></td><td>0.85 ± 0.09</td></tr><tr><td>B108</td><td>...</td><td>...</td><td>697 ± 6</td><td>108 ± 1</td><td>...</td><td>...</td><td>R</td><td>5.5 ± 1.4<sup>a</sup></td><td>...</td></tr><tr><td rowspan="5">HD 142666</td><td>B6</td><td>...</td><td>...</td><td>41 ± 7</td><td>6 ± 1</td><td>...</td><td>...</td><td>R</td><td>5.3 ± 1.4</td><td>...</td></tr><tr><td>D16</td><td>...</td><td>...</td><td>108 ± 7</td><td>16 ± 1</td><td>...</td><td>...</td><td>R</td><td>3.5 ± 1.4</td><td>0.73 ± 0.02</td></tr><tr><td>B20</td><td>...</td><td>...</td><td>135 ± 7</td><td>20 ± 1</td><td>...</td><td>...</td><td>R</td><td>7.8 ± 1.4</td><td>...</td></tr><tr><td>D37</td><td>...</td><td>...</td><td>243 ± 7</td><td>37 ± 1</td><td>...</td><td>...</td><td>R</td><td>&lt;20</td><td>0.99 ± 0.02</td></tr><tr><td>B40</td><td>...</td><td>...</td><td>270 ± 7</td><td>40 ± 1</td><td>...</td><td>...</td><td>R</td><td>&lt;18</td><td>...</td></tr></tbody></table>
<table><caption>Table 1
(Continued)</caption><thead><tr><th>Source<br>(1)</th><th>Feature<br>(2)</th><th>Δx<br>(mas)<br>(3)</th><th>Δy<br>(mas)<br>(4)</th><th>r<sub>0</sub><br>(mas)<br>(5)</th><th>r<sub>0</sub><br>(au)<br>(6)</th><th>Incl.<br>(°)<br>(7)</th><th>P.A.<br>(°)<br>(8)</th><th>Method<br>(9)</th><th>Width<br>(au)<br>(10)</th><th>Depth<br>(11)</th></tr></thead><tbody><tr><td></td><td>D55</td><td>...</td><td>...</td><td>372 ± 7</td><td>55 ± 1</td><td>...</td><td>...</td><td>R</td><td>2.2 ± 1.4<sup>a</sup></td><td>0.96 ± 0.06</td></tr><tr><td></td><td>B58</td><td>...</td><td>...</td><td>392 ± 7</td><td>58 ± 1</td><td>...</td><td>...</td><td>R</td><td>2.1 ± 1.4<sup>a</sup></td><td>...</td></tr><tr><td rowspan="5">HD 143006</td><td>B6<sup>c</sup></td><td>...</td><td>...</td><td>36 ± 6</td><td>6 ± 1</td><td>...</td><td>...</td><td>R</td><td>5 ± 1.4<sup>a</sup></td><td>—</td></tr><tr><td>D22</td><td>...</td><td>...</td><td>133 ± 6</td><td>22 ± 1</td><td>...</td><td>...</td><td>R</td><td>21.7 ± 1.0</td><td>0.04 ± 0.02</td></tr><tr><td>B41</td><td>-6.3 ± 0.6</td><td>21.5 ± 0.6</td><td>247.3<sup>+0.8</sup><sub>-0.7</sub></td><td>40.8 ± 0.1</td><td>18.9<sup>+0.8</sup><sub>-0.9</sub></td><td>169 ± 3</td><td>E</td><td>12.2 ± 1.0</td><td>...</td></tr><tr><td>D51</td><td>...</td><td>...</td><td>309 ± 6</td><td>51 ± 1</td><td>...</td><td>...</td><td>R</td><td>12.8 ± 1.4</td><td>0.53 ± 0.02</td></tr><tr><td>B65</td><td>-1.4 ± 1.9</td><td>24 ± 2</td><td>393 ± 3</td><td>64.9 ± 0.4</td><td>17.2<sup>+1.9</sup><sub>-2.0</sub></td><td>164<sup>+7</sup><sub>-6</sub></td><td>E</td><td>11.5 ± 1.4</td><td>...</td></tr><tr><td rowspan="7">HD 163296</td><td>D10</td><td>...</td><td>...</td><td>100 ± 10</td><td>10 ± 1</td><td>...</td><td>...</td><td>R</td><td>3.2 ± 1.4<sup>a</sup></td><td>0.93 ± 0.03</td></tr><tr><td>B14</td><td>...</td><td>...</td><td>140 ± 10</td><td>14 ± 1</td><td>...</td><td>...</td><td>R</td><td>3.6 ± 1.4<sup>a</sup></td><td>...</td></tr><tr><td>D48</td><td>...</td><td>...</td><td>480 ± 10</td><td>48 ± 1</td><td>...</td><td>...</td><td>R</td><td>20.2 ± 1.0</td><td>0.033 ± 0.006</td></tr><tr><td>B67</td><td>-4.3 ± 1.2</td><td>6.6 ± 1.2</td><td>663.8 ± 1.6</td><td>67.04 ± 0.12</td><td>46.9 ± 0.2</td><td>133.4 ± 0.3</td><td>E</td><td>15.8 ± 1.0</td><td>...</td></tr><tr><td>D86</td><td>1 ± 3</td><td>2 ± 3</td><td>855 ± 4</td><td>86.4 ± 0.3</td><td>46.9 ± 0.4</td><td>132.7 ± 0.5</td><td>E</td><td>16.2 ± 0.3</td><td>0.151 ± 0.008</td></tr><tr><td>B100</td><td>-2.4 ± 0.9</td><td>9.0 ± 0.9</td><td>987.0 ± 1.2</td><td>99.69 ± 0.09</td><td>46.63 ± 0.12</td><td>133.37 ± 0.17</td><td>E</td><td>12.0 ± 0.3</td><td>...</td></tr><tr><td>D145<sup>b</sup></td><td>...</td><td>...</td><td>1450 ± 10</td><td>145 ± 1</td><td>...</td><td>...</td><td>R</td><td>13.4 ± 1.4</td><td>0.76 ± 0.06</td></tr><tr><td rowspan="2">IM Lup</td><td>B155</td><td>...</td><td>...</td><td>1550 ± 10</td><td>155 ± 1</td><td>...</td><td>...</td><td>R</td><td>14.6 ± 1.4</td><td>...</td></tr><tr><td>D117</td><td>-3 ± 3</td><td>-2 ± 3</td><td>743 ± 4</td><td>117.4 ± 0.5</td><td>47.5 ± 0.5</td><td>145.2 ± 0.7</td><td>E</td><td>15.8 ± 0.7</td><td>0.66 ± 0.02</td></tr><tr><td rowspan="4">MY Lup</td><td>B134</td><td>0 ± 3</td><td>3 ± 3</td><td>845 ± 4</td><td>133.5 ± 0.5</td><td>47.5 ± 0.5</td><td>143.9 ± 0.6</td><td>E</td><td>18.4 ± 0.7</td><td>...</td></tr><tr><td>D8</td><td>...</td><td>...</td><td>~50</td><td>~8</td><td>...</td><td>...</td><td>V</td><td>...</td><td>...</td></tr><tr><td>B20</td><td>...</td><td>...</td><td>~130</td><td>~20</td><td>...</td><td>...</td><td>V</td><td>...</td><td>...</td></tr><tr><td>D30</td><td>...</td><td>...</td><td>~190</td><td>~30</td><td>...</td><td>...</td><td>V</td><td>...</td><td>...</td></tr><tr><td rowspan="7">RU Lup</td><td>B40</td><td>...</td><td>...</td><td>~260</td><td>~40</td><td>...</td><td>...</td><td>V</td><td>...</td><td>...</td></tr><tr><td>D14</td><td>...</td><td>...</td><td>~90</td><td>~14</td><td>...</td><td>...</td><td>V</td><td>...</td><td>...</td></tr><tr><td>B17</td><td>...</td><td>...</td><td>~110</td><td>~17</td><td>...</td><td>...</td><td>V</td><td>...</td><td>...</td></tr><tr><td>D21</td><td>...</td><td>...</td><td>132 ± 6</td><td>21 ± 1</td><td>...</td><td>...</td><td>R</td><td>&lt;7</td><td>...</td></tr><tr><td>B24</td><td>...</td><td>...</td><td>151 ± 6</td><td>24 ± 1</td><td>...</td><td>...</td><td>R</td><td>&lt;8</td><td>...</td></tr><tr><td>D29</td><td>-17.0 ± 1.2</td><td>87.8<sup>+1.2</sup><sub>-1.1</sub></td><td>183.0 ± 1.5</td><td>29.1 ± 0.2</td><td>20 ± 2</td><td>117 ± 6</td><td>E</td><td>4.5 ± 0.3</td><td>0.784 ± 0.013</td></tr><tr><td>B34</td><td>-17.2<sup>+1.2</sup><sub>-1.1</sub></td><td>88.4 ± 1.2</td><td>213.7 ± 1.4</td><td>34.0 ± 0.2</td><td>17 ± 2</td><td>126 ± 7</td><td>E</td><td>5.5 ± 0.3</td><td>...</td></tr><tr><td rowspan="2">SR 4</td><td>D42</td><td>...</td><td>...</td><td>~260</td><td>~42</td><td>...</td><td>...</td><td>V</td><td>&lt;16</td><td>...</td></tr><tr><td>B50</td><td>...</td><td>...</td><td>~310</td><td>~50</td><td>...</td><td>...</td><td>V</td><td>...</td><td>...</td></tr><tr><td rowspan="2">Sz 114</td><td>D11</td><td>...</td><td>...</td><td>82 ± 7</td><td>11 ± 1</td><td>...</td><td>...</td><td>R</td><td>6.3 ± 1.4</td><td>0.23 ± 0.02</td></tr><tr><td>B18</td><td>...</td><td>...</td><td>134 ± 3</td><td>18 ± 1</td><td>...</td><td>...</td><td>R</td><td>13.3 ± 1.4</td><td>...</td></tr><tr><td rowspan="2">Sz 129</td><td>D39</td><td>0 ± 5</td><td>3 ± 4</td><td>238 ± 3</td><td>38.6 ± 0.6</td><td>21 ± 2</td><td>165 ± 7</td><td>E</td><td>4.3 ± 0.8<sup>a</sup></td><td>0.94 ± 0.05</td></tr><tr><td>B45</td><td>-2 ± 3</td><td>4 ± 3</td><td>280 ± 3</td><td>45.4 ± 0.5</td><td>21.6<sup>+1.6</sup><sub>-1.7</sub></td><td>165 ± 5</td><td>E</td><td>5.4 ± 0.8</td><td>...</td></tr><tr><td rowspan="5">WaOph 6</td><td>B10</td><td>...</td><td>...</td><td>62 ± 6</td><td>10 ± 1</td><td>...</td><td>...</td><td>R</td><td>17.6 ± 1.1</td><td>...</td></tr><tr><td>D41</td><td>4 ± 3</td><td>4 ± 2</td><td>255 ± 3</td><td>41.0 ± 0.4</td><td>32.4<sup>+1.9</sup><sub>-2</sub></td><td>148 ± 4</td><td>E</td><td>4.1 ± 0.6<sup>a</sup></td><td>0.98 ± 0.02</td></tr><tr><td>B46</td><td>7 ± 3</td><td>2 ± 2</td><td>287 ± 3</td><td>46.1 ± 0.4</td><td>35.3<sup>+1.7</sup><sub>-1.8</sub></td><td>153 ± 3</td><td>E</td><td>3.5 ± 0.6<sup>a</sup></td><td>...</td></tr><tr><td>D64</td><td>...</td><td>...</td><td>~400</td><td>~64</td><td>...</td><td>...</td><td>V</td><td>&lt;23</td><td>...</td></tr><tr><td>B69</td><td>...</td><td>...</td><td>~430</td><td>~69</td><td>...</td><td>...</td><td>V</td><td>...</td><td>...</td></tr><tr><td>WAOph 6</td><td>D79</td><td>...</td><td>...</td><td>642 ± 8</td><td>79 ± 1</td><td>...</td><td>...</td><td>R</td><td>6.7 ± 1.1</td><td>0.98 ± 0.03</td></tr></tbody></table>
<table><caption>**Table 1**
(Continued)</caption><thead><tr><th>Source<br>(1)</th><th>Feature<br>(2)</th><th>Δx<br>(mas)<br>(3)</th><th>Δy<br>(mas)<br>(4)</th><th>r<sub>0</sub><br>(mas)<br>(5)</th><th>r<sub>0</sub><br>(au)<br>(6)</th><th>Incl.<br>(°)<br>(7)</th><th>P.A.<br>(°)<br>(8)</th><th>Method<br>(9)</th><th>Width<br>(au)<br>(10)</th><th>Depth<br>(11)</th></tr></thead><tbody><tr><td></td><td>B88</td><td>-244 ± 3</td><td>-361 ± 3</td><td>712 ± 4</td><td>87.5 ± 0.4</td><td>47.3 ± 0.7</td><td>174.2 ± 0.8</td><td>E</td><td>3.6 ± 1.1<sup>a</sup></td><td></td></tr><tr><td rowspan="2">WSB 52</td><td>D21<sup>d</sup></td><td>...</td><td>...</td><td>~150</td><td>~21</td><td>...</td><td>...</td><td>V</td><td>...</td><td>...</td></tr><tr><td>B25<sup>d</sup></td><td>...</td><td>...</td><td>~180</td><td>~25</td><td>...</td><td>...</td><td>V</td><td>...</td><td>...</td></tr><tr><td rowspan="10">TW Hya<sup>e</sup></td><td>D1</td><td>...</td><td>...</td><td>~17</td><td>~1</td><td>...</td><td>...</td><td>V</td><td>...</td><td>...</td></tr><tr><td>B3</td><td>...</td><td>...</td><td>50 ± 8</td><td>3.0 ± 0.5</td><td>...</td><td>...</td><td>R</td><td>2.2 ± 0.6</td><td>...</td></tr><tr><td>D26</td><td>9 ± 2</td><td>13 ± 2</td><td>427.0 ± 1.6</td><td>25.62 ± 0.14</td><td>...</td><td>...</td><td>E</td><td>4.0 ± 0.5</td><td>0.804 ± 0.008</td></tr><tr><td>B30</td><td>...</td><td>...</td><td>492 ± 8</td><td>29.5 ± 0.5</td><td>...</td><td>...</td><td>R</td><td>&lt;6</td><td>...</td></tr><tr><td>D32</td><td>...</td><td>...</td><td>525 ± 8</td><td>31.5 ± 0.5</td><td>...</td><td>...</td><td>R</td><td>&lt;4</td><td>0.986 ± 0.008</td></tr><tr><td>B33<sup>b</sup></td><td>...</td><td>...</td><td>550 ± 8</td><td>33 ± 0.5</td><td>...</td><td>...</td><td>R</td><td>&lt;10</td><td>...</td></tr><tr><td>D42</td><td>8.2 ± 1.5</td><td>13.7 ± 1.5</td><td>694.0 ± 1.1</td><td>41.64 ± 0.09</td><td>...</td><td>...</td><td>E</td><td>3.40 ± 0.14</td><td>0.763 ± 0.009</td></tr><tr><td>B45</td><td>6.8 ± 1.2</td><td>15.3 ± 1.2</td><td>745.3 ± 0.9</td><td>44.72 ± 0.07</td><td>...</td><td>...</td><td>E</td><td>2.80 ± 0.14</td><td>...</td></tr><tr><td>D48</td><td>...</td><td>...</td><td>~800</td><td>~48</td><td>...</td><td>...</td><td>V</td><td>&lt;7</td><td>...</td></tr><tr><td>B52</td><td>...</td><td>...</td><td>~870</td><td>~52</td><td>...</td><td>...</td><td>V</td><td>...</td><td>...</td></tr><tr><td rowspan="15">HL Tau<sup>f</sup></td><td>D14</td><td>...</td><td>...</td><td>94.6 ± 1.4</td><td>13.9 ± 0.2</td><td>...</td><td>...</td><td>...</td><td>5.7 ± 0.2</td><td>0.41 ± 0.03</td></tr><tr><td>B21</td><td>...</td><td>...</td><td>145.6 ± 0.7</td><td>21.4 ± 0.1</td><td>...</td><td>...</td><td>...</td><td>12.9 ± 0.2</td><td>...</td></tr><tr><td>D34</td><td>...</td><td>...</td><td>230.6 ± 0.7</td><td>33.9 ± 0.1</td><td>...</td><td>...</td><td>...</td><td>4.9 ± 0.1</td><td>0.67 ± 0.02</td></tr><tr><td>B40</td><td>...</td><td>...</td><td>272.1 ± 0.7</td><td>40.0 ± 0.1</td><td>...</td><td>...</td><td>...</td><td>6.4 ± 0.1</td><td>...</td></tr><tr><td>D44</td><td>...</td><td>...</td><td>299 ± 7</td><td>44 ± 1</td><td>...</td><td>...</td><td>...</td><td>3.2 ± 1.4</td><td>0.92 ± 0.02</td></tr><tr><td>B49</td><td>...</td><td>...</td><td>333 ± 7</td><td>49 ± 1</td><td>...</td><td>...</td><td>...</td><td>4.2 ± 1.4</td><td>...</td></tr><tr><td>D53</td><td>...</td><td>...</td><td>~360</td><td>~53</td><td>...</td><td>...</td><td>...</td><td>&lt;10</td><td>...</td></tr><tr><td>B58</td><td>...</td><td>...</td><td>~390</td><td>~58</td><td>...</td><td>...</td><td>...</td><td>&lt;15</td><td>...</td></tr><tr><td>D67</td><td>...</td><td>...</td><td>458.5 ± 0.7</td><td>67.4 ± 0.1</td><td>...</td><td>...</td><td>...</td><td>4.4 ± 0.1</td><td>0.84 ± 0.03</td></tr><tr><td>B72</td><td>...</td><td>...</td><td>491.2 ± 0.7</td><td>72.2 ± 0.1</td><td>...</td><td>...</td><td>...</td><td>5.1 ± 0.1</td><td>...</td></tr><tr><td>D77</td><td>...</td><td>...</td><td>526.5 ± 0.7</td><td>77.4 ± 0.1</td><td>...</td><td>...</td><td>...</td><td>6.1 ± 0.1</td><td>0.62 ± 0.02</td></tr><tr><td>B85</td><td>...</td><td>...</td><td>581.0 ± 0.7</td><td>85.4 ± 0.1</td><td>...</td><td>...</td><td>...</td><td>12.0 ± 0.1</td><td>...</td></tr><tr><td>D96</td><td>...</td><td>...</td><td>~650</td><td>~96</td><td>...</td><td>...</td><td>...</td><td>&lt;20</td><td>...</td></tr><tr><td>B102</td><td>...</td><td>...</td><td>~690</td><td>~102</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td></tr></tbody></table>

**Notes.** Column descriptions: (1) Name of host star. (2) Substructure label. (3) R.A. offset from phase center. (4) Decl. offset from phase center. (5) Radial location of substructure in mas. (6) Radial location of substructure in astronomical units (the uncertainties in mas are simply scaled from the fitting procedure and do not account for the uncertainty in the distance to the source). (7) Inclination of substructure (0° is face-on and 90° is edge-on. (8) Position angle of substructure (east of north). (9) Method used to derive radial location of substructure: “E” indicates ellipse-fitting, “R” indicates identification of local extrema in the radial profile, and “V” indicates identification through visual inspection. (10) Width of substructure (only measured for those identified via “E” and “R” methods). (11) Depth of gap. All uncertainties are $1\sigma$.

a Width of feature is narrower than the minor axis FWHM of the synthesized beam (see Table 2 for the synthesized beam size).

b This feature shows tentative signs of being multiple gaps and rings.

c This ring may be misaligned relative to the outer two rings and therefore could have a larger radius than implied by the radial profile. See Pérez et al. (2018).

d Substructures tentatively classified as annular.

e Fit to the 290 GHz continuum image from Huang et al. (2018a). The P.A. and inclination are fixed at 155° and 7°, respectively, based on values derived from $^{12}$CO observations in Qi et al. (2004).

f Measurements for substructure positions taken from ALMA Partnership et al. (2015) and rescaled to a distance of 147 pc (Galli et al. 2018).

(This table is available in machine-readable form.)
<table><caption>**Table 2**
Disk Geometries</caption><thead><tr><th>Source<br>(1)</th><th>Δx<br>(mas)<br>(2)</th><th>Δy<br>(mas)<br>(3)</th><th>Incl.<br>(deg)<br>(4)</th><th>P.A.<br>(deg)<br>(5)</th><th>Method<br>(6)</th><th>R<sub>dust</sub><br>(au)<br>(7)</th><th>θ<sub>b,min</sub> × d<br>(au)<br>(8)</th></tr></thead><tbody><tr><td>AS 209</td><td>1.9 ± 0.5</td><td>-2.5 ± 0.5</td><td>34.97 ± 0.13</td><td>85.76 ± 0.16</td><td>E</td><td>139 ± 1</td><td>4.4</td></tr><tr><td>DoAr 25</td><td>38 ± 2</td><td>-494 ± 2</td><td>67.4 ± 0.2</td><td>110.6 ± 0.2</td><td>E</td><td>165 ± 1</td><td>3.0</td></tr><tr><td>DoAr 33</td><td>1.5 ± 0.8</td><td>0.6 ± 0.6</td><td>41.8 ± 0.8</td><td>81.1 ± 1.2</td><td>G</td><td>27 ± 1</td><td>3.4</td></tr><tr><td>Elias 20</td><td>-54.5 ± 1.5</td><td>-491.0 ± 1.3</td><td>49 ± 1</td><td>153.2 ± 1.3</td><td>E</td><td>64 ± 1</td><td>3.2</td></tr><tr><td>Elias 24</td><td>110.8 ± 0.8</td><td>-386.8 ± 0.9</td><td>29.0 ± 0.3</td><td>45.7 ± 0.7</td><td>E</td><td>136 ± 1</td><td>4.7</td></tr><tr><td>Elias 27</td><td>-5 ± 5</td><td>-8 ± 3</td><td>56.2 ± 0.8</td><td>118.8 ± 0.7</td><td>E</td><td>254 ± 1</td><td>5.5</td></tr><tr><td>GW Lup</td><td>-2.4 ± 1.1</td><td>0.9 ± 1.2</td><td>38.7 ± 0.3</td><td>37.6 ± 0.5</td><td>E</td><td>105 ± 1</td><td>6.6</td></tr><tr><td>HD 142666</td><td>-43.5 ± 0.3</td><td>28.9 ± 0.6</td><td>62.22 ± 0.14</td><td>162.11 ± 0.15</td><td>G</td><td>59 ± 1</td><td>3.3</td></tr><tr><td>HD 143006</td><td>-5.9 ± 0.6</td><td>21.7 ± 0.6</td><td>18.6 ± 0.8</td><td>169 ± 2</td><td>E</td><td>82 ± 1</td><td>7.4</td></tr><tr><td>HD 163296</td><td>-2.8 ± 0.7</td><td>7.7 ± 0.7</td><td>46.7 ± 0.1</td><td>133.33 ± 0.15</td><td>E</td><td>169 ± 1</td><td>3.9</td></tr><tr><td>IM Lup</td><td>-1.5 ± 2</td><td>1 ± 2</td><td>47.5 ± 0.3</td><td>144.5 ± 0.5</td><td>E</td><td>264 ± 1</td><td>6.9</td></tr><tr><td>MY Lup</td><td>-77.9 ± 0.8</td><td>62.9 ± 0.6</td><td>73.2 ± 0.1</td><td>58.8 ± 0.1</td><td>G</td><td>87 ± 1</td><td>6.7</td></tr><tr><td>RU Lup</td><td>-17.1 ± 0.8</td><td>88.1 ± 0.8</td><td>18.8 ± 1.6</td><td>121 ± 5</td><td>E</td><td>63 ± 1</td><td>3.9</td></tr><tr><td>SR 4</td><td>-56.4 ± 1.2</td><td>-507.4 ± 1.3</td><td>22 ± 2</td><td>18 ± 5</td><td>G</td><td>31 ± 1</td><td>4.6</td></tr><tr><td>Sz 114</td><td>-1 ± 2</td><td>4 ± 2</td><td>21.3 ± 1.3</td><td>165 ± 4</td><td>E</td><td>58 ± 1</td><td>4.6</td></tr><tr><td>Sz 129</td><td>5.4 ± 1.9</td><td>3.0 ± 1.7</td><td>34.1 ± 1.3</td><td>151 ± 2</td><td>E</td><td>76 ± 1</td><td>5.0</td></tr><tr><td>WaOph 6</td><td>-244 ± 3</td><td>-361 ± 3</td><td>47.3 ± 0.7</td><td>174.2 ± 0.8</td><td>E</td><td>103 ± 1</td><td>6.6</td></tr><tr><td>WSB 52</td><td>-119.5 ± 0.4</td><td>-432.8 ± 0.4</td><td>54.4 ± 0.3</td><td>138.4 ± 0.3</td><td>G</td><td>32 ± 1</td><td>3.7</td></tr></tbody></table>

**Note.** Column descriptions: (1) R.A. offset from phase center of disk. (2) Decl. offset from phase center of disk. (3) Disk inclination. (4) Disk position angle. (5) Method used to derive disk P.A., inclination, and phase offset. “E” denotes that the offset from phase center, P.A., and inclination are derived by fitting ellipses to individual annular substructures. “G” denotes that these quantities are estimated by fitting a 2D Gaussian to the image. (6) Radial extent of 1.25 mm continuum emission. (7) FWHM of the synthesized beam along the minor axis multiplied by the distance to the source (see Table 4 of Andrews et al. 2018a for beam dimensions, distance to source, and references).

order in which they appear in the disk, we use the radial location of the substructure because precise distances are now available from Gaia (Gaia Collaboration et al. 2016, 2018), and we favor a system that can flexibly accommodate additional substructures identified in future observations.

The second part of our procedure addresses features that do not lend themselves well to ellipse-fitting due to lower S/N or insufficient angular resolution along the projected disk minor axis. The positions of these substructures are instead identified through examination of the deprojected and azimuthally averaged radial intensity profiles. In order to deproject the disks, P.A.s, inclinations and center positions first need to be calculated. For most disks, these quantities are computed by taking weighted averages of the ellipse fits for the annular substructures. For WSB 52, MY Lup, DoAr 33, HD 142666, and SR 4, P.A.s, inclinations, and center positions are instead computed by fitting two-dimensional Gaussians to the images with the imfit task in CASA. Even when an intensity profile deviates from a Gaussian, a Gaussian fit should still provide a good estimate of the source inclination and P.A. provided that the emission is axisymmetric and much larger than the PSF (e.g., Bernstein & Jarvis 2002). To double-check that Gaussian fitting yields reasonable results, we also fit Gaussians to disks with well-defined rings (and therefore highly non-Gaussian profiles) and find that the results generally agree with those from ellipse-fitting within one degree. The measurements for all disks are listed in Table 2.

Radial intensity profiles are then calculated for each disk by using the previously derived geometries to deproject the disk, binning the pixels in annuli 1 au wide, and averaging the intensities in each bin. For most disks, the average is taken through all azimuthal angles. For HD 143006 and HD 163296, angles where the outer disk emission asymmetry is visible are

excluded from the average (80° < θ < 144° and 50° < θ < 130°, respectively). For disks with higher inclinations (Elias 20, DoAr 25, MY Lup, WSB 52, HD 142666, and HL Tau), only azimuthal angles that are within ±20° (as measured in deprojected coordinates) of the projected disk major axis are included in the average because most of the substructures are poorly resolved along the minor axis. The radial profile is calculated in a similar manner for Sz 114 as well due to the extreme elongation of the synthesized beam. Figure 2 shows the labeled radial intensity profiles. Figure 3 provides logarithmically scaled radial intensity profiles of the disks exceeding 100 au in radius in order to show the substructures in the faint outer regions more clearly.

The radial locations of additional annular substructures are identified by searching for local maxima and minima in the averaged radial profiles and cross-referencing with the original image to verify that they correspond to annular features rather than nonaxisymmetric features such as spiral arms. The uncertainty on the radial location of each extremum is estimated as the width of the radial bin (i.e., 1 au). The positions of features measured through ellipse-fitting are checked for consistency with the positions measured from the averaged radial intensity profiles.

Finally, for certain disks (e.g., MY Lup and DoAr 33), there are regions of the radial profile that are “plateau”-like rather than local maxima or minima. We classify them as annular substructures because they still appear to represent deviations from standard “smooth” disk profiles (e.g., Andrews et al. 2009; Tazzari et al. 2017). This is in part justified by optical depth estimates (Section 3.5) indicating that many of these “plateaus” in the intensity profile may trace local extrema in the surface density profile. The inner and outer edges of the “plateau” regions are identified through visual inspection, with some guidance from checking where $\frac{1}{l_r(r)} \frac{dl_r(r)}{dr}$, the radial intensity profile slope divided by the radial intensity, exceeds
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Deprojected and azimuthally averaged radial intensity profiles ordered by decreasing stellar luminosity from left to right and top to bottom. Light blue ribbons show the $1\sigma$ scatter at each radial bin divided by the number of beams spanning the angles over which the intensities are measured. Solid gray lines mark emission rings and dotted black lines mark gaps listed in Table 1. The Gaussian profiles show the FWHM of the minor axis of the synthesized beams.</figcaption></figure>

−0.05 (i.e., where the radial decrease of the intensity is relatively small). The inner edge of the “plateau” is labeled with the prefix “D” followed by its radial location in astronomical units rounded to the nearest whole number.

Similarly, the outer edge of the “plateau” is labeled with the prefix “B” followed by its radial location. Because of the greater degree of subjectivity involved in identifying the locations of these features compared to identifying local
<figure><img src="image_7.png" /><figcaption>**Figure 3.** Radial intensity profiles plotted on a logarithmic scale for the largest DSHARP disks, ordered by descending stellar luminosity from left to right and top to bottom. Light blue ribbons show the $1\sigma$ scatter at each radial bin divided by the number of beams spanning the angles over which the intensities are measured. The y-axis starts at $10^{-2}$ mJy beam$^{-1}$, corresponding to slightly less than the typical rms level of the continuum images. Note that the Elias 27, IM Lup, and WaOph 6 disks have spiral arms that also create substructures in their radial profiles.</figcaption></figure>

extrema in the emission profiles, approximate locations are listed in Table 1 without formal error estimates. However, the precision should be on the order of the scale of the synthesized beam (i.e., a few astronomical units).

For comparison, the annular substructure locations for the TW Hya and HL Tau disks are appended to Table 1. The small difference in the continuum frequency compared to the DSHARP sources is not expected to significantly affect the substructure positions. The HL Tau substructure positions from Table 2 of ALMA Partnership et al. (2015) are adjusted for the new distance estimate of 147 pc from Galli et al. (2018) and renamed accordingly. The TW Hya substructure positions are measured using the same procedure described for the other sources, except the P.A. is fixed to 155° and the inclination to 7° because nearly face-on disk orientations are more reliably measured from gas observations than from continuum observations (Qi et al. 2004). We do not rely on molecular line observations to measure the inclinations of the DSHARP disks because a large fraction of them are cloud-contaminated (see the $^{12}$CO channel maps in Andrews et al. 2018a). The distance to TW Hya is 60 pc (Gaia Collaboration et al. 2016),

and radial bins are in increments of 0.5 au because of the smaller distance.

Of the DSHARP sample, millimeter continuum substructures have previously been reported for the disks around AS 209 (Fedele et al. 2018; Huang et al. 2017), DoAr 25 (Cox et al. 2017), Elias 24 (Cieza et al. 2017; Cox et al. 2017; Dipierro et al. 2018), Elias 27 (Pérez et al. 2016), HD 142666 (Rubinstein et al. 2018), HD 143006 (Barenfeld et al. 2016), and HD 163296 (Isella et al. 2016; Zhang et al. 2016). They have also been tentatively reported for the disks around IM Lup (Cleeves et al. 2016) and Sz 129 (Tazzari et al. 2017). Although a few studies have classified MY Lup as a transition disk (e.g., Romero et al. 2012; Tazzari et al. 2017; van der Marel et al. 2018), the emission observed at high angular resolution is still centrally peaked. However, there is radial substructure in the vicinity of where Tazzari et al. (2017) inferred a cavity wall from visibility modeling, which may account for its previous classification as a transition disk.

The improved resolution of the DSHARP observations reveals additional structural complexity in all disks with previously known or proposed substructures. Millimeter continuum substructures are
reported for the first time in the other DSHARP disks. However, the substructure in the WSB 52 disk is only tentatively classified as annular because the inclination of the disk makes it difficult to discern the exact nature of the features observed along the major axis of the projected emission.

## 3.2. Widths and Contrasts of Annular Substructures

Defining the sizes of disk substructures is a more nebulous task compared to measuring their radial positions. The widths and amplitudes of emission rings are sometimes measured by modeling them as Gaussian rings, as is done for the DSHARP observations of AS 209, GW Lup, Elias 24, HD 143006, and HD 163296 (Dullemond et al. 2018; Guzmán et al. 2018; Isella et al. 2018; Pérez et al. 2018). These sources lend themselves well to such a procedure because of the high contrast and large separation of many of their substructures. In other cases where the Gaussian components overlap substantially or where the substructures do not closely resemble Gaussians, it is less useful to define ring widths in terms of Gaussian parameters. In disk simulations, the widths and depths of gaps are often defined with respect to an initial surface density profile (e.g., Duffell 2015; Kanagawa et al. 2016), but these definitions do not translate well to observations, where the initial surface density profile is unknown.

Instead, similarly to Zhang et al. (2018), we take a more empirical approach and define width and depth without assuming a functional form for the substructures or an initial surface density. For an adjacent gap-ring pair, we define the gap radius $r_d$ as the value of $r_0$ listed in Table 1 for a substructure labeled with the prefix “D” and the reference ring radius $r_b$ as the $r_0$ value listed for the “B”-labeled substructure directly outside the gap. The gap depth is then defined as

$$ \frac{I_d}{I_b}, \qquad (2) $$

where the gap intensity $I_d$ is the mean intensity value in the radial bin containing $r_d$ and the reference intensity $I_b$ is the mean intensity value in the radial bin containing $r_b$. The uncertainty adopted for the mean intensity value is $\sigma_{bin}/\sqrt{N_{beam}}$, where $\sigma_{bin}$ is the $1\sigma$ scatter of the pixel values in the radial bin and $N_{beam}$ is the number of synthesized beams spanned by the ellipse or arc over which the intensity is being averaged.

For measurements of substructure widths, the radial intensity profiles are then linearly interpolated onto radial gridpoints spaced by 0.1 au. The dividing point between the outer edge of a gap and the inner edge of its immediate exterior ring is defined as the radius $r_{d,o}$ (where $r_d < r_{d,o} < r_b$ at which the intensity is equal to $I_{mean} = 0.5(I_d + I_b)$). The inner edge of the gap $r_{d,i}$ is defined as the largest value $r$ satisfying the criteria that $I(r) = I_{mean}$ and $r < r_d$. The outer edge of a ring $r_{b,o}$ is defined as the smallest value $r$ satisfying the criteria that $I(r) = I_{mean}$ and $r > r_b$. Thus, the width of a gap is $r_{d,o} - r_{d,i}$ and the width of a ring is $r_{b,o} - r_{d,o}$. For rings with emission profiles that can be modeled as isolated Gaussians (e.g., Dullemond et al. 2018), this definition of ring width is essentially the FWHM. A diagram is shown in Appendix B, and a couple special cases are addressed. The measurements are listed in Table 1. These definitions of width and depth are only applicable to the “D”/ “B” pairs that correspond to local extrema in the radial profiles, i.e., those with measurement methods labeled either “E” or “R.” Features can have measured widths smaller than the synthesized

beam because the peak-to-peak or trough-to-trough separations of the gaps/rings on either side of the feature are still larger than the beam.

For shallow emission features that occur within a series of annular substructures, upper bounds for the widths can be estimated by measuring the peak-to-peak or trough-to-trough distances of the adjacent substructures. It is not clear whether the shallowest emission features correspond to surface density features that are shallow and wide or deep and narrow. Observations at higher sensitivity and angular resolution are necessary to measure the widths and contrasts of shallow substructures.

These definitions have a few important limitations. First, the effects of the beam are not accounted for. Thus, the true gap widths should generally be larger and the ring widths should be narrower. In addition, the true ratio of $I_d$ to $I_b$ should generally be smaller because beam smearing will reduce peak intensities and “fill in” gaps. For well-resolved features, the effect of beam convolution on the width and depth estimates should not be too significant, but the effects are quite significant for features near the angular resolution limit. Therefore, the values quoted for such narrow features may be better viewed as upper or lower bounds, as appropriate. Second, the gap contrast is defined with respect to the exterior ring because a local maximum is guaranteed to occur outside a local minimum in the radial intensity profile, whereas there is not necessarily a ring interior to a gap. However, this can sometimes lead to counterintuitive outcomes, such as D90 in the AS 209 disk having a relatively shallow measurement for gap depth even though the gap intensity is at most a few percent of the peak intensity of the immediately interior B74 feature. Similarly, if the local minimum outside a ring has a much lower intensity value than the local minimum interior to a ring, such as B39 in the AS 209 disk or B65 in the HD 143006 disk, it is not clear whether it is more sensible to view the ring as a narrow enhancement on top of a background profile (which our definition implicitly does) or as a single wide structure without additional background emission (as is done when disk emission is modeled as a series of Gaussian rings). Overall, though, the definitions we adopt are still useful for comparing substructures within the sample and for assessing whether the features are resolved.

The substructure widths and contrasts span a wide range. Substructures range from widths of a few astronomical units to more than 20 au, but the majority of the substructures are narrower than 10 au. Variations can be quite wide even within the same sources, as exemplified by AS 209 and HD 163296. Neither substructure widths nor depths in the DSHARP sources increase monotonically with distance from the central star. A majority of the disks have substructures that are marginally resolved (not necessarily exclusively so), suggesting that there could be additional substructures under the resolution limit. The deepest gaps have intensity levels that are at most a few percent of the peak intensity of adjacent rings, but more typically the intensity variations between adjacent gaps and rings are less than 20%. However, many of the low-contrast gaps have widths comparable to the size of the synthesized beam, so their apparent shallowness may result from insufficient resolution.

Figure 4 shows the ratio of substructure widths to their radial positions, $\Delta r/r_0$, for all but the lowest contrast disks (i.e., MY Lup, DoAr 33, and WSB 52). The measured values are generally largest in the inner 20 au of these disks, which is a limitation imposed by angular resolution, as shown by the gray
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Plots of substructure widths divided by radial positions for the 15 disks with sufficiently high-contrast substructures. The positions of substructures that are too shallow to allow a width measurement are marked with blue (gap) or orange (ring) vertical lines. The gray dotted curves correspond to $\frac{\theta_{b,min}}{r} \times d$, separating resolved features above the curve from marginally resolved or unresolved features at or below the curve.</figcaption></figure>

curves in Figure 4. Otherwise, values often fall between 0.05 and 0.2, comparable to the typical values for disk aspect ratios, $h/r$ (e.g., Andrews et al. 2010). Equivalently, the widths are often comparable to expected values for the pressure scale height. A few gaps in the outer disk have notably high values of $\Delta r/r_0 > 0.25$: D48 in the HD 163296 disk, D57 in the Elias 24 disk, D22 in the HD 143006 disk, and D61 in the AS 209 disk. These features are among the widest DSHARP gaps in absolute, not simply relative, terms. They are also among the deepest gaps observed, with intensities inside the gap being at most a few percent of those of adjacent rings.

### 3.3. Dust Disk Radii

To assess how substructure characteristics relate to disk size, we also compute the radius of the millimeter disk, $R_{\text{dust}}$, for each source using a “curve-of-growth”-type method similar to that described in Ansdell et al. (2016). For each disk, continuum fluxes are measured within a series of increasingly large elliptical

apertures with the same P.A. and major axis/minor axis ratio as the projected disk image. The total continuum flux is taken to be the asymptotic value of the flux curve. $R_{\text{dust}}$ is defined to be the radius at which the enclosed flux is equal to 95% of the total flux. The values are listed in Table 2. The uncertainty is estimated as the spacing between elliptical apertures (1 au, limited by the pixel size); the relative uncertainty in the continuum flux does not contribute substantially to the uncertainty in $R_{\text{dust}}$ due to the high S/N of the images.

### 3.4. Additional Substructures

Sz 129, HD 143006, and HD 142666 exhibit decreases in intensity toward the center of the disk. The angular resolution is not sufficient to determine whether these features are cavities or gaps encircling inner disks. However, high resolution millimeter continuum observations that have detected inner disks inside the cavities of a few transition disks suggest that the distinction between a cavity and an annular gap could largely
<figure><img src="image_9.png" /><figcaption>**Figure 5.** Insets of the inner disk brightness asymmetries of HD 142666, Sz 129, and HD 163296. Gray contours are drawn at [70, 80, 90]× rms for HD 142666, [45, 50, 55, 60]× rms for Sz 129, and [120, 140, 160, 180]× rms for HD 163296.</figcaption></figure>

be a matter of angular resolution (e.g., Andrews et al. 2016; Boehler et al. 2017; Dong et al. 2018b).

The innermost emission ring of Sz 129 is brighter in the northeast and that of HD 142666 is brighter on the southwest side. The millimeter continuum intensity peak of the HD 163296 disk is also offset southwest from the estimated center of the disk by ~25 mas and may similarly be tracing an inner emission ring. Insets of these disks are shown in Figure 5. For both HD 163296 and HD 142666, the ¹²CO J = 2–1 emission indicates that the continuum brightness asymmetries in the inner disk are on the far side of the disk (Andrews et al. 2018a; Isella et al. 2018). Thus, the brightness asymmetry may arise from viewing the heated and puffed up interior of a ring. The absolute geometry of Sz 129 is not immediately clear from the ¹²CO observations.

A few disks have tentative annular substructures in addition to what is listed in Table 1. The radial profile of DoAr 25 shows some departures from smoothness inward of D76. However, because the PSF for this disk is poorer than most of the other DSHARP disks due to the abbreviated integration time, better uv coverage will be needed to confirm the structure of the inner disk. The radial profile of Elias 24 has hints of very-low-amplitude variations between D89 and B123. The emission profile in this region is also unusual among the DSHARP disks for overall being quite flat over a radial span of tens of astronomical units, as best seen in Figure 3. The radial profile of Elias 20 shows some low-amplitude variations outside B36, but additional annular substructures are not obvious in the image (note that this is one of the sources for which the radial profile is only averaged through wedges around the major axis). The disks with spiral arms (WaOph 6, IM Lup, and Elias 27) may have additional gaps and rings beyond what is listed in Table 1, but as discussed further in Huang et al. (2018b), these features may instead be branches or tightly wrapped extensions of the spiral arms.

Finally, HD 163296 has a crescent-like asymmetry interior to B67 and HD 143006 has an asymmetry exterior to B65. These are discussed further in Isella et al. (2018) and Pérez et al. (2018), respectively.

### 3.5. Continuum Optical Depths

It has usually been assumed that the millimeter dust continuum of disks is optically thin (e.g., Beckwith et al. 1990). Recent analyses of disk spectral indices and emission sizes have explored whether disks only appear to be optically thin at low resolution because of beam dilution of optically thick substructures (e.g., Ricci et al. 2010; Tripathi et al. 2017), motivating an analysis of whether the DSHARP substructures are optically thick. Multi-wavelength observations are desirable for constraining dust optical depths (e.g., Pérez et al. 2012; Tazzari et al. 2016), but a rough estimate can still be obtained from the single-wavelength DSHARP observations by approximating the midplane temperature profile with the simplified expression for a passively heated, flared disk in radiative equilibrium (e.g., Chiang & Goldreich 1997; D’Alessio et al. 1998; Dullemond et al. 2001):

$$T_{\text{mid}}(r) = \left( \frac{\varphi L_{*}}{8\pi r^2 \sigma_{\text{SB}}} \right)^{0.25}, \quad (3)$$

where $\sigma_{\text{SB}}$ is the Stefan–Boltzmann constant, $L_*$ is the stellar luminosity (taken from Table 1 of Andrews et al. 2018a), and $\varphi$ is the flaring angle. The optical depth $\tau_{\nu}$ is then calculated using the relationship

$$I_{\nu}(r) = B_{\nu}(T_{\text{mid}}(r))(1 - \exp(-\tau_{\nu}(r))), \quad (4)$$

where $I_{\nu}(r)$ is the deprojected, azimuthally averaged radial intensity profile. The full Planck expression is used because the Rayleigh–Jeans approximation loses accuracy at millimeter wavelengths.

Since colder temperature profiles lead to a higher estimate of $\tau$ and we wish to examine how large the optical depths can possibly be, we select a somewhat conservative value of $\varphi = 0.02$. This corresponds to $h/r \approx 0.07$ at $r = 100$ au for a star of solar mass and luminosity. The value of $\varphi$ is also chosen to be consistent with the analysis of a subset of the DSHARP sources in Dullemond et al. (2018), which contains a more detailed discussion of the effect of the assumed temperature profile on other derived disk properties. A higher value of
<figure><img src="image_10.png" /><figcaption>**Figure 6.** Estimated 1.25 mm continuum optical depths profiles for the DSHARP sample, ordered by descending stellar luminosity from left to right and top to bottom. The light orange ribbon shows the range of τ values using the 1σ errors for the stellar luminosity from Andrews et al. (2018a). The inner 5 au of each profile is shaded gray to show where beam dilution artificially lowers the derived optical depth. The gray horizontal dotted lines show where τ = 1. The y-axis starts at 10⁻², commensurate with the optical depth value that the DSHARP observations are expected to be sensitive to. For the disks with radial profiles averaged through narrow azimuthal wedges (i.e., HD 142666, Elias 20, DoAr 25, MY Lup, WSB 52, and Sz 114), the optical depth profiles are only plotted out to where the radial intensity profile drops to ~2× the rms, since the noise fluctuations in the profile are larger than for the disks averaged through all azimuthal angles.</figcaption></figure>

φ = 0.05 is selected for the RU Lup disk because lower values
yield dust temperature estimates below the brightness temper-
ature of the inner disk. Assuming typical values for the beam
size (~50 mas), rms (1.5 × 10⁻⁵ Jy beam⁻¹), and disk
temperature (10 K), the DSHARP observations are sensitive
to dust optical depths on the order of 10⁻². Given uncertainties
in φ and L*, the optical depth estimates may be up to a factor of
two off from the true optical depths (there is also a 10%

absolute flux calibration uncertainty, but this is swamped by the
uncertainties in the thermal profile).

The 1.25 mm continuum optical depth profiles of the DSHARP sample are shown in Figure 6. All of the profiles exhibit an apparent decrease in optical depth in the inner few astronomical units, which is likely primarily due to beam dilution. The central dips in the profiles for HD 142666, Sz 129, and HD 143006 are also in part due to bona fide optical
depth decreases corresponding to the small emission cavities detected. The MY Lup and DoAr 25 optical depth profiles have relatively broad dips toward the disk center, but the extent to which the highly inclined viewing geometry contributes to this is not clear. At high inclinations, deprojection becomes less reliable due to disk vertical structure, so multidimensional radiative transfer modeling would be more appropriate to estimate the optical depths of these two disks.

Factoring in the uncertainty in the stellar host luminosity, Figure 6 indicates that the millimeter continuum optical depths of most of the DSHARP disks are plausibly of order unity out to radii of 5–10 au. This is consistent with the finding from Ansdell et al. (2018) that most of the disks observed in the Lupus star-forming region have cores that are optically thick at millimeter wavelengths. IM Lup, RU Lup, and WSB 52 appear to have quite substantial optically thick cores that extend to radii of 15 au or beyond. Cleeves et al. (2016) previously concluded that the inner disk of IM Lup was optically thick at submillimeter wavelengths based on radiative transfer modeling of the 870 µm continuum and the SED. The low continuum optical depth of the HD 143006 rings (aside from the asymmetric feature on B65 that is ~3x brighter than the rest of the ring) makes it unusual among the DSHARP sources. Otherwise, the estimated optical depths at the bright rings in the DSHARP sources typically range from 0.3 to 0.6. Observations of decrements in $^{12}$CO J = 2–1 emission at the locations of the outer bright continuum emission rings for the AS 209 and HD 163296 disks suggest that these substructures may be optically thick despite the moderate estimates shown in Figure 6 (Guzmán et al. 2018; Isella et al. 2018). Dullemond et al. (2018) propose three possible reasons for this discrepancy: the temperatures are overestimated, scattering is significant, or the dust rings consist of optically thick clumps that only appear to be optically thin due to beam dilution.

While the absolute values of the optical depth profiles are highly uncertain, the relative values within each disk are better-constrained because they only depend on the power-law exponent of the temperature profile. Outside the regions where the continuum is optically thick, the optical depth profile roughly traces the product of the dust surface density $\Sigma_d$ and the dust opacity $\kappa_{\nu}$. In the disks around HD 163296, Elias 24, HD 143006, AS 209, SR 4, and DoAr 25, local variations in the radial optical depth profiles reach one to two orders of magnitude, indicating extremely large local variations in the surface density and/or the dust opacity. Large variations in dust opacity would suggest dramatic changes in grain properties across the gap. To verify this, multiwavelength observations to measure the spectral index would be necessary.

# 4. Aggregate Characteristics of Annular Substructures

## 4.1. The Radial Distribution of Annular Substructures

Histograms of the positions of the radial substructures of the DSHARP sample, along with HL Tau and TW Hya, are shown in Figure 7. We also use the scikit-learn implementation of Gaussian kernel density estimation to estimate the distribution. The bandwidth is optimized through leave-one-out cross-validation to ensure that it is narrow enough to capture significant features but wide enough to smooth over spurious features. The resulting estimates are plotted over the histograms.

Annular substructures occur at virtually all radii probed by millimeter emission, as far inward as $r < 5$ au (i.e., essentially within ALMA's resolution limits) and as far out as ~155 au.

<figure><img src="image_11.png" /><figcaption>**Figure 7.** Top: a histogram of the radial positions of the gaps in the DSHARP sample and the HL Tau and TW Hya disks. The Gaussian kernel density estimate (scaled to match the histogram area) is shown in blue. Center: same as above, but for the bright rings. Bottom: plot of the number of disks where $R_{\text{dust}}$ exceeds radius $R$.</figcaption></figure>

The dropoff in the number of substructures detected beyond a radius of ~40 au appears to be a limitation imposed by the overall radial extent of the disk. The bottom of Figure 7 shows a “survival function” for the disk size, plotting the number of disks with $R_{\text{dust}} > R$ over a range of values of $R$. Among the DSHARP disks, it is quite common for annular substructures to occur near the outer edge of the detected millimeter continuum emission. Elias 27 and IM Lup are notable exceptions, with their millimeter continuum emission extending more than 100 au beyond the farthest detected annular substructures. The presence of spiral arms and the faintness of emission in the outer disk, though, may complicate the identification of more distant substructures. Substructures can occur at even larger radii than what is measured for the DSHARP sample—the
V1094 Sco disk has an emission ring located at $r \sim 220$ au (van Terwisga et al. 2018).

Intuitively, one might expect the distribution of substructures to decrease monotonically with radius because the disks have varying sizes and the more compact disks by definition cannot have substructures at large radii. Thus, the slight dropoff in the occurrence of substructures at smaller radii is also interesting to consider. Selection effects must play at least a partial role, because the DSHARP sample excluded known transition disks, which by definition have substructure in the inner disk.

The optical depth profiles in Figure 6 suggest that surface density substructures may be obscured in the intensity profile due to the steep temperature changes in the inner disk. Some disks, such as DoAr 33 and WSB 52, have substructures that appear to be “plateaus” in the radial intensity profiles but correspond to distinct local maxima and minima in the radial optical depth profiles. In addition, although the inner 30 au of the Elias 24 and GW Lup images look featureless, their optical depth profiles show slight dips suggestive of additional substructure. However, more accurate temperature profiles derived through radiative transfer modeling are needed to investigate this possibility further.

Substructures may also be more difficult to detect in the inner disk due to smaller characteristic scales and higher optical depths. The scale height, which increases with distance from the star, influences the width of gaps opened by zonal flows and planets (e.g., Bryden et al. 1999; Johansen et al. 2009). Even if surface density substructures in the inner disk are sufficiently wide, they may still be obscured in the inner disk due to high optical depth. As shown in Figure 6, multiple disks appear to have 1.25 mm continuum optical depths of order unity inside radii of 10 au. Notably, the WSB 52 disk appears to be optically thick throughout much of its extent and has the most marginal substructure detection.

## 4.2. Relationship with Temperature

Using the radial temperature profiles calculated in Section 3.5, Figure 8 shows the normalized, azimuthally averaged intensity profiles as a function of midplane temperature from 80 to 10 K. This temperature range is chosen to include the freezeout temperatures of key disk volatiles, because several works have proposed that annular substructures trace the locations of various molecular snowlines (e.g., Zhang et al. 2015; Okuzumi et al. 2016; Pinilla et al. 2017). To examine whether the DSHARP substructures are connected to snowlines, we shade in the freezeout temperature ranges for CO$_2$, N$_2$, and CO, following the approach in Hollenbach et al. (2009) and using the binding energies and abundances listed in Table 3. While H$_2$O, NH$_3$, and CH$_3$OH are other volatiles expected to be abundant in disks based on cometary measurements and ice observations toward star-forming regions (e.g., Mumma & Charnley 2011; Öberg et al. 2011a), their binding energies are so high that their snowline locations are expected to be below the resolution limit of the DSHARP observations for all but a couple of the most luminous sources (it is worth noting that HD 163296 and HD 142666, the sources with the highest stellar luminosities, have substructures in their inner 5 au).

There does not appear to be a strong correspondence between temperature and substructure positions. More generally, there is no obvious relationship between stellar luminosity (which controls the temperature estimate) and the substructure positions. With respect to specific snowlines, few substructures are observed within the expected vicinity of the CO$_2$ snowline. About half the

disks have substructures that occur in the expected vicinity of the CO and N$_2$ snowlines. However, the possible range of locations for the CO snowline can be quite wide depending on whether CO is largely binding to water ice or pure ice. Moreover, there is a high level of dissimilarity in substructures that occur near the expected locations of the CO and N$_2$ snowlines; some disks have a series of narrow annuli while others have a single wide, deep gap. Thus, it seems unlikely that all of these substructures are directly associated with the CO and N$_2$ snowlines.

## 4.3. The Spacings of Annular Substructures

Figure 9 plots the distribution of the radius ratios of substructures in the DSHARP sample as well as the HL Tau and TW Hya disks. Because of the large concentration of low ratios and the sparse population of very large ratios, Gaussian kernel density estimates are derived for the distribution of log-ratios rather than for the ratios. Leave-one-out cross-validation is again used to optimize the bandwidth. The estimates are then transformed back to ratio distributions. The plots for pairs of gaps and for pairs of rings (including nonadjacent pairs) only include disks that have multiple gaps or multiple rings (i.e., excluding disks such as SR 4). The plot for pairs of all substructures includes pairs consisting of one gap and one ring, so it is not simply a combination of the two plots above.

Within individual disks, the spread in the locations of substructures can be substantial. The innermost and outermost substructures in disks can be separated by more than 100 au, translating to position ratios exceeding 10. Overall, the distributions are dominated by ratios less than 2. Very low ratios ($<$1.1) are not often observed, which is probably at least partly due to angular resolution limits.

Assuming Keplerian rotation and that the disk mass is small compared to the stellar mass, semimajor axis ratios corresponding to low-order mean-motion resonances are overplotted in Figure 9. There is no significant clustering around mean-motion resonances, although some pairs lie near these values. A subset of pairs that are close to resonant locations are identified in Table 4. We list pairs with ratios within 1% of the exact value, excluding substructures with positions that were estimated visually because the uncertainties are large. The B65:B41 pair from HD 143006 is the most intriguing candidate because it corresponds to two adjacent, narrow, and high-contrast rings that are very close to a 2:1 resonance. Fedele et al. (2018) previously suggested a 2:1 resonance that corresponds to B97 and D61 in the AS 209 disk, although D90-B97-D105 were thought to be a single large gap at the time due to the lower angular resolution of their observations. Multiple candidate resonances are identified for the AS 209 and HD 163296 disks, but they may arise coincidentally due to the large number of annular substructures present in these sources. On the flip side, the radial positions of shallow substructures in disks such as MY Lup and HD 142666 have large uncertainties, so it is difficult to ascertain how well their relative positions correspond to expected resonances. In addition, Tamayo et al. (2015) point out that pericenter precession can significantly shift resonant location ratios in massive disk systems. Thus, Table 4 is not exhaustive.

One might also expect substructure pair radius ratios to cluster around characteristic values if the substructures trace molecular snowlines. Assuming a standard power-law expression for the midplane temperature profile, $T_{\text{mid}}(r) = T_0(r/r_0)^{-q}$ (note that Equation (3) can be reformulated in this way), then the radius ratio
<figure><img src="image_12.png" /><figcaption>**Figure 8.** Normalized, azimuthally averaged intensity profiles (in descending order by stellar luminosity) as a function of midplane temperature. Shaded regions show the freezeout temperature ranges for CO₂ (pink), CO on water ice (light gray), CO on pure ice (dark gray), and N₂ on pure ice (orange). The plotted ranges for CO and N₂ on pure ice overlap due to the range of gas densities explored, but the CO snowline will be interior to the N₂ snowline in a given disk.</figcaption></figure>
<table><caption>**Table 3**
Major Icelines in Disk Regions Probed by ALMA</caption><thead><tr><th>Species<br>(1)</th><th>Abundance<br>(2)</th><th><i>E</i><sub>b</sub><br>(K)<br>(3)</th><th>Reference<br>(4)</th><th><i>T</i><sub>frz</sub><br>(K)<br>(5)</th></tr></thead><tbody><tr><td>CO<sub>2</sub></td><td>3e−5<sup>a</sup></td><td>2605<sup>b</sup></td><td>Martín-Doménech et al. (2014)</td><td>57–72</td></tr><tr><td rowspan="2">CO</td><td>5e−5</td><td>1155<sup>c</sup></td><td>Fayolle et al. (2016)</td><td>26–32</td></tr><tr><td>5e−5</td><td>866<sup>b</sup></td><td>Fayolle et al. (2016)</td><td>19–24</td></tr><tr><td>N<sub>2</sub></td><td>4e−5<sup>d</sup></td><td>770<sup>b</sup></td><td>Fayolle et al. (2016)</td><td>17–21</td></tr></tbody></table>

**Notes.** (1) Species. (2) Abundance with respect to total number of H atoms. (3) Binding energy. (4) Reference for $E_b$ value. (5) Freezeout temperatures computed for gas number densities from $10^8–10^{12}$ cm$^{-3}$ (chosen to be representative of midplane values).

<sup>a</sup> Following Öberg et al. (2011b), the CO$_2$ abundance is estimated from ice measurements toward the CBRR 2422.8-3423 disk from Pontoppidan (2006).

<sup>b</sup> Pure ice substrate.

<sup>c</sup> Compact H$_2$O substrate.

<sup>d</sup> Following Piso et al. (2016), we adopt the proto-Sun elemental nitrogen abundance of $8 \times 10^{-5}$ with respect to elemental hydrogen (Lodders 2003) and assume that N$_2$ is the dominant nitrogen carrier.

corresponding to the snowlines of species *i* and *j* is

$$ \frac{r_i}{r_j} = \left( \frac{T_{\text{frz},j}}{T_{\text{frz},i}} \right)^{1/q}, \qquad (5) $$

where $T_{\text{frz}}$ denotes the freezeout temperature. While the analysis in Section 4.2 is subject to large uncertainties in both the thermal profile and the gas number densities, the ratios are more robust to assumptions about these quantities. Using a representative value of $q = 0.5$ (i.e., the value corresponding to the thermal profiles used throughout this work), then $\frac{r_{\text{CO}}}{r_{\text{CO}_2}} \sim 5$ if we adopt the binding energy for CO on water ice and $\sim 9$ using the binding energy for CO on pure ice. As shown in Figure 9, few substructure pairs have these radius ratios. On the other hand, $\frac{r_{\text{N}_2}}{r_{\text{CO}}}$ $\sim 1.3$ using the binding energy for CO on pure ice and $\sim 2.3$ using the binding energy for CO on water ice. A number of substructure pairs have radius ratios between these two values, but not all of them are likely to be associated with the CO and N$_2$ snowlines. First, many of the DSHARP disks have multiple pairs of close substructures. Furthermore, based on “chemical imaging” techniques, the CO snowline has been estimated to lie at $\sim 20–30$ au (Qi et al. 2013; van’t Hoff et al. 2017) in the TW Hya disk and 75 au in the HD 163296 disk (Qi et al. 2015). These two sources bracket the stellar luminosities of the DSHARP sources, which suggests that substructure pairs at radii much smaller than 20 au or much larger than 75 au are also unlikely to be associated with the CO and N$_2$ snowlines.

## 4.4. Relationship with Stellar Parameters

The radial positions of the substructures of the DSHARP sample and TW Hya are plotted as a function of stellar mass ($M_*$), stellar age ($t_*$), and stellar mass accretion rate ($\dot{M}_*$) in Figure 10. The stellar parameters for the DSHARP sample are taken from Table 1 of Andrews et al. (2018a). The adopted parameters for TW Hya are $\log M_*/M_{\odot} = -0.09^{+0.04}_{-0.15}$, $\log t_*/\text{yr} = 6.8 \pm 0.4$, and $\log \dot{M}_*/M_{\odot} \text{yr}^{-1} = -8.6 \pm 0.5$ (Andrews et al. 2018b). HL Tau is excluded from these plots due to the high uncertainty in its

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

**Figure 9.** Top: histogram of the distribution of the position ratios of all pairs of gaps within the same disk. Inset shows the distribution of position ratios between 1 and 4. The Gaussian kernel density estimate (scaled to match the histogram area) is shown in blue. Position ratios corresponding to low-order mean-motion resonances are marked with vertical lines. Center: same as above, but for pairs of rings. Bottom: same as above, but for all pairs of annular substructures (i.e., pairs of gaps, pairs of rings, and pairs of gaps and rings).

stellar properties (e.g., Robitaille et al. 2007). The uncertainties for the positions of substructures not corresponding to well-defined local maxima or minima in the radial intensity profile are estimated as the standard deviation corresponding to the major
<figure><img src="image_16.png" /><figcaption>**Figure 10.** Scatterplots of the radial positions of substructures as a function of stellar mass ($M_*$), stellar mass accretion rate ($\dot{M}_*$), and stellar age ($t_*$) respectively. $1\sigma$ error bars are shown for the stellar parameters and for the feature locations in the disk. Upper limits are shown as arrows.</figcaption></figure>

<table><caption>**Table 4**
Candidate Substructure Pairs in Resonance</caption><thead><tr><th>m:n</th><th>Theoretical <sup>a</sup><sub>m</sub>/<sub>a</sub><sub>n</sub></th><th>Source</th><th>Features</th><th>Measured Ratio</th></tr></thead><tbody><tr><td rowspan="2">2:1</td><td rowspan="2">1.587</td><td>HD 143006</td><td>B65:B41</td><td>1.591 ± 0.011</td></tr><tr><td>AS 209</td><td>B97:D61</td><td>1.590 ± 0.011</td></tr><tr><td>3:1</td><td>2.080</td><td>HD 163296</td><td>B100:D48</td><td>2.08 ± 0.04</td></tr><tr><td rowspan="2">3:2</td><td rowspan="2">1.310</td><td>AS 209</td><td>D137:D105</td><td>1.303 ± 0.011</td></tr><tr><td>Elias 20</td><td>D33:D25</td><td>1.31 ± 0.04</td></tr><tr><td rowspan="5">4:3</td><td rowspan="5">1.211</td><td>AS 209</td><td>B74:D61</td><td>1.220 ± 0.008</td></tr><tr><td>AS 209</td><td>D90:B74</td><td>1.212 ± 0.006</td></tr><tr><td>GW Lup</td><td>D103:B85</td><td>1.206 ± 0.012</td></tr><tr><td>RU Lup</td><td>D29:B24</td><td>1.22 ± 0.05</td></tr><tr><td>AS 209</td><td>B39:B28</td><td>1.39 ± 0.02</td></tr><tr><td rowspan="5">5:3</td><td rowspan="5">1.406</td><td>AS 209</td><td>D137:B97</td><td>1.417 ± 0.012</td></tr><tr><td>DoAr 25</td><td>B137:D98</td><td>1.40 ± 0.01</td></tr><tr><td>HD 163296</td><td>B14:D10</td><td>1.40 ± 0.17</td></tr><tr><td>HD 163296</td><td>B67:D48</td><td>1.40 ± 0.03</td></tr><tr><td>RU Lup</td><td>B34:B24</td><td>1.42 ± 0.06</td></tr></tbody></table>

axis of the synthesized beam. No obvious trend emerges with these stellar properties, although the plots do demonstrate that annular substructures are present in disks across a large region of parameter space. The absence of an obvious trend, though, could be a consequence of the selection for bright continuum sources swamping out substructure variations with stellar properties.

## 5. Discussion

### 5.1. Millimeter versus Scattered Light Substructures in the DSHARP Sources

Recent high angular resolution scattered light observations suggest that annular substructures are common not just in the disk midplane, but also the surface layers. High resolution, high sensitivity scattered light observations have been published for six of the DSHARP targets. Many come from the DARTTS-S

survey by Avenhaus et al. (2018), which used SPHERE/IRDIS to observe in *J* and *H* bands.

In some cases, the scattered light substructures appear in close proximity to the radial locations of the millimeter continuum substructures. From SPHERE scattered light observations of HD 163296 in *J* and *H* bands, Muro-Arena et al. (2018) resolve a ring that matches the location of B67. SPHERE observations in the *J* band of HD 143006 reveal substructures that appear to correspond to D22, B41, D51, and B65 (Benisty et al. 2018). In SPHERE observations of AS 209, Avenhaus et al. (2018) identify possible annular substructures close to B74 and B120. They also identify narrow, ringlike structures in the inner disk that are ascribed to PSF artifacts, but the presence of narrow ringlike structures in this region in the millimeter continuum raises the question of whether there are also genuine scattered light features blended in with the artifacts. A series of rings are identified in scattered light imaging of IM Lup, although Avenhaus et al. (2018) comment that some of the structure could instead be tightly wrapped spiral arms. The scattered light ring located at 0.′′58 (92 au) is close to the apparent end of the millimeter continuum spiral arms. The second scattered light ring occurs at 0.′′96 (152 au), which lies just outside B134.

In other cases, scattered light substructures have no obvious correspondence to millimeter continuum substructures. Avenhaus et al. (2018) detect a scattered light ring structure in the IM Lup disk at a radius of 1.′′52 (240 au), but there is no nearby substructure seen in the DSHARP observations. However, the S/N is low in the outer disk and the *uv* coverage is not optimal for recovering faint, diffuse emission. The scattered light ring at *r* = 2.′′1 (332 au) lies outside the detected extent of the millimeter continuum emission. Likewise, Avenhaus et al. (2018) identify a ring in scattered light at a radius of 0.′′77 (120 au) in the MY Lup disk, which lies outside the detected extent of millimeter continuum emission. Interestingly, although the RU Lup millimeter continuum is highly structured, the disk appears featureless in scattered light (Avenhaus et al. 2018).
<table><caption>**Table 5**
Other Disks with Millimeter Continuum Annular Substructures</caption><thead><tr><th>Source</th><th>SpT</th><th>References<sup>a</sup></th><th>θ<sub>b,min</sub> × d (au)<sup>b</sup></th><th>Notes on Other Millimeter Features</th></tr></thead><tbody><tr><td>HD 97048</td><td>A0-B9</td><td>(1), (2)</td><td>48</td><td>Classified as transition disk</td></tr><tr><td>MWC 480</td><td>A5</td><td>(3), (4)</td><td>18</td><td></td></tr><tr><td>HD 169142</td><td>A5</td><td>(5), (6)</td><td>16</td><td>Classified as transition disk</td></tr><tr><td>MWC 758</td><td>A8</td><td>(7), (8)</td><td>6.2</td><td>Classified as transition disk; high-contrast asymmetry; spiral arm(s)</td></tr><tr><td>V1247 Ori</td><td>F0</td><td>(9), (10)</td><td>13</td><td>Classified as transition disk; high-contrast asymmetry</td></tr><tr><td>HD 142527</td><td>F6</td><td>(11), (12)</td><td>20</td><td>Classified as transition disk; high-contrast asymmetry; circumbinary disk</td></tr><tr><td>SAO 206462</td><td>F8</td><td>(13), (14)</td><td>8.5</td><td>Classified as transition disk; high-contrast asymmetry</td></tr><tr><td>RY Tau</td><td>G0</td><td>(15), (4)</td><td>14</td><td>Classified as transition disk; high-contrast asymmetry</td></tr><tr><td>CI Tau</td><td>K5.5</td><td>(15), (16)</td><td>4.7</td><td></td></tr><tr><td>DL Tau</td><td>K5.5</td><td>(15), (4)</td><td>17</td><td></td></tr><tr><td>GM Aur</td><td>K6</td><td>(15), (17)</td><td>32</td><td>Classified as transition disk</td></tr><tr><td>V1094 Sco</td><td>K6</td><td>(18), (19)</td><td>28</td><td></td></tr><tr><td>PDS 70</td><td>K7</td><td>(20), (21)</td><td>19</td><td>Classified as transition disk</td></tr><tr><td>Sz 98</td><td>K7</td><td>(18), (22)</td><td>39</td><td></td></tr><tr><td>DN Tau</td><td>M0.3</td><td>(15), (4)</td><td>14</td><td></td></tr><tr><td>DS Tau</td><td>M0.4</td><td>(15), (4)</td><td>16</td><td></td></tr><tr><td>AA Tau</td><td>M0.6</td><td>(15), (23)</td><td>15</td><td></td></tr><tr><td>IQ Tau</td><td>M1.1</td><td>(15), (4)</td><td>14</td><td></td></tr><tr><td>UZ Tau E</td><td>M1.9</td><td>(15), (4)</td><td>14</td><td></td></tr><tr><td>GO Tau</td><td>M2.3</td><td>(15), (4)</td><td>16</td><td></td></tr><tr><td>FT Tau</td><td>M2.8</td><td>(15), (4)</td><td>14</td><td></td></tr><tr><td>DM Tau</td><td>M3.0</td><td>(15), (24)</td><td>4.5</td><td>Classified as transition disk</td></tr><tr><td>WSB 82</td><td>...</td><td>—, (25)</td><td>28</td><td>Classified as transition disk</td></tr></tbody></table>

**Notes.**

**References.** (1) Whittet et al. (1987); (2) van der Plas et al. (2017); (3) Mora et al. (2001); (4) Long et al. (2018a); (5) Dunkin et al. (1997); (6) Fedele et al. (2017); (7) Beskrovnaya et al. (1999); (8) Dong et al. (2018b); (9) Vieira et al. (2003); (10) Kraus et al. (2017); (11) Christiaens et al. (2018); (12) Ohashi et al. (2018); (13) Coulson & Walther (1995); (14) Cazzoletti et al. (2018); (15) Herczeg & Hillenbrand (2014); (16) Clarke et al. (2018); (17) Macías et al. (2018); (18) Alcalá et al. (2017); (19) van Terwisga et al. (2018); (20) Pecaut & Mamajek (2016); (21) Long et al. (2018b); (22) Tazzari et al. (2017); (23) Loomis et al. (2017); (24) Kudo et al. (2018); (24) Cox et al. (2017).

## 5.2. Comparison of the DSHARP Sources to Millimeter Continuum Substructures in Other Disks

Besides HL Tau, TW Hya, and the DSHARP sources, annular substructures have been identified in the millimeter continuum of at least 23 other single-disk systems hosted by Class II stars (see Table 5). The list above only includes transition disks if substructures besides the large central cavity have been identified in the millimeter continuum. Otherwise, we refer the reader to van der Marel et al. (2018) for a compilation of all transition disks imaged at millimeter/submillimeter wavelengths up to the end of 2017 and to Long et al. (2018a) and Cieza et al. (2019) for more recent transition disk discoveries.

All spectral types from M through A are represented among disks with substructures. Among the disks observed at millimeter wavelengths prior to DSHARP, annular substructures occur far more frequently than spiral arms or azimuthal asymmetries. The DSHARP sample conforms to this trend. Three sources (WaOph 6, IM Lup, and Elias 27) have spiral arms (see Huang et al. 2018b for specific discussion of the spirals). Two sources, HD 143006 and HD 163296, have crescent-shaped azimuthal asymmetries near well-defined ring substructures. Outside the DSHARP sample, several other disks have also been shown to exhibit annular substructure in conjunction with spiral arms and azimuthal asymmetries (e.g., Cazzoletti et al. 2018; Dong et al. 2018b). Previous observations and DSHARP together suggest a trend of high-contrast azimuthal asymmetries (i.e., azimuthal contrasts larger than a

factor of a few) accompanying ring structures preferentially in disks hosted by stars with spectral types of G and earlier.

Because most disks outside DSHARP have been observed at lower angular resolution, detections of disk substructures have been biased toward high-contrast gaps and rings with widths of tens of astronomical units. The DSHARP observations, as well as a recent Taurus survey at a spatial resolution of ~15 au by Long et al. (2018a), suggest that narrower substructures are more typical. Higher-resolution observations will likely reveal additional annular substructures in the other disks listed in Table 5.

While DSHARP excluded large-cavity transition disks, many have now been observed at moderate-to-high angular resolution (see Table 5). Large-cavity transition disks have often been analyzed as a group distinct from “full” disks (e.g., Pinilla et al. 2018; van der Marel et al. 2018), but the distinction between the two categories ironically becomes blurrier at higher angular resolutions. Substructures detected in “transition” disks exhibit a level of diversity similar to the DSHARP sample. Observed features also include annular gaps and rings, spiral arms, and azimuthal asymmetries (e.g., Fedele et al. 2017; Dong et al. 2018b). Millimeter continuum observations have also begun to detect small inner disks within transition disk cavities (e.g., Boehler et al. 2017; Dong et al. 2018b; Kudo et al. 2018), suggesting that at least some transition disk cavities could be considered (particularly wide) annular gaps. Meanwhile, several of the DSHARP disks turned out to have small central emission deficits. In addition, while the millimeter continuum emission of the SR 4 disk is centrally peaked, it features a deep annular gap in

<footnote><sup>a</sup> First reference is for spectral type, second reference is for millimeter continuum observations of annular substructures.</footnote>

<footnote><sup>b</sup> Distances computed using Gaia parallaxes (Gaia Collaboration et al. 2018).</footnote>
the inner disk and the SED shows a near-infrared deficit that is reminiscent of transition disk SEDs (Andrews et al. 2018a). It is not clear whether these various types of sources represent different stages along the same evolutionary path, but they at least suggest that the “transition” versus “full” disk dichotomy is too simplistic.

## 5.3. The Prevalence of Millimeter Continuum Substructures

### 5.3.1. Prevalence in Ophiuchus and Lupus

The widespread presence of annular substructures in the DSHARP sample naturally raises the question of whether such structures are ubiquitous. Since DSHARP selected for bright disks, it is not straightforward to derive an occurrence rate for millimeter continuum substructures. Nevertheless, interesting lower bounds can be set for the Lupus and Ophiuchus star-forming regions, from which the bulk of the DSHARP sample originates.

The Lupus star-forming region is ~1–3 Myr old and lies ~150–170 pc away (e.g., Merín et al. 2008; Alcalá et al. 2014; Gaia Collaboration et al. 2016, 2018). Based on optical and IR spectra, 96 Class II or flat spectrum protoplanetary disks hosted by stars with $M_* > 0.1 M_\odot$ were identified in the Lupus I-IV clouds and subsequently targeted at millimeter/submillimeter wavelengths at an angular resolution of ~0.′′3 (Ansdell et al. 2016, 2018). Of the 71 disks with millimeter detections, at least 8 exhibit annular substructure: GW Lup, IM Lup, RU Lup, Sz 114, Sz 129, MY Lup, HK Lup, and V1094 Sco (this work, Tazzari et al. 2017; van Terwisga et al. 2018). Thus the fraction of Class II Lupus disks with annular substructure is ≳10%. This lower bound is extremely conservative, considering that the majority of the substructures observed in the Lupus disks are not evident at the angular resolution of the Lupus survey. For perspective, though, 10 large-cavity transition disks have been identified in Lupus (van der Marel et al. 2018). If considered as distinct categories, disks with annular substructure cannot be much rarer than large-cavity transition disks. Taken together, disks with some kind of substructure comprise at least a quarter of all Class II Lupus disks detected at millimeter wavelengths. Moreover, at least among the bright disks, substructures are almost always present. Of the 15 Lupus disks with 1.3 mm continuum fluxes equal to or brighter than the faintest DSHARP disk (DoAr 33), 12 are known so far to have either cavities or annular substructures. Another of the 15, HT Lup, has spiral arms in its primary disk but does not exhibit annular substructure at the angular resolution of DSHARP (Kurtovic et al. 2018).

The c2d Spitzer Legacy project identified 223 flat spectrum and Class II disks in the Ophiuchus star-forming region, which is similar in age to Lupus but features a higher star formation rate and stellar number density (e.g., Evans et al. 2009). The largest millimeter continuum survey so far of Ophiuchus disks, ODISEA, targeted ~130 Class II and flat spectrum disks (and a small number of Class I disks) with ALMA at an angular resolution of ~0.′′2 (Cieza et al. 2019). At least seven Class II Oph disks are known to have annular substructures: SR 4, Elias 20, DoAr 25, Elias 24, Elias 27, DoAr 33, and WSB 82 (this work, Cox et al. 2017), along with an eighth tentative detection in WSB 52. WaOph 6 and AS 209 also have substructures but are only sometimes classified as Ophiuchus sources. Twenty-seven of the flat spectrum/Class II disks in ODISEA have millimeter continuum fluxes equal to or brighter than the faintest DSHARP disk (including DoAr 33 itself). Thirteen of these are known so far to have some kind of millimeter continuum substructure, including cavities and annular substructures. Again, it should be

emphasized that the other sources were observed at much coarser angular resolution than the DSHARP sample, but even the existing observations suggest that the prevalence of millimeter continuum substructures is high among bright disks at least. However, the angular resolution of ODISEA is sufficient to show that it is not common for the Ophiuchus disks to have the kind of deep, wide gaps observed in disks around sources such as AS 209, Elias 24, HD 163296, and HD 143006.

### 5.3.2. Small versus Large Disks

Although the radial extents of the DSHARP disks span an order of magnitude, the selection criteria created a sample in which large disks are overrepresented (Andrews et al. 2018a). The DSHARP sources are larger than about 40% of the 105 disk targets drawn from various star-forming regions (Tripathi et al. 2017; Andrews et al. 2018b) and about three-quarters of the flat spectrum/Class II disks in Ophiuchus (Cieza et al. 2019). Likewise, other disks with annular substructures reported in millimeter continuum emission (see Section 5.2) have sizes similar to the DSHARP sources.

SR 4, WSB 52, and DoAr 33 are the smallest of the DSHARP sources as traced by millimeter continuum emission. While the gap in SR 4 has one of the highest contrasts measured for the DSHARP sources, DoAr 33 and WSB 52 have very low-contrast substructures compared to the larger DSHARP disks. Furthermore, as shown in Figure 7, a comparatively small number of annular substructures are detected in the inner 40 au of the DSHARP disks. In a separate survey, Long et al. (2018a) find indications that larger disks tend to have a larger number of rings. This raises the question of whether the high prevalence of annular substructures in the DSHARP sample signals ubiquity in the overall disk population, or just among larger disks. It has long been a puzzle how large protoplanetary disks persist for megayear timescales when radial drift of millimeter-size grains is expected to act quickly (e.g., Takeuchi & Lin 2002; Brauer et al. 2007). A possible solution is that disk substructures impede radial drift, allowing large disk sizes to be maintained (e.g., Whipple 1972; Pinilla et al. 2012).

Surveys imaging smaller disks will be essential to determine the extent to which disks of more typical sizes resemble the DSHARP sources. However, if the small disks are also optically thick out to radii of $r \sim 10$ au, then emission substructures may not be observable even if surface density variations are present. Thus, non-detections of substructures in small disks will need to be interpreted carefully.

### 5.3.3. Relationship with Age

The median age of the DSHARP sources is ~1 Myr, and the youngest sources have estimated ages of only a few hundred thousand years (with the caveat that age estimates have large uncertainties). The high prevalence of annular substructures suggests that they form relatively quickly in disks. Given how readily substructures have been detected in Class II disks, a natural question is whether some or all of these structures actually formed while the sources were still embedded. Substructures have been reported in a few Class I (or borderline Class I/Class II) disks (ALMA Partnership et al. 2015; Cieza et al. 2016; Sheehan & Eisner 2018), but a systematic high angular resolution study of Class I disks will be necessary to establish whether these substructures are common at this earlier stage.

Not many disks older than a few megayears have been observed at high angular resolution, but several older disks
(e.g., TW Hya, HD 143006, HD 142666, and HD 163296) are known to have annular substructure. This suggests either that these structures are long-lived in disks or that their formation can take place at very different points in the lifetime of a disk. Other than TW Hya, though, the older sources have spectral types of G or earlier and the younger sources are K or M stars, which makes it difficult to isolate trends with disk age.

## 5.4. Possible Origins of Annular Substructures

### 5.4.1. Planets

The most popular hypothesis for the origin of annular substructures in disks is that they result from gravitational interactions between the disk and one or more planets (e.g., Papaloizou & Lin 1984; Paardekooper & Mellema 2004; Fouchet et al. 2010). Of the DSHARP sources, the companion hypothesis has specifically been explored for the AS 209, Elias 24, Elias 27, and HD 163296 disks via hydrodynamical simulations attempting to reproduce previous millimeter continuum observations (e.g., Isella et al. 2016; Meru et al. 2017; Dipierro et al. 2018; Dong et al. 2018a; Fedele et al. 2018; Forgan et al. 2018; Pinte et al. 2018; Teague et al. 2018a), although models of the Elias 27 disk have focused on the spiral arms rather than the annular substructures. Planet-disk interactions involving protoplanets ranging from Saturn to Jupiter mass can reproduce the lower-resolution observations, but the detection of additional substructures in the DSHARP observations makes these inferences worth revisiting. Zhang et al. (2018) perform a parameter space study to estimate the planet masses compatible with the gaps observed in the DSHARP sample.

Although large masses have been inferred for many candidate protoplanets in disks with substructure, recent works have explored the possibility that these substructures are created instead by lower-mass planets (i.e., super-Earths) in low-viscosity disks (e.g., Duffell & MacFadyen 2013; Zhu et al. 2014; Bae et al. 2017; Dong et al. 2017). Whereas each gap has to be opened by a massive planet in higher-viscosity disks, a single planet can open multiple gaps in low-viscosity disks. In particular, Zhu et al. (2014) and Dong et al. (2017) demonstrate that a super-Earth can create a “double gap” feature wherein annular gaps open interior and exterior to its orbit at a distance of about one pressure scale height. If $h/r \sim 0.1$, then potential “double gap” features among DSHARP sources include D90:D105 in AS 209, D74:D98 in DoAr 25, D25:D33 in Elias 20, and D21:D29 in RU Lup. What is particularly interesting about these pairs is that they create “w-features” in the radial intensity profile—that is, the intensity just interior and exterior to the double gap region is comparable to or brighter than the ring separating the two gaps, forming a section in the radial profile that looks like the letter “w.” This is not in general a trait of an arbitrary pair of consecutive gaps in the DSHARP sample, but it is a predicted characteristic of double gaps opened by a planet in a low-viscosity disk.

Hydrodynamical simulations of planet-disk interactions indicate that the eccentricity of annular substructures increases for more massive planets and for lower disk viscosities (e.g., Lubow 1991; Kley & Dirksen 2006; Zhang et al. 2018). Several disks have been proposed to be eccentric based either on offsets between the measured centers of various features (ALMA Partnership et al. 2015; Dong et al. 2018b) or residuals from subtracting an axisymmetric model (Isella et al. 2016, 2018). For the most part, the DSHARP substructures appear to be concentric within the stated uncertainties of the fitted ellipse center positions. The $1\sigma$ uncertainties in $\Delta x$ and $\Delta y$ are

typically 1–3 mas each, translating to a few tenths of an astronomical unit at a typical distance of 140 au. For a face-on disk with well-resolved substructure, offsets of ~1 au between the centers of substructures should be detectable. If taken to represent the distance between the center of an eccentric substructure and the star at one focus, this offset would correspond to an eccentricity of ~0.05 at $r = 20$ au and ~0.01 at $r = 100$ au. In the AS 209 disk, the measured center of B120 appears offset from B74 by 6.5 ± 1.4 mas (0.79 ± 0.17 au), which might suggest a very small eccentricity of ~0.01 at $r = 120$ au. In practice, it is not straightforward to translate measured offsets to eccentricities because vertical structure and optical depth complicate the apparent offsets of substructures in inclined disks (e.g., Pinte et al. 2016), and substructures with the same eccentricity will not be offset from one another. Forward radiative transfer modeling will be useful to search for more subtle signatures of eccentricity. The inner rings of Sz 129 and HD 142666 may be worthy of further investigation due to the asymmetries noted in Section 3.4. Isella et al. (2018) also discuss possible eccentricity in the HD 163296 disk.

Of the substructures in the DSHARP targets, the ones in the AS 209 and HD 163296 disk may present the most intriguing case so far for originating from planet-disk interactions. Substructures that have been observed in CO isotopologue emission (Huang et al. 2016; Isella et al. 2016) may be indicative of gas surface density substructures created by a perturber. Teague et al. (2018a) also find evidence of non-Keplerian gas motion coinciding with the millimeter continuum substructures in the HD 163296 disk. In addition, Guidi et al. (2018) reported the detection of a point source in $L'$ band within the D48 feature of the HD 163296 disk, although further observations will be necessary to establish whether the source is a protoplanet. Further arguments for the planet-disk interactions scenario for these two disks are discussed in Guzmán et al. (2018), Isella et al. 2018, and Zhang et al. (2018).

### 5.4.2. Snowlines

Previously detected substructures in the HL Tau, TW Hya, and HD 163296 disks have been hypothesized to be due to volatile freezeout altering the coagulation and fragmentation properties of dust grains (e.g., Zhang et al. 2015; Okuzumi et al. 2016; Pinilla et al. 2017). Within the DSHARP sample, the evidence does not seem to point to a snowline origin for most of the substructures. First, because the stellar luminosity plays a central role in setting the thermal profile of the disk, the radii of snowline-related substructures would be expected to scale with luminosity. However, as shown in Figures 2 and 8, there is not an obvious correspondence between the stellar luminosities/thermal profiles and the location of substructures. In addition, the wide variations in the number of substructures identified in each disk would seemingly require that the major volatile species differ from disk to disk. Furthermore, as discussed in Section 4.3, the relative spacings of the substructures do not suggest that they are primarily associated with snowlines. Of course, if only a single substructure in each disk is related to a snowline, the relative spacings would not follow a pattern. However, because most disks exhibit multiple annular substructures, this would still mean that the majority of the substructures are unrelated to snowlines.

On the other hand, ascertaining whether an individual substructure is related to a snowline is challenging. To help distinguish between snowline and planetary origins, Pinilla et al. (2017) suggest comparing the millimeter emission profiles to scattered light emission profiles—snowlines should not alter
the gas surface density and would therefore lead to scattered light substructures that are wider and higher contrast than the millimeter substructures. The emission gaps appear to be deeper and wider at millimeter wavelengths for the AS 209, HD 143006, MY Lup, and RU Lup disks compared to scattered light observations, contrary to expectations for snowline-created structure. In a complementary analysis of several DSHARP sources, Dullemond et al. (2018) find evidence for dust trapping in the high-contrast outer rings in the AS 209, HD 163296, GW Lup, and Elias 24 disks based on comparisons between the ring widths and pressure scale heights. This provides a further argument against the snowline origin, which is not expected to create the requisite gas pressure bumps to trap dust. Finally, while the annular substructures appear relatively shallow and narrow in both scattered light and millimeter continuum emission for the IM Lup disk, Cleeves et al. (2018) predict that the CO snowline is located at $r \sim 40$ au, whereas the annular substructures in the IM Lup disk are detected well beyond $r = 100$ au. Even given the large uncertainties associated with the thermal profile, the IM Lup substructures appear to be much too distant from the central star to be associated with a major volatile snowline.

The lack of an obvious association between substructure and expected snowline locations is in line with several other studies of disks. In a reanalysis of the thermal structure of the HL Tau disk based on modeling multiband high resolution ALMA images, Pinte et al. (2016) find that the disk substructures do not match with predicted snowline locations. van Terwisga et al. (2018) show that the annular substructures in the V1094 Sco disk occur at temperatures much lower than the freezeout temperatures of major disk volatiles. Most recently, in an analysis of 12 disks in Taurus, Long et al. (2018a) find that annular substructures have a large spread in radial locations and do not occur preferentially near snowlines.

### 5.4.3. The Role of Internal Gas Dynamics

Internal gas dynamics arising from the coupling between the magnetic field and the disk dust and gas may also give rise to annular substructure formation. Zonal flows due to MRI turbulence or spontaneous concentrations of magnetic flux can lead to gas density variations up to a few tens of percent across radial extents of several pressure scale heights (e.g., Johansen et al. 2009; Bai & Stone 2014; Simon & Armitage 2014), but there remain large uncertainties associated with the global magnetic flux evolution (e.g., Lubow et al. 1994; Guilet & Ogilvie 2014; Okuzumi et al. 2014; Bai & Stone 2017). An additional source of uncertainty is the extent of particle trapping within such gas substructures, so it is not yet clear whether the extremely high-contrast continuum substructures observed in a few of the DSHARP sources (AS 209, HD 143006, HD 163296, Elias 24) can be created through these “magnetic” mechanisms. In the dead zone scenario, the locations of annular substructures in the dust and gas are associated with a transition in the disk’s radial resistivity profile that should occur well under a radius of 100 au (e.g., Dzyurkevich et al. 2010; Flock et al. 2015; Lyra et al. 2015; Ruge et al. 2016). The majority of the DSHARP substructures occur within $r = 100$ au, but a third of these disks have substructures that lie outside the largest expected radial extent of a dead zone. Moreover, Ruge et al. (2016) comment that vortices appear along with the ring structures in the majority of their simulations, whereas most of the DSHARP disks do not have obvious asymmetric features indicative of vortices.

Absent a direct detection of a planet, it is not observationally straightforward to distinguish between planet-disk interactions and magnetohydrodynamic effects (e.g., Flock et al. 2015; Ruge et al. 2016). Studies of molecular lines at high angular and/or spectral resolution could help to address some key uncertainties, with some caveats. First, the strength and morphology of magnetic fields in disks remain unknown, although future measurements of Zeeman splitting could shed some light (e.g., Brauer et al. 2017). Efforts have also been made to constrain the disk ionization fraction via observations of molecular ions (e.g., Öberg et al. 2011c; Cleeves et al. 2015; Teague et al. 2015), but identifying the boundary of the dead zone is not straightforward due to the relatively low angular resolution data available thus far and uncertainties about the abundances of charge carriers that have not been detected. Studies of turbulence via molecular line observations may also help to characterize the overall disk gas dynamics, including the extent of the dead zone, but degeneracies with the thermal profile still pose a major challenge (e.g., Flaherty et al. 2015; Teague et al. 2018b).

### 5.4.4. Other Hypotheses

The list of hypotheses for the origins of substructures in disks is quite extensive. We briefly address a couple additional hypotheses that are also of interest given the characteristics of some of the DSHARP sources.

X-ray photoevaporation has been proposed as a mechanism for dispersing disks from the inside-out (e.g., Ercolano et al. 2008; Owen et al. 2010; Ercolano et al. 2017). While this mechanism is not expected to generate the annular substructures observed in the DSHARP sample, it could be responsible for the small inner gaps/cavities observed in the HD 143006, HD 142666, and Sz 129 disks.

Secular gravitational instability may create annular substructures in disks with low turbulence, low gas-to-dust ratios, and $\mathcal{Q} \lesssim 3$ (e.g., Takahashi & Inutsuka 2014). Few estimates exist for disk turbulence, but the turbulent line broadening upper limits derived for HD 163296 and TW Hya are less than a tenth of the sound speed, which is considered weak (e.g., Flaherty et al. 2015, 2018). Miotello et al. (2017) estimated that the gas-to-dust ratio is less than 10 for RU Lup, GW Lup, Sz 114, MY Lup, and Sz 129, which may make them interesting candidates for secular GI. However, Miotello et al. (2017) also present the alternative scenario that the disks have normal gas-to-dust ratios but are depleted in CO, yielding artificially low gas mass estimates.

Overall, the planet-disk interaction scenario so far seems best able to accommodate the diversity of substructures in disks. Planets can form at any radius where substructures are observed. The diverse exoplanet masses measured offer a straightforward way to explain the large variations observed in gap widths and contrasts. Exoplanetary systems also contain differing numbers of planets, which can help to explain the different number of substructures observed from disk to disk. Even a given planet mass at a given radius can create vastly differing dust substructures depending on the disk viscosity and dust size distribution (e.g., Bae et al. 2018). While certain individual substructures may be associated with snowlines, the snowline scenario does not seem flexible enough to explain the widely varying number and locations of substructures from disk to disk. More constraints on turbulence and magnetic field properties in disks will be needed to determine whether internal gas dynamics are sufficient to explain the presence of annular substructures
without invoking planets. Direct imaging searches for planets and high angular and spectral resolution molecular line observations will provide additional key pieces of information for determining the origins of disk substructures.

## 6. Summary

We present a systematic analysis of annular substructures detected in high angular resolution millimeter continuum observations of 18 single-disk systems from DSHARP. The DSHARP observations suggest that protoplanetary disk substructures are common and thus must trace fundamental evolutionary processes in disks. Our findings are as follows:

1. Substructures are detected in the millimeter continuum of all 18 disks. Annular substructures are identified in at least 17 disks and tentatively in an eighteenth, making them far more common than spiral arms and crescent-like azimuthal asymmetries. A minority of disks also feature other kinds of substructures—spiral arms are observed in three of these sources and crescent-like azimuthal asymmetries in another two.

2. The dust disk morphologies are diverse, with no two disks having similar emission profiles. Substructures can occur at essentially any disk radius probed by ALMA, with features detected in the inner 5 au of a disk and past radii of 150 au. DSHARP detected both shallow, tightly packed gaps and deep, well-separated gaps, sometimes within the same disk.

3. Annular substructures are observed in disks hosted by stars across a range of luminosities, masses, and accretion rates. No immediately obvious relationships emerge yet between the substructure locations and stellar properties, although underlying relationships may be obscured by the selection bias toward bright disks.

4. The relative spacings of the substructures and the lack of an apparent correspondence to the thermal profiles suggest that most annular substructures do not trace the locations of snowlines, although it is possible that some individual features are associated with snowlines.

5. The observed annular substructures are reminiscent of features observed in simulations of planet-disk interactions. In at least a few disks, some morphological characteristics are particularly suggestive, including substructures in possible mean-motion resonances and “double gap” features resembling hydrodynamical predictions for perturbations by a planet in a low-viscosity disk.

This Letter makes use of ALMA data ADS/JAO. ALMA#2016.1.00484.L, ADS/JAO.ALMA#2012.1.00694.S, ADS/JAO.ALMA#2013.1.00226.S, ADS/JAO.ALMA#2013.1.00366.S, ADS/JAO.ALMA#2013.1.00498.S, ADS/JAO.ALMA#2013.1.00631.S, ADS/JAO.ALMA#2013.1.00798.S, ADS/JAO.ALMA#2015.1.00486.S, and ADS/JAO.ALMA#2015.1.00964.S. We thank the referee for comments improving this Letter and NAASC and JAO staff for their advice on data calibration and reduction. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada) and NSC and ASIAA (Taiwan), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO, and NAOJ. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative

agreement by Associated Universities, Inc. J.H. acknowledges support from the National Science Foundation Graduate Research Fellowship under grant No. DGE-1144152. S.A. and J.H. acknowledge support from the National Aeronautics and Space Administration under grant No. 17-XRP17_2-0012 issued through the Exoplanets Research Program. C.P.D. acknowledges support by the German Science Foundation (DFG) Research Unit FOR 2634, grants DU 414/22-1 and DU 414/23-1. A.I. acknowledges support from the National Aeronautics and Space Administration under grant No. NNX15AB06G issued through the Origins of Solar Systems program, and from the National Science Foundation under grant No. AST-1715719. L.P. acknowledges support from CONICYT project Basal AFB-170002 and from FCFM/U. de Chile Fondo de Instalación Académica. V.V.G. and J.C acknowledge support from the National Aeronautics and Space Administration under grant No. 15XRP15_20140 issued through the Exoplanets Research Program. Z.Z. and S.Z. acknowledge support from the National Aeronautics and Space Administration through the Astrophysics Theory Program with grant No. NNX17AK40G and the Sloan Research Fellowship. M.B. acknowledges funding from ANR of France under contract number ANR-16-CE31-0013 (Planet Forming disks). T.B. acknowledges funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme under grant agreement No. 714769. L.R. acknowledges support from the ngVLA Community Studies program, coordinated by the National Radio Astronomy Observatory, which is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc.

**Facility:** ALMA.

**Software:** analysisUtils (https://casaguides.nrao.edu/index.php/Analysis_Utility), AstroPy (Astropy Collaboration et al. 2013), CASA (McMullin et al. 2007), emcee (Foreman-Mackey et al. 2013), matplotlib (Hunter 2007), scikit-image (van der Walt et al. 2014), scikit-learn (Pedregosa et al. 2011), SciPy (Jones et al. 2001).

## Appendix A
Azimuthal Angle Convention

This section describes the azimuthal angle convention used for the deprojected coordinate system throughout the DSHARP

<figure><img src="image_17.png" /><figcaption>**Figure 11.** Diagram of the azimuthal angle convention used in the DSHARP Letter, demonstrated for the GW Lup disk. North is up and east is left.</figcaption></figure>
<figure><img src="image_18.png" /><figcaption>**Figure 12.** Schematic of how quantities are defined for measurements of substructure widths in Section 3.2. Case 1 is the standard scenario. Case 2 represents scenarios in which a gap is exterior to a ring with a peak intensity value less than *I*<sub>mean</sub>. Case 3 represents scenarios in which a ring is interior to a gap with a minimum intensity value greater than *I*<sub>mean</sub>.</figcaption></figure>

Letters. The positive y-axis of the deprojected coordinate system is rotated east of north by the P.A. of the disk. The polar angle increases in the clockwise direction (i.e., if the P.A. of a disk is 0°, then $\theta = 90°$ is in the direction of increasing decl. and $\theta = 0°$ is in the direction of increasing R.A.). This is diagrammed in Figure 11 for the GW Lup disk.

## Appendix B
Substructure Width Definitions

A schematic for substructure width definitions is presented in Figure 12. Case 1 is the standard scenario described in Section 3.2. In Case 2, the peak of a ring immediately interior to a gap has a lower intensity value than $I_{\text{mean}}$, so $r_{d,i}$ is set equal to $r_{b,o}$ for the interior ring. This applies to the B97-D105-B120 sequence in the AS 209 disk and the B72-D77-B85 sequence for the HL Tau disk. In Case 3, the trough of a gap exterior to a ring has a higher intensity value than $I_{\text{mean}}$, so $r_{b,o}$ is set equal to $r_{d,i}$ for the exterior gap. This applies to the D34-B40-D44 sequence in the HL Tau disk.

## ORCID iDs

Jane Huang https://orcid.org/0000-0001-6947-6072

Sean M. Andrews https://orcid.org/0000-0003-2253-2270

Cornelis P. Dullemond https://orcid.org/0000-0002-7078-5910

Andrea Isella https://orcid.org/0000-0001-8061-2207

Laura M. Pérez https://orcid.org/0000-0002-1199-9564

Viviana V. Guzmán https://orcid.org/0000-0003-4784-3040

Karin I. Öberg https://orcid.org/0000-0001-8798-1347

Zhaohuan Zhu https://orcid.org/0000-0003-3616-6822

Shangjia Zhang https://orcid.org/0000-0002-8537-9114

Xue-Ning Bai https://orcid.org/0000-0001-6906-9549

Myriam Benisty https://orcid.org/0000-0002-7695-7605

Tilman Birnstiel https://orcid.org/0000-0002-1899-8783

John M. Carpenter https://orcid.org/0000-0003-2251-0602

David J. Wilner https://orcid.org/0000-0003-1526-7587

## References

Alcalá, J. M., Manara, C. F., Natta, A., et al. 2017, A&A, 600, A20

Alcalá, J. M., Natta, A., Manara, C. F., et al. 2014, A&A, 561, A2

ALMA Partnership, Brogan, C. L., Pérez, L. M., et al. 2015, ApJL, 808, L3

Andrews, S. M., Huang, J., Pérez, L. M., et al. 2018a, ApJL, 869, L41

Andrews, S. M., Terrell, M., Tripathi, A., et al. 2018b, ApJ, 865, 157

Andrews, S. M., Wilner, D. J., Hughes, A. M., Qi, C., & Dullemond, C. P. 2009, ApJ, 700, 1502

Andrews, S. M., Wilner, D. J., Hughes, A. M., Qi, C., & Dullemond, C. P. 2010, ApJ, 723, 1241

Andrews, S. M., Wilner, D. J., Zhu, Z., et al. 2016, ApJL, 820, L40

Ansdell, M., Williams, J. P., Trapman, L., et al. 2018, ApJ, 859, 21

Ansdell, M., Williams, J. P., van der Marel, N., et al. 2016, ApJ, 828, 46

Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33

Avenhaus, H., Quanz, S. P., Garufi, A., et al. 2018, ApJ, 863, 44

Bae, J., Pinilla, P., & Birnstiel, T. 2018, ApJL, 864, L26

Bae, J., Zhu, Z., & Hartmann, L. 2017, ApJ, 850, 201

Bai, X.-N., & Stone, J. M. 2014, ApJ, 796, 31

Bai, X.-N., & Stone, J. M. 2017, ApJ, 836, 46

Banzatti, A., Pinilla, P., Ricci, L., et al. 2015, ApJL, 815, L15

Barenfeld, S. A., Carpenter, J. M., Ricci, L., & Isella, A. 2016, ApJ, 827, 142

Beckwith, S. V. W., Sargent, A. I., Chini, R. S., & Guesten, R. 1990, AJ, 99, 924

Benisty, M., Juhasz, A., Facchini, S., et al. 2018, A&A, 619, A171

Bernstein, G. M., & Jarvis, M. 2002, AJ, 123, 583

Beskrovnaya, N. G., Pogodin, M. A., Miroshnichenko, A. S., et al. 1999, A&A, 343, 163

Béthune, W., Lesur, G., & Ferreira, J. 2017, A&A, 600, A75

Boehler, Y., Weaver, E., Isella, A., et al. 2017, ApJ, 840, 60

Borucki, W. J., Koch, D. G., Basri, G., et al. 2011, ApJ, 736, 19

Brauer, F., Dullemond, C. P., Johansen, A., et al. 2007, A&A, 469, 1169

Brauer, R., Wolf, S., & Flock, M. 2017, A&A, 607, A104

Bryden, G., Chen, X., Lin, D. N. C., Nelson, R. P., & Papaloizou, J. C. B. 1999, ApJ, 514, 344

Cazzoletti, P., van Dishoeck, E. F., Pinilla, P., et al. 2018, A&A, 619, A161

Chiang, E. I., & Goldreich, P. 1997, ApJ, 490, 368

Christiaens, V., Casassus, S., Absil, O., et al. 2018, A&A, 617, A37

Cieza, L. A., Casassus, S., Pérez, S., et al. 2017, ApJL, 851, L23

Cieza, L. A., Casassus, S., Tobin, J., et al. 2016, Natur, 535, 258

Cieza, L. A., Ruíz-Rodríguez, D., Hales, A., et al. 2019, MNRAS, 482, 698

Clarke, C. J., Tazzari, M., Juhasz, A., et al. 2018, ApJL, 866, L6

Cleeves, L. I., Bergin, E. A., Qi, C., Adams, F. C., & Öberg, K. I. 2015, ApJ, 799, 204

Cleeves, L. I., Öberg, K. I., Wilner, D. J., et al. 2016, ApJ, 832, 110

Cleeves, L. I., Öberg, K. I., Wilner, D. J., et al. 2018, ApJ, 865, 155

Coulson, I. M., & Walther, D. M. 1995, MNRAS, 274, 977

Cox, E. G., Harris, R. J., Looney, L. W., et al. 2017, ApJ, 851, 83

D'Alessio, P., Cantò, J., Calvet, N., & Lizano, S. 1998, ApJ, 500, 411

Dipierro, G., Ricci, L., Pérez, L., et al. 2018, MNRAS, 475, 5296

Dong, R., Li, S., Chiang, E., & Li, H. 2017, ApJ, 843, 127

Dong, R., Li, S., Chiang, E., & Li, H. 2018a, ApJ, 866, 110

Dong, R., Liu, S.-y., Eisner, J., et al. 2018b, ApJ, 860, 124

Duffell, P. C. 2015, ApJL, 807, L11

Duffell, P. C., & MacFadyen, A. I. 2013, ApJ, 769, 41

Dullemond, C. P., Birnstiel, T., Huang, J., et al. 2018, ApJL, 869, L46

Dullemond, C. P., Dominik, C., & Natta, A. 2001, ApJ, 560, 957

Dullemond, C. P., & Penzlin, A. B. T. 2018, A&A, 609, A50

Dunkin, S. K., Barlow, M. J., & Ryan, S. G. 1997, MNRAS, 286, 604

Dzyurkevich, N., Flock, M., Turner, N. J., Klahr, H., & Henning, T. 2010, A&A, 515, A70
Ercolano, B., Drake, J. J., Raymond, J. C., & Clarke, C. C. 2008, ApJ, 688, 398

Ercolano, B., Rosotti, G. P., Picogna, G., & Testi, L. 2017, MNRAS, 464, L95

Evans, Neal J. I., Dunham, M. M., Jørgensen, J. K., et al. 2009, ApJS, 181, 321

Fayolle, E. C., Balfe, J., Loomis, R., et al. 2016, ApJL, 816, L28

Fedele, D., Carney, M., Hogerheijde, M. R., et al. 2017, A&A, 600, A72

Fedele, D., Tazzari, M., Booth, R., et al. 2018, A&A, 610, A24

Flaherty, K. M., Hughes, A. M., Rosenfeld, K. A., et al. 2015, ApJ, 813, 99

Flaherty, K. M., Hughes, A. M., Teague, R., et al. 2018, ApJ, 856, 117

Flock, M., Ruge, J. P., Dzyurkevich, N., et al. 2015, A&A, 574, A68

Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306

Forgan, D. H., Ilee, J. D., & Meru, F. 2018, ApJL, 860, L5

Fouchet, L., Gonzalez, J.-F., & Maddison, S. T. 2010, A&A, 518, A16

Fulton, B. J., Petigura, E. A., Howard, A. W., et al. 2017, AJ, 154, 109

Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2016, A&A, 595, A2

Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2018, A&A, 616, A1

Galli, P. A. B., Loinard, L., Ortiz-Léon, G. N., et al. 2018, ApJ, 859, 33

Goodman, J., & Weare, J. 2010, Communications in Applied Mathematics and Computational Science, 5, 65

Guidi, G., Ruane, G., Williams, J. P., et al. 2018, MNRAS, 479, 1505

Guilet, J., & Ogilvie, G. I. 2014, MNRAS, 441, 852

Guzmán, V. V., Huang, J., Andrews, S. M., et al. 2018, ApJL, 869, L48

Herczeg, G. J., & Hillenbrand, L. A. 2014, ApJ, 786, 97

Hollenbach, D., Kaufman, M. J., Bergin, E. A., & Melnick, G. J. 2009, ApJ, 690, 1497

Howard, A. W., Marcy, G. W., Johnson, J. A., et al. 2010, Sci, 330, 653

Huang, J., Andrews, S. M., Cleeves, L. I., et al. 2018a, ApJ, 852, 122

Huang, J., Andrews, S. M., Pérez, L. M., et al. 2018b, ApJL, 869, L43

Huang, J., Öberg, K. I., & Andrews, S. M. 2016, ApJL, 823, L18

Huang, J., Öberg, K. I., Qi, C., et al. 2017, ApJ, 835, 231

Hunter, J. D. 2007, CSE, 9, 90

Isella, A., Guidi, G., Testi, L., et al. 2016, PhRvL, 117, 251101

Isella, A., Huang, J., Andrews, S. M., et al. 2018, ApJL, 869, L49

Johansen, A., Youdin, A., & Klahr, H. 2009, ApJ, 697, 1269

Jones, E., Oliphant, T., Peterson, P., et al. 2001, SciPy: Open Source Scientific Tools for Python

Kanagawa, K. D., Muto, T., Tanaka, H., et al. 2016, PASJ, 68, 43

Kley, W., & Dirksen, G. 2006, A&A, 447, 369

Kraus, S., Kreplin, A., Fukugawa, M., et al. 2017, ApJL, 848, L11

Kudo, T., Hashimoto, J., Muto, T., et al. 2018, ApJL, 868, L5

Kurtovic, N. T., Pérez, L. M., Benisty, M., et al. 2018, ApJL, 869, L44

Lee, E. J., & Chiang, E. 2016, ApJ, 817, 90

Lodders, K. 2003, ApJ, 591, 1220

Long, F., Pinilla, P., Herczeg, G. J., et al. 2018a, ApJ, 869, 17

Long, Z. C., Akiyama, E., Sitko, M., et al. 2018b, ApJ, 858, L12

Loomis, R. A., Öberg, K. I., Andrews, S. M., & MacGregor, M. A. 2017, ApJ, 840, 23

Lubow, S. H. 1991, ApJ, 381, 259

Lubow, S. H., Papaloizou, J. C. B., & Pringle, J. E. 1994, MNRAS, 267, 235

Lyra, W., Turner, N. J., & McNally, C. P. 2015, A&A, 574, A10

Macías, E., Espaillat, C. C., Ribas, Á, et al. 2018, ApJ, 865, 37

Martín-Doménech, R., Muñoz Caro, G. M., Bueno, J., & Goesmann, F. 2014, A&A, 564, A8

McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, K. 2007, in ASP Conf. Ser. 376, Astronomical Data Analysis Software and Systems XVI, ed. R. A. Shaw, F. Hill, & D. J. Bell (San Francisco, CA: ASP), 127

Merín, B., Jørgensen, J., Spezzi, L., et al. 2008, ApJS, 177, 551

Meru, F., Juhász, A., Ilee, J. D., et al. 2017, ApJL, 839, L24

Miotello, A., van Dishoeck, E. F., Williams, J. P., et al. 2017, A&A, 599, A113

Mora, A., Merín, B., Solano, E., et al. 2001, A&A, 378, 116

Mordasini, C., Alibert, Y., Klahr, H., & Henning, T. 2012, A&A, 547, A111

Mumma, M. J., & Charnley, S. B. 2011, ARA&A, 49, 471

Muro-Arena, G. A., Dominik, C., Waters, L. B. F. M., et al. 2018, A&A, 614, A24

Muto, T., Grady, C. A., Hashimoto, J., et al. 2012, ApJL, 748, L22

Nelson, R. P., Papaloizou, J. C. B., Masset, F., & Kley, W. 2000, MNRAS, 318, 18

Öberg, K. I., Boogert, A. C. A., Pontoppidan, K. M., et al. 2011a, ApJ, 740, 109

Öberg, K. I., Murray-Clay, R., & Bergin, E. A. 2011b, ApJL, 743, L16

Öberg, K. I., Qi, C., Wilner, D. J., & Andrews, S. M. 2011c, ApJ, 743, 152

Ohashi, S., Kataoka, A., Nagai, H., et al. 2018, ApJ, 864, 81

Okuzumi, S., Momose, M., Sirono, S.-i., Kobayashi, H., & Tanaka, H. 2016, ApJ, 821, 82

Okuzumi, S., Takeuchi, T., & Muto, T. 2014, ApJ, 785, 127

Owen, J. E., Ercolano, B., Clarke, C. J., & Alexander, R. D. 2010, MNRAS, 401, 1415

Paardekooper, S.-J., & Mellema, G. 2004, A&A, 425, L9

Papaloizou, J., & Lin, D. N. C. 1984, ApJ, 285, 818

Pecaut, M. J., & Mamajek, E. E. 2016, MNRAS, 461, 794

Pedregosa, F., Varoquaux, G., Gramfort, A., et al. 2011, Journal of Machine Learning Research, 12, 2825

Pérez, L. M., Benisty, M., Andrews, S. M., et al. 2018, ApJL, 869, L50

Pérez, L. M., Carpenter, J. M., Andrews, S. M., et al. 2016, Sci, 353, 1519

Pérez, L. M., Carpenter, J. M., Chandler, C. J., et al. 2012, ApJL, 760, L17

Pinilla, P., Birnstiel, T., Ricci, L., et al. 2012, A&A, 538, A114

Pinilla, P., Pohl, A., Stammler, S. M., & Birnstiel, T. 2017, ApJ, 845, 68

Pinilla, P., Tazzari, M., Pascucci, I., et al. 2018, ApJ, 859, 32

Pinte, C., Dent, W. R. F., Ménard, F., et al. 2016, ApJ, 816, 25

Pinte, C., Price, D. J., Ménard, F., et al. 2018, ApJL, 860, L13

Piso, A.-M. A., Pegues, J., & Öberg, K. I. 2016, ApJ, 833, 203

Pontoppidan, K. M. 2006, A&A, 453, L47

Qi, C., Ho, P. T. P., Wilner, D. J., et al. 2004, ApJL, 616, L11

Qi, C., Öberg, K. I., Andrews, S. M., et al. 2015, ApJ, 813, 128

Qi, C., Öberg, K. I., Wilner, D. J., et al. 2013, Sci, 341, 630

Ricci, L., Testi, L., Natta, A., et al. 2010, A&A, 512, A15

Robitaille, T. P., Whitney, B. A., Indebetouw, R., & Wood, K. 2007, ApJS, 169, 328

Romero, G. A., Schreiber, M. R., Cieza, L. A., et al. 2012, ApJ, 749, 79

Rubinstein, A. E., Macías, E., Espaillat, C. C., et al. 2018, ApJ, 860, 7

Ruge, J. P., Flock, M., Wolf, S., et al. 2016, A&A, 590, A17

Sheehan, P. D., & Eisner, J. A. 2018, ApJ, 857, 18

Simon, J. B., & Armitage, P. J. 2014, ApJ, 784, 15

Suriano, S. S., Li, Z.-Y., Krasnopolsky, R., & Shang, H. 2018, MNRAS, 477, 1239

Takahashi, S. Z., & Inutsuka, S.-i. 2014, ApJ, 794, 55

Takeuchi, T., & Lin, D. N. C. 2002, ApJ, 581, 1344

Tamayo, D., Triaud, A. H. M. J., Menou, K., & Rein, H. 2015, ApJ, 805, 100

Tazzari, M., Testi, L., Ercolano, B., et al. 2016, A&A, 588, A53

Tazzari, M., Testi, L., Natta, A., et al. 2017, A&A, 606, A88

Teague, R., Bae, J., Bergin, E. A., Birnstiel, T., & Foreman-Mackey, D. 2018a, ApJL, 860, L12

Teague, R., Henning, T., Guilloteau, S., et al. 2018b, ApJ, 864, 133

Teague, R., Semenov, D., Gorti, U., et al. 2017, ApJ, 835, 228

Teague, R., Semenov, D., Guilloteau, S., et al. 2015, A&A, 574, A137

Tripathi, A., Andrews, S. M., Birnstiel, T., & Wilner, D. J. 2017, ApJ, 845, 44

Tsukagoshi, T., Nomura, H., Muto, T., et al. 2016, ApJL, 829, L35

van Boekel, R., Henning, T., Menu, J., et al. 2017, ApJ, 837, 132

van der Marel, N., van Dishoeck, E. F., Bruderer, S., et al. 2013, Sci, 340, 1199

van der Marel, N., Williams, J. P., Ansdell, M., et al. 2018, ApJ, 854, 177

van der Plas, G., Wright, C. M., Ménard, F., et al. 2017, A&A, 597, A32

van der Walt, S., Schönberger, J. L., Nunez-Iglesias, J., et al. 2014, PeerJ, 2, e453

van Terwisga, S. E., van Dishoeck, E. F., Ansdell, M., et al. 2018, A&A, 616, A88

van't Hoff, M. L. R., Walsh, C., Kama, M., Facchini, S., & van Dishoeck, E. F. 2017, A&A, 599, A101

Vieira, S. L. A., Corradi, W. J. B., Alencar, S. H. P., et al. 2003, AJ, 126, 2971

Whipple, F. L. 1972, in From Plasma to Planet, ed. A. Elvius (New York: Wiley), 211

Whittet, D. C. B., Kirrane, T. M., Kilkenny, D., et al. 1987, MNRAS, 224, 497

Winn, J. N., & Fabrycky, D. C. 2015, ARA&A, 53, 409

Zhang, K., Bergin, E. A., Blake, G. A., et al. 2016, ApJL, 818, L16

Zhang, K., Blake, G. A., & Bergin, E. A. 2015, ApJL, 806, L7

Zhang, S., Zhu, Z., Huang, J., et al. 2018, ApJL, 869, L47

Zhang, Z. 1997, Parameter Estimation Techniques: A Tutorial with Application to Conic Fitting, 15 (Amsterdam: Elsevier)

Zhu, Z., Stone, J. M., Rafikov, R. R., & Bai, X.-N. 2014, ApJ, 785, 122
<figure><img src="image_19.png" /></figure>

# Erratum: “The Disk Substructures at High Angular Resolution Project (DSHARP). II. Characteristics of Annular Substructures” (2018, ApJL, 869, L42)

Jane Huang<fnref n="1" />, Sean M. Andrews<fnref n="1" />, Cornelis P. Dullemond<fnref n="2" />, Andrea Isella<fnref n="3" />, Laura M. Pérez<fnref n="4" />, Viviana V. Guzmán<fnref n="5" />,<fnref n="6" />, Karin I. Öberg<fnref n="1" />, Zhaohuan Zhu<fnref n="7" />, Shangjia Zhang<fnref n="7" />, Xue-Ning Bai<fnref n="8" />, Myriam Benisty<fnref n="9" />,<fnref n="10" />, Tilman Birnstiel<fnref n="11" />, John M. Carpenter<fnref n="5" />, A. Meredith Hughes<fnref n="12" />, Luca Ricci<fnref n="13" />, Erik Weaver<fnref n="3" />, and David J. Wilner<fnref n="1" />

Received 2020 July 21; published 2020 August 4

The second sentence in each of the captions of Figures 2 and 3 of the published article should be “Light blue ribbons show the 1σ scatter at each radial bin divided by the square root of the number of beams spanning the angles over which the intensities are measured” rather than “Light blue ribbons show the 1σ scatter at each radial bin divided by the number of beams spanning the angles over which the intensities are measured.” The figures themselves are unchanged from the published versions. The results and conclusions of the Letter are not affected by this error.

## ORCID iDs

Jane Huang https://orcid.org/0000-0001-6947-6072

Sean M. Andrews https://orcid.org/0000-0003-2253-2270

Cornelis P. Dullemond https://orcid.org/0000-0002-7078-5910

Andrea Isella https://orcid.org/0000-0001-8061-2207

Laura M. Pérez https://orcid.org/0000-0002-1199-9564

Viviana V. Guzmán https://orcid.org/0000-0003-4784-3040

Karin I. Öberg https://orcid.org/0000-0001-8798-1347

Zhaohuan Zhu https://orcid.org/0000-0003-3616-6822

Shangjia Zhang https://orcid.org/0000-0002-8537-9114

Myriam Benisty https://orcid.org/0000-0002-7695-7605

Tilman Birnstiel https://orcid.org/0000-0002-1899-8783

John M. Carpenter https://orcid.org/0000-0003-2251-0602

A. Meredith Hughes https://orcid.org/0000-0002-4803-6200

David J. Wilner https://orcid.org/0000-0003-1526-7587

<footnote n="1">¹Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA; jane.huang@cfa.harvard.edu</footnote>

<footnote n="2">²Zentrum für Astronomie, Heidelberg University, Albert Ueberle Str. 2, D-69120 Heidelberg, Germany</footnote>

<footnote n="3">³Department of Physics and Astronomy, Rice University, 6100 Main Street, Houston, TX 77005, USA</footnote>

<footnote n="4">⁴Departamento de Astronomía, Universidad de Chile, Camino El Observatorio 1515, Las Condes, Santiago, Chile</footnote>

<footnote n="5">⁵Joint ALMA Observatory, Avenida Alonso de Córdova 3107, Vitacura, Santiago, Chile</footnote>

<footnote n="6">⁶Instituto de Astrofísica, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, 7820436 Macul, Santiago, Chile</footnote>

<footnote n="7">⁷Department of Physics and Astronomy, University of Nevada, Las Vegas, 4505 S. Maryland Pkwy, Las Vegas, NV 89154, USA</footnote>

<footnote n="8">⁸Institute for Advanced Study and Tsinghua Center for Astrophysics, Tsinghua University, Beijing 100084, People's Republic of China</footnote>

<footnote n="9">⁹Unidad Mixta Internacional Franco-Chilena de Astronomía, CNRS/INSU UMI 3386, Departamento de Astronomía, Universidad de Chile, Camino El Observatorio 1515, Las Condes, Santiago, Chile</footnote>

<footnote n="10">¹⁰Univ. Grenoble Alpes, CNRS, IPAG, F-38000 Grenoble, France</footnote>

<footnote n="11">¹¹University Observatory, Faculty of Physics, Ludwig-Maximilians-Universität München, Scheinerstr. 1, D-81679 Munich, Germany</footnote>

<footnote n="12">¹²Department of Astronomy, Van Vleck Observatory, Wesleyan University, 96 Foss Hill Drive, Middletown, CT 06459, USA</footnote>

<footnote n="13">¹³Department of Physics and Astronomy, California State University Northridge, 18111 Nordhoff Street, Northridge, CA 91130, USA</footnote>