The results of a search for individual fast radio bursts with the Large Phased Array of the Lebedev Physical Institute at 111 MHz during July 2012 through August 2018 are presented. The signals were distinguished by convolving the data with a template with a fixed form, followed by convolution with test dispersion measures. Areas of sky containing the galaxies M31 and M33 were chosen for the search. Three radio bursts were detected in the vicinity of M33, five in the vicinity of M31, and one in a region offset from the center of M31 by an hour in right ascension. The dispersion measures of the detected bursts range from 203 to 1262 $ pc \cdot cm^{-3}$.

We present results of simulations with the Max Planck Institute general circulation model (MPI-MGCM) implementing a hydrological cycle scheme. The simulations reveal a seasonal water "pump" mechanism responsible for the upward transport of water vapor. This mechanism occurs in high latitudes above 60$^\circ$ of the southern hemisphere at perihelion, when the upward branch of the meridional circulation is particularly strong. A combination of the mean vertical flux with variations induced by solar tides facilitates penetration of water across the "bottleneck" at approximately 60 km. The meridional circulation then transports water across the globe to the northern hemisphere. Since the intensity of the meridional cell is tightly controlled by airborne dust, the water abundance in the thermosphere strongly increases during dust storms.

Results of a search for individual impulsive signals on the Large Scanning Antenna of the Lebedev Physical Institute at 111 MHz carried out from July 2012 through May 2018 are presented. The data were convolved with a template of a specified form and convolved with a test dispersion measure. A region of sky with central coordinates $\alpha = 05^h 32^m;$ $\delta = +41.72^\circ$ and also a region of sky around the coordinates fixed earlier for FRB 121102 ($\alpha = 05^h 32^m;$ $\delta = +33.1^\circ$) were chosen for the analysis. In all, 355 hours of observations were processed for each beam. Three radio bursts with dispersion measures of $247$ $pc \cdot cm^{-3}$, $570$ $ pc \cdot cm^{-3}$, $1767$ $ pc \cdot cm^{-3}$ were detected in the course of reducing the data.

Timing of highly stable millisecond pulsars provides the possibility of independently verifying terrestrial time scales on intervals longer than a year. An ensemble pulsar time scale is constructed based on pulsar timing data obtained on the 64-m Parkes telescope (Australia) in 1995--2010. Optimal Wiener filters were applied to enhance the accuracy of the ensemble time scale. The run of the time--scale difference PT$_{ens}$ $-$ TT(BIPM2011) does not exceed $0.8 {\pm} 0.4$ $\mu s$ over the entire studied time interval. The fractional instability of the difference $_{ens}$ $-$ TT(BIPM2011) over 15 years is ${\sigma_z} = (0.6\pm 1.6){\cdot}10^{-15}$, which corresponds to an upper limit for the energy density of the gravitational--wave background ${\Omega_g}h^2 {\sim} 10^{-10}$ and variations in the gravitational potential ${\sim} 10^{-15}$ Hz at the frequency $2{\cdot}10^{-9}$.

Results of timing measurements of the pulsar PSR B0329+54 obtained in 1968--2012 using the Big Scanning Antenna of the Pushchino Radio Astronomy Observatory (at 102 and 111 MHz), the DSS 13 and DSS 14 telescopes of the Jet Propulsion Laboratory (2388 MHz), and the 64 m telescope of the Kalyazin Radio Astronomy Observatory (610 MHz) are presented. The astrometric and rotational parameters of the pulsar are derived at a new epoch. Periodic variations in the barycentric timing residuals have been found, which can be explained by the presence of a planet orbiting the pulsar, with an orbital period $P_{1}$ = 27.8 yr, mass \textit{$m_{c}$}sin\textit{i} = 2$M_{\oplus}$, and orbital semi-major axis $a$ = 10.26 AU. The results of this study do not confirm existence of a proposed second planet with orbital period $P_{2}$ = 3 yr.

The algorithm of the ensemble pulsar time scale (PT$_{\rm ens}$) based on the optimal Wiener filtration method has been proposed. This algorithm allows the separation of the contributions to the post-fit pulsar timing residuals of the atomic clock and pulsar itself. Filters were designed with the use of the cross-spectra of the timing residuals. The method has been applied to the timing data of six millisecond pulsars. Direct comparison with the classical method of the weighted average showed that use of the optimal Wiener filters before averaging allows noticeably to improve the fractional instability of the ensemble time scale. Application of the proposed method to the most stable millisecond pulsars with the fractional instability $\sigma_z < 10^{-15}$ may improve the fractional instability of PT$_{\rm ens}$ up to the level $\sim 10^{-16}$.

An algorithm is proposed for constructing a group (ensemble) pulsar time based on the application of optimal Wiener filters. This algorithm makes it possible to separate the contributions of variations of the atomic time scale and of the pulsar rotation to barycentric residual deviations of the pulse arrival times. The method is applied to observations of the pulsars PSR B1855+09 and PSR B1937+21, and is used to obtain corrections to UTC relative to the group pulsar time PT$_{\rm ens}$. Direct comparison of the terrestrial time TT(BIPM06) and the group pulsar time PT$_{\rm ens}$ shows that they disagree by no more than $0.4\pm 0.17\; \mu$s. Based on the fractional instability of the time difference TT(BIPM06) -- PT$_{\rm ens}$, a new limit for the energy density of the gravitational-wave background is established at the level $\Omega_g {h}^2\sim 10^{-9}$.

We suggest a new approach to the detection of gravitational waves using observations of a group of millisecond pulsars. In contrast to the usual method, based on increasing the accuracy of the arrival times of pulses by excluding possible distorting factors, our method supposes that the additive phase noise that is inevitably present even in the most accurate observational data has various spectral components, which have characteristic amplitudes and begin to appear on different time scales. We use the "Caterpillar" (Singular Spectral Analysis, SSA) method to decompose the signal into its components. Our initial data are the residuals of the pulse arrival times for six millisecond pulsars. We constructed the angular correlation function for components of the decomposition of a given number, whose theoretical form for the case of an isotropic and homogeneous gravitational-wave background is known. The individual decomposition components show a statistically significant agreement with the theoretical expectations (correlation coefficient $\rho=0.92\pm 0.10$).

A new approach to the problem of gravitational waves detection based on simultaneous timing of several pulsars and subsequent expansion of the post-fit timing data into components of different spectral kind (with different spectral indices) is proposed. Presence of a signal caused by stochastic gravitational waves in spectral components is tested with the two-point angular correlation function as proposed in the pulsar timing array. This new approach was applied to timing data of a few millisecond pulsars and allowed to detect a signature similar to one predicted for gravitational wave background at relatively high confidence level: correlation coefficient between experimental and theoretical two-point correlations $\rho=0.82\pm 0.07$.

The algorithm of the ensemble pulsar time based on the optimal Wiener filtration method has been constructed. This algorithm allows the separation of the contributions to the post-fit pulsar timing residuals of the atomic clock and pulsar itself. Filters were designed with the use of the cross- and autocovariance functions of the timing residuals. The method has been applied to the timing data of millisecond pulsars PSR B1855+09 and PSR B1937+21 and allowed the filtering out of the atomic scale component from the pulsar data. Direct comparison of the terrestrial time TT(BIPM06) and the ensemble pulsar time PT$_{\rm ens}$ revealed that fractional instability of TT(BIPM06)--PT$_{\rm ens}$ is equal to $\sigma_z=(0.8\pm 1.9)\cdot 10^{-15}$. Based on the $\sigma_z$ statistics of TT(BIPM06)--PT$_{\rm ens}$ a new limit of the energy density of the gravitational wave background was calculated to be equal to $\Omega_g h^2 \sim 3\cdot 10^{-9}$.

An algorithm of the ensemble pulsar time based on the Wiener filtration method has been constructed. This algorithm has allowed the separation of the contributions of an atomic clock and a pulsar itself to the post-fit pulsar timing residuals. The method has been applied to the timing data of the millisecond pulsars PSR B1855+09 and PSR B1937+21 and allowed the filtering out of the atomic scale component from the pulsar phase variations. Direct comparison of the terrestrial time TT(BIPM96) and the ensemble pulsar time PT$_{\rm ens}$ has displayed that the difference TT(BIPM96) -- PT$_{\rm ens}$ is within $\pm0.4 \mu$s range. A new limit of gravitational wave background based on the difference TT(BIPM96) -- PT$_{\rm ens}$ was established to be $\Omega_g {\rm h}^2\sim 10^{-10}$.