The paper presents the results of the longitudinal development of extensive air showers (X$_{max}$) of ultra-high energies and mass composition of cosmic rays. The measurements of X$_{max}$ are based on data from observations of the Cherenkov radiation at the Yakutsk array for the period 1974-2014. The cascade curves of individual showers and the depth of maximum X$_{max}$ were reconstructed over the energy range 10$^{16}$-5.7$\cdot$10$^{19}$ eV. It is shown that the displacement rate of the parameter dX$_{max}$ / dE in the atmosphere is nonlinear and depends on the energy. Such a feature indicates a change in mass composition, which is confirmed by fluctuations of X$_{max}$ in this energy region. The composition of cosmic rays was determined by interpolation using the QGSJetII-04 model.

The Yakutsk Array is designed to study cosmic rays at energy 10$^{15}$-10$^{20}$ eV. It consists several independent arrays that register charged particles, muons with energy E$\geq$1 GeV, Cherenkov light and radio emission. The paper presents a technical description of the Yakutsk Radio Array and some preliminary results obtained from measurements of radio emission at 30-35 MHz frequency induced by air shower particles with energy $\varepsilon$ $\geq$ 1$\cdot$10$^{17}$ eV. The data obtained at the Yakutsk array in 1986-1989 (first set of measurements) and 2009-2014 (new set of measurements). Based on the obtained results we determined: Lateral distribution function (LDF) of air showers radio emission with energy $\geq$ 10$^{17}$ eV. Radio emission amplitude empirical connection with air shower energy. Determination of depth of maximum by the ratio of amplitude at different distances from the shower axis. For the first time, at the Yakutsk array, radio emission from the air shower with energy $>$ 10$^{19}$ eV was registered including the shower with the highest energy ever registered at the Yakutsk array with energy $\sim$2 $\cdot$ 10$^{20}$ eV.

Long-term series of measurement of spectral transparency of the atmosphere (\lambda = 430 nm) and atmospheric optical thickness (AOT) measured by multimode photometer CE 318 in the region of Yakutsk array are analyzed. Correlation of AOT with intensity of air showers with small energies 10^{15} - 10^{16} eV is found. The variability of aerosol composition of the atmosphere during the registration period of the Cherenkov light should be taken into account since it may affect the quality of determining characteristics of air showers.

This paper presents data on the optical thickness AOT in the context of long-term observations, including in times of large-scale forest fires in Yakutia. Sudden changes in such features as the AOT point directly to the impact of fires on aerosol structure of the atmosphere during the summer.

Neutrino air showers can be formed in any part of the atmosphere passing a long was in the matter due to its physical properties. In general, air showers produced by neutrinos are highly inclined and formed near the ground level, i.e. young showers. Therefore, one should expect a large number of peaks in the signal of such air showers [1, 2, 3]. The goal of our work is to search for air shower candidates produced by neutrino. For this purposes, we analyzed large amount of data from scintillation detectors with different area and energy thresholds [1, 2]. Preliminary analysis of Yakutsk array data indicated the absence of air shower produced by neutrino, but it does not mean that such air showers does not exist. It is going to need a further analysis of highly inclined showers. In order to do that improved methodology for recording and processing of air showers required.

In this paper, we present results obtained from the measurements of radio emission at frequency of 32 MHz with energy more than 10$^{19}$ eV. Generalized formula that describe lateral distribution and depends on main characteristic of the air showers: energy E$_0$ and depth of maximum X$_{max}$ was derived. The formula has a good agreement with data at average and large distances from shower axis. Employing the ratio of radio emission amplitude at distances 175 m and 725 m we determined the depth of maximum X$_{max}$ for air shower with energy 3.7$\cdot$10$^{19}$ eV, which in our case is equal to X$_{max}$ = 769$\pm$34g$\cdot$cm$^{-2}$.

The paper presents short technical description of Yakutsk Radio Array and some preliminary results obtained from measurements of radio emission at 32 MHz frequency induced by air shower particles with energy {\epsilon}$\geq$1$\cdot$10$^{19}$ eV. The data obtained at Yakutsk array in 1987-1989 (first set of measurements) and 2009-2014 (new set of measurements). For the first time, at Yakutsk array radio emission from air shower with energy > 10$^{19}$ eV was registered including the shower with highest energy ever registered at Yakutsk array with energy $\sim$2$\cdot$10$^{20}$ eV.

The paper presents long-term observations of the atmosphere in Yakutsk region. Analysis of the data for 40 year period indicates a gradual strengthening of cyclonic activity in the region and hence the increase of the average winter temperature, increase variations of the rest atmosphere, which greatly softens the continental climate of Central Yakutia.

Calibration of a PMT matrix is crucial for the treatment of the data obtained with Cherenkov tracking detector. Furthermore, due to high variability of the aerosol abundance in the atmosphere depending on season, weather etc. A constant monitoring of the atmospheric transparency is required during the measurements. For this purpose, besides traditional methods, a station for laser atmospheric probing is used.

The paper present correlation of radio signal with air shower parameters: shower energy E0 and depth of maximum Xmax. It is shown that from radio emission measurements of air showers one can obtain individual showers parameters and mass composition of cosmic rays. We also derived generalized formula for calculating energy of the air showers.

The paper describes the techniques and method of registration of air shower radio emission at the Yakutsk array of extensive air showers at a frequency of 32 MHz. At this stage, emission registration involves two set of antennas, the distance between them is 500m. One set involves 8 antennas, second - 4 antennas. The antennas are perpendicularly crossed dipoles with radiation pattern North South,West East and raised 1.5 m above the ground. Each set of antennas connected to an industrial PC. The registration requires one of two triggers. First trigger are generated by scintillation detectors of Yakutsk array. Scintillation detectors cover area of 12 km^2 and registers air showers with energy more than 10^17 eV. The second trigger is generated by Small Cherenkov Array that covers area of 1 km^2 and registers air showers with energy 10^15 - 5*10^17 eV. Small Cherenkov Array is part of Yakutsk array and involve Cherenkov detectors located at a distance of 50, 100, 250 m. For further selection we are using an additional criterion the radio pulse must be localized in the area corresponding to the delay time on first and second triggers. In addition, descriptions of the algorithm and the flowcharts of the program for the air shower selection and further analysis are given. This method registers EAS radio emission with energy 10^16 - 10^19 eV. With the absolute calibration, the amplitudes of all antennas converted to a single value. Air shower radio emission dependences from zenith angle and shower energy are plotted.

This paper presents the set of measurements of ultra-high energy air shower radio emission at frequency 32 MHz in period of 2008-2012. The showers are selected by geomagnetic and azimuth angles and then by the energy in three intervals: 3*10^16 3*10^17 eV, 3*10^17 6*10^17 eV and 6*10^17 5*10^18 eV. In each energy interval average lateral distribution function using mathematically averaged data from antennas with different directions are plotted. In the paper, using experimental data the dependence of radio signal averaged amplitude from geomagnetic angle, the shower axis distance and the energy are determined. Depth of maximum of cosmic ray showers Xmax for the given energy range is evaluated. The evaluation is made according QGSJET model calculations and average lateral distribution function shape.

We present of measurements of radio emission from EAS at the frequency 32 MHz performed in 2008-2012. Showers were selected by zenith and azimuth directions, sampled by energy. A mean lateral distribution function was fitted to the data each interval. The paper presents unied approximation formula for all three energy intervals of lateral distribution, using connection between equation (1) and (2) coefficients with energy and geomagnetic angle. Here we analyze the connection between the shape of lateral distribution and the depth of the maximum shower development Xmax too.