In 2019 September, a sudden flare of the 6.7 GHz methanol maser was observed toward the high-mass young stellar object (HMYSO) G24.33+0.14. This may represent the fourth detection of a transient mass accretion event in an HMYSO after S255IR NIRS3, NGC 6334I-MM1, and G358.93−0.03-MM1. G24.33+0.14 is unique among these sources as it clearly shows a repeating flare with an 8 yr interval. Using the Atacama Large Millimeter/submillimeter Array (ALMA), we observed the millimeter continuum and molecular lines toward G24.33+0.14 in the pre-flare phase in 2016 August (ALMA Cycle 3) and the mid-flare phase in 2019 September (ALMA Cycle 6). We identified three continuum sources in G24.33+0.14, and the brightest source, C1, which is closely associated with the 6.7 GHz maser emission, shows only a marginal increase in flux density with a flux ratio (Cycle 6$/$Cycle 3) of 1.16 ± 0.01, considering an additional absolute flux calibration uncertainty of $10\%$. We identified 26 transitions from 13 molecular species other than methanol, and they exhibit similar levels of flux differences with an average flux ratio of 1.12 ± 0.15. In contrast, eight methanol lines observed in Cycle 6 are brighter than those in Cycle 3 with an average flux ratio of 1.23 ± 0.13, and the higher excitation lines tend to show a larger flux increase. If this systematic increasing trend is real, it would suggest radiative heating close to the central HMYSO due to an accretion event which could expand the size of the emission region and/or change the excitation conditions. Given the low brightness temperatures and small flux changes, most of the methanol emission is likely to be predominantly thermal, except for the 229.759 GHz (8−1–70 E) line known as a class I methanol maser. The flux change in the millimeter continuum of G24.33+0.14 is smaller than in S255IR NIRS3 and NGC 6334I-MM1 but is comparable with that in G358.93−0.03-MM1, suggesting different amounts of accreted mass in these events.

Context. Recent studies have shown that 6.7 GHz methanol maser flares can be a powerful tool for verifying the mechanisms of maser production and even the specific signatures of accretion rate changes in the early stages of high-mass star formation. Aims. We characterize the spatial structure and evolution of methanol and water masers during a flare of methanol maser emission at 6.7 GHz in the high-mass young stellar object (HMYSO) G24.33+0.14. Methods. Very Long Baseline Array (VLBA) was used to image the 6.7 and 12.2 GHz methanol and 22.2 GHz water vapor masers at three epochs guided by monitoring the methanol line with the Torun 32m telescope. The 6.7 GHz maser maps were also obtained with the European VLBI Network (EVN) and Long Baseline Array (LBA) during the flare. The Wide-field Infrared Survey Explorer (WISE) data were used to find correlations between the 6.7 GHz maser and infrared (IR) fluxes. Results. The 6.7 GHz methanol maser cloudlets are distributed over ~3500 au, and the morphology of most of them is stable although their brightness varies following the course of the total flux density on a timescale of two months. The 12.2 GHz methanol maser cloudlets cover an area an order of magnitude smaller than that of 6.7 GHz emission, and both transitions emerge from the same masing gas. The 22.2 GHz maser cloudlets lie in the central region and show a systematic increase in brightness and moderate changes in size and orientation, together with the velocity drift of the strongest cloudlet during two months of the Very Long Baseline Interferometry (VLBI) observing period. Time lag estimates imply the propagation of changes in the physical conditions of the masing region with a subluminal speed (~0.3c). A tight correlation of IR (4.6 μm) and 6.7 GHz flux densities is found, supporting the radiative pumping model. Proper motion analysis does not reveal any signs of expansion or inflow of the methanol cloudlets within ~6 mas over ~10 yr. Comparison with the 230 GHz Atacama Large Millimeter Array (ALMA) data indicates that the methanol masers are distributed in the inner part of the rotating disk, whereas the 22.2 GHz emission traces the compact inner component of the bipolar outflow or a jet structure. Conclusions. The maser morphology in the target is remarkably stable over the course of the flare and is similar to the quiescent state, possibly due to less energetic accretion events that can repeat on a timescale of ~8 yr.

The variability study of 6.7 GHz methanol masers has become a useful way to improve our understanding of the physical conditions in high-mass star-forming regions. Based on the single-dish monitoring using the Irbene telescopes, we selected three sources with close sky positions.We imaged them using the European Very Long Baseline Interferometer Network and searched available data on VLBI archives to follow detailed changes in their structures and single maser spot variability. All three targets show a few groups of maser cloudlets of a typical size of 3.5 mas and the majority of them show linear or arched structures with velocity gradients of order 0.22 km s−1 mas−1. The cloudlets and overall source morphologies are remarkably stable on time scales of 7-15 yr supporting a scenario of variability due to changes in the maser pumping rate.