Rapidly Varying Anisotropic Methanol (CH$_3$OH) Production in the Inner Coma of Comet 46P/Wirtanen as Revealed by the ALMA Atacama Compact Array

We report the detection of CH$_3$OH emission in comet 46P/Wirtanen on UT 2018 December 8 and 9 using the Atacama Compact Array (ACA), part of the Atacama Large Millimeter/Submillimeter Array (ALMA). These interferometric measurements of CH$_3$OH along with continuum emission from dust probed the inner coma ($<$2000 km from the nucleus) of 46P/Wirtanen approximately one week before its closest approach to Earth ($\Delta$ = 0.089 -- 0.092 au), revealing rapidly varying and anisotropic CH$_3$OH outgassing during five separate ACA executions between UT 23:57 December 7 and UT 04:55 December 9, with a clear progression in the spectral line profiles over a timescale of minutes. We present spectrally integrated flux maps, production rates, rotational temperatures, and spectral line profiles of CH$_3$OH during each ACA execution. The variations in CH$_3$OH outgassing are consistent with Wirtanen's 9 hr nucleus rotational period derived from optical and millimeter wavelength measurements and thus are likely coupled to the changing illumination of active sites on the nucleus. The consistent blue offset of the line center indicates enhanced CH$_3$OH sublimation from the sunward hemisphere of the comet, perhaps from icy grains. These results demonstrate the exceptional capabilities of the ACA for time-resolved measurements of comets such as 46P/Wirtanen.


INTRODUCTION
Comets, kilometer-sized bodies of ice and dust, were assembled in the solar nebula during the era of planet formation and subsequently scattered to the Kuiper disk or the Oort cloud. Stored with minimal processing for the last ∼4.5 Gyr in the deep freeze of the cold outer solar system, they may serve as "fossils" of solar system formation, with the volatile composition of their nuclei reflecting the chemistry and prevailing conditions present in the midplane of the solar nebula (Bockelée- Morvan et al. 2004;Mumma & Charnley 2011;Bockelée-Morvan & Biver 2017). Volatiles are released as they enter the inner solar system (heliocentric distances, r H < 5 au) and encounter increasing solar radiation, activating sublimation and leading to the formation of a coma (a freely expanding atmosphere). Remote sensing of coma volatiles over the past three decades has revealed that comets display significant compositional diversity (e.g., A'Hearn et al. 1995;Crovisier et al. 2009;Cochran et al. 2015;Dello Russo et al. 2016a), providing a window into the diversity of their initial conditions and the degree of chemical processing during their lifetime.
The composition and structure of the coma can be measured in exceptional detail through pure rotational transitions of coma molecules, such as carbon monoxide (CO), carbon monosulfide (CS), hydrogen cyanide (HCN), formaldehyde (H 2 CO), and methanol (CH 3 OH). Single-dish and interferometric observations at millimeter wavelengths (e.g., Bockelée-Morvan et al. 1994;Biver et al. 1999;Fray et al. 2006;Milam et al. 2006;Boissier et al. 2007;Lis et al. 2008;Cordiner et al. 2014) have shown that the distributions of HCN and CH 3 OH in some comets are consistent with direct release from the nucleus, whereas CS, H 2 CO, and HNC are produced by "distributed" sources in the coma. More recently, the high sensitivity and angular resolution of the Atacama Large Millimeter/Submillimeter Array (ALMA) have confirmed these results, and enabled spatially and temporally resolved measurements capable of discerning variations in outgassing on small timescales (Cordiner et al. 2017b). Note-T int is the total on-source integration time. r H, ∆, and φ are the heliocentric distance, geocentric distance, and phase angle (Sun-Comet-Earth), respectively, of Wirtanen at the time of observations. ν is the mean frequency of each instrumental setting. N ants is the number of antennas utilized during each observation, with the range of baseline lengths indicated for each. θmin is the angular resolution (synthesized beam) at ν, and PWV is the mean precipitable water vapor at zenith during the observations.
The 2018 perihelion passage of the Jupiter-family comet (JFC) 46P/Wirtanen afforded an opportunity to couple the high spatial and spectral resolution of ALMA with the most favorable apparition of a JFC in the previous decade. Here we report time-resolved pre-perihelion observations of CH 3 OH in 46P/Wirtanen taken with the ALMA Atacama Compact Array (ACA). The 12 × 7 m antennas of the ACA provide short baselines between 9 and ∼50 m (corresponding to an angular resolution ranging from 5. 45 to 29. 0 at 230 GHz), resulting in greater sensitivity to extended flux than the main 12 m array (which was in configuration C43-4 at the time of our observations with an angular resolution ranging from 0. 39 to 4. 89 at 230 GHz). We note that these observations were most sensitive to the inner coma of 46P/Wirtanen, and that the largest recoverable scale of the ACA is small compared to the size of 46P/Wirtanen's extended coma (up to several arcminutes). We sampled CH 3 OH emission in the coma of 46P/Wirtanen with multiple executions of the ACA spread over two dates, searching for changes in molecular production rates, outgassing, and spatial distribution of the emission. Significant variations in molecular production and anisotropic outgassing of CH 3 OH were measured throughout these data. We report spectrally integrated flux maps, production rates, rotational temperatures, and spectral line profiles of CH 3 OH emission. Section 2 presents the observations and data analysis. Results are found in Section 3. Section 4 discusses the variation in CH 3 OH during these measurements, identifies potential mechanisms, and places these results in the context of other comets characterized to date. We conducted pre-perihelion observations toward Wirtanen on UT 2018 December 8 and 9 during Cycle 6 using the ALMA ACA with the Band 6 receiver, covering frequencies between 241.5 and 242.0 GHz (λ = 1.23 -1.24 mm) in one 500 MHz wide spectral window. The observing log is shown in Table 1. We tracked the comet position using JPL Horizons ephemerides (JPL #K181/21). One correlator setting was employed in each execution to simultaneously sample multiple CH 3 OH transitions and continuum. Weather conditions were good (mean precipitable water vapor at zenith, zenith PWV, = 0.89 -1.09 mm on December 8; zenith PWV = 2.15 -3.23 mm on December 9). Quasar observations were used for bandpass and phase calibration, as well as calibrating Wirtanen's flux scale. The spatial scale (the range in semi-minor and semi-major axes of the synthesized beam) was 3. 98 -7. 99 and the channel spacing was 122 kHz, leading to a spectral resolution of 0.15 km s −1 . The data were flagged, calibrated, and imaged using standard routines in Common Astronomy Software Applications (CASA) package version 5.6.1 (McMullin et al. 2007). We deconvolved the point-spread function with the Högbom algorithm, using natural visibility weighting and a flux threshold of twice the rms noise in each image. The deconvolved images were then convolved with the synthesized beam and corrected for the (Gaussian) response of the ALMA primary beam. We transformed the images from astrometric coordinates to projected cometocentric distances, with the location of the peak continuum flux chosen as the origin, which was in good agreement with the comet's predicted ephemeris position. The topocentric frequency of each spectral channel was converted to cometocentric velocity as where c is the speed of light (km s −1 ), f rest is the rest frequency of a given transition (GHz), f is the frequency of the channel (GHz), and d∆/dt is the topocentric velocity of the comet at the time of the observations (km s −1 ).

RESULTS
We detected molecular emission from multiple strong CH 3 OH transitions near 241.7 GHz as well as continuum emission from dust in the coma of Wirtanen. Molecular line emission was modeled using a three-dimensional radiative transfer method based on the Line Modeling Engine (LIME; Brinch & Hogerheijde 2010) adapted for cometary atmospheres, including a full non-LTE treatment of coma gases, collisions with H 2 O and electrons, and pumping by solar radiation (for further details, see Cordiner et al. 2019a). We calculated the number density of molecules released from the nucleus as a function of distance following the Haser formalism (Haser 1957) as where Q is the production rate (molecules s −1 ), v exp is the expansion velocity (km s −1 ) and β is the photodissociation rate adopted from Huebner et al. (1992), with Q and v exp allowed to vary as free parameters. We assumed that the collision rate between CH 3 OH and H 2 O is the same as that with H 2 (taken from the LAMDA database; Schöier et al. 2005). This assumption has negligible impact on the derived parameters as the coma contained within the ACA beam was close to LTE due to the proximity of Wirtanen and the relatively small beam size (Section 3.2). Our models were run through the CASA "simobserve" tool to mimic the effects of the ACA using the exact antenna configuration and hour angle ranges of our observations. We found that interferometric filtering by the ACA resulted in a 24% loss in the peak flux of our models; therefore, an efficiency factor η = 0.76 was applied to the radiative transfer models when comparing with observed spectra and calculating production rates. Detected spectral lines, including upper-state energies (E u ) and integrated fluxes, are listed in Table 2. December 9, with traces and labels as in Figure 1. Contour intervals are 10% of the projected peak of the gas emission in each map. The rms noise (σ, mJy beam −1 km s −1 ) and contour spacing, δ, for each map are (A): σ = 16, δ = 3.3σ, (B): σ = 18, δ = 3.3σ, and (C): σ = 21, δ = 2.7σ. Continuum emission is shown as a color map behind the contours. Each execution is numbered in the upper left as denoted in Table 1.

Molecular Maps and Extracted Spectra
Figures 1 and 2 show spectrally integrated flux maps for the CH 3 OH 5 0 -4 0 A + transition near 241.791 GHz (the strongest CH 3 OH transition in these data) for each ACA execution in Wirtanen. The shape and spatial extent of the emissions are relatively consistent across each execution.
We extracted nucleus-centered spectra of the 5 0 -4 0 A + transition in each ACA execution. Figure 3 shows the significant variation in the velocity profile (and to an extent, amplitude) measured in each execution. We found that the velocity profiles of the remaining five detected CH 3 OH transitions in each execution (Table 2) were consistent with those of the strongest 5 0 -4 0 A + transition. Taken together, these spectra indicate that CH 3 OH outgassing in Wirtanen was anisotropic, and furthermore, that this anisotropy varied with time (Section 4).  Table 3. Each execution is numbered as denoted in Table 1.

Rotational and Kinetic Temperatures
We detected emission from six strong CH 3 OH transitions in each ACA execution on December 8 and 9. We constrained the coma rotational temperature along the line of sight during each execution using the rotational diagram method (Bockelée-Morvan et al. 1994), and list them in Table 2. Retrieved rotational temperatures are marginally higher on December 8 than on December 9. However, rotational temperatures for executions on a given date are consistent within 1σ uncertainty.
Given the very small geocentric distance (∆ ∼ 0.09 au) of Wirtanen during our observations, the half-width at half-maximum of the ACA synthesized central beam extended to nucleocentric distances of ∼132 × 265 km. This corresponds to the very inner coma of Wirtanen, where the rotational states should be in thermal equilibrium due to collisions at the kinetic temperature of the gas; therefore, the rotational temperature is expected to be similar to the kinetic temperature for the CH 3 OH sampled by our measurements.
We tested this by running full non-LTE models for all six detected CH 3 OH transitions in Execution 4 (the execution with the strongest transitions and highest signal-to-noise ratio, S/N), setting the model kinetic temperature equal to the best-fit rotational temperature for the execution (T rot = 60 ± 5 K; Table 2). We then performed a rotational diagram analysis on the model-generated lines. We find T rot (model) = 55 K, in formal agreement with the rotational temperature calculated from the observed spectra ( Figure 4). Thus, the use of LTE models (i.e., setting the kinetic temperature of the model to the derived rotational temperature for each execution) is justified and results in negligible changes in production rate (8% difference for ±5 K changes in T ). Given that these conditions were satisfied for Execution 4, it is reasonable to assume that they were satisfied for the other executions. These results are in excellent agreement with data obtained with the IRAM 30 m (Biver et al. 2021 Tables 1 and 2), along with the best-fit line (black, solid), corresponding to T rot = 60 ± 5 K, and the 1σ uncertainty in the fit (blue shaded area). Also shown is the rotational diagram for a non-LTE model with T kin = 60 K (red squares) and the best-fit line (red, dashed) corresponding to T rot = 55 K, illustrating formal agreement between the observations and the model.

Expansion Velocities and Anisotropic Outgassing
We derived the expansion velocities and outgassing asymmetry factors for the CH 3 OH 5 0 -4 0 A + transition in each execution by performing nonlinear least-squares fits of our radiative transfer models to extracted spectral line profiles. Hemispheric asymmetry in outgassing along the Sun-comet vector was considered when fitting the observations for each execution. Consistent with previous ALMA studies of comets, we modeled CH 3 OH as a parent species (Cordiner et al. 2017a(Cordiner et al. , 2019b. These models produced excellent fits to the data (reduced chi-square, χ 2 red = 0.97-1.35) (see Figure 3), allowing the spectral line asymmetry to be fully reproduced. Table 3 lists our results, and Figure 3 shows our extracted spectra with best-fit models overlain. Table 3 and Figure 3 show that the line velocity profile for CH 3 OH emission in Wirtanen varied over the course of our observations. As Table 1 2021) reported P = 8.91 ± 0.03 hr, P = 9.19 hr, P = 9 hr, and P = 9.13 hr, respectively. Coupled with this ∼9 hr rotational period, the line morphologies in Figure 3 and the production rates in Table 3 suggest that these effects could be periodic and related to the rotation of Wirtanen's nucleus, as Executions 1 and 4 as well as Executions 2 and 5 would be separated by approximately three complete rotations. These variations may be consistent with changes in coma outgassing as active regions rotated from the day side to the night side of Wirtanen's nucleus. Examining trends in integrated flux (Table 2) and Q(CH 3 OH) (Table 3) compared to trends in weather conditions (Table 1) reveals no discernible link between the weather and the observed variability. We examine potential mechanisms for variability further below.

Variation on Timescales of Minutes
To further address these potential rotational effects, we searched for changes in outgassing on smaller timescales. Each ACA execution consisted of six scans, with each scan lasting ∼8 minutes. The S/N in the CH 3 OH 5 0 -4 0 A + transition in individual scans was not high enough to draw definitive conclusions, so we averaged the line profiles of the three strongest transitions in each scan (namely the 5 0 -4 0 E, 5 −1 -4 −1 E, and 5 0 -4 0 A + transitions) to improve S/N. Figure 5 shows our results for December 8, and Figure 6 shows the same for December 9. Figure 5 shows a clear progression from an asymmetric profile with a blueshifted peak ( Figure 3A) to a redshifted peak ( Figure 3B) on December 8. The progression from scan to scan in Figure 6 on December 9 is less dramatic than on December 8, perhaps owing to the poorer weather (zenith PWV was higher overall on December 9 compared to December 8, and increased over the course of our observations on December 9; see Table 1). Regardless, the progression from a roughly symmetric profile in Execution 3, to an asymmetric profile with a blueshifted peak and higher amplitude in Execution 4, followed by a decrease in amplitude and return to a more symmetric profile in Execution 5, is clearly visible.
To investigate the potential periodicity of the variation in line asymmetry, we examined the mean Doppler shift (velocity first moment; i.e., the average velocity weighted by the emission intensity in each channel) as well as the Figure 5. Averaged spectra of the CH3OH 50-40 E, 5−1-4−1 E, and 50-40 A + transitions in each scan for Execution 1 (Ex. 1, left column) and Execution 2 (Ex. 2, right column) on December 8 (see Table 1), showing the evolution in the velocity profile of CH3OH outgassing over each ∼8 minute scan. maximum amplitude of the blueward peak (taken as the maximum value for v < 0 km s −1 ) in each scan on both dates and plotted them as a function of time (Figure 7). We calculated least-squares fits for sine waves with periods near 3, 6, 9, 12, and 27 hr for both plots and calculated the goodness-of-fit (reduced χ 2 ). The constraint Max[T A ] (K) ≥ 0 excludes P = 13.6 hr and 27.15 hr. Based on the reduced χ 2 , the amplitude of the blue peaks is best fit by a sine wave with P = 6.83 ± 0.04 hr, amplitude A = -0.07 ± 0.01 K, and A 0 = 0.20 ± 0.01 K. Similarly, the mean Doppler shifts are best fit by a sine wave with P = 9.16 ± 0.08 hr, A = -0.15 ± 0.02 km s −1 , and v 0 = -0.16 ± 0.01 km s −1 (when excluding P = 13 or 27 hr). Thus, within the limitations of the data, our results suggest P ∼ 7-9 hr, consistent with measurements from optical and millimeter wavelengths (P ∼ 9 hr; Farnham et al. 2018;Jehin et al. 2018;Handzlik et al. 2019;Farnham et al. 2021). The consistent periodicity between our observations and those from optical and millimeter wavelengths suggests that the variation in CH 3 OH outgassing in Wirtanen reported here is indeed tied to rotational effects. The rotational variation in the velocity structure in Wirtanen is also clearly evident in the position-velocity diagrams for each scan (Figures 8 and 9). . Averaged spectra of the CH3OH 50-40 E, 5−1-4−1 E, and 50-40 A + transitions in each scan for Execution 3 (Ex. 3, left column), Execution 4 (Ex. 4, middle column), and Execution 5 (Ex. 5, right column) on December 9 (see Table 1), showing the evolution in the velocity profile of CH3OH outgassing over each ∼8 minute scan.
In addition to its measured period, the nucleus of Wirtanen is well studied due to its status as the original target for the Rosetta mission. It is a small (r eff ∼ 0.5 km) prolate ellipsoid (axial ratio > 1.4; Boehnhardt et al. 2002). Although previous observations were unable to discern whether Wirtanen is in a state of non-principal axis rotation (Samarasinha et al. 1996), measurements during the 2018 perihelion passage demonstrated that Wirtanen is in a state of simple rotation (Farnham et al. 2021). Given Wirtanen's prolate elliptical nucleus shape and rotational period, coupled with the periodicity measured in our observations, it is possible that the variations in CH 3 OH outgassing are due to differences in illumination on the nucleus, along with active sites rotating from the day side to the night side. Such rotational effects have been observed in a comet with a prolate elliptical nucleus before, namely 103P/Hartley 2. At 103P, the EPOXI mission, coupled with ground-based observations with IRAM and CSO, found that varying illumination of the small lobe of the nucleus due to rotation (along with the nucleus shape) caused significant changes in volatile release (A'Hearn et al. 2011;Drahus et al. 2012;Boissier et al. 2014). Furthermore, at comet 67P/Churyumov-Gerasimenko, the Rosetta mission found diurnal variations in CH 4 , CO, and C 2 H 6 outgassing (Luspay-Kuti et al. 2015;Fink et al. 2016)  Overlaid are best-fit sine waves with periods near 3, 6, 9, 12, and 27 hr, as well as the reduced χ 2 for each fit. Bottom: mean doppler shift (velocity first moment) in each scan on December 8 and 9 as a function of UT time. Overlaid are best-fit sine waves with periods near 3, 6, 9, 12, and 27 hr, as well as the reduced χ 2 for each fit.  reported strong variations in the spectral line profile of CO consistent with a rotating jet. This rotating jet spiraled along with the rotating nucleus and was active during both the day and night. Clearly, nucleus rotation and complex outgassing are strongly linked in multiple comets.

Asymmetric Line Center and Hyperactivity
In light of the rapidly varying asymmetry in Wirtanen's CH 3 OH spectral line profile, an equally striking feature is the consistent blueward offset of the line center observed in every spectrum across the two dates of our measurements, including in individual 8 minute scans. Coupled with the higher expansion velocity in the sunward hemisphere compared to the anti-sunward hemisphere (Table 3), this suggests that there is an excess of CH 3 OH outflow on the sunward side of the nucleus. Furthermore, the redward side of the profile is remarkably stable across all of our observations, even including Execution 2 on December 8 (Figure 5), where it retains a relatively constant amplitude as the blueward peak continually diminishes. This implies that the overall variability in CH 3 OH production may be largely due to activity occurring on the dayward side. This is consistent with optical observations of CN in Wirtanen that indicated two sources of CN production: one which remained active throughout a rotation, and a second which turned on and off as the nucleus rotated (Farnham et al. 2021).
Wirtanen is considered a hyperactive comet (e.g., Lis et al. 2019, and references therein), with an unusually high Q(H 2 O) for its small nucleus size attributed to a significant release of H 2 O from icy grains in the coma. Observations of H 2 O hot bands near 2.9 µm with NIRSPEC-2 at the W. M. Keck Observatory (Martin et al. 2016(Martin et al. , 2018 on December 17 and 18 ) and iSHELL at the NASA-IRTF (Rayner et al. 2012(Rayner et al. , 2016 on December 6, December 18 (Roth et al. 2021), and December 21 (Khan et al. 2021) found significant enhancements in the spatial profiles of Cometocentric Velocity (km s -1 ) Figure 8. Position-velocity diagrams of CH3OH in Wirtanen for each scan on December 8, taken parallel to the solar vector and with the spatial origin at the peak emission position for each scan. Each diagram shows the cometocentric velocity as a function of distance from the image center, with each pixel color-coded by the line flux, thereby showing the evolution in velocity structure across each scan. Scans are labeled and ordered as in Figure 5. Negative offsets are in the anti-sunward direction.
H 2 O emission, supporting the icy grain contribution to Wirtanen's H 2 O production (as much as 40% from extended sources; Bonev et al. 2021).
Given its hyperactivity, an intriguing possible explanation for the consistent blueward offset in the CH 3 OH line center is enhanced sublimation on the sunward side of the nucleus, perhaps from icy grains in the coma. Although the presence of icy grains would change the absolute Q(CH 3 OH) for each execution (compared to the parent species models used here), the trends in production rate and line morphology would remain the same. This is consistent with measurements of Wirtanen at other wavelengths during its 2018 perihelion passage. Cordiner et al. (2019b) reported extended CH 3 OH release in Wirtanen on UT 2018 December 7, most likely from an icy grain source, based on ALMA observations with the main 12 m array. Wang et al. (2020) found that the HCN (J = 1-0) transition near 3.4 mm displayed a strong blueward offset in its line center on December 14 and 15, consistent with anisotropic outgassing on its sunward side. Furthermore, simultaneous iSHELL measurements of CH 3 OH, OH* (prompt emission, a direct tracer of H 2 O production; see Bonev et al. 2006), and dust continuum taken on December 6 found that CH 3 OH and H 2 O displayed significant peaks in emission in the sunward-facing hemisphere ( Figure 10). Post-perihelion measurements Cometocentric Velocity (km s -1 ) Figure 9. Position-velocity diagrams of CH3OH in Wirtanen for each scan on December 9, taken parallel to the solar vector and with the spatial origin at the peak emission position for each scan. Each diagram shows the cometocentric velocity as a function of distance from the image center, with each pixel color-coded by the line flux, thereby showing the evolution in velocity structure across each scan. Scans are labeled and ordered as in Figure 6. Negative offsets are in the anti-sunward direction.
with NIRSPEC-2 ) and iSHELL (Roth et al. 2021;Khan et al. 2021) showed that the enhancements in CH 3 OH and H 2 O persisted, although now in the projected anti-sunward direction (but at a considerably smaller phase angle, φ ∼ 18 • ). If the CH 3 OH enhancement in Wirtanen is indeed owing to icy grain sublimation, it is not the first time that icy grain contributions to multiple species have been found in a hyperactive comet. In comet 103P/Hartley 2, Drahus et al. (2012) found that HCN displayed a varying blueshifted peak corresponding to production from a jet, while the remaining HCN production was attributed to icy grains. In contrast, CH 3 OH did not display a blueshifted peak, suggesting the CH 3 OH was produced solely from icy grains. Measurements by Boissier et al. (2014) also indicated the production of HCN and CH 3 OH from icy grains in 103P/Hartley 2, implied by the phase offset between the CO 2 production curve (produced purely from nucleus sources) and the HCN, CH 3 OH, and H 2 O production curves. The similarities between these two comets suggest that their outgassing behaviors and sources may have been similar, and emphasize the complex comae of hyperactive comets. Future work analyzing our full ALMA data set of observations toward Wirtanen will enable us to clarify the extent to which icy grain sublimation contributed to CH 3 OH production in its coma.

CONCLUSION
The historic 2018 perihelion passage of Wirtanen together with the capabilities of ALMA enabled us to measure its inner coma with high spatial, spectral, and temporal resolution. Our ALMA ACA observations presented here revealed rapidly varying and anisotropic CH 3 OH production in the inner coma of Wirtanen on two pre-perihelion dates shortly before its closest approach to Earth, with the evolution of the spectral line asymmetry clearly visible on a timescale of minutes. We found that the variations in CH 3 OH were periodic, in good agreement with the nucleus rotational period derived from optical and millimeter wavelength measurements (P ∼ 9 hr), and consistent with hemispheric outgassing along the Sun-comet vector, suggesting variations in outgassing owing to changing illumination of active sites on the nucleus. Furthermore, the consistent blueward offset of the CH 3 OH line center, coupled with the higher expansion velocity in the sunward hemisphere, suggested enhanced CH 3 OH sublimation in this direction, perhaps from icy grains. The work presented here demonstrates that time-resolved studies of the composition, spatial extent, and outgassing profiles of coma species are possible with the high sensitivity of ALMA on timescales as small as minutes, and lays the groundwork for future detailed cometary coma studies using ALMA ACA.