Water maser variability in a high-mass YSO outburst: VERA and ALMA observations of S255 NIRS 3

Context. Clarifying the relationship between mass accretion and ejection history is one of the key issues in understanding high-mass star formation processes.Aims. We aim to investigate the possible relationship between the mass accretion burst event in mid-June 2015 and the jet ejection in the high-mass protostar S255 NIRS 3.Methods. The Very Long Baseline Interferometer (VLBI) monitoring observations of the 22 GHz H2O masers were carried out using VLBI Exploration of Radio Astrometry (VERA) to reveal the 3D velocity and spatial structure of the outflow/jet traced by the H2O masers in S255 NIRS 3. In addition, we conducted follow-up observations of the submillimeter continuum and the 321 GHz H2O masers with the Atacama Large Millimeter/submillimeter Array (ALMA) at Band 7.Results. We successfully measured the proper motions of the 22 GHz H2O masers associated with a bipolar outflow. The structure is almost the same as was observed in 2005 and 2010. The expansion velocity of the blueshifted bow shock traced by the 22 GHz H2O masers was measured to be 28 km s(-1), corresponding to a dynamical timescale of 60 yr. The direction of the maser outflow is slightly tilted compared with the radio jet, which could suggest a more recent ejection episode during the accretion burst event. The total flux density of the 22 GHz H2O masers gradually increases from the beginning of the VLBI monitoring in early 2017 and becomes almost constant in subsequent single-dish monitoring in 2018. The brightening of the H2O masers is more prominent in the northeast outflow lobe. For the first time, we revealed extended H2O maser emission at 22 GHz in a star-forming region, which is partly resolved out by VERA and even by the most extended Very Large Array (VLA) configurations. We find that the flux variation of such an extended component is similar to that of the unresolved maser emission. The ALMA Band 7 continuum emission did not show significant variations compared with the previous observations performed five months before. We mapped the 321 GHz H2O masers in S255 NIRS 3 providing the fourth example, for this maser, of the spatial distribution in a high-mass star-forming region.Conclusions. We conclude that the bow shock structure traced by the 22 GHz H2O maser features is unlikely to originate at the interface between the radio jet powered by the recent accretion outburst and the surrounding medium. The brightening of the 22 GHz H2O masers could be due to radiative excitation by photons form the (declining) infrared (IR) outburst escaping along the cavity created by the newly ejected material. The lower ratio of the 22 GHz/321 GHz maser luminosity in the blueshifted bow shock suggests a temperature ( >1000 K), higher than for the other maser features in this region.