We present the first astrometry catalog from the Japanese VLBI (very long baseline interferometer) project VERA (VLBI Exploration of Radio Astrometry). We have compiled all the astrometry results from VERA, providing accurate trigonometric annual parallax and proper motion measurements. In total, 99 maser sources are listed in the VERA catalog. Among them, 21 maser sources are newly reported while the rest of 78 sources are referred to previously published results or those in preparation for forthcoming papers. The accuracy in the VERA astrometry are revisited and compared with those from the other VLBI astrometry projects such as BeSSeL (The Bar and Spiral Structure Legacy) Survey and GOBELINS (the Gould's Belt Distances Survey) with the VLBA (Very Long Baseline Array). We have confirmed that most of the astrometry results are consistent with each other, and the largest error sources are due to source structure of the maser features and their rapid variation, along with the systematic calibration errors and different analysis methods. Combined with the BeSSeL results, we estimate the up-to-date fundamental Galactic parameter of $R_{0}=7.92\pm0.16_{\rm{stat.}}\pm0.3_{\rm{sys.}}$~kpc and $\Omega_{\odot}=30.17\pm0.27_{\rm{stat.}}\pm0.3_{\rm{sys.}}$~km~s$^{-1}$~kpc$^{-1}$, where $R_{0}$ and $\Omega_{\odot}$ are the distance from the Sun to the Galactic center and the Sun's angular velocity of the Galactic circular rotation, respectively.

We present searches for 22 GHz H$_2$O masers toward 36 narrow-line Seyfert 1 galaxies (NLS1s), selected from known NLS1s with $v_{sys}$ $\lesssim$ 41000 km/s. Out of the 36 NLS1s in our sample, 11 NLS1s have been first surveyed in our observations, while the observations of other NLS1s were previously reported in literature. In our survey, no new water maser source from NLS1s was detected at 3 $\sigma$ rms level of 8.4 mJy to 144 mJy, which depends on different observing conditions or inhomogeneous sensitivities of each observation using three different telescopes. It is likely that non-detection of new masers in our NLS1 sample is primarily due to insufficient sensitivities of our observations. Including the five known NLS1 masers, the total detection rate of the H$_2$O maser in NLS1s is not remarkably different from that of type 2 Seyfert galaxies or LINERs. However, more extensive and systematic searches of NLS1 would be required for statistical discussion of the detection rate of the NLS1 maser, compared with that of type 2 Seyferts or LINERs.

We report on interferometric observations of a face-on accretion system around the High-Mass young stellar object, G353.273+0.641. The innermost accretion system of 100 au radius was resolved in a 45 GHz continuum image taken with the Jansky-Very Large Array. Our spectral energy distribution analysis indicated that the continuum could be explained by optically thick dust emission. The total mass of the dusty system is $\sim$ 0.2 $M_{\sun}$ at minimum and up to a few $M_{\sun}$ depending on the dust parameters. 6.7 GHz CH$_{3}$OH masers associated with the same system were also observed with the Australia Telescope Compact Array. The masers showed a spiral-like, non-axisymmetric distribution with a systematic velocity gradient. The line-of-sight velocity field is explained by an infall motion along a parabolic streamline that falls onto the equatorial plane of the face-on system. The streamline is quasi-radial and reaches the equatorial plane at a radius of 16 au. This is clearly smaller than that of typical accretion disks in High-Mass star formation, indicating that the initial angular momentum was very small, or the CH$_{3}$OH masers selectively trace accreting material that has small angular momentum. In the former case, the initial specific angular momentum is estimated to be 8 $\times$ 10$^{20}$ ($M_{*}$$/$10 $M_{\sun}$)$^{0.5}$ cm$^{2}$ s$^{-1}$, or a significant fraction of the initial angular momentum was removed outside of 100 au. The physical origin of such a streamline is still an open question and will be constrained by the higher-resolution ($\sim$ 10 mas) thermal continuum and line observations with ALMA long baselines.

We detected flaring flux variability that regularly occurred with the period of 23.9 days on a 6.7 GHz methanol maser emission at Vlsr = 25.30 km/s in G 014.23-00.50 through highly frequent monitoring using the Hitachi 32-m radio telescope. By analyzing data from 05 January 2013 to 21 January 2016, the periodic variability has persisted in at least 47 cycles, corresponding to approximately 1,100 days. The period of 23.9 days is the shortest one observed in masers at around high-mass young stellar objects so far. The flaring component normally falls below the detection limit (3 sigma) of 0.9 Jy. In the flaring periods, the component rises above the detection limit with the ratio of the peak flux density more than 180 in comparison with a quiescent phase, showing intermittent periodic variability. The time-scale of the flux rise was typically two days or shorter, and both symmetric and asymmetric profiles of flux variability were observed through intraday monitoring. These characteristics might be explained by a change in the flux of seed photons by a colliding-wind binary (CWB) system or a variation of the dust temperature by an extra heating source of a shock formed by the CWB system within a gap region in a circumbinary disk, in which the orbital semi-major axes of the binary are 0.26-0.34 au.

We have conducted the first parallax and proper motion measurements of 6.7 GHz methanol maser emission using the Australian Long Baseline Array (LBA). The parallax of G339.884$-$1.259 measured from five epochs of observations is 0.48$\pm $0.08 mas, corresponding to a distance of $2.1^{+0.4}_{-0.3}$ kpc, placing it in the Scutum spiral arm. This is consistent (within the combined uncertainty) with the kinematic distance estimate for this source at 2.5$\pm $0.5 kpc using the latest Solar and Galactic rotation parameters. We find from the Lyman continuum photon flux that the embedded core of the young star is of spectral type B1, demonstrating that luminous 6.7 GHz methanol masers can be associated with high-mass stars towards the lower end of the mass range.

We report on a long-term single-dish and VLBI monitoring for intermittent flare activities of a Dominant Blue-Shifted H$_{2}$O Maser (DBSM) associated with a southern high mass young stellar object, G353.273+0.641. Bi-weekly single-dish monitoring using Hokkaido University Tomakomai 11-m radio telescope has shown that a systematic acceleration continues over four years beyond a lifetime of individual maser features. This fact suggests that the H$_{2}$O maser traces a region where molecular gas is steadily accelerated. There were five maser flares during five-years monitoring, and maser distributions in four of them were densely monitored by the VLBI Exploration of Radio Astrometry (VERA). The overall distribution of the maser features suggests the presence of a bipolar jet, with the 3D kinematics indicating that it is almost face-on (inclination angle of $\sim$ 8$^{\fdg}$--17$^{\fdg}$ from the line-of-sight). Most of maser features were recurrently excited within a region of 100$\times$100 AU$^{2}$ around the radio continuum peak, while their spatial distributions significantly varied between each flare. This confirms that episodic propagations of outflow shocks recurrently invoke intermittent flare activities. We also measured annual parallax, deriving the source distance of 1.70 $^{+0.19}_{-0.16}$ kpc that is consistent with the commonly-used photometric distance.

We report parallaxes and proper motions of three water maser sources in high-mass star-forming regions in the Outer Spiral Arm of the Milky Way. The observations were conducted with the Very Long Baseline Array as part of Bar and Spiral Structure Legacy Survey and double the number of such measurements in the literature. The Outer Arm has a pitch angle of 14.9 +/- 2.7 deg and a Galactocentric distance of 14.1 +/- 0.6 kpc toward the Galactic anticenter. The average motion of these sources toward the Galactic center is 10.7 +/- 2.1 km/s and we see no sign of a significant fall in the rotation curve out to 15 kpc from the Galactic center. The three-dimensional locations of these star-forming regions are consistent with a Galactic warp of several hundred parsecs from the plane.

We report trigonometric parallaxes and proper motions of water masers for 12 massive star forming regions in the Perseus spiral arm of the Milky Way as part of the Bar and Spiral Structure Legacy (BeSSeL) Survey. Combining our results with 14 parallax measurements in the literature, we estimate a pitch angle of 9.9 +/- 1.5 degrees for a section of the Perseus arm. The three-dimensional Galactic peculiar motions of these sources indicate that on average they are moving toward the Galactic center and slower than the Galactic rotation.

Over 100 trigonometric parallaxes and proper motions for masers associated with young, high-mass stars have been measured with the BeSSeL Survey, a VLBA key science project, the EVN, and the Japanese VERA project. These measurements provide strong evidence for the existence of spiral arms in the Milky Way, accurately locating many arm segments and yielding spiral pitch angles ranging from 7 to 20 degrees. The widths of spiral arms increase with distance from the Galactic center. Fitting axially symmetric models of the Milky Way with the 3-D position and velocity information and conservative priors for the solar and average source peculiar motions, we estimate the distance to the Galactic center, Ro, to be 8.34 +/- 0.16 kpc, a circular rotation speed at the Sun, To, to be 240 +/- 8 km/s, and a rotation curve that is nearly flat (a slope of -0.2 +/- 0.4 km/s/kpc) between Galactocentric radii of 5 and 16 kpc. Assuming a "universal" spiral galaxy form for the rotation curve, we estimate the thin disk scale length to be 2.44 +/- 0.16 kpc. The parameters Ro and To are not highly correlated and are relatively insensitive to different forms of the rotation curve. Adopting a theoretically motivated prior that high-mass star forming regions are in nearly circular Galactic orbits, we estimate a global solar motion component in the direction of Galactic rotation, Vsun = 14.6 +/- 5.0 km/s. While To and Vsun are significantly correlated, the sum of these parameters is well constrained, To + Vsun = 255.2 +/- 5.1 km/s, as is the angular speed of the Sun in its orbit about the Galactic center, (To + Vsun)/Ro = 30.57 +/- 0.43 km/s/kpc. These parameters improve the accuracy of estimates of the accelerations of the Sun and the Hulse-Taylor binary pulsar in their Galactic orbits, significantly reducing the uncertainty in tests of gravitational radiation predicted by general relativity.

Very long baseline interferometry (VLBI) monitoring of the 6.7 GHz methanol maser allows us to measure the internal proper motions of the maser spots and therefore study the gas motion around high-mass young stellar objects. To this end, we have begun monitoring observations with the East-Asian VLBI Network. In this paper we present the results of the first epoch observation for 36 sources, including 35 VLBI images of the methanol maser. Since two independent sources were found in three images, respectively, images of 38 sources were obtained. In 34 sources, more than or equal to 10 spots were detected. The observed spatial scale of the maser distribution was from 9 to 4900 astronomical units, and the following morphological categories were observed: elliptical, arched, linear, paired, and complex. The position of the maser spot was determined to an accuracy of approximately 0.1 mas, sufficiently high to measure the internal proper motion from two years of monitoring observations. The VLBI observation, however, detected only approximately 20% of all maser emission, suggesting that the remaining 80% of the total flux was spread into an undetectable extended distribution. Therefore, in addition to high-resolution observations, it is important to observe the whole structure of the maser emission including extended low-brightness structures, to reveal the associated site of the maser and gas motion.

Astrometric Very Long Baseline Interferometry (VLBI) observations of maser sources in the Milky Way are used to map the spiral structure of our Galaxy and to determine fundamental parameters such as the rotation velocity ($\Theta_0$) and curve and the distance to the Galactic center (R$_0$). Here, we present an update on our first results, implementing a recent change in the knowledge about the Solar motion. It seems unavoidable that the IAU recommended values for R$_0$ and $\Theta_0$ need a substantial revision. In particular the combination of 8.5 kpc and 220 \kms\, can be ruled out with high confidence. Combining the maser data with the distance to the Galactic center from stellar orbits and the proper motion of Sgr\,A* gives best values of R$_0$ = 8.3 $\pm$ 0.23 kpc and $\Theta_0$ = 239 or 246 $\pm$ 7 \kms, for Solar motions of V$_ \odot$ = 12.23 and 5.25 \kms, respectively. Finally, we give an outlook to future observations in the Bar and Spiral Structure Legacy (BeSSeL) Survey.

We have conducted VLBA phase-referencing monitoring of H2O masers around the red supergiant, S Persei, for six years. We have fitted maser motions to a simple expanding-shell model with a common annual parallax and stellar proper motion, and obtained the annual parallax as 0.413+/-0.017 mas, and the stellar proper motion as (-0.49+/-0.23 mas/yr, -1.19+/-0.20 mas/yr) in right ascension and declination, respectively. The obtained annual parallax corresponds to the trigonometric distance of 2.42+0.11-0.09 kpc. Assuming the Galactocentric distance of the Sun of 8.5 kpc, the circular rotational velocity of the LSR at the distance of the Sun of 220 km/s, and a flat Galactic rotation curve, S Persei is suggested to have a non-circular motion deviating from the Galactic circular rotation for 15 km/s, which is mainly dominated by the anti rotation direction component of 12.9+/-2.9 km/s. This red supergiant is thought to belong to the OB association, Per OB1, so that this non-circular motion is representative of a motion of the OB association in the Milky Way. This non-circular motion is somewhat larger than that explained by the standard density-wave theory for a spiral galaxy, and is attributed to either a cluster shuffling of the OB association, or to non-linear interactions between non-stationary spiral arms and multi-phase interstellar media. The latter comes from a new view of a spiral arm formation in the Milky Way suggested by recent large N-body/smoothed particle hydrodynamics numerical simulations.

We are using the VLBA and the Japanese VERA project to measure trigonometric parallaxes and proper motions of masers found in high-mass star-forming regions across the Milky Way. Early results from 18 sources locate several spiral arms. The Perseus spiral arm has a pitch angle of 16 +/- 3 degrees, which favors four rather than two spiral arms for the Galaxy. Combining positions, distances, proper motions, and radial velocities yields complete 3-dimensional kinematic information. We find that star forming regions on average are orbiting the Galaxy ~15 km/s slower than expected for circular orbits. By fitting the measurements to a model of the Galaxy, we estimate the distance to the Galactic center R_o = 8.4 +/- 0.6 kpc and a circular rotation speed Theta_o = 254 +/- 16 km/s. The ratio Theta_o/R_o can be determined to higher accuracy than either parameter individually, and we find it to be 30.3 +/- 0.9 km/s/kpc, in good agreement with the angular rotation rate determined from the proper motion of Sgr A*. The data favor a rotation curve for the Galaxy that is nearly flat or slightly rising with Galactocentric distance. Kinematic distances are generally too large, sometimes by factors greater than two; they can be brought into better agreement with the trigonometric parallaxes by increasing Theta_o/R_o from the IAU recommended value of 25.9 km/s/kpc to a value near 30 km/s/kpc. We offer a "revised" prescription for calculating kinematic distances and their uncertainties, as well as a new approach for defining Galactic coordinates. Finally, our estimates of Theta_o and To/R_o, when coupled with direct estimates of R_o, provide evidence that the rotation curve of the Milky Way is similar to that of the Andromeda galaxy, suggesting that the dark matter halos of these two dominant Local Group galaxy are comparably massive.

We performed astrometric observations with the VLBA of WB89-437, an H2O maser source in the Outer spiral arm of the Galaxy. We measure an annual parallax of 0.167 +/- 0.006 mas, corresponding to a heliocentric distance of 6.0 +/- 0.2 kpc or a Galactocentric distance of 13.4 +/- 0.2 kpc. This value for the heliocentric distance is considerably smaller than the kinematic distance of 8.6 kpc. This confirms the presence of a faint Outer arm toward l = 135 degrees. We also measured the full space motion of the object and find a large peculiar motion of ~20 km/s toward the Galactic center. This peculiar motion explains the large error in the kinematic distance estimate. We also find that WB89-437 has the same rotation speed as the LSR, providing more evidence for a flat rotation curve and thus the presence of dark matter in the outer Galaxy.

We present the results of a high sensitivity survey for 6.7 GHz methanol masers towards 22 GHz water maser using the 100 m Efflesberg telescope. A total of 89 sources were observed and 10 new methanol masers were detected. The new detections are relatively faint with peak flux densities between 0.5 and 4.0 Jy. A nil detection rate from low-mass star forming regions enhances the conclusion that the masers are only associated with massive star formation. Even the faintest methanol maser in our survey, with a luminosity of 1.1 $10^{-9} L_\odot$ is associated with massive stars as inferred from its infrared luminosity.

We report phase-referencing VLBA observations of H2O masers near the star-forming region W3(OH) to measure their parallax and absolute proper motions. The measured annual parallax is 0.489 +/- 0.017 milli-arcseconds (2.04 +/- 0.07 kpc), where the error is dominated by a systematic atmospheric contribution. This distance is consistent with photometric distances from previous observations and with the distance determined from CH3OH maser astrometry presented in a related paper. We also find that the source driving the H2O outflow, the ``TW-object'', moves with a 3-dimensional velocity of > 7 km/s relative to the ultracompact HII region W3(OH).

The most powerful tool for measuring distances within our Galaxy is the annual parallax. We carried out phase-referencing VLBI observations of H$_{2}$O masers in the star forming region W3(OH) with respect to the extragalactic continuum source ICRF 0244+624 to measure their absolute proper motions. The measured annual parallax is 0.484 $\pm$ 0.004 milli-arcseconds which corresponds to a distance of 2.07^{+0.01}_{-0.02}$ kpc from the sun. This distance is consistent with photometric and kinematic distances from previous observations.