Comparison of the H₂O maser features detected toward VLA 1 (left panel) in the four EVN epochs (2014.46, 2016.45, 2018.44, and 2020.82) and superimposed to the uniform-weighted continuum map at Q-band (central frequency 44 GHz) of the thermal radio jet driven by VLA 1 obtained with the VLA by Rodríguez-Kamenetzky et al. (2020). The light gray contours are 5, 6, 7, 8, 10, 12, 14 times σ = 100 μJy beam⁻¹ and the VLA beam is shown on the bottom right corner. The positions of the maser features are corrected assuming the proper motion of the region W75N(B) equal to the median proper motion measured for the 6.7 GHz CH₃OH maser features by Rygl et al. (2012), ⟨μ_α⟩ = (−1.97 ± 0.10) mas yr⁻¹ and ⟨μ_δ⟩ = (−4.16 ± 0.15) mas yr⁻¹. In addition, the positions of the maser features associated with VLA 1 and detected in the last three epochs have been further corrected considering the proper motion of VLA 1 within W75N(B). This proper motion is ⟨μ_α⟩^VLA 1 = ( − 6.3 ± 0.4) mas yr⁻¹ and ⟨μ_δ⟩^VLA 1 = ( + 5.0 ± 0.4) mas yr⁻¹ (see Appendix B). The size of the symbols are ten times the uncertainties of the absolute positions of the maser features (see Table 2). For clarity we show the linear fit of the maser features for each epoch (right panel), along with the length of the segments correspond to the position of the maser features used in the fitting. The reference position is α₂₀₀₀ = 20^h38^m36^s​.43399 and δ₂₀₀₀ = +42° 37′34″​​.8710. The parameters of the linear fit are reported in Table 4. The solid light-gray line indicates the orientation of the thermal radio jet (PA = +42°±5°; Rodríguez-Kamenetzky et al. 2020).