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.

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 Yakutsk array group is developing a wide field-of-view Cherenkov telescope to be operated in coincidence with the surface detectors of the extensive air shower array. Currently, the engineering prototype of the reflecting telescope with the front-end electronics is designed, assembled, and tested to demonstrate the feasibility of the conceived instrument. The status and specifications of the prototype telescope are presented, as well as the modernization program of the already existing Cherenkov light detectors subset of the array measuring ultra-high energy cosmic rays.

From the data on Cherenkov light from EASs collected over the period of more than 15 years, the spectrum of cosmic rays was obtained in the energy range 1.0E15-1.0E18 eV. This spectrum has two features: at ~3.0E15 eV (first knee) and at ~1.0E17 eV (second knee). The first knee is characterized by the indexes \gamma_1=2.7+/-0.03 and \gamma_2 = 3.12+/-0.03 and the second knee - by \gamma_1 = 2.92+/-0.03 and \gamma_2=3.24+/-0.04. In the second case the difference amounts to \Delta\gamma_23=0.32{\pm}0.03{\pm}0.05 which is less than the case of the first knee \Delta\gamma_12=0.42+/-0.03+/-0.05 (here dual errors represent statistics and systematics correspondingly). The lesser difference \Delta\gamma_23 can be explained with the influx of cosmic rays from Meta-galaxy and, hence, with some increase of the cosmic rays intensity in the energy range 5.0E16-1.0E18 eV, which compensates the escape of heavier nuclei from the Galaxy. The presence of the second knee could be confirmed by abrupt change in mass composition from from lnA~3 at ~1.0E17 eV to lnA~1.5 at 1.0E18eV found from the analysis of longitudinal development of EAS.

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.

The aim of the Yakutsk array enhancement project is to create an instrument to study the highest-energy galactic cosmic rays (CRs) -- their sources, energy spectrum, and mass composition. Additionally, there will be unique capabilities for investigations in the transition region between galactic and extragalactic components of CRs. Using the well-developed imaging atmospheric Cherenkov telescope technique adapted to the energy region $E>10^{15}$ eV, we plan to measure the longitudinal structure parameters of the shower, e.g., angular and temporal distributions of the Cherenkov signal related to $X_{max}$ and the mass composition of CRs. The main advantages of the Yakutsk array, such as its multi-component measurements of extensive air showers, and model-independent CR energy estimation based on Cherenkov light measurements, will be inherited by the instrument to be created.

Signal from scintillation detectors in showers with energies 10**17-10**20eV has been analyzed. It was determined, that for core distances >=200m the time scale of a particle arrival time at the detector constitutes 40-5000ns and longer. Events with fine structure of the signal have been observed. The pulse consists of several peaks (for vertical showers) and a single pulse (for inclined showers). There are some anomalous events with two pulses wherein such a pulse shape is presented in two, three and more detectors. A period between these pulses makes 60-350ns.

The experimental distributions of $\Xmax$ obtained with the Yakutsk EAS array at fixed energies of $5 \times10^{17}$, $1\times10^{18}$ and $5\times10^{18}$ eV are analysed. A recent version of the QGSJET model is used as a tool of our analysis. In the framework of this model, the most adequate mass composition of primary particles satisfying the experimental data on $\Xmax$ is selected.

At the Yakutsk EAS array E_em is determined by using measurements of EAS Cherenkov light flux and charged particle flux. It is known from calculations that these characteristics depend on a sort of primary particle and, therefore, the estimation of E_em depends on a primary particle mass. In the work the dependence of the E_em/E_0 ratio on the energy is given and experimental data are compared with calculations by the QGSJET model. The calculations have been carried out for the primary proton and iron nucleus. The average calculated meaning of the value of E_em/E_0 ratio (between the proton and iron nucleus) within experimental errors is in agreement with experimental data that doesnt contradict to the mixed mass composition of primary cosmic radiation.

In the work the data of the Yakutsk complex EAS array and their comparison with calculation in the case of primary nuclei of different chemical elements are presented. The calculation by QGSJET model have been used interpreting experimental data.