The knots in solar prominences are often observed to fall with nearly constant velocity, but the associated physical mechanism is currently not well understood. In this letter, we presented a prominence observed by New Vacuum Solar Telescope (NVST) in H-alpha wavelength. Knots that rose within the prominence appear to have been preferentially located at higher altitude, whereas those that fell were found throughout the entire prominence structure. The descending speed of the knots near the solar surface was higher than that far away from the solar surface. We noted that the knots near the solar surface may run along a set of coronal loops observed from the Atmospheric Imaging Assembly. Elsewhere, the majority of knots are interpreted to have descended across more horizontal magnetic field with a nearly constant speed. This lack of acceleration indicates that the liberated gravitational potential energy may not manifest as an increase in kinetic energy. Assuming instead that the descending knots were capable of exciting Alfven waves that could then dissipate within the local corona, the gravitational potential energy of the knots may have been converted into thermal energy. Assuming a perfectly elastic system, we therefore estimate that the gravitational energy loss rate of these observed knots amounts to 1/2000 of that required to heat the entire quiet-Sun, increasing to 1/320 when considering possibly further downward motions of the knots having disappeared in the H-alpha observations. This result suggests such a mechanism may contribute to the heating of the corona local to these prominences.

For filter-based magnetographs, the linear calibration method under the weak-field assumption is usually adopted; this leads to magnetic saturation effect in the regions with strong magnetic field. This article explores a new method to overcome the above disadvantage using a multilayer perceptron network, which we call MagMLP, based on a back-propagation algorithm with one input layer, five hidden layers, and one output layer. We use the data from the \textit{Spectropolarimeter} (SP) on board \textit{Hinode} to simulate single-wavelength observations for the model training, and take into account the influence of the Doppler velocity field and the filling factor. The training results show that the linear fitting coefficient (LFC) of the transverse field reaches above 0.91, and that of the longitudinal field is above 0.98. The generalization of the models is good because the corresponding LFCs are above 0.9 for the test subsets. Compared with the linear calibration method, the MagMLP is much more effective on dealing with the magnetic saturation effect. Analyzing an active region, the results of the linear calibration present an evident magnetic saturation effect in the umbra regions; the corresponding systematic error reaches values greater than 1000 G in most areas, or even exceeds 2000 G at some pixels. However, the results of MagMLP at these locations are very close to the inversion results, and the systematic errors are basically within 300 G. In addition, we find that there are many "bright spots" and "dark spots" on the inclination angle images from the inversion results of \textit{Hinode}/SP with values of 180 and 0 degrees, respectively, where the inversion is not reliable and does not produce a good result; the MagMLP handles these points well.

This paper is devoted to the primary spectro-polarimetric observation performed at the New Vacuum Solar Telescope of China since 2017, and our aim is to precisely evaluate the real polarimetric accuracy and sensitivity of this polarimetry by using full Stokes spectro-polarimetric observations of the photospheric line Fe I 532.4 nm. In the work, we briefly describe the salient characteristic of the NVST as a polarimeter in technology and then characterize its instrumental polarization based on the operation in 2017 and 2019. It is verified that the calibration method making use of the instrumental polarization calibration unit (ICU) is stable and credible. The calibration accuracy can reach up to 3$\times 10^{-3}$ . Based on the scientific observation of the NOAA 12645 on April 5th, 2017, we estimate that the residual cross-talk from Stokes $I$ to Stokes $Q$, $U$ and $V$, after the instrumental polarization calibration, is about 4$\times10^{-3}$ on average, which is consistent with the calibration accuracy and close to the photon noise. The polarimetric sensitivity (i.e., the detection limit) for polarized light is of the order of $10^{-3}$ with an integration time over 20 seconds. Slow modulation rate is indeed an issue for the present system. The present NVST polarimeter is expected to be integrated with an high-order adaptive optics system and a field scanner to realize 2D magnetic field vector measurements in the following instrumentation update.

We developed a spectroscopic monitoring project to investigate the kinematics of the broad-line region (BLR) in active galactic nuclei (AGN) with ultra-fast outflows (UFOs). Mrk~79 is a radio-quiet AGN with UFOs and warm absorbers, had been monitored by three reverberation mapping (RM) campaigns, but its BLR kinematics is not understood yet. In this paper, we report the results from a new RM-campaign of Mrk~79, which was undertaken by Lijiang 2.4-m telescope. Mrk~79 is seeming to come out the faint state, the mean flux approximates a magnitude fainter than historical record. We successfully measured the lags of the broad emission lines including H$\beta~\lambda4861$, H$\gamma~\lambda4340$, He II $\lambda4686$ and He I $\lambda5876$ with respect to the varying AGN continuum. Based on the broad H$\beta~\lambda4861$ line, we measured black hole (BH) mass of $M_{\bullet}=5.13^{+1.57}_{-1.55}\times10^{7}M_{\odot}$, estimated accretion rates of ${\dot{M}_{\bullet}}=(0.05\pm0.02)~L_{\rm Edd}~c^{-2}$, indicating that Mrk~79 is a sub-Eddington accretor. We found that Mrk~79 deviates from the canonical Radius$-$Luminosity relationship. The marginal blueshift of the broad He II $\lambda4686$ line detected from rms spectrum indicates outflow of high-ionization gas. The velocity-resolved lag profiles of the broad H$\gamma~\lambda4340$, H$\beta~\lambda4861$, and He I $\lambda5876$ lines show similar signatures that the largest lag occurs in the red wing of the lines then the lag decreases to both sides. These signatures should suggest that the BLR of Keplerian motion probably exists the outflow gas motion. All findings including UFOs, warm absorbers, and the kinematics of high- and low-ionization BLR, may provide an indirect evidence that the BLR of Mrk~79 probably originates from disk wind.

We study the early evolution of a hot-channel-like magnetic flux rope (MFR) toward eruption. Combining with imaging observation and magnetic field extrapolation, we find that the hot channel possibly originated from a pre-existing seed MFR with a hyperbolic flux tube (HFT). In the precursor phase, three-dimensional tether-cutting reconnection at the HFT is most likely resulting in the heating and buildup of the hot channel. In this process, the forming hot channel was rapidly enlarged at its spatial size and slipped its feet to two remote positions. Afterward, it instantly erupted outwards with an exponential acceleration, leaving two core dimmings near its feet. We suggest that pre-flare reconnection at the HFT played a crucial role in enlarging the seed MFR and facilitating the onset of its final solar eruption. Moreover, a recently predicted drifting of MFR's footpoints was detected at both core dimmings. In particular, we find that MFR's west footpoint drift was induced by a new reconnection geometry among the erupting MFR's leg and thereby inclined arcades. As MFR's west footpoints gradually drifted to a new position, a set of newborn atypical flare loops connected into the west core dimming, causing a rapid decrease of dimmed area inside this core dimming and also generating a secondary flare ribbon at their remote feet. This reveals that core dimmings may suffer a pronounced diminishment due to the eruptive MFR's footpoint drift, implying that mapping the real footpoints of the erupting MFR down to the Sun's surface is more difficult than previously thought.

Lijiang 2.4-meter Telescope(LJT), the largest common-purpose optical telescope in China, has been applied to the world-wide astronomers since 2008. It is located at Gaomeigu site, Lijiang Observatory(LJO), the southwest of China. The site has very good observational conditions. Since 10-year operation, several instruments have been equipped on the LJT. Astronomers can perform both photometric and spectral observations. The main scientific goals of LJT include photometric and spectral evolution of supernova, reverberation mapping of active galactic nucleus, physical properties of binary star and near-earth object(comet and asteroid), identification of exoplanet, and all kinds of transients. Until now, the masses of 41 high accretion rate black holes have been measured, and more than 168 supernova have been identified by the LJT. More than 190 papers related to the LJT have been published. In this paper, the general observation condition of the Gaomeigu site is introduced at first. Then, the LJT structure is described in detail, including the optical, mechanical, motion and control system. The specification of all the instruments, and some detailed parameters of the YFOSC is also presented. Finally, some important scientific results and future expectations are summarized.

Dynamical jets are generally found on Light bridges (LBs), which are key to studying sunspots decays. So far, their formation mechanism is not fully understood. In this paper, we used state-of-the-art observations from the Goode Solar Telescope, the Interface Region Imaging Spectrograph, the Spectro-Polarimeter on board Hinode and the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory to analyze the fan shape jets on LBs in detail. Continuous upward motion of the jets in ascending phase is found from the H$\alpha$ velocity, which lasts for 12 minutes and is associated with the H$\alpha$ line wing enhancements. Two mini jets appear upon the bright fronts of the fan shape jets visible in the AIA 171 {\AA} and 193 {\AA} channels, with a time interval as short as 1 minute. Two kinds of small scale convective motions are identified in the photospheric images, along with the H$\alpha$ line wing enhancements. One seems to be associated with the formation of a new convection cell and the other manifests as the motion of a dark lane passing through the convection cell. The finding of three lobes Stokes V profiles and their inversion with NICOLE code indicates that there is magnetic field lines with opposite polarities in LBs. From the H$\alpha$ -0.8 {\AA} images, we found ribbon like brightenings propagating along the LBs, possibly indicating slipping reconnection. Our observation supports that the fan shape jets under study are caused by the magnetic reconnection and photospheric convective motions play an important role in triggering the magnetic reconnection.

Using the 135-second cadence of the photospheric vector data provided by the Helioseismic and Magnetic Imager telescope on board the Solar Dynamic Observatory, we examined the time-evolution of magnetic helicity fluxes across the photosphere during 16 flares with the energy class lower than M5.0. During the flare in 4 out of 16 events, we found impulsive changes in the helicity fluxes. This indicates that even the flare with less energy could be associated with anomalistic transportation of the magnetic helicity across the photosphere. Accompanying the impulsive helicity fluxes, the poynting fluxes across the photosphere evolved from positive to negative. As such, the transportations of magnetic energy across the photosphere were toward solar interior during these flares. In each of the 4 events, the impulsive change in the helicity flux was always mainly contributed by abrupt change in horizontal velocity field on a sunspot located near the flaring polarity inversion line. The velocity field on each sunspot shows either an obvious vortex patten or an shearing patten relative to the another magnetic polarity, which tended to relax the magnetic twist or shear in the corona. During these flares, abrupt change in the Lorentz force acting on these sunspots were found. The rotational motions and shearing motions of these sunspots always had the same directions with the resultant Lorentz forces. These results support the view that the impulsive helicity transportation during the flare could be driven by the change in the Lorentz force applied on the photosphere.

We present observations of the formation process of a small-scale filament on the quiet Sun during 5-6 February 2016 and investigate its formation cause. Initially, a small dipole emerged and its associated arch filament system was found to reconnect with overlying coronal fields accompanied by numerous EUV bright points. When bright points faded out, many elongated dark threads formed bridging the positive magnetic element of dipole and external negative network fields. Interestingly, an anti-clockwise photospheric rotational motion (PRM) set in within the positive endpoint region of newborn dark threads following the flux emergence and lasted for more than 10 hours. Under the drive of the PRM, these dispersive dark threads gradually aligned along the north-south direction and finally coalesced into an inverse S-shaped filament. Consistent with the dextral chirality of the filament, magnetic helicity calculations show that an amount of negative helicity was persistently injected from the rotational positive magnetic element and accumulated during the formation of the filament. These observations suggest that twisted emerging fields may lead to the formation of the filament via reconnection with pre-existing fields and release of its inner magnetic twist. The persistent PRM might trace a covert twist relaxation from below photosphere to the low corona.

This paper presents a fast algorithm for obtaining high-accuracy subpixel translation of low PSNR images. Instead of locating the maximum point on the upsampled images or fitting the peak of correlation surface, the proposed algorithm is based on the measurement of centroid on the cross correlation surface by Modified Moment method. Synthetic images, real solar images and standard testing images with white Gaussian noise added were tested, and the results show that the accuracies of our algorithm are comparable with other subpixel registration techniques and the processing speed is higher. The drawback is also discussed at the end of this paper.

Solid evidence of magnetic reconnection is rarely reported within sunspots, the darkest regions with the strongest magnetic fields and lowest temperatures in the solar atmosphere. Using the world's largest solar telescope, the 1.6-meter Goode Solar Telescope, we detect prevalent reconnection through frequently occurring fine-scale jets in the H${\alpha}$ line wings at light bridges, the bright lanes that may divide the dark sunspot core into multiple parts. Many jets have an inverted Y-shape, shown by models to be typical of reconnection in a unipolar field environment. Simultaneous spectral imaging data from the Interface Region Imaging Spectrograph show that the reconnection drives bidirectional flows up to 200~km~s$^{-1}$, and that the weakly ionized plasma is heated by at least an order of magnitude up to $\sim$80,000 K. Such highly dynamic reconnection jets and efficient heating should be properly accounted for in future modeling efforts of sunspots. Our observations also reveal that the surge-like activity previously reported above light bridges in some chromospheric passbands such as the H${\alpha}$ core has two components: the ever-present short surges likely to be related to the upward leakage of magnetoacoustic waves from the photosphere, and the occasionally occurring long and fast surges that are obviously caused by the intermittent reconnection jets.

In this paper we illustrate the technique used by the New Vacuum Solar Telescope to increase the spatial resolution of two-dimensional (2D) solar spectroscopy observation involving two dimensions of space and one of wavelength. Without an image stabilizer at the NVST, a large scale wobble motion is present during the spatial scanning, whose instantaneous amplitude could reach up to 1.3" due to the earth's atmosphere and the precision of the telescope guiding system, and seriously decreases the spatial resolution of 2D spatial maps composed with the scanning spectra. We make the following effort to resolve this problem: the imaging system (e.g., the TiO-band) is used to record and detect the displacement vectors of solar image motion during the raster scan, in both the slit and scanning directions. The spectral data (e.g., the Ha line) which are originally obtained in time sequence are corrected and re-arranged in space according to those displacement vectors. Raster scans are carried out in several active regions with different seeing conditions (two rasters are illustrated in this paper). Given a certain spatial sample and temporal resolution, the spatial resolution of the composed 2D map could be close to that of the slit-jaw image. The resulting quality after correction is quantitatively evaluated with two methods. Two-dimensional physical quantity, such as the line-of-sight velocities in multi-layer of the solar atmosphere, is also inferred demonstrating the effect of this technique.

We present a precise and complete procedure for processing spectral data observed by the 1-meter New Vacuum Solar Telescope (NVST). The procedure is suitable for both the sit-and-stare and raster-scan spectra. In this work, the geometric distortions of the spectra are firstly corrected for subsequent processes. Then, considering the temporal changes and the remnants of spectral lines in the flat-field, the original flat-field matrix is split into four independent components for ensuring a high precision flat-fielding correction, consisting of the continuum gradient matrix, slit non-uniform matrix, CCD dust matrix, and interference fringe matrix. Subsequently, the spectral line drifts and intensity fluctuations of the science data are further corrected. After precise reduction with this procedure, the measuring accuracies of the Doppler velocities for different spectral lines and of the oscillation curves of the chromosphere and photosphere are measured. The results show that the highest measuring accuracy of the Doppler velocity is within 100 ms-1, which indicates that the characteristics of the photosphere and chromosphere can be studied co-spatially and co-temporally with the reduced spectra of NVST.

All next generation ground-based and space-based solar telescopes require a good quality assessment metric in order to evaluate their imaging performance. In this paper, a new image quality metric, the median filter gradient similarity (MFGS) is proposed for photospheric images. MFGS is a no-reference/blind objective image quality metric (IQM) by a measurement result between 0 and 1 and has been performed on short-exposure photospheric images captured by the New Vacuum Solar Telescope (NVST) of the Fuxian Solar Observatory and by the Solar Optical Telescope (SOT) onboard the Hinode satellite, respectively. The results show that: (1)the measured value of MFGS changes monotonically from 1 to 0 with degradation of image quality; (2)there exists a linear correlation between the measured values of MFGS and root-mean-square-contrast (RMS-contrast) of granulation; (3)MFGS is less affected by the image contents than the granular RMS-contrast. Overall, MFGS is a good alternative for the quality assessment of photospheric images.

The New Vacuum Solar Telescope (NVST) is a 1-m solar telescope that aims to observe the fine structures in both the photosphere and the chromosphere of the Sun. The observational data acquired simultaneously from one channel for the chromosphere and two channels for the photosphere bring great challenges to the data storage of NVST. The multi-channel instruments of NVST, including scientific cameras and multi-band spectrometers, generate at least 3 terabytes data per day and require high access performance while storing massive short-exposure images. It is worth studying and implementing a storage system for NVST which would balance the data availability, access performance and the cost of development. In this paper, we build a distributed data storage system (DDSS) for NVST and then deeply evaluate the availability of real-time data storage on a distributed computing environment. The experimental results show that two factors, i.e., the number of concurrent read/write and the file size, are critically important for improving the performance of data access on a distributed environment. Referring to these two factors, three strategies for storing FITS files are presented and implemented to ensure the access performance of the DDSS under conditions of multi-host write and read simultaneously. The real applications of the DDSS proves that the system is capable of meeting the requirements of NVST real-time high performance observational data storage. Our study on the DDSS is the first attempt for modern astronomical telescope systems to store real-time observational data on a low-cost distributed system. The research results and corresponding techniques of the DDSS provide a new option for designing real-time massive astronomical data storage system and will be a reference for future astronomical data storage.

The New Vacuum Solar Telescope (NVST) is a 1-meter vacuum solar telescope that aims to observe the fine structures of active regions on the Sun. The main tasks of the NVST are high resolution imaging and spectral observations, including the measurements of the solar magnetic field. The NVST has been collecting more than 20 million FITS files since it began routine observations in 2012 and produces a maximum observational records of 120 thousand files in a day. Given the large amount of files, the effective archiving and retrieval of files becomes a critical and urgent problem. In this study, we implement a new data archiving system for the NVST based on the Fastbit Not Only Structured Query Language (NoSQL) database. Comparing to the relational database (i.e., MySQL; My Structured Query Language), the Fastbit database manifests distinctive advantages on indexing and querying performance. In a large scale database of 40 million records, the multi-field combined query response time of Fastbit database is about 15 times faster and fully meets the requirements of the NVST. Our study brings a new idea for massive astronomical data archiving and would contribute to the design of data management systems for other astronomical telescopes.

The Chinese Spectral RadioHeliograph (CSRH) is a synthetic aperture radio interferometer built in Inner Mongolia, China. As a solar-dedicated interferometric array, CSRH is capable of producing high quality radio images at frequency range from 400 MHz to 15 GHz with high temporal, spatial, and spectral resolution.To implement high cadence imaging at wide-band and obtain more than 2 order higher multiple frequencies, the implementation of the data processing system for CSRH is a great challenge. It is urgent to build a pipeline for processing massive data of CSRH generated every day. In this paper, we develop a high performance distributed data processing pipeline (DDPP) built on the OpenCluster infrastructure for processing CSRH observational data including data storage, archiving, preprocessing, image reconstruction, deconvolution, and real-time monitoring. We comprehensively elaborate the system architecture of the pipeline and the implementation of each subsystem. The DDPP is automatic, robust, scalable and manageable. The processing performance under multi computers parallel and GPU hybrid system meets the requirements of CSRH data processing. The study presents an valuable reference for other radio telescopes especially aperture synthesis telescopes, and also gives an valuable contribution to the current and/or future data intensive astronomical observations.

In this paper, we report our first-step results of high resolution He\,\textsc{i} 10830 \AA\ narrow-band imaging (bandpass: 0.5 {\AA}) of an M1.8 class two-ribbon flare on July 5, 2012. The flare was observed with the 1.6 meter aperture New Solar Telescope at Big Bear Solar Observatory. For this unique data set, sunspot dynamics during flaring were analyzed for the first time. By directly imaging the upper chromosphere, running penumbral waves are clearly seen as an outward extension of umbral flashes, both take the form of absorption in the 10830 \AA\ narrow-band images. From a space-time image made of a slit cutting across a flare ribbon and the sunspot, we find that the dark lanes for umbral flashes and penumbral waves are obviously broadened after the flare. The most prominent feature is the sudden appearance of an oscillating absorption strip inside the ribbon when it sweeps into the sunspot's penumbral and umbral regions. During each oscillation, outwardly propagating umbral flashes and subsequent penumbral waves rush out into the inwardly sweeping ribbon, followed by a returning of the absorption strip with similar speed. We tentatively explain the phenomena as the result of a sudden increase in the density of ortho-Helium atoms in the area of the sunspot being excited by the flare's EUV illumination. This explanation is based on the observation that 10830 \AA\ absorption around the sunspot area gets enhanced during the flare. Nevertheless, questions are still open and we need further well-devised observations to investigate the behavior of sunspot dynamics during flares.

G-band bright points (GBPs) are regarded as good manifestations of magnetic flux concentrations. We aim to investigate the relationship between the dynamic properties of GBPs and their longitudinal magnetic field strengths. High spatial and temporal resolution observations were recorded simultaneously with G-band filtergrams and Narrow-band Filter Imager (NFI) Stokes I and V images with Hinode /Solar Optical Telescope. The GBPs are identified and tracked in the G-band images automatically, and the corresponding longitudinal magnetic field strength of each GBP is extracted from the calibrated NFI magnetograms by a point-to-point method. After categorizing the GBPs into five groups by their longitudinal magnetic field strengths, we analyze the dynamics of GBPs of each group. The results suggest that with increasing longitudinal magnetic field strengths of GBPs correspond to a decrease in their horizontal velocities and motion ranges as well as by showing more complicated motion paths. This suggests that magnetic elements showing weaker magnetic field strengths prefer to move faster and farther along straighter paths, while stronger ones move more slowly in more erratic paths within a smaller region. The dynamic behaviors of GBPs with different longitudinal magnetic field strengths can be explained by that the stronger flux concentrations withstand the convective flows much better than weaker ones.

Umbral dots (UDs) are small isolated brightenings observed in sunspot umbrae. They are convective phenomena existing inside umbrae. UDs are usually divided into central UDs (CUDs) and peripheral UDs (PUDs) according to their positions inside an umbra. Our purpose is to investigate UD properties and analyze their relationships, and further to find whether or not the properties depend on umbral magnetic field strengths. Thus, we selected high-resolution TiO images of four active regions (ARs) taken under the best seeing conditions with the \textit{New Vacuum Solar Telescope} in the Fuxian Solar Observatory of the Yunnan Astronomical Observatory, China. The four ARs (NOAA 11598, 11801, 12158, and 12178) include six sunspots. A total of 1220 CUDs and 603 PUDs were identified. Meanwhile, the radial component of the vector magnetic field of the sunspots taken with the \textit{Helioseismic and Magnetic Imager} on-board the \textit{Solar Dynamics Observatory} was used to analyze relationships between UD properties and umbral magnetic field strengths. We find that diameters and lifetimes of UDs exhibit an increasing trend with the brightness, but velocities do not. Moreover, diameters, intensities, lifetimes and velocities depend on the surrounding magnetic field. A CUD diameter was found larger, the CUD brighter, its lifetime longer, and its motion slower in a weak umbral magnetic field environment than in a strong one.

Six high-resolution TiO-band image sequences from the New Vacuum Solar Telescope (NVST) are used to investigate the properties of intergranular bright points (igBPs). We detect the igBPs using a Laplacian and morphological dilation algorithm (LMD) and track them using a three-dimensional segmentation algorithm automatically, and then investigate the morphologic, photometric and dynamic properties of igBPs, in terms of equivalent diameter, the intensity contrast, lifetime, horizontal velocity, diffusion index, motion range and motion type. The statistical results confirm the previous studies based on G-band or TiO-band igBPs from the other telescopes. It illustrates that the TiO data from the NVST have a stable and reliable quality, which are suitable for studying the igBPs. In addition, our method is feasible to detect and track the igBPs in the TiO data from the NVST. With the aid of the vector magnetograms obtained from the Solar Dynamics Observatory /Helioseismic and Magnetic Imager, the properties of igBPs are found to be influenced by their embedded magnetic environments strongly. The area coverage, the size and the intensity contrast values of igBPs are generally larger in the regions with higher magnetic flux. However, the dynamics of igBPs, including the horizontal velocity, the diffusion index, the ratio of motion range and the index of motion type are generally larger in the regions with lower magnetic flux. It suggests that the absence of strong magnetic fields in the medium makes it possible for the igBPs to look smaller and weaker, diffuse faster, move faster and further in a straighter path.

We study the motions of G band bright points (GBPs) in the quiet Sun to obtain the characteristics of different motion types. A high resolution image sequence taken with the Hinode/Solar Optical Telescope (SOT) is used, and GBPs are automatically tracked by segmenting 3D evolutional structures in a space time cube. After putting the GBPs that do not move during their lifetimes aside, the non stationary GBPs are categorized into three types based on an index of motion type. Most GBPs that move in straight or nearly straight lines are categorized into a straight motion type, a few moving in rotary paths into a rotary motion, and the others fall into a motion type we called erratic. The mean horizontal velocity is 2.18 km/s, 1.63 km/s and 1.33 km/s for straight, erratic and rotary motion type, respectively. We find that a GBP drifts at a higher and constant velocity during its whole life if it moves in a straight line. However, it has a lower and variational velocity if it moves in a rotary path. The diffusive process is ballistic, super and sub diffusion for straight, erratic and rotary motion type, respectively. The corresponding diffusion index and coefficients are 2.13 and 850 km2/s, 1.82 and 331 km2/s, 0.73 and 13 km2/s. In terms of direction of motion, it is homogeneous and isotropical, and usually persists between neighbouring frames, no matter what motion type a GBP belongs to.